基于随机子空间法的海洋平台模态参数识别技术研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
海洋平台结构复杂、造价昂贵,且服役环境恶劣。为了避免重大恶性事故的发生,需要对海洋平台结构在服役期内进行定期/不定期的检测和安全评估。特别是服役中后期的海洋平台及超期服役的平台,实施有效的健康检测具有更加重要的意义。基于海洋平台振动测试响应的结构健康检测技术能够对结构的健康状况进行全面的检测,是一种全局检测方法,其中模态参数的识别是该技术的基础性和关键性环节。因此,准确掌握海洋平台模态参数显得尤为重要。
     早期的模态参数识别方法主要是基于频域响应的方法,但是由于频域方法存在频谱估计的局限性、频谱泄漏等问题,所以目前模态参数识别方法的发展趋势是使用基于结构时域振动响应的输入-输出或者输出的模态参数识别方法。其中复指数法(PRONY)、特征系统实现法(ERA)、随机子空间法(SSI)是在其他领域中模态识别效果较好的几种方法。
     本文通过对现有典型时域方法(复指数法(PRONY)、特征系统实现法(ERA)、随机子空间法(SSI))的性能比较,旨在发展一套有效的海洋平台模态参数精确识别技术。研究了有效获取结构的模态参数的方法,通过数据驱动随机子空间法(SSI/data)和协方差驱动随机子空间法(SSI/cov)差异性以及数据驱动随机子空间法(SSI/data)Hankel矩阵与噪声关系的深入研究完善了现有的模态参数识别技术,为进一步的结构健康检测和安全评价提供了理论基础和技术支持。主要研究工作如下:
     1.在不同激励方式(冲击激励、张拉激励以及白噪声激励)下分别对比复指数法(PRONY)、特征系统实现法(ERA)和随机子空间法(SSI/data和SSI/cov),证明了数据驱动随机子空间法(SSI/data)能够适用于不同激励下的模态参数识别,而且模态识别能力最强。
     2.传统观点认为数据驱动随机子空间法(SSI/data)和协方差驱动随机子空间法(SSI/cov)无论在理论上还是实际应用中都是一致的。通过蒙特卡洛试验、函数试验分别证明了协方差驱动随机子空间法(SSI/cov)识别模态数目少和模态参数估计精度不高等缺点,数据驱动随机子空间法(SSI/data)有更强的模态识别能力:消噪能力强、准确识别临近模态、模态估计结果精度高等。继而从理论上分析了两种方法存在差异的原因:数据有限性和QR分解算法,并通过数值试验进行了验证。
     3.深入研究了数据驱动随机子空间法(SSI/data),推导了Hankel矩阵的构建方式与数据驱动随机子空间法的消噪能力的关系式。提出了一套数据驱动随机子空间法(SSI/data) Hankel矩阵维数的评估方法,即归一化奇异值曲线、稳定图和有限元模态结果相结合的方式。分别利用数值算例和模型试验检验了该方法的有效性,研究表明:数据驱动随机子空间法(SSI/data)的Hankel矩阵应设置成非方阵。为今后数据驱动随机子空间法(SSI/data)的有效使用提供了参考和指导。
     4.基于海洋平台物理模型多激励(冲击激励、张拉激励以及白噪声激励)试验和渤海某平台现场实测等振动响应数据,系统地证明了数据驱动随机子空间法(SSI/data)在海洋平台模态参数识别中的适用性和优越性。同时验证了激励方式与海洋平台模态之间的关系。
To avoid major malignant accident, it is necessary to conduct periodic/aperiodic detection and safety assessment for offshore platforms with the action of various environmental forces, complex structure and high cost.So effective health detection for offshore platform, particularly those in old age or in overterm service, is more significant. The structural health monitoring technology based on the response data of the vibration test which is a global detection method is able to construct a comprehensive testing for offshore platforms, of which modal parameter identification is fundamental and critical. Therefore, the accurate modal parameter for offshore platform is particularly important.
     Early methods of modal identification were developed for the frequency domain. Because of limitations in the frequency resolution of spectral estimates and leakage errors in the estimates for the frequency domain methods, the new trend is to employ either input-output or output-only time-domain modal identification methods, of which the PRONY, Eigensystem Realization Algorithm(ERA) and SSI are more successful methods in other areas.
     With the aim at develop a effective identification technique for modal parameters of offshore platforms more exactly, the thesis studies on how to identify modal parameters more efficiently and improves existing modal parameters identification technologies by investigating the difference between SSI/data and SSI/cov and relation between Hankel matrix and noise. The main contributions are as follows:
     1. Based on the data associated with different excitations (impact loading, step relaxation, white noise loading), Prony, ERA and SSI/data are compared. It is proved that SSI/data has a better capacity of identifying the modal parameters of structures under different excitations (impact loading, step relaxation, white noise loading).
     2. The data-driven and covariance-driven stochastic subspace identification traditionally is thought to be consistent with each other theoretically and practically for modal identification. The paper investigates the reason that probably produce the difference,then confirms the reason of difference, the numerical study demonstrates that data-driven stochastic subspace identification method outperforms the covariance driven subspace identification method not only on accuracy of identification parameter but also on capacity of identifying weaker mode.
     3. The paper deduces theoretically the relation formula between noise and Hankel matrix of data-driven subspace identification method. And the paper also proposes a verification procedure to justify the noise can be eliminated properly by data-driven subspace identification method with selected Hankel matrix, the procedure includes SVD, stability diagram and finite element result (FE). Finally based on data associated with numerical study and jacket-type platform vibration test separately, we demonstrate systematically that data-driven stochastic subspace identification method with non-square Hankel matrix has better capacity of denoising and estimating the modal parameters with higher accuracy.
     4. With the data collected from a jacket-type physical model under different excitations(including impact loading, step-relaxation loading, white noise loading) and test data of realistic offshore jacket-type platform in Bohai sea in China, this paper demonstrates that SSI/data can be applied to modal identification of offshore platforms and be able to finish it better, and different strong modes can be excited by different excitation.
引文
[1]杨和振,李华军,黄维平.环境载荷激励下海洋平台的模态参数识别.第八届全国振动理论及应用学术会议,上海,2003.
    [2]欧进萍,段忠东,肖仪清等.海洋平台结构安全评定--理论、方法与应用.科学出版社,北京,2003.
    [3]张敏.固定式海洋结构模型转换与动力特性研究.博士论文,中国海洋大学,2010.
    [4]萧龙翔,曹树谦,张文德.振动结构模态分析——理论、实验与应用,天津大学出版社,天津,2001.
    [5]陈兆能,续绣忠,华宏星.基于环境激励模态参数辨识方法综述.振动与冲击,21(3):1-6,2002.
    [6]杨和振.环境激励下海洋平台结构模态参数识别与损伤诊断研究.博士论文,中国海洋大学,2004.
    [7]顾培英,邓昌.结构模态分析及其损伤诊断.东南大学出版社,南京,2008.
    [8]周晶.基于环境振动模态参数识别随机子空间方法与应用.硕士论文,兰州理工大学,2008.
    [9]Shuqing Wang, Huajun Li, and Xingxian Bao. Structural condition assessment and engineering treatment of an offshore platform with excessive vibration.In Proceedings of the International Offshore and Polar Engineering Conference}, pages 3813-3817, Lisbon, Portugal, July,2007.
    [10]李华军,杨和振.海洋平台结构环境激励的实验模态分析.振动与冲击,5:129-133,2005.
    [11]林循泓.振动模态参数识别及其应用,东南大学出版社,南京,1994.
    [12]冯志强.南黄海的油气勘查工作急需加强.第六届世界华人地质科学研讨会,2005.
    [13]N.Moller. R.Brincker, P.Andersen.Output-only modal testing of a car body subject to engine excitation.In Proc.of the 18th International Modal Analysis Conference,San Antonio,Texas,USA, page 786-792,2000.
    [14]陆秋海,李德葆.实验模态分析及其应用.科学出版社,北京,2001.
    [15]包兴先.基于模型定阶和信号消噪的海洋平台结构模态参数识别研究.博士论文,中国海洋大学,2010.
    [16]Andersen P, Brincker R, et al. Comparison of system identification methods using ambientbridge test data[C].17th IMAC,1032-1041,1999.
    [17]Bendat J S, Piersol A G. Engineering applications of correlation and spectral analysis[M], 2ndedition.New York:John Wiley-Sons,1993.
    [18]Brincker R, et al.Modal identification from ambient responses using frequency domain decomposition[C].18th IMAC,625-630,2000.
    [19]Mercer C.A. Clarkson, B.L. Use of cross correlation in studying the response of lightly damped structures to random forces. AIAA,3 (12):2287-2291,1965.
    [20]林裕裕.基于输出的海洋平台结构模态参数识别技术研究.硕士论文,中国海洋大学,2011.
    [21]Akaike H.Power spectrum estimation though autoregressive model fitting. Annals of the Institute of Statistical Mathematics,21:407-419,1969.
    [22]J.Lardies, N.Larbi.A new method for model order selection and modal parameter estimation in time domain. Journal of Sound and Vibration},245 (2):187-203,2001.
    [23]Ibrahim S.R.Application of random time domain analysis to dynamic flight measurements. Shock and Vibration Bullatin, Vol.49/2:165-170.,1979.
    [24]姜大正.环境激励下船舶结构模态分析实验与理论研究.硕士论文,大连理工大学,2009.
    [25]Cole H.A.on line failure detection and damping measurement of aerospace structure by random decrement signature. NASACR,22:(05),1973.
    [26]孟庆丰,何正嘉.随机减量技术中周期激励的影响及消除方法[J].振动与冲击,22(01):100-102,2003.
    [27]D.J. Ewins. Modal Testing:Theory, Practice and Applications.Research Studies Press, Baldock, Hertfordshire, England,2nd edition,2000.
    [28]G.H. James, T.G. Carner, and J.P. Lauffer.The natural excitation technique (next) for modal parameter extraction from operation structures. The International Journal of Analytical and Experimental ModalAnalysis,10 (4):260-277,1995.
    [29]Laufer J.P James G.H., Came T.G. The natural excitation technique (next) for modal parameter extraction from operating structures. Modal Analysis, VI 0 (14):260-277,1995.
    [30]James P. Lauffer George H. James 111, Thomas G. Carrie.The natural excitation technique (next) for modal parameter extraction from operating wind turbines.Technical report, Sandia National Laboratories,1993.
    [31]P.V.Oversehee, DeRoeck G, Pollet.T., Schueremans L. Stochastic subspace techniques applied to parameter identifieation of civil engineering struetures. In Proeeedings of the International Conferenee MVZ on New Advances in Modal Synthesis of Large Struetures, Non — Linear, Damped and Non-Deterministie Case, pp:151-162,Lyon,Franee,1995.
    [32]张笑华.结构环境振动模态参数识别随机子空间方法与应用.硕士论文,福州大学,2005.
    [33]姚志远.大型工程结构模态识别的理论和方法研究.博士论文,东南大学,2004.
    [34]P.R. Hammond, J.K. White. The analysis of the non-stationary signals using time-frequency methods. Journal of Sound and Vibration},190(3):419-447,1996.
    [35]禹丹江.土木工程结构模态参数识别.博士论文,福州大学,2006.
    [36]李萍.海洋平台结构动力特性参数识别的理论与实验研究.博士论文,中国海洋大学,2011.
    [37]丁志凯.基于人工神经网络的海洋平台损伤诊断技术研究.硕士论文,武汉理工大学,2007.
    [38]吴雅,杨叔子.时间序列分析的工程应用.华中理工大学出版社,1991.
    [39]朱东生.曲线桥地震反应研究.中国公路学报,15:43-48,2002.
    [40]史东峰.组合梁的物理非线性分析.南京建筑工程学院学报,2:12-19,1995.
    [41]郑敏.利用互相关函数进行环境激励下的模态分析.航空学报,6:535-537,2000.
    [42]李惠彬.青马悬索桥的时域模态识别.土木工程学报,34(5):12-19,2001.
    [43]欧进萍,何林,肖仪清.利用互相关函数进行环境激励下的模态分析.航空学报,6:535-537,2000.
    [44]何鹍.时域模态分析方法在桥梁动态无损检测中的应用研究.硕士论文,吉林大学,2008.
    [45]Huajun Li, Ping Li, Sau-Lon James Hu.Modal parameter estimation for jacket-type platforms using noise free-vibration data. Applied Ocean Research 37 (2012) 45-53,2011.
    [46]常军.随机子空间方法在桥梁模态参识别中的应用研究.博士论文,同济大学,2006.
    [47]F.R. Spitznogle, A.H. Quazi.Representation and analysis of time-limited signals using a complex exponential algorithm. Journal of the Acoustical Society of America, 47:1150-1155,1970.
    [48]华宏星,傅志方.模态分析理论与应用.上海交通大学出版社,上海,2002.
    [49]D.J.Ewins. Modal Testing:Theory, Practice and Applications[M],2nd. Baldoek Hertfordshire, England:Research Studies Press,2000.
    [50]N.Maia, J.Silva, J.He, N.Lieven, R-M.Lin, G.Skingle, W.TO, A.Urgueira. Theoretical and Experimental Modal Analysis[M].Taunton, Somerset, England:Research Studies Press,1997.
    [51]丁屹峰,古勇,程浩忠,孙毅斌,陆融.基于Prony算法的暂态电能质量扰动信号分析[J].上海交通大学学报,2005,39:36-38.
    [52]Johnson, M.A, Zarafonitis.I.P, Calligaris.M. Prony analysis and Power System stability-some recent theoretical and applications research[CI.Power Engineering Society Summer Meeting, 2000.IEEE, Volume.3:16-20, July 2000:1918--1923.
    [53]魏伟,赵书强,马燕峰.基于Prony算法的模糊电力系统稳定器设计[J].电力自动化设备,2005,25(13):54-56.
    [54]M.Amano, M.Watanabe, M.Banjo. Self-testing and Self-tuning of Power System Stabilizers Using PronyAnalysis[J].1998, IEEE.
    [55]刘红超,李兴源.基于PRONY辨识的交直流并联输电系统直流阻尼控制的研究.中国电机工程学报,2002,22(7):54-57.
    [56]刘红超,李兴源,邱晓燕,王建.基于Prony算法的HVDC输电系统模型辨识[J].四川大学学报(工程科学版),2004,36(2):86-89.
    [57]Sau-Lon James Hu,Xingxian Bao, Huajun Li. Model order determination and noise removal for modal parameter estimation. Mechanical systems and signal proeessing, 24(4):1605-1620,2010.
    [58]宫晨.提高海洋平台模态识别精度的研究.硕士论文,中国海洋大学,2008.
    [59]Juang J N, Pappa R S.An eigensystem realization algorithm for modal paremeter identification and modal reduction. Control, and Dynamics, pages 620-627,1985.
    [60]Jer-NanJuang, Riehard S.Pappa. Effects of Noise on Modal Parameters Identifiedby the Eigensystem Realization Algorithm[J]. Journal of Guidance,Control and Dynamies,1986, (3):194-303.
    [61]Darby J.Luscher,et al..parameter extraction of Z24 Bridge data[A].Proc.Of Modal IMAC 19[C],Kissimmee,Florida,USA,2001,836-841.
    [62]李惠彬,秦权,钱良忠.青马悬索桥的时域模态识别.土木工程学报,2001,34(5):52-56
    [63]王刚,姚谦峰.环境激励下框架结构体系的模态识别方法研究.工业建筑.2004,34(12):45-47.
    [64]Caicedo,juan M.Dyke,Shirley J.. Implementation of a SHM method on a numerical model of a cable-stayed bridge.Proceedings of the American Control Conference,v5,Proceedings of the 2004 American Control Conference (AAC),2004:4213-4218.
    [65]Dionysius M.Siringoringo and Yozo Fujino.Experimental study of laser Doppler vibrometer and ambient vibration for vibration-based damage detection Engineering Structures,2006,28 (13):1803-1815.
    [66]杨和振,李华军等.基于环境激励下海洋平台的试验模态分析.振动与冲击,2005,24(2):129-133.
    [67]纪晓东,钱稼如,徐龙河.模拟环境激励下结构模态参数识别试验研究.清华大学学报(自然科学版).2006,46(6):769-772.
    [68]叶锡钧,颜全胜,王卫锋,李健,基于多参考点稳定图的斜拉桥模态参数识别.华南理工大学学报(自然科学版).2011 Vol.39 No.9,41-50.
    [69]Ward Heylen, Stefan Lammens, Paul Sas. Modal Analysis Theory and Testing.Katholieke Universiteit Leuven, Leuven, Belgium,2nd edition,1997.
    [70]P.V.Oversehee, B.DeMoor. Subspace Algorithms for the Stoehastic Identifieation Problem. Proeeedings of the 30th IEEE Conferenee on Decision and Conirol, Brighion,649-660,1991.
    [71]P.V.Oversehee, B.DeMoor.Subspace algoritllms for the stochastic identification problem. Automatiea,29 (3):649-660,1993.
    [72]Lju K.Modal Parameter estimation using the state space method[J]. Journal of Sound and Vibration, 197 (4):387-402,1996.
    [73]L.Ljung.System identification-theory for the user. Prentice Hall, page 197-247,2002.
    [74]吴福生,顾培英,邓昌.结构模态分析及其损伤诊断.东南大学出版社,2007.
    [75]Bart Peeters, Guido DeRoeck. Stochastic system identification for operational modal analysis:A review. Journal of Dynamic Systems, Measurement, and Control, Vol.123:659-667,2001.
    [76]Edwin Reyndersa,Rik Pintelon,Guido DeRoeek.Uncertainty bounds onmodal Parameters obtained from stoehastic subspace identifieation. Mechanical systems and signal proeessing,22:948-969, 2008.
    [77]Virote Boonyapinyo, Tharach Janesupasaeree Data-driven stochastic subspace identification of flutter derivatives of bridge decks. Journal of Wind Engineering and Industrial Aerodynamics. 98(12). Page 784-799.
    [78]任伟新.环境振动系统识别方法的比较分析.福州大学学报(自然科学版),9(4):123-126,2001.
    [79]徐良,江见鲸,过静裙.随机子空间识别在悬索桥实验模态分析中的应用.工程力学,19(4):221-225,2002.
    [80]姚志远,汪凤泉,刘艳.工程结构模态的连续型随机子空间分解识别方法.东南大学学报(自然科学版),34(3):465-470,2004.
    [81]樊可清,倪一清,高赞明.改进随机子空间系统辨识方法及其桥梁状态监测中的应用.中国公路学报,30(3):342-346,2004.
    [82]禹丹江.土木工程结构模态参数识别.博士论文,福州大学,2006.
    [83]张笑华.结构环境振动模态参数识别随机子空间方法与应用.硕士论文,福州大学,2005.
    [84]徐士代.环境激励下工程结构模态参数识别.博士论文,东南大学,2006.
    [85]林友勤.基于随机子空间的结构损伤识别方法的研究.博士论文,福州大学博士学位论文,2007.
    [86]彭细荣,路新瀛,陈肇元.结构应变模态识别的随机子空间方法[J].振动与冲击,27(06):125-129,2008.
    [87]F.H.Miller.Dynamic properties of offshore structures with different founcation types by the random decrement technique[D].Master thesis,University of Maryland,1981.
    [88]Thoft-Christensen P, MUROT SU Y. Applicat ion of structural systems reliability theory [M]. Berlin: Spr inger-Verlag,1986.
    [89]要明伦,李大成,黄黎忠.箱腿深卧式海洋平台-土的动力相互作用[J].岩土工程学报.1993,15(2) :40-45.
    [90]刘方成.土—结构动力相互作用非线性分析及基于SSI效应的结构隔震研究[D].湖南大学,2008
    [91]金栋平.冰、风激励下海洋平台结构的非线性动力学和弛振分析[D].博士论文,天津大学,1995.
    [92]曹树谦.高维复杂转子系统非线性动力学的若干现代问题研究[D].天津大学,2003
    [93]L.Mangal,V.G.Idichandy,C.Canapathy.Structural monitoring of offshore plat-forms using impluse and relaxation response[J].Ocean Engineering.2001,28:689-705.
    [94]王俊荣.海洋平台结构物损伤检测与模型修正方法研究[D].中国海洋大学,2009.
    [95]陆建辉,李玉辉,付方.海洋平台动力响应及模态参数辨识[J]. Journal of Vibration and Shock 2006,25 (4):119-122
    [96]李玉辉.导管架平台结构优化设计及疲劳分析[D].中国海洋大学,2006
    [97]何晓宇,李宏男.地震与波浪联合作用下海洋平台动力特性分析[J].海洋工程.2007,25(3)
    [98]全伟.大跨桥梁多维多点地震反应分析研究[D].大连理工大学,2008
    [99]S.L.J.Hu,H.J.Li,S.Q.Wang.Cross-model cross-mode method for model updating[J].Mechanical Systems and Signal Processing.2007,21:1690-1703.
    [100]王熙堃.基于遗传算法的海洋平台水下构件的损伤识别研究[D].中国海洋大学,2008
    [1]D.J.Ewins.ModalTesting:Theory,PracticeandApplications[M],2nd.Baldoek.Hertfordshire,England:R esearch Studies Press,2000.
    [2]Juang J N, Pappa R S.An eigensystem realization algorithm for modal paremeter identification and modal reduction. Control, and Dynamics, pages 620-627,1985.
    [3]P.V.Oversehee, B.DeMoor.Subspace algoritllms for the stochastic identification problem. Automatiea,29(3):649-660,1993.
    [4]姚志远.大型工程结构模态识别的理论和方法研究.博士论文,东南大学,2004.
    [5]包兴先.基于模型定阶和信号消噪的海洋平台结构模态参数识别研究.博士论文,中国海洋大学,2010.
    [6]N.Maia,J.Silva,J.He,N.Lieven,R-M.Lin,G.Skingle,W.TO,A.Urgueira.TheoreticalandExperi mental Modal Analysis [M]. Taunton,Somerset,England:Researeh.Studies Press,1997.
    [7]林裕裕.基于输出的海洋平台结构模态参数识别技术研究.硕士论文,中国海洋大学,2011.
    [8]P. Overschee, B. DeMoor, Subspace Identification for Linear System:Theory-Implementation Applications, Kluwer Academic Publishers Boston/London/Dordrecht, (1996).
    [1]李萍.海洋平台结构动力特性参数识别的理论与实验研究.博士论文,中国海洋大学,2011.
    [2]包兴先.基于模型定阶和信号消噪的海洋平台结构模态参数识别研究.博士论文,中国海洋大学,2010.
    [3]The MathWorks, Control system toolbox user's guide, The MathWorks Inc., Natic, MA (2004).
    [4]B. Peeters, G. DeRoeck, Reference-based stochastic subspace identification for output-only modal analysis, Mechanical Systems and Signal Processing 13(6) (1999) 855-878.
    [5]ANSYS. Inc, ed.,2004. ANSYS Academic Research.Release 9.0.Help System. Coupled Field Analysis Guide 10.
    [1]B. Peeters, G. DeRoeck, Stochastic system identification for operational modal analysis:A review, Journal of Dynamic Systems, Measurement, and Control Vol.123 (2001) 659-667.
    [2]P. Overschee, B. DeMoor, Subspace Identification for Linear System:Theory-Implementation Applications, Kluwer Academic Publishers Boston/London/Dordrecht, (1996).
    [1]包兴先.基于模型定阶和信号消噪的海洋平台结构模态参数识别研究.博士论文,中国海洋大学,2010.
    [2]S.-L.J. Hu, P. Li, H. Vincent, H. Li, Modal parameter estimation for jacket-type platforms using free-vibration data,, Journal of Waterway, Port, Coastal, and Ocean Engineering, (2011) 1943-5460.
    [3]A. Cunha, E. Caetano, Experimental modal analysis of civil engineering structures, Sound and Vibration (2006) 12-20.
    [4]刘进明,应怀樵等时域模态分析方法的研究及软件研发振动与冲击2004年第23卷第4期Pg.123-126.
    [5]B. Peeters, G. DeRoeck, Reference-based stochastic subspace identification for output-only modal analysis, Mechanical Systems and Signal Processing 13(6) (1999) 855-878.
    [6]S.-L.J. Hu, X. Bao, H. Li, Model order determination and noise removal for modal parameter estimation, Mechanical Systems and Signal Processing 24 (2010) 1605-1620.
    [7]The Math Works, Control system toolbox user's guide, The Math Works Inc., Natic, MA (2004).
    [1]The MathWorks, Control system toolbox user's guide, The MathWorks Inc., Natic, MA (2004).
    [2]ANSYS. Inc, ed.,2004. ANSYS Academic Research.Release 9.0.Help System. Coupled Field Analysis Guide10.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700