首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 78 毫秒
1.
Seismic data from the MVT-SLN sesmic station located 7 km from the summit area of Mt Etna volcano, which has been operating steadily for the last two decades, have been analysed together with the volcanic activity during the same period. Cross-correlation techniques are used to investigate possible relationships between seismic and volcanic data and to evaluate the statistical significance of the results. A number of significant correlations have been identified, showing that there is an evident relation between seismic events and flank eruptions, and a less clear relation with summit activity, which appears more linked to tremor rather than to the low-frequency events. Particularly interesting are the low-frequency events whose rate of occurrence increases, starting from 17 to 108 days, prior to the onset of the flank activity and are candidates for a useful precursor. On the other hand, a tendency towards the increase in both the duration and the occurrence rate of transients in the volcanic tremor was observed before the onset of summit eruptions. As a result of this study different stages in the volcanic activity of Mt Etna, represented by changes in the characteristics of the recorded seismic phenomena, are identified.  相似文献   

2.
The volcano Hekla in south Iceland had its latest eruption in January–March 1991. The eruption was accompanied and followed by considerable seismic activity. This study examines the seismicity in the Hekla region (63°42′–64°18′N, 18°30′–20°12′W) during a period when the high activity related to the eruption had ceased, from July 1991 to October 1995. The aim is to define the level of the normal background seismicity of the area that can be compared to the eruption-related activity. The Hekla Volcano proper was generally aseismic during the study period. The most prominent earthquake cluster is in the neighbouring Torfajökull Volcano. The epicentres are concentrated in the western part of the caldera and west of it. The hypocentres are located at all depths from the surface down to 14?km, with highest activity at 5–12?km. Inside this cluster, in the northwest part of the caldera, is a spherical volume void of earthquakes, approximately 4?km in diameter and centred at 8?km depth. This is interpreted as a cooling magma body. Small, low-frequency events of volcanic origin were occasionally recorded at Torfajökull. This activity has mainly occurred in swarms and was most abundant during the first year of the study period, presumably reflecting some kind of connection to the 1991 Hekla eruption. Our study area also includes the easternmost section of the South Iceland seismic zone, a transform zone characterized by bookshelf faulting on transverse faults. Two lineaments of epicentres were identified, roughly corresponding to mapped faults of the South Iceland seismic zone. The hypocentres are relatively deep, mainly at 6–12?km, matching the general trend of hypocentral depth increasing toward the east. The seismicity is highest in the area of the mapped faults. However, the epicentres extend beyond them and indicate greater width of the South Iceland seismic zone, or 20–30?km rather than approximately 10?km as indicated by the length of the surface faults. The seismicity in the volcanic systems of Hekla and Vatnafjöll shows some characteristics of the South Iceland seismic zone. Epicentres are concentrated into two N–S lineaments, one of which coincides with the location of the 1987 Vatnafjöll earthquake (Mw=5.9), a strike-slip event on a N- to S-trending fault. The hypocentres of the Hekla–Vatnafjöll events are mainly at 8–13?km depth, which indicates a continuation of the depth trend of the earthquakes of the South Iceland seismic zone. The events located at Hekla proper and immediately north of it are all of low-frequency character, which can be held as an indication of volcanic origin. On the other hand, they show clear S arrivals at observing stations like normal high-frequency tectonic earthquakes.  相似文献   

3.
Various xenoliths have been found in lavas of the 1763 (“La Montagnola”), 2001, and 2002–03 eruptions at Mt. Etna whose petrographic evidence and mineral chemistry exclude a mantle origin and clearly point to a cognate nature. Consequently, cognate xenoliths might represent a proxy to infer the nature of the high-velocity body (HVB) imaged beneath the volcano by seismic tomography. Petrography allows us to group the cognate xenoliths as follows: i) gabbros with amphibole and amphibole-bearing mela-gabbros, ii) olivine-bearing leuco-gabbros, iii) leuco-gabbros with amphibole, and iv) Plg-rich leuco gabbros. Geobarometry estimates the crystallization pressure of the cognate xenoliths between 1.9 and 4.1 kbar. The bulk density of the cognate xenoliths varies from 2.6 to 3.0 g/cm3. P wave velocities (V P ), calculated in relation to xenolith density, range from 4.9 to 6.1 km/s. The integration of mineralogical, compositional, geobarometric data, and density-dependent V P with recent literature data on 3D V P seismic tomography enabled us to formulate the first hypothesis about the nature of the HVB which, in the depth range of 3–13 km b.s.l., is likely made of intrusive gabbroic rocks. These are believed to have formed at the “solidification front”, a marginal zone that encompasses a deep region (>5 km b.s.l.) of Mt. Etna’s plumbing system, within which magma crystallization takes place. The intrusive rocks were afterwards fragmented and transported as cognate xenoliths by the volatile-rich and fast-ascending magmas of the 1763 “La Montagnola”, 2001 and 2002–03 eruptions.  相似文献   

4.
Anomalous changes in the diffuse emission of carbon dioxide within the Masaya caldera have been observed before two seismic events that occurred at 10 and 30 km from the observation site. Their epicenters are located, respectively, south of Managua in Las Colinas (4.3 magnitude) and the Xiloa caldera (3.6 magnitude), in 2002 and 2003, recorded by the geochemical station located at El Comalito, Masaya volcano (Nicaragua). Anomalous increases were observed, which occurred around 50 and 8 days before the main seismic event that took place in Las Colinas, and 4 days before the seismic swarm at the Xiloa caldera, with a maximum CO2 efflux of 9.3 and 10.7 kg m?2 day?1, respectively. The anomalous CO2 efflux increases remained after filtering with multiple regression analysis was applied to the CO2 efflux time series, which indicated that atmospheric variables, during the first 4 months, explained 23 % CO2 variability, whereas, during the rest of the time series, CO2 efflux values are poorly controlled with only 6 %. The observed anomalies of the diffuse CO2 emission rate might be related to pressure changes within the volcanic–hydrothermal system and/or to geostructural changes in the crust due to stress/strain changes caused before and during the earthquakes’ formation, and seem not to be related to the activity of the main crater of Masaya volcano.  相似文献   

5.
 Approximately 20 km south of Mt. Etna craters, at the contact between volcanic and sedimentary formations, three mud volcanoes discharge CO2-rich gases and Na–Cl brines. The compositions of gas and liquid phases indicate that they are fed by a hydrothermal system for which temperatures of 100–150  °C were estimated by means of both gas and solute geothermometry. The hydrothermal system may be associated with CO2-rich groundwaters over a large area extending from the central part of Etna to the mud volcanoes. Numerous data on the He, CH4, CO2 composition of the gases of the three manifestations, sampled over the past 5 years, indicate clearly that variations are due to separation processes of a CO2-rich gas phase from the liquid. The effects of these processes have to be taken into account in the interpretation of the monitoring data collected for the geochemical surveillance of Etna volcano. Received: 4 September 1995 / Accepted: 14 February 1996  相似文献   

6.
In the Central Mediterranean region, the production of chemically diverse volcanic products (e.g., those from Mt. Etna and the Aeolian Islands archipelago) testifies to the complexity of the tectonic and geodynamic setting. Despite the large number of studies that have focused on this area, the relationships among volcanism, tectonics, magma ascent, and geodynamic processes remain poorly understood. We present a tomographic inversion of P-wave velocity using active and passive sources. Seismic signals were recorded using both temporary on-land and ocean bottom seismometers and data from a permanent local seismic network consisting of 267 seismic stations. Active seismic signals were generated using air gun shots mounted on the Spanish Oceanographic Vessel ‘Sarmiento de Gamboa’. Passive seismic sources were obtained from 452 local earthquakes recorded over a 4-month period. In total, 184,797 active P-phase and 11,802 passive P-phase first arrivals were inverted to provide three different velocity models. Our results include the first crustal seismic active tomography for the northern Sicily area, including the Peloritan–southern Calabria region and both the Mt. Etna and Aeolian volcanic environments. The tomographic images provide a detailed and complete regional seismotectonic framework and highlight a spatially heterogeneous tectonic regime, which is consistent with and extends the findings of previous models. One of our most significant results was a tomographic map extending to 14 km depth showing a discontinuity striking roughly NW–SE, extending from the Gulf of Patti to the Ionian Sea, south-east of Capo Taormina, corresponding to the Aeolian–Tindari–Letojanni fault system, a regional deformation belt. Moreover, for the first time, we observed a high-velocity anomaly located in the south-eastern sector of the Mt. Etna region, offshore of the Timpe area, which is compatible with the plumbing system of an ancient shield volcano located offshore of Mt. Etna.  相似文献   

7.
The July 17 – August 9, 2001 flank eruption of Mt. Etna was preceded and accompanied by remarkable changes in volcanic tremor. Based on the records of stations belonging to the permanent seismic network deployed on the volcano, we analyze amplitude and frequency content of the seismic signal. We find considerable changes in the volcanic tremor which mark the transition to different styles of eruptive activity, e.g., lava fountains, phreatomagmatic activity, Strombolian explosions. In particular, the frequency content of the signal decreases from 5 Hz to 3 Hz at our reference station ETF during episodes of lava fountains, and further decreases at about 2 Hz throughout phases of intense lava emission. The frequency content and the ratios of the signal amplitude allow us to distinguish three seismic sources, i.e., the peripheral dike which fed the eruption, the reservoir which fed the lava fountains, and the central conduit. Based on the analysis of the amplitude decay of the signal, we highlight the migration of the dike from a depth of ca. 5 km to about 1 km between July 10 and 12. After the onset of the effusive phase, the distribution of the amplitude decay at our stations can be interpreted as the overall result of sources located within the first half kilometer from the surface. Although on a qualitative basis, our findings shed some light on the complex feeding system of Mt. Etna, and integrate other volcanological and geophysical studies which tackle the problem of magma replenishment for the July–August, 2001 flank eruption. We conclude that volcanic tremor is fundamental in monitoring Mt. Etna, not only as a marker of the different sources which act within the volcano edifice, but also of the diverse styles of eruptive activity. An erratum to this article is available at .  相似文献   

8.
Volcanic tremor on Etna seems to have its origin within the main magma feeding system. On the basis of both spectral analyses at two permanent seismic stations and periodical measurements along the slopes of the volcano, two distinct sources are proposed. The former, characterized by low frequency contents (f<1.5 Hz), is located in a 2 km deep flat magma chamber, whereas the latter source seems to be linked to the upper part of the active vents.Turbulent motions in the magma-gas mixture, induced by escaping gases within the conduits, is one proposed cause of volcanic tremor on Etna (Seidl et al., 1981).From spectral analyses we propose approximate models of the feeding system of the main summit craters.Time variations of tremor energy were also investigated, and no regular patterns have been observed for the studied eruptions.More systematic information seems to be needed for a better knowledge of both the source model and location, and correlation between tremor features and volcanic activity.  相似文献   

9.
On January 30, 1974, an explosive eruption began on the western side of Etna. The activity evolved into two eruptive periods (January 30–February 17 and March 11–29). Two spatter cones (Mount De Fiore I and Mount De Fiore II) were formed at a height of about 1650 m a.s.l. and a distance of 6 km from the summit area. The effusive activity was very irregular with viscous lava flows of modest length.A seismic network of four stations was established around the upper part of the volcano on February 3. Moreover additional mobile stations were set up at several different sites in order to obtain more detailed informations on epicenter locations and spectral content of volcanic tremor.The volcanic activity is discussed in relation to the distribution of epicenters and the time-space distribution of the spectral characteristics of volcanic earthquakes and tremor. The characteristics of the seismic activity suggest that the flank eruption of Mount Etna was probably feed by a lateral branch of the main conduit yielding the activity at the Central Crater.  相似文献   

10.
A new set of three-dimensional velocity models beneath Mt. Etna volcano is derived in the present work. We have used P- and S-wave arrivals from local earthquakes recorded at permanent and temporary seismic networks installed since 1980. A set of 1249 earthquakes recorded at more than four seismic stations was selected for traveltime inversion. The velocity models obtained by using different data selection criteria and parametrization display similar basic features, showing a high P-wave velocity at shallow depth in the SE quadrant, in close connection with a high gravimetric Bouguer anomaly. This area shares a low Vp/Vs ratio. High P-wave velocities and high Vp/Vs ratios are obtained along the central conduits, suggesting the presence of dense, intrusive magmatic bodies extending to a depth of about 20 km. The central intrusive core is surrounded by lower P-wave velocities. The relocated earthquake hypocenters also display the presence of an outward dipping brittle region, away from the central conduits, surrounding a ductile zone spatially related to the high P-wave velocity anomalies located in proximity to the central craters.  相似文献   

11.
龙岗火山区尾波Q值的初步研究   总被引：2，自引：0，他引：2       下载免费PDF全文
基于sato单次散射模型,利用2007年5月至2010年2月龙岗火山测震台网记录到的发生在火山区及邻近地区内的41次Mt.≥1.6地震的波形资料,计算了龙岗火山区的尾波Q值,得出尾波Q值与频率的关系为Q(f)=(42.65±7.53)f(0.845 9+0.164 2),具有以低Q值高η值为特征的火山构造活跃地区的尾波...  相似文献   

12.
2002年夏季长白山天池火山区的地震活动研究   总被引：32，自引：8，他引：24       下载免费PDF全文
2002年6月以来,长白山天池火山区的地震活动明显增加. 本文利用2002年夏季布设在长白山天池火山区15套宽频带流动地震台站的记录资料,对天池火山区的地震活动进行了研究. 地震观测结果表明,2002年夏季长白山天池火山日平均地震发生频次超过30次. 地震主要位于长白山天池西南部和东北部两个区域,震源深度较浅,离地表的深度一般小于5km. 天池西南部和东北部的地震,b值存在较大的差异. 火山区地震记录的频谱分析和时频分析结果表明,这些地震主要为火山构造型地震. HSZ和DZD等台站地震记录中丰富的低频成分,可能与台站附近的局部介质或断层带有关. 我们认为2002年夏季频繁发生的地震和小震震群活动是由火山深部活动诱发的局部断裂活动引起.  相似文献   

13.
Volcanoes generate a broad range of seismo-volcanic and infrasonic signals, whose features and variations are often closely related to volcanic activity. The study of these signals is hence very useful in the monitoring and investigation of volcano dynamics. The analysis of seismo-volcanic and infrasonic signals requires specifically developed techniques due to their unique characteristics, which are generally quite distinct compared with tectonic and volcano-tectonic earthquakes. In this work, we describe analysis methods used to detect and locate seismo-volcanic and infrasonic signals at Mt. Etna. Volcanic tremor sources are located using a method based on spatial seismic amplitude distribution, assuming propagation in a homogeneous medium. The tremor source is found by calculating the goodness of the linear regression fit (R 2) of the log-linearized equation of the seismic amplitude decay with distance. The location method for long-period events is based on the joint computation of semblance and R 2 values, and the location method of very long-period events is based on the application of radial semblance. Infrasonic events and tremor are located by semblance–brightness- and semblance-based methods, respectively. The techniques described here can also be applied to other volcanoes and do not require particular network geometries (such as arrays) but rather simple sparse networks. Using the source locations of all the considered signals, we were able to reconstruct the shallow plumbing system (above sea level) during 2011.  相似文献   

14.
Two explosive eruptions occurred on 2 January 1996 at Karymsky Volcanic Center (KVC) in Kamchatka, Russia: the first, dacitic, from the central vent of Karymsky volcano, and the second, several hours later, from Karymskoye lake in the caldera of Akademia Nauk volcano. The main significance of the 1996 volcanic events in KVC was the phreatomagmatic eruption in Karymskoye lake, which was the first eruption in this lake in historical time, and was a basaltic eruption at the acidic volcanic center. The volcanic events were associated with the 1 January Ms 6.7 (Mw 7.1) earthquake that occurred at a distance of about 9–17 km southeast from the volcanoes just before the eruptions. We study the long-term (1972–1995) and short-term (1–2 January 1996) characteristics of crustal deformations and seismicity before the double eruptive event in KVC. The 1972–1995 crustal deformation was homogeneous and characterized by a gradual extension with a steady velocity. The seismic activity in 1972–1995 developed at the depth interval from 0 to 20 km below the Akademia Nauk volcano and spread to the southeast along a regional fault. The seismic activity in January 1996 began with a short sequence of very shallow microearthquakes (M ~0) beneath Karymsky volcano. Then seismic events sharply increased in magnitude (up to mb 4.9) and moved along the regional fault to the southeast, culminating in the Ms 6.7 earthquake. Its aftershocks were located to the southeast and northwest from the main shock, filling the space between the two active volcanoes and the ancient basaltic volcano of Zhupanovsky Vostryaki. The eruption in Karymskoye lake began during the aftershock sequence. We consider that the Ms 6.7 earthquake opened the passageway for basic magma located below Zhupanovsky Vostryaki volcano that fed the eruption in Karymskoye lake.  相似文献   

15.
长白山天池火山区的震群活动研究   总被引：4，自引：3，他引：4       下载免费PDF全文
2002和2003年夏季流动地震观测揭示，天池火山口附近存在大量的微震活动和一系列震群活动.地震定位结果表明地震主要发生在火山口附近，以震群形式发生的地震全部集中在天池火山口西南部，东北部地震密集区没有观测到震群活动.在夏季以外的其他季节，天池火山区只有一个固定地震台站（CBS）用于地震监测.利用CBS台不同时间的观测纪录，通过波形相关分析发现其他季节的主要震群活动仍然集中在天池西南部.震群的高精度相对定位揭示震源位置沿北西－南东向分布，倾向西南，倾角约80°. 2003年7月13日的震群发生期间，地震震源位置出现从深到浅的迁移现象，同时震源深度较大的地震在不同台站的地震波初动方向几乎全部向上，表明震源具有明显的膨胀分量.考虑到长白山天池火山2002年以来出现明显的地表形变、地球化学异常和谐频地震等现象，我们认为震群活动可能与5 km深度附近存在岩浆热液活动和岩浆增压有关.  相似文献   

16.
Mount Cameroon is an active volcano located in the Gulf of Guinea, west of Central Africa. After the March–April 1999 eruption on the SW flank, another eruption of the volcano occurred in 2000. It took place from three sites on the southwest flank and near the summit. The first eruptive site was located 500 m to the southwest of the summit, at 3900 m altitude. Activity on this site was mainly explosive with no lava flow. The second site was located between 3220 and 3470 m altitude. Lava was emitted along NNE–SSE fissures from this site and flew towards Buea, the main city of the area, stopping ~ 4 km from the first houses. The last site was located in the south western flank at 2750 m altitude. The lava ejected from an old cone near the first 1999 eruptive site was divided into two branches, for a total length of around 1 km. The location of active volcanic cones in 1999 and 2000 seems to be linked to the local tectonics. The pre-eruptive period was characterized by a seismic swarm which may be a precursor recorded in March 2000 by an analogue seismic station. The main shock was a magnitude 3.2 event, and was felt by the population in Ekona town located on the eastern flank. It had a Modified Mercalli intensity of III–IV. When the eruption started, a temporary network of short period 3-component seismic stations was set up around the volcano to improve the monitoring of seismic activity. The co-eruptive period from late May to September was characterized by sequences of earthquake swarms, volcanic tremor and a family of earthquakes having similar waveform and appearing regularly in August and early September. Some of the earthquakes were felt by the population in Buea and its environments. The largest seismic event recorded had a magnitude of 4. During the post-eruptive period from mid-September to December, seismicity returned to its background level of 1–3 earthquakes per 3 days. Hypocenter locations reveal a linear narrow structure under the summit zone which could represent the magmatic conduit of the volcano. The frequency/magnitude relationship revealed a b-value of 1.43 higher than those previously determined, but more representative of volcanic media. Seismic energy release was gradual after the 2000 eruption started.  相似文献   

17.
We describe the seismicity at Iwate volcano, northeastern Japan, during the volcanic unrest of 1998 with reference to a three-dimensional P and S wave velocity model from tomographic analysis. The abnormal seismic activity beneath Iwate volcano started under the caldera in February, 1998 and migrated westward in the period February to August, 1998. Previous geodetic modeling [Sato and Hamaguchi, Chikyu Monthly 21 (1999) 312–317] suggested the growth of a dike in the time of the seismic activity. Comparing the seismicity and dike extension with the tomographic images of the P and S wave velocity structure, we find that the trace of the growing dike coincides with the region of the high Vp and high Vp/Vs ratio beneath the volcano. The seismic and geodetic data are consistent with an intrusion of magma or other fluid under the caldera in 1998. Another pressure source causing the predominant crustal deformation at Iwate volcano was detected from geodetic data, which was located in the region with high Vp/Vs ratio under the western end of the volcano through the period from February to August. It is suggested that the activation of the point pressure source probably associated with the inflation of a hot fluid reservoir relate to a geothermal region adjacent to the western edge of the volcano.  相似文献   

18.
Systematic investigation of discrete gravity measurements has continued at Mount Etna since 1986. The network now covers an area of 400 km2 with about 70 stations 0.5–3 km apart. Mass redistributions occurring at depths ranging between about 8 km below sea level and a few hundred metres below the surface (magma level changes within the shallower parts of the feeding conduits) have been identified from these data. Conventional (discrete) microgravity monitoring on a network of stations furnishes only instantaneous states of the mass distribution at continuously active systems. In order to obtain information on the rate at which the volcanic processes (and thus mass transfers) occur, three stations for continuously recording gravity where installed on Mount Etna in 1998. A 16-month long sequence from one of the continuously running stations (PDN, located 2 km from the active northeast crater at the summit of Etna volcano) is presented. After removing the effects of Earth Tide and tilt, the correlation of the residual gravity sequence with simultaneous recordings of meteorological parameters acquired at the same station was analysed. Once the meteorological effects have also been removed, continuous gravity changes are within 10 μGal of gravity changes measured using conventional microgravity observations at sites very close to the continuous station. This example shows how discrete and continuous gravity observations can be used together at active volcanoes to get a fuller and more accurate picture of the spatial and temporal characteristics of volcanic processes.  相似文献   

19.
21 earthquakes recorded by a temporary seismic network in the Changbaishan Tianchi volcanic area in Northeast China operated during the summer of 2002 and 2003 were analyzed to estimate the S coda attenuation. The attenuation quality factor Qc was estimated using the single scattering attenuation model of Sato (1977) in the frequency band from 4 to 24 Hz. All the events studied in this paper occurred at depths from 2 to 6 km with ML of 1.4–2.8. The epicentral distances are less than 25 km. For all events which occurred near the Tianchi Lake (caldera), the Qc patterns obtained at the stations near the lake are similar, and the Qc values are relatively small. At the stations located about 15 km east of the Tianchi Lake, however, the average Qc is significantly higher. For an event which occurred 25km from the lake to the west, Qc patterns derived at the stations near the lake are quite similar to the above mentioned Qc for stations located in the east. Further study shows that Qc value in the north and central areas of the volcano is relatively lower than that in the surrounding area. Compared to other volcanic areas in the world, the average Qc of the Changbaishan Tianchi volcanic area is obviously lower. The deep seismic sounding and teleseismic receiver function studies indicated more than one lower velocity layer in the crust. The MT studies suggested the presence of high conductive bodies beneath the area. We interpret the strong attenuation of coda waves near the Changbaishan Tianchi volcano as being possibly related to high temperature medium caused by shallow magma chambers.  相似文献   

20.
This paper reports a study of the Tolud earthquake sequence; the sequence was a burst of shallow seismicity between November 28 and December 7, 2012; it accompanied the initial phase in the Tolbachik Fissure Eruption of 2012?2013. The largest earthquake (the Tolud earthquake of November 30, 2012, to be referred to as the Tolud Earthquake in what follows, with KS = 11.3, ML = 4.9, MC = 5.4, and MW = 4.8) is one of the five larger seismic events that have been recorded at depths shallower than 10 km beneath the entire Klyuchevskoi Volcanic Cluster in 1961?2015. It was found that the Tolud earthquake sequence was the foreshock–aftershock process of the Tolud Earthquake. This is one of the larger seismicity episodes ever to have occurred in the volcanic areas of Kamchatka. Data of the Kamchatka seismic stations were used to compute some parameters for the Tolud Earthquake and its largest (ML = 4.3) aftershock; the parameters include the source parameters and mechanisms, and the moment magnitudes, since no information on these is available at the world seismological data centers. The focal mechanisms for the Tolud Earthquake and for its aftershock are consistent with seismic ruptures at a tension fault in the rift zone. Instrumental data were used to estimate the intensity of shaking due to the Tolud Earthquake. We discuss the sequence of events that was a signature of the time-dependent seismic and volcanic activity that took place in the Tolbachik zone in late November 2012 and terminated in the Tolud burst of seismicity. Based on the current ideas of the tectonics and magma sources for the Tolbachik volcanic zone, we discuss possible causes of these earthquakes.  相似文献   

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

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