1、第 卷第期 年月有色金属(矿山部分)N o n f e r r o u sM e t a l s(M i n i n gS e c t i o n)V o l ,N o M a y d o i:/j i s s n 无人机摄影测量在尾矿坝变形监测分析中的应用廖文景,戴军,唐鹏飞(长沙矿山研究院有限责任公司,长沙 ;湖南省矿山地质灾害防治与环境再造工程技术研究中心,长沙 ;金川集团股份有限公司,甘肃 金昌 )摘要:由于尾矿坝安全管理范围大,山地环境特殊,传统测绘手段难以开展,尾矿坝变形位移数据仅通过单点或多点的监测方式获得,因此本文提出一种基于无人机(UAV)摄影测量的尾矿坝变形监测分析方法.该
2、方法通过无人机对尾矿坝两期空间三维地理信息数据的采集,将采集到的两期数据分别进行三维重建并生成点云模型,再利用地理信息三维时空对比分析技术对两期点云模型进行对比分析.结果表明:该方法实现了对尾矿坝由点到面的全方位监测,监测精度可达厘米级,为尾矿坝宏观变形监测及异常变形区域的快速核查及安全巡检等提供了依据.关键词:尾矿坝;无人机;摄影测量;点云模型;变形分析中图分类号:T D 文献标志码:A文章编号:()A p p l i c a t i o no fU A Vp h o t o g r a mm e t r y i nm o n i t o r i n ga n da n a l y s i
3、so ft a i l i n g sd a md e f o r m a t i o nL I AO W e n j i n g,D A I J u n,T ANGP e n g f e i(C h a n g s h aI n s t i t u t eo fM i n i n gR e s e a r c hC o,L t d,C h a n g s h a ,C h i n a;H u n a nM i n eG e o l o g i c a lD i s a s t e rP r e v e n t i o na n dE n v i r o n m e n t a lR e c
4、o n s t r u c t i o nE n g i n e e r i n gT e c h n o l o g yR e s e a r c hC e n t e r,C h a n g s h a ,C h i n a;J i n c h u a nG r o u pC o,L t d,J i n c h a n gG a n s u ,C h i n a)A b s t r a c t:D u e t o t h e l a r g e s c o p eo f s a f e t ym a n a g e m e n t o f t a i l i n g sd a ma n d
5、 t h e s p e c i a lm o u n t a i n o u s e n v i r o n m e n t,t r a d i t i o n a l s u r v e y i n ga n dm a p p i n gm e a n sa r ed i f f i c u l tt oc a r r yo u t,a n dt h ed e f o r m a t i o na n dd i s p l a c e m e n td a t ao ft a i l i n g sd a mc a no n l yb eo b t a i n e dt h r o u g
6、 hs i n g l eo rm u l t i p o i n tm o n i t o r i n g T h e r e f o r e,t h i sp a p e rp r o p o s e sat a i l i n g sd a m d e f o r m a t i o n m o n i t o r i n ga n da n a l y s i s m e t h o d b a s e do n UAV p h o t o g r a mm e t r y,w h i c hc o l l e c t st h et h r e e d i m e n s i o n
7、 a lg e o g r a p h i ci n f o r m a t i o nd a t ao ft w op h a s e so fat a i l i n g sd a mb yUAV T h ec o l l e c t e dd a t ao ft h et w op e r i o d sw e r er e c o n s t r u c t e d i nt h r e ed i m e n s i o n sa n dt h ep o i n tc l o u dm o d e lw a sg e n e r a t e d T h ep o i n tc l o
8、u dm o d e lo ft h e t w op e r i o d sw a sc o m p a r e da n da n a l y z e db yu s i n gt h et h r e e d i m e n s i o n a ls p a t i o t e m p o r a lc o m p a r a t i v ea n a l y s i st e c h n o l o g yo fg e o g r a p h i c i n f o r m a t i o n T h er e s u l t s s h o wt h a t t h em e t h
9、 o dc a nr e a l i z e t h e c o m p r e h e n s i v em o n i t o r i n go ft a i l i n g sd a mf r o mp o i n t t os u r f a c e,a n dt h em o n i t o r i n ga c c u r a c yc a nr e a c hc ml e v e l,w h i c hp r o v i d e sab a s i sf o rt h em a c r od e f o r m a t i o n m o n i t o r i n go ft
10、h et a i l i n g sd a m a n dt h er a p i dv e r i f i c a t i o na n ds a f e t yi n s p e c t i o no fa b n o r m a ld e f o r m a t i o na r e a K e yw o r d s:t a i l i n g sd a m;UAV;p h o t o g r a mm e t r y;p o i n t c l o u dm o d e l;d e f o r m a t i o na n a l y s i s收稿日期:基金项目:湖南省应急管理厅安全
11、生产预防及应急专项项目(H u n a n )作者简介:廖文景(),男,高级工程师,硕士,安全工程专业,主要从事矿山安全技术咨询与研究工作.E m a i l:q q c o m.尾矿库是由筑坝拦截谷口或围地构成的,用于储存选矿厂废弃尾砂的重要生产设施 .尾矿库作为矿山最重要最危险的部分,是矿山事故发生的高频区域,它的安全对矿山生产具有至关重要的作用,因此需要实时监测尾矿库坝体变形情况,找出变第期廖文景等:无人机摄影测量在尾矿坝变形监测分析中的应用形的 位 置,对 于 尾 矿 库 安 全 稳 定 运 行 具 有 重 要意义 .目前,针对尾矿坝坝体变形监测分析,已有不少专家学者做了大量研究.李宗
12、伟等基于M o h r C o u l o m b破坏准则的弹塑性模型原理,结合某尾矿坝介绍了其建模分析过程,分析了其应力及变形计算结果的合理性.赵小稚采用灰色GM(,)模型对尾矿坝的变形进行预测,并结合某金矿尾矿坝变形监测实际数据进行预测,结果表明,模型用于尾矿坝变形趋势预测具有很好的适用性.邱宇等提出了基于V e r h u l s t灰色模型的坝体变形定量分析的思路,得出坝体的变形具有空间和时间效应.张一帆等 利用三维激光扫描技术设计一种新的尾矿坝变形监测方法,通过控制网测量、数据采集、数据预处理、点云建模、变形分析等操作完成对尾矿坝变形监测.张卫中等 针对某尾矿坝的沉降变形监测问题,采
13、用岩土数值模拟软件S I GMA/W,对其闭库后多年的沉降量进行了模拟计算分析.王乾等建立了基于P S O B P组合的尾矿坝变形预测模型,并验证了P S O B P组合模型在尾矿坝变形监测数据处理中具有较高的预测精度和拟合精度.L I U等 以尾矿坝变形为研究重点,以在线监测为行数据,以变形因素为主线,采用小波神经网络的方法,建立预测模型并进行分析.近年来,随着无人机摄影测量技术和影像三维重建技术的发展,为地理信息的获取带来了极大的便利.通过使用无人机摄影测量技术对同一地形进行多期数据的采集,然后进行三维重建生成带有坐标信息的点云模型,利用点云对比分析软件得出地理空间信息随着时间发生的变化
14、.本文以湖南省某尾矿坝为工程背景,采用无人机摄影测量技术,将该尾矿坝无人机摄影测量获取的两期原始影像数据进行三维重建,然后基于三维重建生成的点云模型对坝坡变形量进行分析.尾矿库概况该尾矿库为山谷型尾矿库,设计库容 万m.初期坝设计为透水坝,在构筑过程中,所选原料为堆积石块和选厂废石,其两侧设计构筑坡比均为,在内侧铺设土工布反滤层,坝底标高 m,坝高 m,坝顶宽度 m.堆积坝初期设计所采用的方法为向上游堆积尾砂筑坝,后因尾矿粒度太细,坝前难以形成干滩,因此在原有筑坝方案上进行改进,联合采用废石、土工布、旋流器沉砂筑坝,其中,堆积坝设计外坡比为,设计堆积高度 m,坝顶宽度 m,尾矿堆积至最终标高
15、m时,总高为 m.根 据 尾 矿 库 安 全 规 程(G B ),该尾矿库为三等库,但其总坝高接近二等库的高度,因此提高管理标准,将其按二等库管理.由于该尾矿坝安全管理范围大,山地环境特殊,传统测绘手段难以开展,尾矿坝变形位移数据仅通过单点或多点的监测方式获得,难以实现尾矿坝由点到面的全方位监测,并且监测精度不高.尾矿坝无人机影像获取及三维点云数据建立本研究选取湖南省某尾矿坝为研究区域,利用大疆精灵P h a n t o mR T K进行倾斜摄影测量,相机型号为:F C R;相机焦距:mm;像元大小:m.尾矿坝区域飞行相对高度设为 m,分辨率为 c m/p i x e l,相机倾角设置为 ,航
16、向重叠度和旁向重叠度均为,航线之间采用交叉飞行的方式进行航线设计,数据采集周期为 个月,共获取了 年月和 年月两期无人机航测影像集,然后利用D J IT e r r a软件分别对两期航拍影像集进行三维重建生成点云模型.如图所示为无人机摄影测量技术流程图.图摄影测量技术流程图F i g F l o wc h a r t o fp h o t o g r a m e t r i c t e c h n i q u e s在进行点云变形分析之前,对点云模型进行预处理,去除噪声点,并对其进行点云配准,预处理后两期点云模型如图所示.在两期点云模型中可以明显的看出,在第一次数据采集时尾矿坝部分区域进行了覆
17、绿工作,而在第二次采集数据时由于当地气候连续多天干旱,导致覆绿的植被大多已经干枯,难免会产生误差,本次研究为了尽可能降低覆绿对尾矿坝变形分析的影响,选取坝体中间包含道路和部分区域未覆绿的部分分析其变形量.有色金属(矿山部分)第 卷图 年月(左)和 年月(右)尾矿坝点云模型F i g T a i l i n g sd a mp o i n t c l o u dm o d e l s i nA u g u s t (l e f t)a n dJ u n e (r i g h t)尾矿坝变形监测分析 两期点云模型对比分析将两期摄影测量生成的点云模型导入到点云分析软件中,首先使用软件中的裁剪功能对其
18、分别进行裁剪,使得两期点云模型对比范围一致,图为两期点云模型裁剪后进行对比分析的范围.两期点云模型对比范围确定后,依据点云分析软件中提供的对齐功能,将第一期重建点云模型设置为参考模型,第二期点云模型设置为测试模型,然后建立特征点,使得两期模型对齐,对齐后不同位置的变形量通过色谱图来渲染,并且软件会自动计算出对比后的 D偏差、最大和最小偏差以及平均偏差和标准偏差,图为两期点云模型对比色谱图,从图中左下角可以看出两期尾矿坝模型对比区域中的最大正偏差为 c m,最大负偏差为 c m,平均正偏差为 c m,平均负偏差为 c m,标准偏差为 c m.特征点变形量对比分析点云分析软件提供了色谱注释功能,通
19、过对待分析的点添加色谱注释,能够在色谱图中直观地分析出每一点的色值分别在X、Y、Z轴方向以及 D(模型整体)的变形量,并对每个点设置允许偏差值,软件将依据设置的偏差允许范围,自动计算出当前点的状态,若当前点的偏差大于设定值,会显示失败,当偏差在设定的范围内时,则该点的状态显示为通过,因此利用该功能可以快速地识别出异常点的位置.图两期点云模型对比范围F i g C o m p a r i s o nr a n g eo f t w op h a s ep o i n t c l o u dm o d e l s图两期点云模型对比色谱图F i g C o m p a r i s o nc h r
20、o m a t o g r a mo f t w op h a s ep o i n tc l o u dm o d e l本次试验选取两期点云模型对比色谱图中具有代表性的个点进行注释,分别为A、A、A、A、A、A、A、A,如图点云模型色谱注释图所示.从图中可以直观地看出所注释的色谱分别在X、Y、第期廖文景等:无人机摄影测量在尾矿坝变形监测分析中的应用Z轴方向以 D(模型整体)的变形量,如图中A 点处D X的偏差为 c m、D Y的偏差为 c m、D Z的偏差为 c m,D偏差为 c m,可以看出该点在Z轴方向的偏差大,而在X和Y轴方向的偏差很小,再对比其它点在三个坐标轴方向的偏差也会发现其主
21、要偏差在D Z,也就是高程方向的偏差较大,而水平方向的偏差很小,结合实际情况中坝坡的台阶上通常有工作人员进行巡检以及覆绿植被等均会造成影响,因此在高程方向偏差较大,和实际情况相符.此外,从两期点云模型对比色谱图中也可以看出,色谱图明显成条带状分布,结合实际情况分析可知这是由于坝坡台阶覆绿造成的影响.图点云模型色谱注释图F i g C h r o m a t o g r a p h i ca n n o t a t i o no fp o i n t c l o u dm o d e l 计算结果综合分析点云模型对比完成后,根据对比结果,软件将自动生成点云模型对比分析报告,点云模型分析报告主要为
22、两期点云对比中点云数量在不同变形量区间的分布情况以及最大偏差和平均偏差,表为两期点云对比分析报告中的偏差分布表,图为点云对表点云模型对比偏差分布表T a b l eC o m p a r i s o nd e v i a t i o nd i s t r i b u t i o nt a b l eo fp o i n t c l o u dm o d e lM i n/c mM a x/c m点云数量/个占比/有色金属(矿山部分)第 卷图点云模型对比偏差分布图F i g C o n t r a s td e v i a t i o nd i s t r i b u t i o no fp o
23、 i n t c l o u dm o d e l比分析中的偏差分布柱状图.从点云模型对比偏差分布表和分布图中可以看出,本次点云模型中点云的总数为 个,其中偏差在 c m区间的点云数量占比最大,达到 ,绝对偏差大于 c m的点云数量总共占比为 ,可以看出偏差大于 c m的点云数量很少,在图两期点云模型对比色谱图找出变形量较大的点,结合原始点云模型进行对比,未发现异常,综上可以得出两期点云模型变形量基本在厘米级范围内,标准偏差为 c m,考虑到实际情况中坝坡表面通常有工作人员进行巡检以及坝坡覆绿植被等均会造成影响.该方法实现了对尾矿坝由点到面的全方位监测,为尾矿坝宏观变形监测及异常变形区域的快速
24、核查及安全巡检等提供了依据.讨论本文提出的基于无人机摄影测量的尾矿坝表面变形监测分析方法目前仍有需要继续改进与完善的方面,无人机摄影测量作为一种新型的测量方法,虽然有很多的优势,但和使用高精度便携仪器测量的坝面位移监测点数据相比,坝面位移监测点的精度达到毫米级,而摄影测量数据的精度为厘米级,摄影测量的精度不及坝面位移监测点的精度.造成无人机摄影测量精度低的主要原因是无人机在航飞过程中拍照时,由于网络等原因,难免在拍摄时会存在延迟的情况;此外,无人机在进行航拍作业时由于电机转动和风速的影响,也会影响每张影像的定位精度.随着无人机低空遥感技术的发展,将无人机搭载L i D A R、可见光、热红外等
25、传感器在同一研究区域进行多次飞行获取不同种类的数据,实现多类型遥感载荷集成,可以获取更高精度的数据 .将 G技术应用到无人机摄影测量领域,G网络由于具有超低时延的特性,数据传输可以实现毫秒级的传输时延,尤其对于在执行航拍任务时的无人机,在拍摄影像时,响应速度更快,地面控制人员可以准确地操控无人机,从而使每张影像的位置信息也更加准确.结论)通过对该尾矿坝两期点云模型对比分析可得:对比区域中的最大正偏差为 c m,最大负偏差为 c m,平均正偏差为 c m,平均负偏差为 c m,标准偏差为 c m.)通过对两期点云模型个特征点X、Y、Z轴方向以及 D(模型整体)的变形量进行分析可得,所有特征点在Z
26、轴方向的偏差大,而在X和Y轴方向的偏差很小,是由于坝坡台阶覆绿造成的影响.)相较于目前尾矿坝应用较多的单点或多点监测方式,基于无人机摄影测量的尾矿坝边坡表面变形监测分析方法能够快速准确地获取尾矿坝边坡的整体变形和细节变化,监测精度可达厘米级.)该方法实现了对尾矿坝由点到面的全方位监测,为尾矿坝宏观变形监测及异常变形区域的快速核查及安全巡检等提供了依据.参考文献田文旗,薛剑光尾矿库安全技术与管理M北京:煤炭工业出版社,T I ANW e n q i,XU E J i a n g u a n g S a f e t y t e c h n o l o g y a n dm a n a g e m
27、e n to ft a i l i n g s p o n d sM B e i j i n g:C o a lI n d u s t r yP r e s s,沈楼燕,李连通,张超,等我国尾矿库安全监测技术发展综述J有色金属工程,():S HE NL o u y a n,L IL i a n t o n g,Z HANGC h a o,e ta l O v e r v i e wo ft h ed e v e l o p m e n to fs a f e t y m o n i t o r i n gt e c h n o l o g yf o rt a i l i n g sp o n
28、d si n C h i n aJ N o n f e r r o u s M e t a l s E n g i n e e r i n g,():侯建斌尾矿库工程地质特性和稳定性分析 评 尾矿库工程地质 特 性 与 稳 定 性 研 究 J矿 业 研 究 与 开 发,():HOU J i a n b i n A n a l y s i s o f g e o l o g i c a l c h a r a c t e r i s t i c s a n ds t a b i l i t yo ft a i l i n g sp o n de n g i n e e r i n g:r e v
29、 i e wo f r e s e a r c ho ng e o l o g i c a l c h a r a c t e r i s t i c s a n d s t a b i l i t y o f t a i l i n g s p o n de n g i n e e r i n g J M i n i n g R e s e a r c ha n d D e v e l o p m e n t,():王乾,赵杰,王要善基于P S O B P组合模型的尾矿坝变形预测J地理空间信息,():,WAN G Q i a n,Z HA O J i e,WAN G Y a o s h a
30、n P r e d i c t i o n o ft a i l i n g s d a m d e f o r m a t i o n b a s e d o n P S O B P c o m b i n a t i o nm o d e lJ G e o s p a t i a l I n f o r m a t i o n,():,第期廖文景等:无人机摄影测量在尾矿坝变形监测分析中的应用袁利伟,李素敏尾矿库溃坝灾害风险分析理论与实践M北京:电子工业出版社,YUANL i w e i,L IS u m i n T h e o r ya n dp r a c t i c eo fr i s
31、 ka n a l y s i sf o rt a i l i n g s d a m f a i l u r e d i s a s t e r sM B e i j i n g:E l e c t r o n i cI n d u s t r yP r e s s,詹京晶,肖田元基于G P S数据的尾矿坝监测系统的研究与应用J中国安全生产科学技术,():Z HANJ i n g j i n g,X I AO T i a n y u a n R e s e a r c ha n da p p l i c a t i o no ft a i l i n g sd a m m o n i t o
32、 r i n gs y s t e mb a s e do nG P Sd a t aJ C h i n aS a f e t yP r o d u c t i o nS c i e n c ea n dT e c h n o l o g y,():李宗伟,许灵基于弹塑性理论的尾矿坝应力及变形研究J东北水利水电,():,L IZ o n g w e i,XU L i n g R e s e a r c ho ns t r e s sa n dd e f o r m a t i o no ft a i l i n g sd a m sb a s e do ne l a s t i c p l a
33、 s t i ct h e o r yJ N o r t h e a s tW a t e rR e s o u r c e sa n dH y d r o p o w e r,():,赵小稚基于GM()模型的尾矿坝变形趋势预测J山东理工大学学报(自然科学版),():Z HA O X i a o z h i P r e d i c t i o n o fd e f o r m a t i o nt r e n d o ft a i l i n g sd a m sb a s e do n GM()m o d e lJ J o u r n a lo fS h a n d o n gU n i v
34、 e r s i t yo f T e c h n o l o g y(N a t u r a lS c i e n c e E d i t i o n),():邱宇,项宏海,汪斌基于V e r h u l s t模型的尾矿坝变形监测数据分析J现代矿业,():Q I UY u,X I AN GH o n g h a i,WAN GB i n A n a l y s i s o fd e f o r m a t i o n m o n i t o r i n gd a t af o rt a i l i n g sd a m sb a s e do nt h eV e r h u l s tm
35、o d e lJ M o d e r nM i n i n g,():张一帆,张硕三维激光扫描技术在尾矿坝变形监测中的应用研究J世界有色金属,():Z HAN GY i f a n,Z HAN GS h u o A p p l i c a t i o nr e s e a r c ho ft h r e e d i m e n s i o n a l l a s e r s c a n n i n g t e c h n o l o g y i n d e f o r m a t i o nm o n i t o r i n go ft a i l i n g sd a m sJ W o r
36、l d N o n f e r r o u s M e t a l s,():张卫中,陈赞宇,康钦容,等某铜尾矿库子坝坝体变形与沉降量数值模拟J矿业研究与开发,():Z HAN G W e i z h o n g,C HE N Z a n y u,KAN G Q i n r o n g,e ta l N u m e r i c a l s i m u l a t i o no fd e f o r m a t i o na n ds e t t l e m e n to f t h es u bd a mb o d yo fac o p p e rt a i l i n g sp o n dJ
37、 M i n i n gR e s e a r c ha n dD e v e l o p m e n t,():L I UPP,HA OYJ T h er e s e a r c ho f t a i l i n g sp o n dd e f o r m a t i o nd i s p l a c e m e n tp r e d i c t i o nb a s e do nw a v e l e tn e u r a ln e t w o r kJA p p l i e dM e c h a n i c sa n dM a t e r i a l s,:杨超,杨鹏,吕文生,等基于无人
38、机摄影测量的尾矿坝边坡表面变形监测J中国安全生产科学技术,():YAN G C h a o,YAN G P e n g,L YU W e n s h e n g,e ta l S u r f a c ed e f o r m a t i o nm o n i t o r i n go ft a i l i n g sd a ms l o p e sb a s e do nd r o n ep h o t o g r a mm e t r yJ C h i n a S a f e t y P r o d u c t i o n S c i e n c e a n dT e c h n o l o
39、 g y,():闫石,林荡,李胜宣,等基于无人机倾斜摄影的三维重建技术在崩岸监测中的应用J成都工业学院学报,():YANS h i,L I N D a n g,L IS h e n g x u a n,e ta l T h ea p p l i c a t i o no ft h r e e d i m e n s i o n a lr e c o n s t r u c t i o nt e c h n o l o g yb a s e do n d r o n eo b l i q u ep h o t o g r a p h yi nl a n d s l i d e m o n i t
40、 o r i n gJ J o u r n a lo fC h e n g d uI n s t i t u t eo fT e c h n o l o g y,():毛亚纯,伏雨文,曹旺,等基于无人机点云数据的露天采场矿车提取方法J东北大学学报(自然科学版),():,MAOY a c h u n,F U Y u w e n,C AO W a n g,e ta l A m e t h o df o re x t r a c t i n go p e n p i tm i n i n gv e h i c l e sb a s e do nu n m a n n e da e r i a lv
41、e h i c l ep o i n t c l o u dd a t aJ J o u r n a l o fN o r t h e a s t e r nU n i v e r s i t y(N a t u r a lS c i e n c eE d i t i o n),():,侯方国,刘欣,任秀波无人机倾斜摄影与L i D A R融合监测技术J测绘通报,():HOUF a n g g u o,L I U X i n,R E N X i u b o D r o n et i l tp h o t o g r a p h ya n dL i D A Rf u s i o n m o n
42、i t o r i n gt e c h n o l o g yJS u r v e y i n ga n dM a p p i n gB u l l e t i n,():廉旭刚,韩雨,刘晓宇,等无人机低空遥感矿山地质灾害监测研究进展及发展趋势J金属矿山,():L I ANX u g a n g,HANY u,L I UX i a o y u,e t a l R e s e a r c hp r o g r e s sa n dd e v e l o p m e n t t r e n do fu n m a n n e da e r i a l v e h i c l e l o wa l
43、 t i t u d er e m o t es e n s i n gf o rm o n i t o r i n gg e o l o g i c a lh a z a r d s i nm i n e sJM e t a lM i n e,():樊冰,胡晓宁,马良,等无人机L i D A R及倾斜摄影测量技术在复杂河道工程测 绘中的融合研 究J水运工程,():,F A NB i n g,HUX i a o n i n g,MAL i a n g,e ta l T h e f u s i o nr e s e a r c ho fu n m a n n e da e r i a lv e
44、h i c l eL i D A Ra n do b l i q u ep h o t o g r a mm e t r yt e c h n o l o g y i nc o m p l e xr i v e re n g i n e e r i n gs u r v e y i n gJ W a t e rT r a n s p o r tE n g i n e e r i n g,():,宣润通,杨浩,刘成,等基于 G技术的无人机侦查系统设计J价值工程,():XUAN R u n t o n g,YAN G H a o,L I U C h e n g,e ta l D e s i g no fu n m a n n e da e r i a lv e h i c l er e c o n n a i s s a n c es y s t e m b a s e do n Gt e c h n o l o g yJ V a l u eE n g i n e e r i n g,():
©2010-2024 宁波自信网络信息技术有限公司 版权所有
客服电话:4008-655-100 投诉/维权电话:4009-655-100