首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
2017年9月8日发生了一次强磁暴,Kp指数最大值达到8.利用区域电离层格网模型(Regional Ionosphere Map,RIM)和区域ROTI(Rate of TEC Index)地图,分析了磁暴期间中国及其周边地区电离层TEC扰动特征和低纬地区电离层不规则体的产生与发展情况,同时利用不同纬度IGS(International GNSS Service)测站BJFS(39.6°N,115.9°E),JFNG(30.5°N,114.5°E)和HKWS(22.4°N,114.3°E)的GPS双频观测值,获取各测站的ROTI和DROT(Standard Deviation of Differential ROT)指数变化趋势.结果表明:此次磁暴发生期间电离层扰动先以正相扰动为主,主要发生在中低纬区域,dTEC(differential TEC)最大值达到14.9TECU,随后电离层正相扰动逐渐衰减,在低纬区域发生电离层负相扰动,dTEC最小值达到-7.2TECU;在12:30UT-13:30UT时段,中国南部低纬地区发生明显的电离层不规则体事件;相比BJFS和JFNG两个测站,位于低纬的HKWS测站的ROTI和DROT指数变化更为剧烈,这表明电离层不规则体结构存在纬度差异.   相似文献   

2.
基于加拿大地区高纬度电离层观测网的电离层闪烁观测数据,分析了2018年8月26日地磁暴事件引发的北半球高纬度地区电离层总电子含量(TEC)异常变化、TEC变化率指数(ROTI)及电离层相位闪烁的变化特征.结果表明:加拿大地区最大异常值约6 TECU,磁暴引发全球电离层TEC异常峰值高达20 TECU;加拿大地区电离层相位闪烁发生率最大增至12.6%,而磁静日期间约为1%;强电离层闪烁期间,电离层相位闪烁指数与ROTI之间具有较强的一致性.对GPS双频精密单点定位(Precise Point Positioning,PPP)结果进行分析发现:无闪烁期间定位误差随测站纬度的增高呈现出增大趋势,但均方根误差小于0.4m;闪烁发生期间各测站的定位误差均显著增大,水平和垂直方向定位误差分别增至约0.9m及1.7m.   相似文献   

3.
The Earth's ionosphere and especially its equatorial part is a highly dynamical medium. Geostationary satellites are known to be a powerful tool for ionospheric studies. Recent developments in BDS-GEO satellites allow such studies on the new level due to the best noise pattern in TEC estimations, which corresponds to those of GPS/GLONASS systems. Here we used BDS-GEO satellites to demonstrate their capability for studying equatorial ionosphere variability on different time scales. Analyzing data from the equatorial SIN1 IGS station we present seasonal variations in geostationary slant TEC for the periods of high (October 2013 - October 2014) and low (January 2017 - January 2018) solar activity, which show semi-annual periodicity with amplitudes about 10 TECU during solar maximum and about 5 TECU during the solar minimum. The 27-day variations are also prominent in geostationary slant TEC variations, which correlates quite well with the variations in solar extreme UV radiation. We found semi-annual pattern in small scale ionospheric disturbances evaluated based on geostationary ROTI index: maximal values correspond to spring and fall equinoxes and minimum values correspond to summer and winter solstices. The seasonal asymmetry in ROTI values was observed: spring equinox values were almost twice as higher than fall equinox ones. We also present results on the 2017 May 28–29 G3 geomagnetic storm, when ~30 TECU positive anomaly was recorded, minor and final major sudden stratospheric warmings in February and March 2016, with positive daytime TEC anomalies up to 15–20 TECU, as well as the 2017 September 6 X9.3 solar flare with 2 TECU/min TEC rate. Our results show the large potential of geostationary TEC estimations with BDS-GEO signals for continuous monitoring of space weather effects in low-latitude and equatorial ionosphere.  相似文献   

4.
Earthquake prediction stimulates the searches for a correlation between seismic activity and ionospheric anomalies. Contrary to common focuses on strong earthquakes, we report the ionospheric disturbances, 2 days before a moderate Ms = 4.7 Chongqing earthquake (29.4°N, 105.5°E, depth = 7.0 km, occurred at 21:21 LT, 10 September, 2010) with the data of ground-based ionosondes and IGS receivers. The data covering the period under the quiet geomagnetic conditions and a geomagnetic storm was analyzed with upper and lower bounds. It is found that there were significant enhancements of foF2 and total electron content (TEC) on the afternoon of 8 September, 2010, with a limited area close to the epicentre, which was different from the feature of ionospheric perturbations triggered by the geomagnetic storm on 15 September. Taking into account the heliogeomagnetical condition, we conclude that the observed ionospheric enhancements were very likely associated with the forthcoming moderate Chongqing earthquake, which implies that the relationship between the amplitudes of ionospheric disturbances and earthquakes is very complicated.  相似文献   

5.
The present study reveals the features of ionospheric parameters variations during the geomagnetic storm of September 7–8, 2017. In particular, parameters of vertical (foF2, foEs) and oblique ionospheric sounding (MOF, modes), absorption level, Total Electron Content (TEC) and particle fluxes at high altitudes were under analysis. The storm was characterized by two Dst-index mimima and can be considered as a sequence of two storms: first - with Dstmin?=??142?nT at 02 UT on September 8th and second - with Dstmin?=??122?nT and at 15 UT on September 8th. It was found that these two storms had different impacts on the ionosphere and HF propagation at mid- and high-latitudes of Northern Hemisphere. The signals of vertical and oblique ionospheric sounding were present in all ionograms before the first storm. Further, at the maximum of the first storm these signals were totally absorbed. Then, before the second storm and during its maximum the signals were detected again in the ionograms due to the low absorption. GOES satellite data showed the significant burst of electrons and protons only during the first storm and small particle fluxes - during the second storm. This feature was also confirmed with GPS data: TEC increased during the first storm and decreased during the second storm.  相似文献   

6.
全球电离层对2000年4月6-7日磁暴事件的响应   总被引:1,自引:0,他引:1  
利用分布于全球的电离层台站的测高仪观测数据,对扰动期间,foF2值与其宁静期间参考值进行比较,研究了2000年4月6—7日磁暴期间全球不同区域电离层的响应形态,并通过对比磁扰期间NmF2的变化与由MSISR90经验模式估算的中性大气浓度比(no/nN2)的变化,探讨了本次事件期间的电离层暴扰动机制.结果表明,在磁暴主相和恢复相早期,出现了全球性的电离层F2层负相暴效应.最大负相暴效应出现在磁暴恢复相早期,即电离层暴恢复相开始时间滞后于磁暴恢复相开始时间.在磁暴恢复相后期,一些台站出现正相扰动.研究结果表明,本次事件期间的电离层暴主要是由磁暴活动而诱发的热层暴环流引起的.  相似文献   

7.
The dual-frequency satellite-based measurements from Global Positioning System (GPS) may provide feasible ways of studying and potentially detecting of earthquake (EQ) related anomalies in the ionosphere. In this paper, GPS based Total Electron Content (TEC) data are studied for three major M?>?7.0 EQs in Nepal and Iran-Iraq border during 2015–2017 by implementing statistical procedures on temporal and spatial scale. Previous studies presented different time interval of pre-seismic ionospheric anomalies, however, this study showed that EQs ionospheric precursors may occur within 10?days. Furthermore, the ionospheric anomalies on the suspected day occurred during UT?=?08:00–12:00?h before the main shock. The Global Ionospheric Map TEC (GIM-TEC) data retrieved over the epicenter of M7.8 (Nepal EQ) showed a significant increase of 6 TECU on April 24, 2015 (one day before the main shock), which is recorded by the ground GPS station data of Islamabad (station lies within the EQ preparation zone). Furthermore, the spatial GIM-TEC result imply significant anomalies over the epicenter during the time interval between UT?=?08:00–12:00?h (LT?=?13:00–17:00). For M7.3 (Nepal EQ), the TEC anomalies were detected on May 10, 2015 (2?days before the event) in the temporal data. The spatial TEC data imply the huge clouds over the epicenter at about UT?=?08:00–12:00?h on May 10, 2015, that may be associated with this EQ in the quiet geomagnetic storm conditions. Similarly, temporal and spatial TEC showed anomaly on November 3, 2017, during UT?=?08:00–12:00 (9?days before the Iran-Iraq border EQ) after implementing the statistical method on it. Conversely, there exists a short-term but low magnitude TEC anomaly synchronized with a geomagnetic storm on November 7–8, 2017 (4 to 5?days prior to M7.3 Iran-Iraq border EQ). The diurnal and hourly GIM-TEC and VTEC data also imply the execution of ionospheric anomalies within 10?days prior to all events. All these positive anomalies in TEC may be due to the existence of a huge energy from the epicenter during the EQ preparation period.  相似文献   

8.
The responses of the ionospheric F region using GPS–TEC measurements during two moderate geomagnetic storms at equatorial, low-, and mid-latitude regions over the South American and African sectors in May 2010, during the ascending phase of solar cycle 24, are investigated. The first moderate geomagnetic storm studied reached a minimum Dst value of −64 nT at 1500 UT on 02 May 2010 and the second moderate geomagnetic storm reached a minimum Dst value of −85 nT at 1400 UT on 29 May 2010. In this paper, we present vertical total electron content (VTEC) and phase fluctuations (in TECU/min) from Global Positioning System (GPS) observations from the equatorial to mid-latitude regions in the South American and African sectors. Our results obtained during these two moderate geomagnetic storms from both sectors show significant positive ionospheric storms during daytime hours at the equatorial, low-, and mid-latitude regions during the main and recovery phases of the storms. The thermospheric wind circulation change towards the equator is a strong indicator that suggests an important mechanism is responsible for these positive phases at these regions. A pre-storm event that was observed in the African sector from low- to the mid-latitude regions on 01 May 2010 was absent in the South American sector. This study also showed that there was no generation or suppression of ionospheric irregularities by storm events. Therefore, knowledge about the suppression and generation of ionospheric irregularities during moderate geomagnetic storms is still unclear.  相似文献   

9.
This paper presents the response of the ionosphere during the intense geomagnetic storms of October 12–20, 2016 and May 26–31, 2017 which occurred during the declining phase of the solar cycle 24. Total Electron Content (TEC) from GPS measured at Indore, Calcutta and Siliguri having geomagnetic dips varying from 32.23°N, 32°N and 39.49°N respectively and at the International GNSS Service (IGS) stations at Lucknow (beyond anomaly crest), Hyderabad (between geomagnetic equator and northern crest of EIA) and Bangalore (near magnetic equator) in the Indian longitude zone have been used for the storms. Prominent peaks in diurnal maximum in excess of 20–45 TECU over the quiet time values were observed during the October 2016 storm at Lucknow, Indore, Hyderabad, Bangalore and 10–20 TECU for the May 2017 storm at Siliguri, Indore, Calcutta and Hyderabad. The GUVI images onboard TIMED spacecraft that measures the thermospheric O/N2 ratio, showed high values (O/N2 ratio of about 0.7) on October 16 when positive storm effects were observed compared to the other days during the storm period. The observed features have been explained in terms of the O/N2 ratio increase in the equatorial thermosphere, CIR-induced High Speed Solar Wind (HSSW) event for the October 2016 storm. The TEC enhancement has also been explained in terms of the Auroral Electrojet (AE), neutral wind values obtained from the Horizontal Wind Model (HWM14) and equatorial electrojet strength from magnetometer data for both October 2016 and May 2017 storms. These results are one of the first to be reported from the Indian longitude sector on influence of CME- and CIR-driven geomagnetic storms on TEC during the declining phase of solar cycle 24.  相似文献   

10.
利用行星际太阳风参数与太阳活动指数、地磁活动指数、电离层总电子含量格点化地图数据,首次基于一种能处理时间序列的深度学习递归神经网络(Recurrent Neural Network,RNN),建立提前24h的单站电离层TEC预报模型.对北京站(40°N,115°E)的预测结果显示,RNN对扰动电离层的预测误差低于反向传播神经网络(Back Propagation Neural Network,BPNN)0.49~1.46TECU,将太阳风参数加入预报因子模型后对电离层正暴预测准确率的提升可达16.8%.RNN对2001和2015年31个强电离层暴预报的均方根误差比BPNN低0.2TECU,将太阳风参数加入RNN模型可使31个事件的平均预报误差降低0.36~0.47TECU.研究结果表明深度递归神经网络比BPNN更适用于电离层TEC的短期预报,且在预报因子中加入太阳风数据对电离层正暴的预报效果有明显改善.   相似文献   

11.
磁暴期间电离层扰动的GPS台网观测分析   总被引:1,自引:3,他引:1  
给出了一种利用GPS台网观测获取TEC快速变化的计算方法,并将该方法用于东亚一澳大利亚扇区的GPS台网观测数据,分析了2000年7月14—18H大磁暴期间的电离层响应,揭示出电离层暴期间赤道异常峰的压缩和移动等特性.计算结果表明,在站点分布不均匀、原始观测数据不足且随时间跳变等多种不利因素的影响下,这种新的算法仍能保持很好的计算稳定性,并能快速地提取给定时空范围内的三维TEC短时变化的特征,适用于研究电离层暴等情况下引起的TEC扰动.  相似文献   

12.
The effects of the 15 May 2005 severe geomagnetic storm on the South African ionosphere are studied using ground-based and satellite observations. Ionospheric disturbances have less frequently been investigated over mid-latitude regions. Recently, a number of studies investigated their evolution and generation over these regions. This paper reports on the first investigation of travelling ionospheric disturbances (TIDs) over mid-latitude South Africa. Using global positioning system (GPS)-derived total electron content (TEC) variations from the South African network of dual frequency GPS receivers, we were able to examine the effects of the disturbance on the TEC. During this storm, two TEC enhancements were observed at low- and mid-latitudes: the first enhancement was observed between 30–45°S geomagnetic latitudes associated with equatorward neutral winds and the passage of a TID, while the second TEC enhancement is associated with a second TID. In addition, the F-region critical frequency (foF2) values observed at two ionosonde stations show response features that differ from those of the TEC during the disturbance period. The dissimilarity between the TEC and the foF2 suggests that two competing drivers may have existed, i.e., the westward electric field and equatorward neutral wind effects.  相似文献   

13.
This paper presents traveling ionospheric disturbances (TIDs) observations from GPS measurements over the South African region during the geomagnetically disturbed period of 29–31 October 2003. Two receiver arrays, which were along two distinct longitudinal sectors of about 18°-20° and 27°-28° were used in order to investigate the amplitude, periods and virtual propagation characteristics of the storm induced ionospheric disturbances. The study revealed a large sudden TEC increase on 28 October 2003, the day before the first of the two major storms studied here, that was recorded simultaneously by all the receivers used. This pre-storm enhancement was linked to an X-class solar flare, auroral/magnetospheric activities and vertical plasma drift, based on the behaviour of the geomagnetic storm and auroral indices as well as strong equatorial electrojet. Diurnal trends of the TEC and foF2 measurements revealed that the geomagnetic storm caused a negative ionospheric storm; these parameters were depleted between 29 and 31 October 2003. Large scale traveling ionospheric disturbances were observed on the days of the geomagnetic storms (29 and 31 October 2003), using line-of-sight vertical TEC (vTEC) measurements from individual satellites. Amplitude and dominant periods of these structures varied between 0.08–2.16 TECU, and 1.07–2.13 h respectively. The wave structures were observed to propagate towards the equator with velocities between 587.04 and 1635.09 m/s.  相似文献   

14.
Moderate geomagnetic storms occurred during January 22–25, 2012 period. The geomagnetic storms are characterized by different indices and parameters. The SYM-H value on January 22 increased abruptly to 67 nT at sudden storm commencement (SSC), followed by a sharp decrease to −87 nT. A second SSC on January 24 followed by a shock on January 25 was also observed. These SSCs before the main storms and the short recovery periods imply the geomagnetic storms are CME  -driven. The sudden jump of solar wind dynamic pressure and IMF BzBz are also consistent with occurrence of CMEs. This is also reflected in the change in total electron content (TEC) during the storm relative to quiet days globally. The response of the ionospheric to geomagnetic storms can also be detected from wave components that account for the majority of TEC variance during the period. The dominant coherent modes of TEC variability are diurnal and semidiurnal signals which account upto 83% and 30% of the total TEC variance over fairly exclusive ionospheric regions respectively. Comparison of TEC anomalies attributed to diurnal (DW1) and semidiurnal (SW2) tides, as well as stationary planetary waves (SPW1) at 12 UTC shows enhancement in the positive anomalies following the storm. Moreover, the impact of the geomagnetic storms are distinctly marked in the daily time series of amplitudes of DW1, SW2 and SPW1. The abrupt changes in amplitudes of DW1 (5 TECU) and SW2 (2 TECU) are observed within 20°S–20°N latitude band and along 20°N respectively while that of SPW1 is about 3 TECU. Coherent oscillation with a period of 2.4 days between interplanetary magnetic field and TEC was detected during the storm. This oscillation is also detected in the amplitudes of DW1 over EIA regions in both hemispheres. Eventhough upward coupling of quasi two day wave (QTDWs) of the same periodicity, known to have caused such oscillation, are detected in both ionosphere and upper stratosphere, this one can likely be attributed to the geomagnetic storm as it happens after the storm commencement. Moreover, further analysis has indicated that QTDWs in the ionosphere are strengthened as a result of coherent oscillation of interplanetary magnetic field with the same frequency as QTDWs. On the otherhand, occurrences of minor SSW and geomagnetic storms in quick succession complicated clear demarcation of attribution of the respective events to variability of QTDWs amplitudes over upper stratosphere.  相似文献   

15.
We investigated the ionospheric anomalies observed before the Tohoku earthquake, which occurred near the northeast coast of Honshu, Japan on 11 March, 2011. Based on data from a ground-based Global Positioning System (GPS) network on the Korean Peninsula, ionospheric anomalies were detected in the total electron content (TEC) during the daytime a few days before earthquake. Ionospheric TEC anomalies appeared on 5, 8 and 11 March. In particular, the ionospheric disturbances on 8 March evidenced a remarkable increase in TEC. The GPS TEC variation associated with the Tohoku earthquake was an increase of approximately 20 total electron content units (TECU), observed simultaneously in local and global TEC measurements. To investigate these pre-earthquake ionospheric anomalies, space weather conditions such as the solar activity index (F10.7) and geomagnetic activity indices (the Kp and Dst indices) were examined. We also created two-dimensional TEC maps to visual the spatial variations in the ionospheric anomalies preceding the earthquake.  相似文献   

16.
This paper presents results pertaining to the response of the mid-latitude ionosphere to strong geomagnetic storms that occurred from 31 March to 02 April 2001 and 07–09 September 2002. The results are based on (i) Global Positioning Systems (GPSs) derived total electron content (TEC) variations accompanying the storm, (ii) ionosonde measurements of the ionospheric electrodynamic response towards the storms and (iii) effect of storm induced travelling ionospheric disturbances (TIDs) on GPS derived TEC. Ionospheric data comprising of ionospheric TEC obtained from GPS measurements, ionograms, solar wind data obtained from Advanced Composition Explorer (ACE) and magnetic data from ground based magnetometers were used in this study. Storm induced features in vertical TEC (VTEC) have been obtained and compared with the mean VTEC of quiet days. The response of the mid-latitude ionosphere during the two storm periods examined may be characterised in terms of increased or decreased level of VTEC, wave-like structures in VTEC perturbation and sudden enhancement in hmF2 and h′F. The study reveals both positive and negative ionospheric storm effects on the ionosphere over South Africa during the two strong storm conditions. These ionospheric features have been mainly attributed to the travelling ionospheric disturbances (TIDs) as the driving mechanism for the irregularities causing the perturbations observed. TEC perturbations due to the irregularities encountered by the satellites were observed on satellites with pseudo random numbers (PRNs) 15, 17, 18 and 23 between 17:00 and 23:00 UT on 07 September 2002.  相似文献   

17.
The response of the ionospheric F-region in the equatorial and low latitude regions in the Brazilian sector during the super geomagnetic storm on 06–07 April 2000 has been studied in the present investigation. The geomagnetic storm reached a minimum Dst of −288 nT at 0100 UT on 07 April. In this paper, we present vertical total electron content (VTEC) and phase fluctuations (in TECU/min) from GPS observations obtained at Imperatriz (5.5°S, 47.5°W; IMPZ), Brasília (15.9°S, 47.9°W; BRAZ), Presidente Prudente (22.12°S, 51.4°W; UEPP), and Porto Alegre (30.1°S, 51.1°W; POAL) during the period 05–08 April. Also, several GPS-based TEC maps are presented from the global GPS network, showing widespread and drastic TEC changes during the different phases of the geomagnetic storm. In addition, ion density measurements on-board the satellite Defense Meteorological Satellite Program (DMSP) F15 orbiting at an altitude of 840 km and the first Republic of China satellite (ROCSAT-1) orbiting at an altitude of 600 km are presented. The observations indicate that one of the orbits of the DMSP satellite is fairly close to the 4 GPS stations and both the DMSP F15 ion-density plots and the phase fluctuations from GPS observations show no ionospheric irregularities in the Brazilian sector before 2358 UT on the night of 06–07 April 2000. During the fast decrease of Dst on 06 April, there is a prompt penetration of electric field of magnetospheric origin resulting in decrease of VTEC at IMPZ, an equatorial station and large increase in VTEC at POAL, a low latitude station. This resulted in strong phase fluctuations on the night of 06–07 April, up to POAL. During the daytime on 07 April during the recovery phase, the VTEC observations show positive ionospheric storm at all the GPS stations, from IMPZ to POAL, and the effect increasing from IMPZ to POAL. This is possibly linked to the equatorward directed meridional wind. During the daytime on 08 April (the recovery phase continues), the VTEC observations show very small negative ionospheric storm at IMPZ but the positive ionospheric storm effect is observed from BRAZ to POAL possibly linked to enhancement of the equatorial ionospheric anomaly.  相似文献   

18.
为研究中国陆态网区域电离层TEC在空间小尺度、高分辨率情况下的变化特性及适用精度范围,利用陆态网260个GNSS连续运行观测站数据,解算并生成2016-2017年731天陆态网区域电离层RIM格网,并进行精度验证.在同一RIM格网中,分别在经度和纬度方向上对间隔不同经纬度的TEC格网点作差分析.结果表明:陆态网区域内经度方向上TEC最大变化率和平均变化率分别为0.30TECU·(°)-1和0.11TECU·(°)-1;经度间隔1°时,TEC差值小于2TECU,且随着经度间隔的增大,其TEC差值也随之增大,并表现出一定的半年和周年变化规律;纬度方向上TEC最大变化率和平均变化率分别为1.7TECU·(°)-1和0.46TECU·(°)-1;陆态网区域内电离层TEC随纬度减小而增大,纬度间隔1°时,99.4%的TEC差值小于4TECU,且随着纬度间隔的增大,其TEC差值也随之增大,并表现出一定的半年和周年变化规律;间隔相同情况下,纬度方向上TEC的变化比经度方向上大.   相似文献   

19.
The F2-region reaction to geomagnetic storms usually called as an ionospheric storm is a rather complicated event. It consists of so called positive and negative phases, which have very complicated spatial and temporal behavior. The main morphological features of ionospheric storms and the main processes governing their behavior were understood at the end of the 1900s and described in a series of review papers. During the recent decade there were many publications dedicated to the problem of ionospheric storms. In this paper a concept of ionospheric storm morphology and physics formulated at the end of the 1990s is briefly summarized and the most interesting results obtained in the 2000s are described. It is shown that the main features of the studies of the previous decade were: the use of GPS TEC data for analyzing the ionospheric storm morphology, attraction of sophisticated theoretical models for studying the processes governing ionospheric behavior in disturbed conditions, and accent to analysis of ionospheric behavior during prominent events (very strong and great geomagnetic storms). Also a special attention was paid to the pre-storm enhancements in foF2 and TEC.  相似文献   

20.
利用两个中纬度台站GPS观测数据提取的GPS卫星硬件延迟,分析了不同太阳活动情况下估算的硬件延迟稳定性和统计特征,结合同期电离层观测数据,研究了电离层状态对硬件延迟估算结果的影响.研究结果表明,基于太阳活动高年(2001年)GPS观测数据估算的硬件延迟稳定性要低于太阳活动低年GPS观测数据的估算结果,利用2001年GPS数据估算的卫星硬件延迟年标准偏差(RMS)平均值约为1TECU,而2009年GPS数据估算的卫星硬件延迟年标准偏差平均值约为0.8TECU.通过对2001年和2009年北京地区电离层F2层最大电子密度(NmF2)变化性分析,结合GPS硬件延迟估算方法对电离层时空变化条件的要求,认为硬件延迟稳定性与太阳活动强度的联系是由不同太阳活动条件下电离层变化的强度差异引起的.   相似文献   

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

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