首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 26 毫秒
1.
Using a new type of oscillation map, made from the radio spectra by the wavelet technique, we study the 18 April 2014 M7.3 flare (SOL2014-04-18T13:03:00L245C017). We find a quasi-periodic character of this flare with periods in the range 65?–?115 seconds. At the very beginning of this flare, in connection with the drifting pulsation structure (plasmoid ejection), we find that the 65?–?115 s oscillation phase slowly drifts towards lower frequencies, which indicates an upward propagating wave initiated at the start of the magnetic reconnection. Many periods (1?–?200 seconds) are found in the drifting pulsation structure, which documents multi-scale and multi-periodic processes. On this drifting structure, fiber bursts with a characteristic period of about one second are superimposed, whose frequency drift is similar to that of the drifting 65?–?115 s oscillation phase. We also checked periods found in this flare by the EUV Imaging Spectrometer (EIS)/Hinode and Interface Region Imaging Spectrograph (IRIS) observations. We recognize the type III bursts (electron beams) as proposed, but their time coincidence with the EIS and IRIS peaks is not very good. The reason probably is that the radio spectrum is a whole-disk record consisting of all bursts from any location, while the EIS and IRIS peaks are emitted only from locations of slits in the EIS and IRIS observations.  相似文献   

2.
In the paper by Kliem, Karlický, and Benz (Astron. Astrophys. 360, 715, 2000) it was suggested, that plasmoids formed during the bursty regime of solar flare reconnection can be “visualised” in the radio spectra as drifting pulsating structures via accelerated particles trapped inside the plasmoid. In the present paper we investigate this idea in detail. First, simple statistical analysis supporting this hypothesis is presented. Then, by using the 2.5-D MHD (including gravity) model solar flare reconnection in the inhomogeneous, stratified atmosphere is simulated and the formation and subsequent ejection of the plasmoid is demonstrated. The ejected plasmoid, which is considered to be a trap for accelerated electrons, is traced and its plasma parameters are computed. To estimate the associated plasma radio emission we need to know locations of accelerated electrons and corresponding plasma frequencies. General considerations predict that these electrons should be distributed mainly along the magnetic separatrix surfaces and this was confirmed by using a particle-in-cell simulation. Finally, under some simplifying assumptions the model dynamic radio spectrum is constructed. The relation between the global frequency drift and the plasmoid motion in the inhomogeneous ambient atmosphere is studied. The results are discussed with respect to the observed drifting pulsation structures and their possible utilisation for flare magnetic field diagnostics.  相似文献   

3.
Solar radio emission features a large number of fine structures demonstrating great variability in frequency and time. We present spatially resolved spectral radio observations of type IIIb bursts in the 30?–?80 MHz range made by the Low Frequency Array (LOFAR). The bursts show well-defined fine frequency structuring called “stria” bursts. The spatial characteristics of the stria sources are determined by the propagation effects of radio waves; their movement and expansion speeds are in the range of \((0.1\,\mbox{--}\,0.6)c\). Analysis of the dynamic spectra reveals that both the spectral bandwidth and the frequency drift rate of the striae increase with an increase of their central frequency. The striae bandwidths are in the range of \({\approx}\,(20\,\mbox{--}\,100)\) kHz and the striae drift rates vary from zero to \({\approx}\,0.3~\mbox{MHz}\,\mbox{s}^{-1}\). The observed spectral characteristics of the stria bursts are consistent with the model involving modulation of the type III burst emission mechanism by small-amplitude fluctuations of the plasma density along the electron beam path. We estimate that the relative amplitude of the density fluctuations is of \(\Delta n/n\sim10^{-3}\), their characteristic length scale is less than 1000 km, and the characteristic propagation speed is in the range of \(400\,\mbox{--}\,800~\mbox{km}\,\mbox{s}^{-1}\). These parameters indicate that the observed fine spectral structures could be produced by propagating magnetohydrodynamic waves.  相似文献   

4.
We report on the results of observations of a type IV burst made by the Ukrainian Radio interferometer of the Academy of Sciences (URAN-2) in the frequency range 22?–?33 MHz. The burst is associated with a coronal mass ejection (CME) initiated by a behind-the-limb active region (N05E151) and was also observed by the Nançay Decameter Array (NDA) radio telescope in the frequency band 30?–?60 MHz. The purpose of the article is the determination of the source of this type IV burst. After analysis of the observational data obtained with the URAN-2, the NDA, the Solar-Terrestrial Relations Observatory (STEREO) A and B spacecraft, and the Solar and Heliospheric Observatory (SOHO) spacecraft, we come to the conclusion that the source of the burst is the core of a behind-the-limb CME. We conclude that the radio emission can escape the center of the CME core at a frequency of 60 MHz and originates from the periphery of the core at a frequency of 30 MHz that is due to occultation by the solar corona at the corresponding frequencies. We find plasma densities in these regions assuming the plasma mechanism of radio emission. We show that the frequency drift of the start of the type IV burst is governed by an expansion of the CME core. The type III bursts that were observed against this type IV burst are shown to be generated by fast electrons propagating through the CME core plasma. A type II burst was registered at frequencies of 44?–?64 MHz and 3?–?16 MHz and was radiated by a shock with velocities of about \(1000~\mbox{km}\,\mbox{s}^{-1}\) and \(800~\mbox{km}\,\mbox{s}^{-1}\), respectively.  相似文献   

5.
As a coronal mass ejection (CME) passes, the flank and wake regions are typically strongly disturbed. Various instruments, including the Large Angle and Spectroscopic Coronagraph (LASCO), the Atmospheric Imaging Assembly (AIA), and the Coronal Multi-channel Polarimeter (CoMP), observed a CME close to the east limb on 26 October 2013. A hot (\({\approx}\,10~\mbox{MK}\)) rising blob was detected on the east limb, with an initial ejection flow speed of \({\approx}\, 330~\mbox{km}\,\mbox{s}^{-1}\). The magnetic structures on both sides and in the wake of the CME were strongly distorted, showing initiation of turbulent motions with Doppler-shift oscillations enhanced from \({\approx}\, \pm 3~\mbox{km}\,\mbox{s}^{-1}\) to \({\approx}\, \pm 15~\mbox{km}\,\mbox{s}^{-1}\) and effective thermal velocities from \({\approx}\,30~\mbox{km}\,\mbox{s}^{-1}\) to \({\approx}\,60~\mbox{km}\,\mbox{s}^{-1}\), according to the CoMP observations at the Fe?xiii line. The CoMP Doppler-shift maps suggest that the turbulence behaved differently at various heights; it showed clear wave-like torsional oscillations at lower altitudes, which are interpreted as the antiphase oscillation of an alternating red/blue Doppler shift across the strands at the flank. The turbulence seems to appear differently in the channels of different temperatures. Its turnover time was \({\approx}\,1000\) seconds for the Fe 171 Å channel, while it was \({\approx}\,500\) seconds for the Fe 193 Å channel. Mainly horizontal swaying rotations were observed in the Fe 171 Å channel, while more vertical vortices were seen in the Fe 193 Å channel. The differential-emission-measure profiles in the flank and wake regions have two components that evolve differently: the cool component decreased over time, evidently indicating a drop-out of cool materials due to ejection, while the hot component increased dramatically, probably because of the heating process, which is suspected to be a result of magnetic reconnection and turbulence dissipation. These results suggest a new turbulence-heating scenario of the solar corona and solar wind.  相似文献   

6.
A large set of coronal mass ejections (CMEs, 3463) has been selected to study their periodic oscillations in speed in the Solar and Heliospheric Observatory (SOHO) mission’s Large Angle and Spectrometric Coronagraph (LASCO) field of view. These events, reported in the SOHO/LASCO catalog in the period of time 1996?–?2004, were selected based on having at least 11 height–time measurements. This selection criterion allows us to construct at least ten-point speed–distance profiles and evaluate kinematic properties of CMEs with a reasonable accuracy. To identify quasi-periodic oscillations in the speed of the CMEs a sinusoidal function was fitted to speed–distance profiles and the speed–time profiles. Of the considered events 22 % revealed periodic velocity fluctuations. These speed oscillations have on average amplitude equal to \(87~\mbox{km}\,\mbox{s}^{-1}\) and period \(7.8 R _{\odot}/241~\mbox{min}\) (in distance/time). The study shows that speed oscillations are a common phenomenon associated with CME propagation implying that all the CMEs have a similar magnetic flux-rope structure. The nature of oscillations can be explained in terms of magnetohydrodynamic (MHD) waves excited during the eruption process. More accurate detection of these modes could, in the future, enable us to characterize magnetic structures in space (space seismology).  相似文献   

7.
We present the results of solar observations at 20 and 25 MHz with the Ukrainian T-shaped Radio telescope of the second modification (UTR-2) in the interferometric session from 27 May to 2 June 2014. In this case, the different baselines 225, 450, and 675 m between the sections of the east–west and north–south arms of UTR-2 were used. On 29 May 2014, strong sporadic radio emission consisting of Type III, Type II, and Type IV bursts was observed. On other days, there was no solar radio activity in the decameter range. We discuss the observation results of the quiet Sun. Fluxes and sizes of the Sun in east–west and north–south directions were measured. The average fluxes were 1050?–?1100 Jy and 1480?–?1570 Jy at 20 and 25 MHz, respectively. The angular sizes of the quiet Sun in equatorial and polar directions were \(55'\) and \(49'\) at 20 MHz and \(50'\) and \(42'\) at 25 MHz. The brightness temperatures of the radio emission were \({T_{\mathrm{b}}} = 5.1 \times{10^{5}}~\mbox{K}\) and \({T_{\mathrm{b}}} = 5.7 \times{10^{5}}~\mbox{K}\) at 20 and 25 MHz, respectively.  相似文献   

8.
The photometry of eclipse white-light (W-L) images showing a moving blob is interpreted for the first time together with observations from space with the PRoject for On Board Autonomy (PROBA-2) mission (ESA). An off-limb event seen with great details in W-L was analyzed with the SWAP imager (Sun Watcher using Active pixel system detector and image Processing) working in the EUV near 174 Å. It is an elongated plasma blob structure of 25 Mm diameter moving above the east limb with coronal loops under. Summed and co-aligned SWAP images are evaluated using a 20-h sequence, in addition to the 11 July, 2010 eclipse W-L images taken from several sites. The Atmospheric Imaging Assembly (AIA) instrument on board the Solar Dynamics Observatory (SDO) recorded the event suggesting a magnetic reconnection near a high neutral point; accordingly, we also call it a magnetic plasmoid. The measured proper motion of the blob shows a velocity up to \(12~\mbox{km}\,\mbox{s}^{-1}\). Electron densities of the isolated condensation (cloud or blob or plasmoid) are photometrically evaluated. The typical value is \(10^{8}~\mbox{cm}^{-3}\) at \(r=1.7~\mathrm{R}_{\odot}\), superposed on a background corona of \(10^{7}~\mbox{cm}^{-3}\) density. The mass of the cloud near its maximum brightness is found to be \(1.6\times10^{13}\) g, which is typically \(0.6\times10^{-4}\) of the overall mass of the corona. From the extrapolated magnetic field the cloud evolves inside a rather broad open region but decelerates, after reaching its maximum brightness. The influence of such small events for supplying material to the ubiquitous slow wind is noticed. A precise evaluation of the EUV photometric data, after accurately removing the stray light, suggests an interpretation of the weak 174 Å radiation of the cloud as due to resonance scattering in the Fe IX/X lines.  相似文献   

9.
An experiment carried out at the Brookhaven National Laboratory over a period of almost 8 years acquired 364 measurements of the beta-decay rates of a sample of \({}^{32}\mbox{Si}\) and, for comparison, of a sample of \({}^{36}\mbox{Cl}\). The experimenters reported finding “small periodic annual deviations of the data points from an exponential decay?…?of uncertain origin”. We find that power-spectrum and spectrogram analyses of these datasets show evidence not only of the annual oscillations, but also of transient oscillations with frequencies near 11 year?1 and 12.5 year?1. Similar analyses of 358 measurements of the solar neutrino flux acquired by the Super-Kamiokande neutrino observatory over a period of about 5 years yield evidence of an oscillation near 12.5 year?1 and another near 9.5 year?1. An oscillation near 12.5 year?1 is compatible with the influence of rotation of the radiative zone. We suggest that an oscillation near 9.5 year?1 may be indicative of rotation of the solar core, and that an oscillation near 11 year?1 may have its origin in a tachocline between the core and the radiative zone. Modulation of the solar neutrino flux may be attributed to an influence of the Sun’s internal magnetic field by the Resonant Spin Flavor Precession (RSFP) mechanism, suggesting that neutrinos and neutrino-induced beta decays can provide information about the deep solar interior.  相似文献   

10.
In this article, we present a multi-wavelength and multi-instrument investigation of a halo coronal mass ejection (CME) from active region NOAA 12371 on 21 June 2015 that led to a major geomagnetic storm of minimum \(\mathrm{Dst} = -204\) nT. The observations from the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory in the hot EUV channel of 94 Å confirm the CME to be associated with a coronal sigmoid that displayed an intense emission (\(T \sim6\) MK) from its core before the onset of the eruption. Multi-wavelength observations of the source active region suggest tether-cutting reconnection to be the primary triggering mechanism of the flux rope eruption. Interestingly, the flux rope eruption exhibited a two-phase evolution during which the “standard” large-scale flare reconnection process originated two composite M-class flares. The eruption of the flux rope is followed by the coronagraphic observation of a fast, halo CME with linear projected speed of 1366 km?s?1. The dynamic radio spectrum in the decameter-hectometer frequency range reveals multiple continuum-like enhancements in type II radio emission which imply the interaction of the CME with other preceding slow speed CMEs in the corona within \(\approx10\)?–?\(90~\mbox{R} _{\odot}\). The scenario of CME–CME interaction in the corona and interplanetary medium is further confirmed by the height–time plots of the CMEs occurring during 19?–?21 June. In situ measurements of solar wind magnetic field and plasma parameters at 1 AU exhibit two distinct magnetic clouds, separated by a magnetic hole. Synthesis of near-Sun observations, interplanetary radio emissions, and in situ measurements at 1 AU reveal complex processes of CME–CME interactions right from the source active region to the corona and interplanetary medium that have played a crucial role towards the large enhancement of the geoeffectiveness of the halo CME on 21 June 2015.  相似文献   

11.
We present a multi-wavelength study of a solar eruption event on 20 July 2004, comprising observations in H??, EUV, soft X-rays, and in radio waves with a wide frequency range. The analyzed data show both oscillatory patterns and shock wave signatures during the impulsive phase of the flare. At the same time, large-scale EUV loops located above the active region were observed to contract. Quasi-periodic pulsations with ???10 and ???15 s oscillation periods were detected both in microwave??C?millimeter waves and in decimeter??C?meter waves. Our calculations show that MHD oscillations in the large EUV loops ?C but not likely in the largest contracting loops ?C could have produced the observed periodicity in radio emission, by triggering periodic magnetic reconnection and accelerating particles. As the plasma emission in decimeter??C?meter waves traces the accelerated particle beams and the microwave emission shows a typical gyrosynchrotron flux spectrum (emission created by trapped electrons within the flare loop), we find that the particles responsible for the two different types of emission could have been accelerated in the same process. Radio imaging of the pulsed decimetric??C?metric emission and the shock-generated radio type II burst in the same wavelength range suggest a rather complex scenario for the emission processes and locations. The observed locations cannot be explained by the standard model of flare loops with an erupting plasmoid located above them, driving a shock wave at the CME front.  相似文献   

12.
We analyze the high-frequency drift radio structures observed by the spectrometer at Purple Mountain Observatory (PMO) over the frequency range of 4.5 – 7.5 GHz during the 18 March 2003 solar flare. The drifting structures take place before the soft X-ray maximum, almost at the maximum of hard X-ray flux at 25 – 50 keV. For the first time, the positive drift in this kind of radio structures is detected in such a high frequency range. Their global drifting rate is roughly estimated as 3.6 GHz s−1. They appear in four groups, lasting in total for less than 6 s, and have a broad bandwidth of more than 2 GHz but a smaller ratio of the bandwidth of the drifting structures to mean frequency than that of the lower frequency range. The lifetime of each individual burst in this event can be derived by using the high temporal resolution of the spectrometer at PMO and has an average value of 36.3 ms. Since the negative drifting structures observed in the 0.6 – 4.5 GHz frequency range were interpreted to be a radio signature of a plasmoid ejected upward (moving out of the Sun), the present observation may imply that it is possible for a plasmoid to move downward during a solar flare. However, for a confirmation of this suggestion direct radio imaging observation would be needed.  相似文献   

13.
We briefly discuss the observed features including the high flux density, short duration, narrow emission band, fast frequency drift, quasi-periodic oscillation and fast variation of polarized components, of 51 spike emission events observed at 2545/2645 MHz in the solar activity peak year, 1991 January–December, and carry out correlation analysis between these events and optical flares, magnetic field intensity and configuration of flare regions, and sunspot evolution types of active regions. In view of the fact that the observed and statistical characteristics of the spike emissions are very different from those of known types of solar radio burst and known solar radio components, we think that the spike emission in the peak years is probably a new type of radio burst excited by electron cyclotron maser instability under wave-particle resonance, or a new solar radio component.  相似文献   

14.
We analyze the time series of Ca?ii H-line obtained from Hinode/SOT on the solar limb. We follow three cases of upwardly propagating kink waves along a spicule and inverted Y-shaped structures at the cusp of it. The time-distance analysis shows that the axis of spicule undergos quasi-periodic transverse displacement at different heights from the photosphere. The mean period of transverse displacement is ~175 s and the mean amplitude is 1 arc?sec. The oscillation periods are increasing linearly with height which may be counted as the signature that the spicule is working as a low pass filter and allows only the low frequencies to propagate towards higher heights. The oscillations amplitude is increasing with height due to decrease in density. The phase speeds are increasing until some heights and then decreasing which may be related to the small scale reconnection at the spicule basis. We conclude that transversal displacement of spicules axis can be related to the propagation of kink waves along them. Moreover, we observe signatures of small-scale magnetic reconnection at the cusp of spicules which may excite kink waves.  相似文献   

15.
本文介绍一组短分米波(1.42GHz)太阳微耀斑的射电和光学辐射特征,它们包含53个叠加在连续辐射背景上的射电快速精细结构(FFS),即准周期快速脉冲链(称微耀斑),它们的形态相似,强度大约在150-200sfu范围内,其寿命(半功率宽)大多为15-50ms,有18个分离的双峰结构,该事件产生的7646活动区中出现两处新浮现的几个小黑子,呈现复杂极性,可能存在多重交叉小磁流环多次重联的复杂状况.本文定性地讨论了其产生机制,支持电流环爆炸性聚合模型的理论。  相似文献   

16.
Images of an east-limb flare on 3 November 2010 taken in the 131 Å channel of the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory provide a convincing example of a long current sheet below an erupting plasmoid, as predicted by the standard magnetic reconnection model of eruptive flares. However, the 171 Å and 193 Å channel images hint at an alternative scenario. These images reveal that large-scale waves with velocity greater than 1000 km?s?1 propagated alongside and ahead of the erupting plasmoid. Just south of the plasmoid, the waves coincided with type-II radio emission, and to the north, where the waves propagated along plume-like structures, there was increased decimetric emission. Initially, the cavity around the hot plasmoid expanded. Later, when the erupting plasmoid reached the height of an overlying arcade system, the plasmoid structure changed, and the lower parts of the cavity collapsed inwards. Hot loops appeared alongside and below the erupting plasmoid. We consider a scenario in which the fast waves and the type-II emission were a consequence of a flare blast wave, and the cavity collapse and the hot loops resulted from the break-out of the flux rope through an overlying coronal arcade.  相似文献   

17.
On 29 March 2014, NOAA Active Region (AR) 12017 produced an X1 flare that was simultaneously observed by an unprecedented number of observatories. We have investigated the pre-flare period of this flare from 14:00 UT until 19:00 UT using joint observations made by the Interface Region Imaging Spectrometer (IRIS) and the Hinode Extreme Ultraviolet Imaging Spectrometer (EIS). Spectral lines providing coverage of the solar atmosphere from the chromosphere to the corona were analysed to investigate pre-flare activity within the AR. The results of the investigation have revealed evidence of strongly blue-shifted plasma flows, with velocities up to \(200~\mbox{km}\,\mbox{s}^{-1}\), being observed 40 minutes prior to flaring. These flows are located along the filament present in the active region and are both spatially discrete and transient. In order to constrain the possible explanations for this activity, we undertake non-potential magnetic field modelling of the active region. This modelling indicates the existence of a weakly twisted flux rope along the polarity inversion line in the region where a filament and the strong pre-flare flows are observed. We then discuss how these observations relate to the current models of flare triggering. We conclude that the most likely drivers of the observed activity are internal reconnection in the flux rope, early onset of the flare reconnection, or tether-cutting reconnection along the filament.  相似文献   

18.
The oscillations of a gaseous polytrope with a magnetic field having both a toroidal and a poloidal component are examined using the second-order tensor virial equations on the assumption that the magnetic energy is small compared with the gravitational energy. The frequencies of oscillation of the transverse shear, the toroidal and the coupled pulsation modes are tabulated for polytropic indicesn=1, 1.5, 2, 3 and 3.5. It is found that the magnetic field decreases the frequency of oscillation of (i) the transverse shear mode and (ii) the mode which starts as a radial pulsation in the absence of a magnetic field while it increases the frequency of oscillation of (i) the toroidal mode and (ii) the Kelvin mode. In all cases the shift in frequency decreases with increasingn.  相似文献   

19.
We investigate the parameters of global solar p-mode oscillations, namely damping width \(\Gamma\), amplitude \(A\), mean squared velocity \(\langle v^{2}\rangle\), energy \(E\), and energy supply rate \(\mathrm{d}E/\mathrm{d}t\), derived from two solar cycles’ worth (1996?–?2018) of Global Oscillation Network Group (GONG) time series for harmonic degrees \(l=0\,\mbox{--}\,150\). We correct for the effect of fill factor, apparent solar radius, and spurious jumps in the mode amplitudes. We find that the amplitude of the activity-related changes of \(\Gamma\) and \(A\) depends on both frequency and harmonic degree of the modes, with the largest variations of \(\Gamma\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le3300~\upmu\mbox{Hz}\) and \(31\le l \le60\) with a minimum-to-maximum variation of \(26.6\pm0.3\%\) and of \(A\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le 3300~\upmu\mbox{Hz}\) and \(61\le l \le100\) with a minimum-to-maximum variation of \(27.4\pm0.4\%\). The level of correlation between the solar radio flux \(F_{10.7}\) and mode parameters also depends on mode frequency and harmonic degree. As a function of mode frequency, the mode amplitudes are found to follow an asymmetric Voigt profile with \(\nu_{\text{max}}=3073.59\pm0.18~\upmu\mbox{Hz}\). From the mode parameters, we calculate physical mode quantities and average them over specific mode frequency ranges. In this way, we find that the mean squared velocities \(\langle v^{2}\rangle\) and energies \(E\) of p modes are anticorrelated with the level of activity, varying by \(14.7\pm0.3\%\) and \(18.4\pm0.3\%\), respectively, and that the mode energy supply rates show no significant correlation with activity. With this study we expand previously published results on the temporal variation of solar p-mode parameters. Our results will be helpful to future studies of the excitation and damping of p modes, i.e., the interplay between convection, magnetic field, and resonant acoustic oscillations.  相似文献   

20.
This review summarizes new trends in studies of magnetic reconnection in solar flares. It is shown that plasmoids play a very important role in this primary flare process. Using the results of magnetohydrodynamic and particle-in-cell simulations, we describe how the plasmoids are formed, how they move and interact, and how a flare current sheet is fragmented into a cascade of plasmoids. Furthermore, it is shown that during the interactions of these plasmoids electrons are not only very efficiently accelerated and heated, but electromagnetic(radio) emission is also produced.We also describe possible mechanisms for the triggering of magnetic reconnection.The relevant X-ray and radio signatures of these processes(such as radio drifting pulsation structures, narrowband dm-spikes, and the loop-top and above-the-loop-top X-ray sources) are then described. It is shown that plasmoids can also be formed in kinked magnetic ropes. A mapping of X-points of the magnetic reconnection on the chromosphere(as e.g. a splitting of flare ribbons) is mentioned. Supporting EUV and white-light observations of plasmoids are added. The significance of all these processes for the fast magnetic reconnection and electron acceleration is outlined. Their role in fusion experiments is briefly mentioned.  相似文献   

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

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