首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 588 毫秒
1.
CONT08 was a 15 days campaign of continuous Very Long Baseline Interferometry (VLBI) sessions during the second half of August 2008 carried out by the International VLBI Service for Geodesy and Astrometry (IVS). In this study, VLBI estimates of troposphere zenith total delays (ZTD) and gradients during CONT08 were compared with those derived from observations with the Global Positioning System (GPS), Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS), and water vapor radiometers (WVR) co-located with the VLBI radio telescopes. Similar geophysical models were used for the analysis of the space geodetic data, whereas the parameterization for the least-squares adjustment of the space geodetic techniques was optimized for each technique. In addition to space geodetic techniques and WVR, ZTD and gradients from numerical weather models (NWM) were used from the European Centre for Medium-Range Weather Forecasts (ECMWF) (all sites), the Japan Meteorological Agency (JMA) and Cloud Resolving Storm Simulator (CReSS) (Tsukuba), and the High Resolution Limited Area Model (HIRLAM) (European sites). Biases, standard deviations, and correlation coefficients were computed between the troposphere estimates of the various techniques for all eleven CONT08 co-located sites. ZTD from space geodetic techniques generally agree at the sub-centimetre level during CONT08, and??as expected??the best agreement is found for intra-technique comparisons: between the Vienna VLBI Software and the combined IVS solutions as well as between the Center for Orbit Determination (CODE) solution and an IGS PPP time series; both intra-technique comparisons are with standard deviations of about 3?C6?mm. The best inter space geodetic technique agreement of ZTD during CONT08 is found between the combined IVS and the IGS solutions with a mean standard deviation of about 6?mm over all sites, whereas the agreement with numerical weather models is between 6 and 20?mm. The standard deviations are generally larger at low latitude sites because of higher humidity, and the latter is also the reason why the standard deviations are larger at northern hemisphere stations during CONT08 in comparison to CONT02 which was observed in October 2002. The assessment of the troposphere gradients from the different techniques is not as clear because of different time intervals, different estimation properties, or different observables. However, the best inter-technique agreement is found between the IVS combined gradients and the GPS solutions with standard deviations between 0.2 and 0.7?mm.  相似文献   

2.
In October 2002, 15 continuous days of Very Long Baseline Interferometry (VLBI) data were observed in the Continuous VLBI 2002 (CONT02) campaign. All eight radio telescopes involved in CONT02 were co-located with at least one other space-geodetic technique, and three of them also with a Water Vapor Radiometer (WVR). The goal of this paper is to compare the tropospheric zenith delays observed during CONT02 by VLBI, Global Positioning System (GPS), Doppler Orbitography Radiopositioning Integrated by Satellite (DORIS) and WVR and to compare them also with operational pressure level data from the European Centre for Medium-Range Weather Forecasts (ECMWF). We show that the tropospheric zenith delays from VLBI and GPS are in good agreement at the 3–7 mm level. However, while only small biases can be found for most of the stations, at Kokee Park (Hawaii, USA) and Westford (Massachusetts, USA) the zenith delays derived by GPS are larger by more than 5 mm than those from VLBI. At three of the four DORIS stations, there is also a fairly good agreement with GPS and VLBI (about 10 mm), but at Kokee Park the agreement is only at about 30 mm standard deviation, probably due to the much older installation and type of DORIS equipment. This comparison also allows testing of different DORIS analysis strategies with respect to their real impact on the precision of the derived tropospheric parameters. Ground truth information about the zenith delays can also be obtained from the ECMWF numerical weather model and at three sites using WVR measurements, allowing for comparisons with results from the space-geodetic techniques. While there is a good agreement (with some problems mentioned above about DORIS) among the space-geodetic techniques, the comparison with WVR and ECMWF is at a lower accuracy level. The complete CONT02 data set is sufficient to derive a good estimate of the actual precision and accuracy of each geodetic technique for applications in meteorology.  相似文献   

3.
CONT campaigns are 2-week campaigns of continuous VLBI observations. The IERS working group on combination at the observation level uses these campaigns to study such combinations. In this work, combinations of DORIS, GPS, SLR, and VLBI technique measurements are studied during CONT08. We present different results concerning the use of common zenith tropospheric delay (ZTD) during the combination. We compare the ZTD obtained separately using each individual technique data processing, the combined ZTD, and the ZTD derived from a meteorological model. This resulted in a high level of consistency between each of these ZTD at a sub-centimeter level, a consistency which especially depends on the number of observations per estimated ZTD and the humidity level in the troposphere. We noted that GPS provides the main information about the combined ZTD, the other techniques providing complementary information when a lack of GPS observations occurs.  相似文献   

4.
Modeling path delays in the neutral atmosphere for the analysis of Very Long Baseline Interferometry (VLBI) observations has been improved significantly in recent years by the use of elevation-dependent mapping functions based on data from numerical weather models. In this paper, we present a fast way of extracting both, hydrostatic and wet, linear horizontal gradients for the troposphere from data of the European Centre for Medium-range Weather Forecasts (ECMWF) model, as it is realized at the Vienna University of Technology on a routine basis for all stations of the International GNSS (Global Navigation Satellite Systems) Service (IGS) and International VLBI Service for Geodesy and Astrometry (IVS) stations. This approach only uses information about the refractivity gradients at the site vertical, but no information from the line-of-sight. VLBI analysis of the CONT02 and CONT05 campaigns, as well as all IVS-R1 and IVS-R4 sessions in the first half of 2006, shows that fixing these a priori gradients improves the repeatability for 74% (40 out of 54) of the VLBI baseline lengths compared to fixing zero or constant a priori gradients, and improves the repeatability for the majority of baselines compared to estimating 24-h offsets for the gradients. Only if 6-h offsets are estimated, the baseline length repeatabilities significantly improve, no matter which a priori gradients are used.  相似文献   

5.
The combination of tropospheric parameters derived from different space-geodetic techniques has not been of large interest in geodesy so far. However, due to the high correlation between station coordinates and tropospheric parameters, the latter should not be neglected in combinations. This paper deals with the comparison and combination of tropospheric parameters derived from global positioning system (GPS) and very long baseline interferometry (VLBI) observations stemming from a 15-day campaign of continuous VLBI observations in 2002 (CONT02). The observation data of both techniques were processed homogeneously to avoid systematic differences between the solutions. We compared the tropospheric estimates of GPS and VLBI at eight co-location sites and found a very good agreement in the temporal behavior of the tropospheric zenith path delays (ZPD), reflected by correlation factors up to 0.98. Following this, a combination of the tropospheric parameters was performed. We demonstrate that the combination of tropospheric parameters leads to a stabilization of combined station networks. This becomes visible in the improvement of the repeatabilities of the station height components. Furthermore, the potential use of independent data from water vapor radiometers (WVRs) to validate space-technique-derived tropospheric parameters was investigated. Correlation coefficients of 0.95 or better were estimated between the tropospheric parameters of WVR and GPS or VLBI. Additionally, the utility of the tropospheric parameters for validation of local tie vectors was investigated. Both tropospheric zenith delays and tropospheric gradients were found to be very suitable to validate the height component and the horizontal components of the local tie, respectively.  相似文献   

6.
Precise positioning requires an accurate a priori troposphere model to enhance the solution quality. Several empirical models are available, but they may not properly characterize the state of troposphere, especially in severe weather conditions. Another possible solution is to use regional troposphere models based on real-time or near-real time measurements. In this study, we present the total refractivity and zenith total delay (ZTD) models based on a numerical weather prediction (NWP) model, Global Navigation Satellite System (GNSS) data and ground-based meteorological observations. We reconstruct the total refractivity profiles over the western part of Switzerland and the total refractivity profiles as well as ZTDs over Poland using the least-squares collocation software COMEDIE (Collocation of Meteorological Data for Interpretation and Estimation of Tropospheric Pathdelays) developed at ETH Zürich. In these two case studies, profiles of the total refractivity and ZTDs are calculated from different data sets. For Switzerland, the data set with the best agreement with the reference radiosonde (RS) measurements is the combination of ground-based meteorological observations and GNSS ZTDs. Introducing the horizontal gradients does not improve the vertical interpolation, and results in slightly larger biases and standard deviations. For Poland, the data set based on meteorological parameters from the NWP Weather Research and Forecasting (WRF) model and from a combination of the NWP model and GNSS ZTDs shows the best agreement with the reference RS data. In terms of ZTD, the combined NWP-GNSS observations and GNSS-only data set exhibit the best accuracy with an average bias (from all stations) of 3.7 mm and average standard deviations of 17.0 mm w.r.t. the reference GNSS stations.  相似文献   

7.
Continuous (CONT) VLBI campaigns have been carried out about every 3 years since 2002. The basic idea of these campaigns is to acquire state-of-the-art VLBI data over a continuous time period of about 2 weeks to demonstrate the highest accuracy of which the current VLBI system is capable. In addition, these campaigns support scientific studies such as investigations of high-resolution Earth rotation, reference frame stability, and daily to sub-daily site motions. The size of the CONT networks and the observing data rate have increased steadily since 1994. Performance of these networks based on reference frame scale precision and polar motion/LOD comparison with global navigation satellite system (GNSS) earth orientation parameters (EOP) has been substantially better than the weekly operational R1 and R4 series. The precisions of CONT EOP and scale have improved by more than a factor of two since 2002. Polar motion precision based on the WRMS difference between VLBI and GNSS for the most recent CONT campaigns is at the 30 \(\upmu \)as level, which is comparable to that of GNSS. The CONT campaigns are a natural precursor to the planned future VLBI observing networks, which are expected to observe continuously. We compare the performance of the most recent CONT campaigns in 2011 and 2014 with the expected performance of the future VLBI global observing system network using simulations. These simulations indicate that the expected future precision of scale and EOP will be at least 3 times better than the current CONT precision.  相似文献   

8.
Within the International Very Long Baseline Interferometry (VLBI) Service for Geodesy and Astrometry (IVS), long time-series of zenith wet and total troposphere delays have been combined at the level of parameter estimates. The data sets were submitted by eight IVS Analysis Centers (ACs) and cover January 1984 to December 2004. In this paper, the combination method is presented and the time-series submitted by the eight IVS ACs are compared with each other. The combined zenith delays are compared with time-series provided by the International Global Navigation Satellite System (GNSS) Service (IGS), and with zenith delays derived from the European Centre for Medium-Range Weather Forecasts (ECMWF). Before the combination, outliers are eliminated from the individual time-series using the robust BIBER (bounded influence by standardized residuals) estimator. For each station and AC, relative weight factors are obtained by variance component estimation. The mean bias of the IVS ACs’ time-series with respect to the IVS combined time-series is 0.89 mm and the mean root mean square is 7.67 mm. Small differences between stations and ACs can be found, which are due to the inhomogeneous analysis options, different parameterizations, and different treatment of missing in-situ pressure records. Compared to the IGS zenith total delays, the combined IVS series show small positive mean biases and different long-term trends. Zenith wet delays from the ECMWF are used to validate the IVS combined series. Inconsistencies, e.g., long-term inhomogeneity of the in-situ pressure data used for the determination of VLBI zenith delays, are identified.  相似文献   

9.
对陆态网223个全球导航卫星系统(Global Navigation Satellite System,GNSS)测站6 a实测对流层天顶延迟(zenith total delay,ZTD)的时空特性进行了分析,结果表明,各测站ZTD平均值随大地高指数递减,衰减因子与纬度近似线性关系,其时域变化呈现年周期和半年周期,周期、振幅、初相位与地域分布有关。综合采用周期函数及格网函数,建立了中国大陆区域ZTD经验模型SHAtrop。模型提供区域内分辨率为2.5°×2.0°的格网,用户使用时,先在相应格网内插得到对应参数,再利用三角函数得到椭球面ZTD,最后利用指数函数计算ZTD。实测ZTD的数据验证结果表明,SHAtrop的均方根误差(root mean square,RMS)为3.4 cm,优于常见经验模型。SHAtrop采用较多测站,对高程改正更精细;使用时只需要输入经纬度与时间,使用方便,能满足中国区域GNSS用户实时定位导航的ZTD改正需求。  相似文献   

10.
The Vienna Mapping Functions 1 (VMF1) as provided by the Institute of Geodesy and Geophysics (IGG) at the Vienna University of Technology are the most accurate mapping functions for the troposphere delays that are available globally and for the entire history of space geodetic observations. So far, the VMF1 coefficients have been released with a time delay of almost two days; however, many scientific applications require their availability in near real-time, e.g. the Ultra Rapid solutions of the International GNSS Service (IGS) or the analysis of the Intensive sessions of the International VLBI Service (IVS). Here we present coefficients of the VMF1 as well as the hydrostatic and wet zenith delays that have been determined from forecasting data of the European Centre for Medium-Range Weather Forecasts (ECMWF) and provided on global grids. The comparison with parameters derived from ECMWF analysis data shows that the agreement is at the 1 mm level in terms of station height, and that the differences are larger for the wet mapping functions than for the hydrostatic mapping functions and the hydrostatic zenith delays. These new products (VMF1-FC and hydrostatic zenith delays from forecast data) can be used in real-time analysis of geodetic data without significant loss of accuracy.  相似文献   

11.
The revitalized Russian GLONASS system provides new potential for real-time retrieval of zenith tropospheric delays (ZTD) and precipitable water vapor (PWV) in order to support time-critical meteorological applications such as nowcasting or severe weather event monitoring. In this study, we develop a method of real-time ZTD/PWV retrieval based on GLONASS and/or GPS observations. The performance of ZTD and PWV derived from GLONASS data using real-time precise point positioning (PPP) technique is carefully investigated and evaluated. The potential of combining GLONASS and GPS data for ZTD/PWV retrieving is assessed as well. The GLONASS and GPS observations of about half a year for 80 globally distributed stations from the IGS (International GNSS Service) network are processed. The results show that the real-time GLONASS ZTD series agree quite well with the GPS ZTD series in general: the RMS of ZTD differences is about 8 mm (about 1.2 mm in PWV). Furthermore, for an inter-technique validation, the real-time ZTD estimated from GLONASS-only, GPS-only, and the GPS/GLONASS combined solutions are compared with those derived from very long baseline interferometry (VLBI) at colocated GNSS/VLBI stations. The comparison shows that GLONASS can contribute to real-time meteorological applications, with almost the same accuracy as GPS. More accurate and reliable water vapor values, about 1.5–2.3 mm in PWV, can be achieved when GLONASS observations are combined with the GPS ones in the real-time PPP data processing. The comparison with radiosonde data further confirms the performance of GLONASS-derived real-time PWV and the benefit of adding GLONASS to stand-alone GPS processing.  相似文献   

12.
Troposphere parameters estimated from space-geodetic techniques, like the Global Positioning System (GPS) or Very Long Baseline Interferometry (VLBI), can be used to monitor the atmospheric water vapor content. Although the troposphere can only be monitored at discrete locations, the distribution of the instruments, at least the GPS antennas, can be assumed to be quasi-global. Critical in the data analysis are systematic effects within each single technique that significantly degrade the accuracy and especially the long-term stability of the zenith delay determination. In this paper, consistent time-series of troposphere zenith delays and gradients from homogeneously reprocessed GPS and VLBI solutions are compared for a time period of 11 years. The homogeneity of these completely reprocessed time-series is essential to avoid misinterpretations due to individual model changes. Co-located sites are used to investigate systematic effects and the long-term behavior of the two space-geodetic techniques. Both techniques show common signals in the troposphere parameters at a very high level of precision. The biases between the troposphere zenith delays are at the level of a few millimeters. On the other hand, long-term trends significantly differ for the two techniques, preventing climatological interpretations at present. Tests assume these differences to be due to mathematical artifacts such as different sampling rates and unmodeled semi-annual signals with varying amplitudes.  相似文献   

13.
从全球国际地球参考框架(International Terrestrial Reference Frame,ITRF)的建立、维护与发展,卫星测高、卫星重力等的发展及应用,全球卫星导航系统(Global Navigation Satellite System,GNSS)、卫星激光测距(Satellite Laser Ranging,SLR)、甚长基线干涉测量(very Long Baseline Interferometry,VLBI)、卫星多普勒定轨定位(Doppler Orbitography by Radiopositioning Integrated on Satellite,DORIS)的融合应用,海洋测绘和室内定位的发展等几个方面综述了大地测量学及卫星导航定位技术的最新进展,并提出中国2000国家大地坐标系与自主卫星导航系统的主要应用及发展目标。  相似文献   

14.
We present comparisons of 10-year-long time series of the atmospheric zenith wet delay (ZWD), estimated using the global positioning system (GPS), geodetic very long baseline interferometry (VLBI), a water vapour radiometer (WVR), radiosonde (RS) observations, and the reanalysis product of the European Centre for Medium-Range Weather Forecasts (ECMWF). To compare the data sets with each other, a Gaussian filter is applied. The results from 10 GPS–RS comparisons using sites in Sweden and Finland show that the full width at half maximum at which the standard deviation (SD) is a minimum increases with the distance between each pair. Comparisons between three co-located techniques (GPS, VLBI, and WVR) result in mean values of the ZWD differences at a level of a few millimetres and SD of less than 7?mm. The best agreement is seen in the GPS–VLBI comparison with a mean difference of ?3.4?mm and an SD of 5.1?mm over the 10-year period. With respect to the ZWD derived from other techniques, a positive bias of up to ~7?mm is obtained for the ECMWF reanalysis product. Performing the comparisons on a monthly basis, we find that the SD including RS or ECMWF varies with the season, between 3 and 15?mm. The monthly SD between GPS and WVR does not have a seasonal signature and varies from 3 to 7?mm.  相似文献   

15.
The very long baseline interferometry (VLBI) Intensive sessions are typically 1-h and single-baseline VLBI sessions, specifically designed to yield low-latency estimates of UT1-UTC. In this work, we investigate what accuracy is obtained from these sessions and how it can be improved. In particular, we study the modeling of the troposphere in the data analysis. The impact of including external information on the zenith wet delays (ZWD) and tropospheric gradients from GPS or numerical weather prediction models is studied. Additionally, we test estimating tropospheric gradients in the data analysis, which is normally not done. To evaluate the results, we compared the UT1-UTC values from the Intensives to those from simultaneous 24-h VLBI session. Furthermore, we calculated length of day (LOD) estimates using the UT1-UTC values from consecutive Intensives and compared these to the LOD estimated by GPS. We find that there is not much benefit in using external ZWD; however, including external information on the gradients improves the agreement with the reference data. If gradients are estimated in the data analysis, and appropriate constraints are applied, the WRMS difference w.r.t. UT1-UTC from 24-h sessions is reduced by 5% and the WRMS difference w.r.t. the LOD from GPS by up to 12%. The best agreement between Intensives and the reference time series is obtained when using both external gradients from GPS and additionally estimating gradients in the data analysis.  相似文献   

16.
We compare precipitable water vapor (PWV) time series measured by water vapor radiometers (WVRs) to PWV time series estimated using global positioning system (GPS) observations in a regional network of stations in western Europe. Inside this network, we focus on the baseline Brussels – Wettzell which presents the advantage to have the collocation of a GPS receiver and a WVR at both endpoints. The comparison between our GPS and WVR estimations of precipitable water vapor shows an agreement at the millimeter level. In addition, we show that the zenith total delay (ZTD) estimations computed with our GPS processing strategy agrees with the GPS estimations of ZTD done by the CODE analysis center at the millimeter level. Electronic Publication  相似文献   

17.
The empirical model GPT (Global Pressure and Temperature), which is based on spherical harmonics up to degree and order nine, provides pressure and temperature at any site in the vicinity of the Earth’s surface. It can be used for geodetic applications such as the determination of a priori hydrostatic zenith delays, reference pressure values for atmospheric loading, or thermal deformation of Very Long Baseline Interferometry (VLBI) radio telescopes. Input parameters of GPT are the station coordinates and the day of the year, thus also allowing one to model the annual variations of the parameters. As an improvement compared with previous models, it reproduces the large pressure anomaly over Antarctica, which can cause station height errors in the analysis of space-geodetic data of up to 1 cm if not considered properly in troposphere modelling. First tests at selected geodetic observing stations show that the pressure biases considerably decrease when using GPT instead of the very simple approaches applied to various Global Navigation Satellite Systems (GNSS) software packages so far. GPT also provides an appropriate model for the annual variability of global temperature. Electronic supplementary material The online version of this article (doi: contains supplementary material, which is available to authorized users.  相似文献   

18.
对流层延迟差异影响合成孔径雷达干涉测量技术(InSAR)形变测量精度;水汽的变化影响天气变化.对流层延迟与水汽具有较好的对应,因此有必要开展全球导航卫星系统(GNSS)对流层延迟的插值研究.以京津冀地区为例,针对GNSS对流层延迟,开展对流层延迟的空间插值研究.首先开展了GNSS对流层延迟与水汽的比较分析,两者存在显著正相关特性,相关性超过91.7%,论证了对流层延迟取代水汽的可行性.然后利用反距离权重法对京津冀地区2016年9月至2017年8月的12组GNSS测站对流层延迟进行空间插值,通过提取插值点对流层延迟与GNSS站点对流层延迟比较验证空间插值精度.全年数据平均偏差最大为1.12 cm,均方根误差最大为0.89 cm;未发生降水过程平均偏差最大为1.25 cm,均方根误差最大为0.82 cm;发生降水过程平均偏差最大为1.08 cm,均方根误差最大为1.38 cm.京津冀平原区域的GNSS对流层延迟空间插值结果精度满足气象等应用要求,可为气象预报和InSAR大气校正提供参考.   相似文献   

19.
High accurate global navigation satellite systems (GNSS) require to correct a signal delay caused by the troposphere. The delay can be estimated along with other unknowns or introduced from external models. We assess the impact of the recently developed augmentation tropospheric model on real-time kinematic precise point positioning (PPP). The model is based on numerical weather forecast and thus reflects the actual state of weather conditions. Using the G-Nut/Geb software, we processed GNSS and meteorological data collected during the experiment using a hot-air balloon flying up to an altitude of 2000 m. We studied the impacts of random walk noise setting of zenith total delay (ZTD) on estimated parameters and the mutual correlations, the use of external tropospheric corrections, the use of data from a single or dual GNSS constellation and the use of Kalman filter and backward smoothing processing methods. We observed a significant negative correlation of the estimated rover height and ZTD which depends on constraining ZTD estimates. Such correlation caused a degraded performance of both parameters when estimated simultaneously, in particular for a single GNSS constellation. The impact of ZTD constraining reached up to 50-cm differences in the rover height. Introducing external tropospheric corrections improved the PPP solution regarding: (1) shortened convergence, (2) better overall robustness, particularly, in case of degraded satellite geometry, (3) less adjusted parameters with lower correlations. The numerical weather model-driven PPP resulted in 9–12- and 5–6-cm uncertainties in the rover altitude using the Kalman filter and the backward smoothing, respectively. Compared to standard PPP, it indicates better performance by a factor of 1–2 depending on the availability of GNSS constellations, the troposphere constraining and the processing strategy.  相似文献   

20.
VLBI-derived troposphere parameters during CONT08   总被引:2,自引:2,他引:0  
Time-series of zenith wet and total troposphere delays as well as north and east gradients are compared, and zenith total delays (ZTD) are combined on the level of parameter estimates. Input data sets are provided by ten Analysis Centers (ACs) of the International VLBI Service for Geodesy and Astrometry (IVS) for the CONT08 campaign (12?C26 August 2008). The inconsistent usage of meteorological data and models, such as mapping functions, causes systematics among the ACs, and differing parameterizations and constraints add noise to the troposphere parameter estimates. The empirical standard deviation of ZTD among the ACs with regard to an unweighted mean is 4.6?mm. The ratio of the analysis noise to the observation noise assessed by the operator/software impact (OSI) model is about 2.5. These and other effects have to be accounted for to improve the intra-technique combination of VLBI-derived troposphere parameters. While the largest systematics caused by inconsistent usage of meteorological data can be avoided and the application of different mapping functions can be considered by applying empirical corrections, the noise has to be modeled in the stochastic model of intra-technique combination. The application of different stochastic models shows no significant effects on the combined parameters but results in different mean formal errors: the mean formal errors of the combined ZTD are 2.3?mm (unweighted), 4.4?mm (diagonal), 8.6?mm [variance component (VC) estimation], and 8.6?mm (operator/software impact, OSI). On the one hand, the OSI model, i.e. the inclusion of off-diagonal elements in the cofactor-matrix, considers the reapplication of observations yielding a factor of about two for mean formal errors as compared to the diagonal approach. On the other hand, the combination based on VC estimation shows large differences among the VCs and exhibits a comparable scaling of formal errors. Thus, for the combination of troposphere parameters a combination of the two extensions of the stochastic model is recommended.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号