首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 281 毫秒
1.
The maximum likelihood estimation method is applied to study the geographical distribution of earthquake hazard parameters and seismicity in 28 seismogenic source zones of NW Himalaya and the adjoining regions. For this purpose, we have prepared a reliable, homogeneous and complete earthquake catalogue during the period 1500–2010. The technique used here allows the data to contain either historical or instrumental era or even a combination of the both. In this study, the earthquake hazard parameters, which include maximum regional magnitude (M max), mean seismic activity rate (λ), the parameter b (or β?=?b/log e) of Gutenberg–Richter (G–R) frequency-magnitude relationship, the return periods of earthquakes with a certain threshold magnitude along with their probabilities of occurrences have been calculated using only instrumental earthquake data during the period 1900–2010. The uncertainties in magnitude have been also taken into consideration during the calculation of hazard parameters. The earthquake hazard in the whole NW Himalaya region has been calculated in 28 seismogenic source zones delineated on the basis of seismicity level, tectonics and focal mechanism. The annual probability of exceedance of earthquake (activity rate) of certain magnitude is also calculated for all seismogenic source zones. The obtained earthquake hazard parameters were geographically distributed in all 28 seismogenic source zones to analyze the spatial variation of localized seismicity parameters. It is observed that seismic hazard level is high in Quetta-Kirthar-Sulaiman region in Pakistan, Hindukush-Pamir Himalaya region and Uttarkashi-Chamoli region in Himalayan Frontal Thrust belt. The source zones that are expected to have maximum regional magnitude (M max) of more than 8.0 are Quetta, southern Pamir, Caucasus and Kashmir-Himanchal Pradesh which have experienced such magnitude of earthquakes in the past. It is observed that seismic hazard level varies spatially from one zone to another which suggests that the examined regions have high crustal heterogeneity and seismotectonic complexity.  相似文献   

2.
A straightforward Bayesian statistic is applied in five broad seismogenic source zones of the northwest frontier of the Himalayas to estimate the earthquake hazard parameters (maximum regional magnitude M max, β value of G–R relationship and seismic activity rate or intensity λ). For this purpose, a reliable earthquake catalogue which is homogeneous for M W ≥ 5.0 and complete during the period 1900 to 2010 is compiled. The Hindukush–Pamir Himalaya zone has been further divided into two seismic zones of shallow (h ≤ 70 km) and intermediate depth (h > 70 km) according to the variation of seismicity with depth in the subduction zone. The estimated earthquake hazard parameters by Bayesian approach are more stable and reliable with low standard deviations than other approaches, but the technique is more time consuming. In this study, quantiles of functions of distributions of true and apparent magnitudes for future time intervals of 5, 10, 20, 50 and 100 years are calculated with confidence limits for probability levels of 50, 70 and 90 % in all seismogenic source zones. The zones of estimated M max greater than 8.0 are related to the Sulaiman–Kirthar ranges, Hindukush–Pamir Himalaya and Himalayan Frontal Thrusts belt; suggesting more seismically hazardous regions in the examined area. The lowest value of M max (6.44) has been calculated in Northern-Pakistan and Hazara syntaxis zone which have estimated lowest activity rate 0.0023 events/day as compared to other zones. The Himalayan Frontal Thrusts belt exhibits higher earthquake magnitude (8.01) in next 100-years with 90 % probability level as compared to other zones, which reveals that this zone is more vulnerable to occurrence of a great earthquake. The obtained results in this study are directly useful for the probabilistic seismic hazard assessment in the examined region of Himalaya.  相似文献   

3.
We performed a tectonophysical analysis of earthquake frequency–size relationship types for large Central Asian earthquakes in the regions of dynamical influence due to major earthquake-generating faults based on data for the last 100 years. We identified four types of frequency–size curves, depending on the presence/absence of characteristic earthquakes and the presence or absence of a downward bend in the tail of the curve. This classification by the shape of the tail in frequency–size relationships correlates well with the values of the maximum observed magnitude. Thus, faults of the first type (there are characteristic earthquakes, but no downward bend) with Mmax ≥ 8.0 are classified as posing the highest seismic hazard; faults with characteristic earthquakes and a bend, and with Mmax = 7.5–7.9, are treated as rather hazardous; faults of the third type with Mmax = 7.1–7.5 are treated as posing potential hazard; and lastly, faults with a bend, without characteristic earthquakes, and with a typical magnitude Mmax ≤ 7.0, are classified as involving little hazard. The tail types in frequency–size curves are interpreted using the model of a nonlinear multiplicative cascade. The model can be used to treat different tail types as corresponding to the occurrence/nonoccurrence of nonlinear positive and negative feedback in earthquake rupture zones, with this feedback being responsible for the occurrence of earthquakes with different magnitudes. This interpretation and clustering of earthquake-generating faults by the behavior the tail of the relevant frequency–size plot shows raises the question about the physical mechanisms that underlie this behavior. We think that the occurrence of great earthquakes is related to a decrease in effective strength (viscosity) in the interblock space of faults at a scale appropriate to the rupture zone size.  相似文献   

4.
Probabilistic seismic hazard analysis (PSHA) has been carried out for Iraq. The earthquake catalogue used in the present study covers an area between latitude 29°–38.5° N and longitude 39°–50° E containing more than a thousand events for the period 1905–2000. The entire Iraq region has been divided into thirteen seismogenic sources based on their seismic characteristics, geological setting and tectonic framework. The completeness of the seismicity catalogue has been checked using the method proposed by Stepp (1972). The analysis of completeness shows that the earthquake catalogue is not complete below Ms=4.8 for all of Iraq and seismic source zones S1, S4, S5, and S8, while it varies for the other seismic zones. A statistical treatment of completeness of the data file was carried out in each of the magnitude classes. The Frequency Magnitude Distributions (FMD) for the study area including all seismic source zones were established and the minimum magnitude of complete reporting (Mc) were then estimated. For the entire Iraq the Mc was estimated to be about Ms=4.0 while S11 shows the lowest Mc to be about Ms=3.5 and the highest Mc of about Ms=4.2 was observed for S4. The earthquake activity parameters (activity rate , b value, maximum regional magnitude mmax) as well as the mean return period (R) with a certain lower magnitude mmin m along with their probability of occurrence have been determined for all thirteen seismic source zones of Iraq. The maximum regional magnitude mmax was estimated as 7.87 ± 0.86 for entire Iraq. The return period for magnitude 6.0 is largest for source zone S3 which is estimated to be 705 years while the smallest value is estimated as 9.9 years for all of Iraq.The large variation of the b parameter and the hazard level from zone to zone reflects crustal heterogeneity and the high seismotectonic complexity. The seismic hazard near the source boundaries is directly and strongly affected by the change in the delineation of these boundaries. The forces, through which the geological structure along the plate boundary in Eastern and Northeastern Iraq are evolved, are still active causing stress-strain accumulation, deformation and in turn producing higher probabilities of earthquake activity. Thus, relatively large destructive earthquakes are expected in this region. The study is intended to serve as a reference for more advanced approaches and to pave the path for the probabilistic assessment of seismic hazard in this region.  相似文献   

5.
The Iranian Plateau does not appear to be a single crustal block, but an assemblage of zones comprising the Alborz—Azerbaijan, Zagros, Kopeh—Dagh, Makran, and Central and East Iran. The Gumbel’s III asymptotic distribution method (GIII) and maximum magnitude expected by Kijko—Sellevoll method is applied in order to check the potentiality of the each seismogenic zone in the Iranian Plateau for the future occurrence of maximum magnitude (Mmax). For this purpose, a homogeneous and complete seismicity database of the instrumental period during 1900–2012 is used in 29 seismogenic zones of the examined region. The spatial mapping of hazard parameters (upper bound magnitude (ω), most probable earthquake magnitude in next 100 years (M100) and maximum magnitude expected by maximum magnitude estimated by Kijko—Sellevoll method (max MK ? Smax) reveals that Central and East Iran, Alborz and Azerbaijan, Kopeh—Dagh and SE Zagros are a dangerous place for the next occurrence of a large earthquake.  相似文献   

6.
In this paper, we calculated the seismic pattern of instrumental recorded small and moderate earthquakes near the epicenter of the 1303 Hongtong M=8 earthquake, Shanxi Province. According to the spatial distribution of small and moderate earthquakes, 6 seismic dense zones are delineated. Temporal distribution of M L≥2 earthquakes since 1970 in each seismic dense zone has been analyzed. Based on temporal distribution characteristics and historical earthquake activity, three types of seismicities are proposed. The relationship between seismic types and crustal medium is analyzed. The mechanism of three types is discussed. Finity of strong earthquake recurrence is proposed. Seismic hazard in mid-long term and diversity of earthquake disaster in Shanxi seismic belt are discussed.  相似文献   

7.
—?Earthquake hazard parameters are estimated by the application of the maximum likelihood method. The technique is based on a procedure which utilizes data of different quality, e.g., those in which the uncertainty in the assessment of the magnitudes is great and those in which the magnitudes are computed with great precision. In other words the data were extracted from both historical (incomplete) and recorded (complete) files. The historical part of the catalogue contains only the strongest events, whereas the complete part can be divided into several sub-catalogues; each one assumed to be complete above a specified magnitude threshold. Uncertainty in the determination of magnitudes has also been taken into account. The method allows us to estimate the earthquake hazard parameters which are the maximum regional magnitude, M max, the activity rate, λ, of the seismic events and the well known value β (b=β?log?e), which is the slope of the magnitude-frequency relationship. All these parameters are of physical significance. The mean return periods, RP, of earthquakes with a certain lower magnitude M?≥?m are also determined. The method is applied in the Island of Crete and the adjacent area, where catastrophic earthquakes are known from the historical era. The earthquake hazard of the whole area is divided in a cellular manner which allow the analysis of the localized hazard parameters and the representation of their regional variation. The seismic hazard analysis, which is expressed by: (a) The annual probability of exceedance of a specified value of magnitude and (b) the return periods (in years) that are expected for given magnitudes, for shallow events is finally performed for shallow events. This hazard analysis is useful for both theoretical and practical reasons and provides a tool for earthquake resistant design in both areas of low and high seismicity.  相似文献   

8.
—?The procedure developed by Kijko and Sellevoll (1989, 1992) and Kijko and Graham (1998, 1999) is used to estimate seismic hazard parameters in north Algeria. The area-specific seismic hazard parameters that were calculated consist of the b value of the Gutenberg–Richter frequency–magnitude relation, the activity rate λ(M) for events above the magnitude M, and the maximum regional magnitude M max. These parameters were calculated for each of the six seismogenic zones of north Algeria. The site-specific seismic hazard was calculated in terms of the maximum possible PGA at hypothetical engineering structures (HES), situated in each of the six seismogenic zones with coordinates corresponding with those of the six most industrial and populated cities in Algeria.  相似文献   

9.
Maximum earthquake size varies considerably amongst the subduction zones. This has been interpreted as a variation in the seismic coupling, which is presumably related to the mechanical conditions of the fault zone. The rupture process of a great earthquake indicates the distribution of strong (asperities) and weak regions of the fault. The rupture process of three great earthquakes (1963 Kurile Islands, MW = 8.5; 1965 Rat Islands, MW = 8.7; 1964 Alaska, MW = 9.2) are studied by using WWSSN stations in the core shadow zone. Diffraction around the core attenuates the P-wave amplitudes such that on-scale long-period P-waves are recorded. There are striking differences between the seismograms of the great earthquakes; the Alaskan earthquake has the largest amplitude and a very long-period nature, while the Kurile Islands earthquake appears to be a sequence of magnitude 7.5 events.The source time functions are deconvolved from the observed records. The Kurile Islands rupture process is characterized by the breaking of asperities with a length scale of 40–60 km, and for the Alaskan earthquake the dominant length scale in the epicentral region is 140–200 km. The variation of length scale and MW suggests that larger asperities cause larger earthquakes. The source time function of the 1979 Colombia earthquake (MW = 8.3) is also deconvolved. This earthquake is characterized by a single asperity of length scale 100–120 km, which is consistent with the above pattern, as the Colombia subduction zone was previously ruptured by a great (MW = 8.8) earthquake in 1906.The main result is that maximum earthquake size is related to the asperity distribution on the fault. The subduction zones with the largest earthquakes have very large asperities (e.g. the Alaskan earthquake), while the zones with the smaller great earthquakes (e.g. Kurile Islands) have smaller scattered asperities.  相似文献   

10.
-- The study addresses the evaluation of earthquake hazard parameters such as maximum regional magnitude (Mmax) and the slope of Gutenberg-Richter law # (where b=# log e) for the Hellenic Wadati-Benioff zone and the overriding lithospheric plate in the area of Crete and its surroundings. The seismicity of the area is divided in a cellular (1.0° 2 1.0°) manner allowing analysis of the localized earthquake hazard parameters and graphical representation of their spatial variation. Our approach incorporates the recently updated earthquake catalogue for Greece and the adjacent areas, the consideration of the morphology of the deep seismically active structures in the studied area and use of a probabilistic procedure for estimating the earthquake hazard parameters.¶One of the main inconsistencies in the earthquake hazard assessment is the estimation of the maximum magnitude and the related uncertaint y. The Bayesian approach, applied in the present, is a straightforward technique for evaluating the earthquake hazard parameters and is based on the following assumptions: Poissonian character of seismic events flow, a frequency-magnitude law of Gutenberg-Richter's type with cutoff maximal value for estimated parameter and a seismic catalogue, having a rather sizeable number of events (i.e., 50 events at least per cell). For five cells in which the number of events is less than 50, an effort is made to produce synthetic data. The re-assessed parameters obtained from the synthetic data show no significant difference and the real data (of the five cells) are finally taken into account although the estimated uncertainty is high.¶For four random cells we constructed hazard curves showing the probabilities that a certain magnitude M will be exceeded in one year and the return periods (in years) that are expected for a given magnitude. These are particularly useful for the mapping of earthquake hazard in regions of either low or high seismic activity, as is Crete and the adjacent area.¶The obtained results show that the W and E parts of both subducting and overriding plates differ in the spatial distribution of all the estimated earthquake hazard parameters. The Mmax distribution indicates strong coupling between the western portions of the interacting plates (Mmax > 6.3) to the south of 36°N. The smaller values of Mmax (Mmax < 6.3) estimated in the SE part of the studied area indicate weak coupling between the eastern portions of the subducting and overriding plates.¶Values of b > 1.0 are found to the south and east of Crete for the Wadati-Benioff zone, and over the central part of the island and the area to the northeast of it (cell 11) for the continental wedge, which suggests nonuniform stress field and/or heterogeneous material.  相似文献   

11.
The seismic microzonation of the Bengal Basin, Haldia region, India is carried out using the Analytical Hierarchy Process (AHP) on the Geographic Information System (GIS). Three themes are used for the seismic microzonation, namely Peak Ground Acceleration (PGA), predominant frequency and elevation map. An analysis of the maximum magnitude (m max) and the b value is carried out after preparing the earthquake catalogue from various sources. On the basis of the tectonic set up and seismicity of the region, five seismic zones are delineated which can be a threat to Haldia. They are broadly classified as Zone 1: Arakan-Yoma Zone (AYZ), Zone 2: Himalayan Zone (HZ), Zone 3: Shillong Plateau Zone (SPZ), Zone 4: Bay of Bengal Zone (BBZ) and Zone 5: Shield Zone (SZ). The m max for Zones 1, 2, 3, 4 and 5 are 8.30 ± 0.51, 9.09 ± 0.58, 9.20 ± 0.51, 6.62 ± 0.43 and 6.61 ± 0.43, respectively. The PGA value is computed for Haldia following the attenuation relationship taking the m max of each source zone. The expected PGA at Haldia varies from 0.09–0.19 g. The predominant frequency of Haldia is also calculated using the H/V ratio with a frequency ranging from 0.1–3.0 Hz. The elevation map of Haldia is also generated using the Shuttle Radar Topography Mission (STRM) data. A first-order seismic microzonation map of Haldia is prepared in which four zones of hazard have been broadly classified for Haldia as very high seismic hazard zone, high seismic hazard zone, moderate seismic hazard zone and less seismic hazard zone. The very high seismic hazard zone is observed along the southern part of Haldia where there are major industrial and port facilities. The PGA for the four hazard zones are: 0.09–0.13 g for low hazard zone, > 0.13–0.15 g for moderate hazard zone, > 0.15–0.16 g for high hazard zone and > 0.16–0.19 g for very high hazard zone.  相似文献   

12.
This paper presents the review of the experience in applying the approach based on the limiting distributions of the extreme value theory (the generalized Pareto distribution, GPS, and generalized extreme value distribution, GEV) for deriving the distributions of maximal magnitudes and related ground accelerations from the earthquakes on the future time intervals of a given duration. The results of analyzing the global and regional earthquake catalogs and the ground peak accelerations during the earthquakes are described. It is shown that the magnitude of the strongest possible earthquake M max (and analogous characteristics for other types of data), which is often used in seismic risk assessment, is potentially unstable. We suggest a stable alternative for M max in the form of quantiles Q q (τ) of the maximal possible earthquake, which could occur during the future time interval of length τ. The quantity of the characteristic maximal event M c, which has been introduced in our previous publications, is another helpful robust scalar parameter. All the cases of approximation of the tails of empirical distributions, which were studied in our works, turned out to be finite (bounded); however, the rightmost point of these distributions, M max, is often poorly detectable and unstable. Therefore, the M max parameter has a low practical value.  相似文献   

13.
—A moderately strong earthquake (M w = 6.2) occurred in the town of Dinar at 17.57 UT on October 1, 1995, taking the lives of 90 people and damaging about 4500 buildings. Its epicenter is located near the Dinar-Çivril fault and its focal mechanism is linked to a northeast-southwesterly tensional stress field arising from the interaction between the subducting African plate and the overriding Aegean-Anatolian plate in the eastern Mediterranean.¶Surface cracks of the October 1 earthquake have been observed 10 km continuously along the Dinar-Çivril fault. The cracks have displayed a mode of dip-slip; however, some have also indicated lateral slip. The different modes of slip are generally in agreement with the fault plane solution and are indicators of the complex nature of the rupture process.¶In investigating the earthquake hazard of the Dinar-Çivril fault and proximity, the maximum likelihood method was used to estimate seismic hazard parameters of b-value, seismicity activity rate λ m and the expected maximum magnitude M max?. The data consisted of the historical data covering the period between 1800–1900 and instrumental data between 1900 and 1992. This method, allowing use of the mixed earthquake catalogue containing both historical and instrumental earthquake data, yielded values of 0.70, 1.92 and 7.14 for b, λ m and M max?, respectively. The recurrence time estimated for an earthquake of a magnitude of M w = 6.2 is 123 years. The non-occurrence probabilities of such an earthquake in 1 and 50 years are 0.21 and 0.04, respectively.  相似文献   

14.
Evaluating Tsunami Hazard in the Northwestern Indian Ocean   总被引：1，自引：0，他引：1  
We evaluate here the tsunami hazard in the northwestern Indian Ocean. The maximum regional earthquake calculated from seismic hazard analysis, was used as the characteristic earthquake for our tsunami hazard assessment. This earthquake, with a moment magnitude of M w 8.3 and a return period of about 1000 years, was moved along the Makran subduction zone (MSZ) and its possible tsunami wave height along various coasts was calculated via numerical simulation. Both seismic hazard analysis and numerical modeling of the tsunami were validated using historical observations of the Makran earthquake and tsunami of the 1945. Results showed that the possible tsunami may reach a maximum height of 9.6 m in the region. The distribution of tsunami wave height along various coasts is presented. We recommend the development of a tsunami warning system in the region, and emphasize the value of education as a measure to mitigate the death toll of a possible tsunami in this region.  相似文献   

15.
Deterministic Seismic Zoning of Eastern Cuba   总被引：1，自引：0，他引：1  
—A deterministic seismic zoning of Cuba is performed by modelling, with modal summation, the complete P-SV and SH waves fields generated by point-source earthquakes buried in flat-layered anelastic media. The results of the computation, performed for periods greater than 1 second, are presented in two sets of maps of maximum displacement (d max), maximum velocity (v max) and design ground acceleration (DGA), obtained by using two different criteria in the definition of the input magnitude: (1) values reported in the earthquake catalogue (M obs) and (2) values determined from seismotectonic considerations (M max). A comparison with the results of a previous probabilistic seismic zoning is made to test the possibility of making intensity — ground motion conversion with the aid of log-linear regressions.  相似文献   

16.
In this work we review earthquakes that happened in Southern Siberia and Mongolia within the coordinates of 42°–62° N and 80°–124° E and first propose relationships between earthquake parameters (a surface-wave earthquake magnitude M s and an epicentral intensity(I 0) based on the MSK-64 scale) and maximal distances from an earthquake epicenter (R e max), hypocenter (R h max), and a seismogenic fault (R f max) to the localities of secondary coseismic effects. Special attention was paid to the study of these relationships for the effects of soil liquefaction. Hence, it was shown that secondary deformations from an earthquake were distributed in space away from an earthquake epicenter, than from an associating seismogenic fault. The effects of soil liquefaction are manifested by several times closer to a seismogenic fault, than all other effects, regardless of the type of tectonic movement in a seismic focus. Within the 40 km zone from an earthquake epicenter 44% of the known manifestations of liquefaction process occurred; within the 40 km zone from a seismogenic fault—90%. We propose the next relationship for effects of soil liquefaction: M s = 0.007 × R e max + 5.168 that increases the limits of the maximum epicentral distance at an earthquake magnitude of 5.2 ≤ M s ≤ 8.1 as compared to the corresponding relationships for different regions of the world.  相似文献   

17.
Estimation of Maximum Earthquakes in Northeast India   总被引：1，自引：0，他引：1  
We attempt to estimate possible maximum earthquakes in the northeast Indian region for four seismic source zones, namely EHZ, MBZ, EBZ, and SHZ, which encapsulates the various seismogenic structures of the region and also for combined source zones taken as a single seismic source regime. The latter case exhibits a high maximum earthquake estimate of MW 9.4 (±0.85) through Bayesian interpretation of frequency magnitude distribution with Gamma function implicating a moderate deviation from the standard Gutenberg Richter model at the higher magnitudes. However, tapering Gutenberg Richter models with corner magnitudes at MW 8.01, 8.7 and 9.1, respectively indicated maximum values corresponding to MW 8.4, 9.0, and 9.3. The former approach was applied to each of the source zones wherein the data are presented in parts according to the data completeness, thereof. EHZ, MBZ, EBZ and SHZ are seen with maximum earthquakes of MW 8.35 (±0.59), 8.79 (±0.31), 8.20 (±0.50), and 8.73 (±0.70), respectively. The maximum possible earthquakes estimated for each individual zone are seen to be lower than that estimated for the single regime. However, the pertaining return periods estimated for the combined zone are far less than those estimated for the demarcated ones.  相似文献   

18.
In the present study, the level of the largest earthquake hazard is assessed in 28 seismic zones of the NW Himalaya and its vicinity, which is a highly seismically active region of the world. Gumbel’s third asymptotic distribution (hereafter as GIII) is adopted for the evaluation of the largest earthquake magnitudes in these seismic zones. Instead of taking in account any type of Mmax, in the present study we consider the ω value which is the largest earthquake magnitude that a region can experience according to the GIII statistics. A function of the form Θ(ω, RP6.0) is providing in this way a relatively largest earthquake hazard scale defined by the letter K (K index). The return periods for the ω values (earthquake magnitudes) 6 or larger (RP6.0) are also calculated. According to this index, the investigated seismic zones are classified into five groups and it is shown that seismic zones 3 (Quetta of Pakistan), 11 (Hindukush), 15 (northern Pamirs), and 23 (Kangra, Himachal Pradesh of India) correspond to a “very high” K index which is 6.  相似文献   

19.
总结2015年以来、2010年以来和1980年以来祁连山地震带活动增强指标,得到祁连山地震带5年、10年以及40年尺度的小震频度预报效能。祁连山地震带ML2.0以上地震月频度持续增强的异常出现时间在3个月以上,异常出现后1~3个月,祁连山地震带及边邻地区发生5级以上地震以及青藏高原东北缘发生6级以上地震的概率较高。作为短期预报指标,祁连山地震带的地震活动增强对甘肃及边邻地区的地震预报工作具有重要的指示意义。  相似文献   

20.
ＰｒｏｂａｂｉｌｉｔｙｆｏｒｅｃａｓｔｏｆｅａｒｔｈｑｕａｋｅｍａｇｎｉｔｕｄｅｉｎＣｈｉｎｅｓｅｍａｉｎｌａｎｄｂｅｆｏｒｅＡ．Ｄ．２００５ＸＩＡＯ－ＱＩＮＧＷＡＮＧ（王晓青），ＺＨＥＮＧ－ＸＩＡＮＧＦＵ（傅征祥）ａｎｄＭＩＮＧＪＩＡＮＧ（蒋铭）...  相似文献   

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

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