首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The 3 strongest earthquakes,M7.0, which have occurred since 1973 in the area of Greece were preceded by a specific increase of the earthquake activity in the lower magnitude range. This activation is depicted by algorithm M8. This algorithm of intermediate term earthquake prediction was originally designed for diagnosis by Times of Increased Probability (TIPs) of the strongest earthquake,M8.0 worldwide (Keilis-Borok andKossobokov, 1984). At present the algorithm is retrospectively tested for smaller magnitudes in different seismic regions (Keilis-Borok andKossobokov, 1986, 1988). A TIP refers to a time period of 5 years and an area whose linear size is proportional and several times larger than that of the incipient earthquake source. Altogether the TIPs diagnosed by the algorithm M8 in the area of Greece occupy less than 20% and the Times of Expectation (TEs) about 10% of the total space-time domain considered. Also there is a current TIP for the southeastern Aegean sea and 1988–1992. It may specify the long-term prediction given inWyss andBaer (1981a,b).The results of this study are further evidence favoring applicability of algorithm M8 in diverse seismotectonic environment and magnitude ranges and support indirectly the hypothesis of self-similarity of the earthquake activity. It also implies the possibility of intermediate term prediction of the strongest earthquakes in the area of Greece.  相似文献   

2.
A sequence of moderate shallow earthquakes (3.5M L5.3) was located within the Vercors massif (France) in the period 1961–1984. This subalpine massif has been a low seismic area for at least 5 centuries. During the period 1962–1963, 12 shallow earthquakes occurred in the neighborhood (10 km) of the Monteynard reservoir, 30 km south of the city of Grenoble. The latest fourM L4.0 earthquakes occurred in 1979–1984 either at larger distance (35 km) or greater depth (10 km) from the reservoir. Two triggering mechanisms are suggested for this sequence: (i) the direct effect of elastic loading through either increased shear stress or strength reducing by increased pore pressure at depth; (ii) the pore pressure diffusion induced by poroelastic stress change due to the reservoir filling.The weekly water levels, local balanced geological cross sections, and focal mechanisms argue for two types of mechanical connection between the earthquake sequence and the filling cycles of the Monteynard reservoir. The seismic sequence started with the 1962–1963 shallow earthquakes that occurred during the first filling of the reservoir and are typical of the direct effect of elastic loading. The 1979 deeper earthquake is located at a 10 km depth below the reservoir. This event occurred 16 years after the initial reservoir impoundment, but one month after the previous 1963 maximum water level was exceeded. Moreover the yearly reservoir level increased gradually in the period 1962–1979 and has decreased since 1980. Accordingly we suggest that the gradual diffusion of water from reservoir to hypocentral depths decreases the strength of the rock matrices through increased pore pressure. The transition between the two types of seismic response is supported by the analysis ofM L3.5 earthquakes which all occurred in the period 1964–1971, ranging between 10 and 30 km distance from the reservoir. The three other delayed earthquakes of the 1961–1984 seismic sequence (M L4 during the 1979–1984 period) are all located 35 km away from the reservoir. Based on the seismic activity, the estimates for the hydraulic diffusivities range between 0.2–10 m2/s, except for the first event that occurred 30 km north of the reservoir, the filling just started. The lack ofin situ measurements of crustal hydrological properties in the area, shared by most of the Reservoir-Induced-Seismicity cases, prevents us from obtaining absolute evidence for the triggering processes. These observations and conceptual models attest that previous recurrence times for moderate natural shocks (4.5M L5.5) estimated within this area using historical data, could be modified by 0.1–1 MPa stress changes. These small changes in deviatoric stress suggest that the upper crust is in this area nearly everywhere at a state of stress near failure. Although the paucity of both number and size of earthquakes in the French subalpine massif shows that aseismic displacements prevail, our study demonstrates that triggered earthquakes are important tools for assessing local seismic risk through mapping fault zones and identifying their possible seismic behavior.  相似文献   

3.
— Analysis of the Koyna-Warna earthquake catalog (1968–1996) shows that on an average there is a positive correlation between the b value (decrease) and fractal dimensions (decrease in both D2s and D2t) of earthquake epicenters 0.5 and 2.5 years prior to 1973 (M5.2) and 1980 (M5.5) events, respectively, except a negative correlation for about five years (1988–1993) prior to the 1993/1994 sequence (M5.4). This positive correlation indicates a weaker clustering, or that the epicenters tend to fill the two-dimensional plane. While the origin of the negative correlation seems to be that during periods of large events (low b value), there is strong clustering around the main shock epicenter (high fractal dimension). Interestingly, during the last year (1995–1996) of the studied period both the b value and correlation dimensions rose significantly, suggesting that stress release occurs through increased levels of low magnitude and increasingly scattered seismicity, suggesting an increased risk of larger magnitude events. Incidentally, during 2000 three earthquakes of magnitude M 5.0, one earthquake of M 4.0, 45 earthquakes of magnitude M 3.0–3.9, and several thousand earthquakes of M < 3 have occurred in the region. Thus it can be inferred that at local scales the relationship yields both positive and negative correlation that appears to be controlled by different modes of failure within the active fault complex.Acknowledgement. The authors are grateful to Dr. B.K. Rastogi of NGRI for providing the catalog of Koyna earthquakes and for useful scientific discussions. The comments of Dr. I. G. Main have improved the quality of paper for which we extend to him our sincere thanks. One of the authors (AOM) thanks the Third World Academy of Science and the Council of Scientific and Industrial Research, India for the Postdoctoral Fellowship award under which this work was carried out.  相似文献   

4.
Reservoir induced earthquakes began to occur in the vicinity of Shivajisagar Lake formed by Koyna Dam in Maharashtra state, western India, soon after its filling started in 1962. Induced earthquakes have continued to occur for the past 34 years in the vicinity of this reservoir, and so far a total of 10 earthquakes of M 5.0, over 100 of M 4 and about 100,000 of M 0.0 have occurred. Every year, following the rainy season, the water level in the reservoir rises and induced earthquakes occur. Seismic activity during 1967–68 was most intense when globally, the largest reservoir induced earthquake occurred on 10 December, 1967. Other years of intense seismic activity are 1973 and 1980. During 1986 another reservoir, Warna, some 20 km south of Koyna, began to be filled. The recent burst of seismic activity in Koyna-Warna region began in August, 1993, and was monitored with a close network of digital and analog seismographs. During August, 1993–December, 1995, 1,272 shocks of magnitude 2 were located, including two earthquakes of M 5.0 and M 5.4 on 8 December, 1993 and 1 February, 1994, respectively. Two parallel epicentral trends in NNE-SSW direction, one passing through Koyna and the other through Warna reservoir are delineated. The 1993 increase in seismicity has followed a loading of 44.15 m in Warna reservoir during 11 June 11, 1993 through August 4, 1993, with a maximum rate of filling being 16 m/week. The larger shocks have been found to be preceded by a precursory nucleation process.  相似文献   

5.
The most complete and reliable data of strong (M s6.5), shallow (h<70 km) earthquakes which occurred in the inner Aegean seismic zone have been utilized to describe its seismicity time variation during 1800–1986 by two independent statistical models. The first is a sequentially stationary model of seismicity rates which shows that intervals of low seismicity rate, lasting for some 37 years, alternate with high rate intervals of 8–12 years duration. The second model is a statistical model according which seismic energy released within 5-year time windows approximates a harmonic curve within a period of about 50 years. This model is in agreement with the notion that the time series of strong earthquake occurrences in the inner Aegean seismic zone consists of a random (shocks withM s=6.5–6.8) and a nonrandom component (M s6.9). Maxima and minima of the harmonic curve coincide with the high and low rate intervals, respectively. A model of regional stationary accumulation of thermal stresses along certain seismic belts and their cyclic relaxation may explain this periodicity.  相似文献   

6.
In this paper we evaluate the present state of the seismic regime in Southern California using the concentration parameter of seismogenic faults (K sf ,Sobolev andZavyalov, 1981). The purpose of this work is to identify potential sites for large earthquakes during the next five or ten years. The data for this study derived from the California Institute of Technology's catalog of southern California earthquakes, and spanned the period between 1932 to June 1982. We examined events as small asM L 1.8 but used a magnitude cutoff atM L =3.3 for a detailed analysis. The size of the target earthquakes (M M ) was chosen as 5.3 and 5.8.The algorithm for calculatingK sf used here was improved over the algorithm described bySobolev andZavyalov (1981) in that it considered the seismic history of each elementary seismoactive volume. The dimensions of the elementary seismoactive volumes were 50 km×50 km and 20 km deep. We found that the mean value ofK sf within 6 months prior to the target events was 6.1±2.0 for target events withM L 5.3 and 5.41.8 for targets withM L 5.8. Seventy-three percent of the targets withM L 5.8 occurred in areas whereK sf was less than 6.1. The variance of the time between the appearance of areas with lowK sf values and the following main shocks was quite large (from a few months to ten years) so this parameter cannot be used here for accurate predictions of occurrence time.Regions where the value ofK sf was below 6.1 at the end of our data set (June, 1982) are proposed as the sites of target earthquakes during the next five to ten years. The most dangerous area is the area east of San Bernardino whereK sf values are presently between 2.9 and 3.7 and where there has been no earthquake withM L 5.3 since 1948.  相似文献   

7.
Summary Mean travel-time residuals of P waves for the period 1964–1970 at stations in North America and Europe are well separated into large domains of positive and negative values. The spatial distribution of residuals in North America is in good agreement with the structure of the upper mantle obtained from Rayleigh wave dispersion and is in accord with magnitude anomalies and heat flow data. A systematic variation of residuals in time is found to be sensitive to major changes in the nature of the earthquake source distribution, as for example aftershocks, as well as to changes in methods of detection of first arrivals. Residuals at neighboring stations are correlated up to distances of the order of 2°. Half year mean residuals at stations within 2° distance from the epicenters of strong earthquakes have a sudden jump 1.5 years before the occurrence of an earthquake. A formal algorithm based on this pattern permits prediction of the times of all 16 earthquakes with magnitudeM7.5 occurring during the interval 1966–1972. Although the average duration of the alarm periods is about half the total time interval of the catalog, so that the algorithm has no practical importance, these results are statistically significant on a 99% level of confidence.  相似文献   

8.
Aftershocks or swarms indicate increase of the flow intensity in the vicinity of the initial earthquakes. By normalizing their number according to the dynamic range of the standard frequency magnitude distribution the increase or positive aftereffect property of the initial earthquakes can be compared for different magnitude intervals, periods of time or regions. After applying accurate formal algorithm of aftershock identification it is possible to study negative aftereffect of the main events (nonaftershocks) in the catalog.Negative aftereffect means decrease of the probability of successive events in a time-space vicinity of the main event, when the aftershocks are over. The negative effect is the most important part of the seismic cycle and seismic gaps approach. Global statistical test give high confidence level for the relative decrease in intensity of the flow of the events withM7 in the first 20–25 years after the events withM8 in their 1o-vicinities in the total time period under study of approximately 60 years. The decrease approximates 32% of the undisturbed intensity of the flow ofM>7 events in the vicinities.Self-similar negative aftereffect was observed 3–7 years after 6M<7 events, it totals approximately 18% of the undisturbed intensity. Another type of self-similarity of seismic regime, with respect to the negative aftereffect, is the decrease of probabilities of aftershocks with large magnitudes in aftershock sequences. When we have adequate dynamic range in the catalog for the study of this property, for example, for main events withM7 in the catalog with low cut-off limitM=4, the statistical significance of the negative aftereffect is clear. However, the absolute value of the effect is also rather small, about 10%, which means that in 90% of the cases the aftershock sequences do not experience lack of energy due to the main shock energy release and follow a standard magnitude distribution for earthquakes in the entire catalog.The small values of the negative aftereffect apparently indicate partial stress relase by earthquakes and may explain short recurrence time intervals after major earthquakes observed periodically in different places.  相似文献   

9.
The algorithm CN makes use of normalized functions. Therefore the original algorithm, developed for the California-Nevada region, can be directly applied, without adjustment of the parameters, to the determination of the Time of Increased Probability (TIP) of strong earthquakes for Central Italy. The prediction is applied to the events with magnitudeMM 0=5.6, which in Central Italy have a return period of about six years. The routinely available digital earthquake bulletins of the Istituto Nazionale di Geofisica (ING), Rome, permits continuous monitoring. Here we extend to November 1994 the first study made by Keilis-Boroket al. (1990b). On the basis of the combined analysis of seismicity and seismotectonic, we formulate a new regionalization, which reduces the total alarm time and the failures to predict, and narrows the spatial uncertainty of the prediction with respect to the results ofKeilis-Borok et al. (1990b).The premonitory pattern is stable when the key parameters of the CN algorithm and the duration of the learning period are changed, and when different earthquake catalogues are used.The anlysis of the period 1904–1940, for whichM 0=6, allows us to identify self-similar properties between the two periods, in spite of the considerably higher seismicity level of the earlier time interval compared with the recent one.  相似文献   

10.
Earthquake Triggering along the Xianshuihe Fault Zone of Western Sichuan,China   总被引：19，自引：0，他引：19  
Western Sichuan is among the most seismically active regions in southwestern China and is characterized by frequent strong (M 6.5) earthquakes, mainly along the Xianshuihe fault zone. Historical and instrumental seismicity show a temporal pattern of active periods separated by inactive ones, while in space a remarkable epicenter migration has been observed. During the last active period starting in 1893, the sinistral strike–slip Xianshuihe fault of 350 km total length, was entirely broken with the epicenters of successive strong earthquakes migrating along its strike. This pattern is investigated by resolving changes of Coulomb failure function (CFF) since 1893 and hence the evolution of the stress field in the area during the last 110 years. Coulomb stress changes were calculated assuming that earthquakes can be modeled as static dislocations in an elastic halfspace, and taking into account both the coseismic slip in strong (M 6.5) earthquakes and the slow tectonic stress buildup associated with major fault segments. The stress change calculations were performed for faults of strike, dip, and rake appropriate to the strong events. We evaluate whether these stress changes brought a given strong earthquake closer to, or sent it farther from, failure. It was found that all strong earthquakes, and moreover, the majority of smaller events for which reliable fault plane solutions are available, have occurred on stress–enhanced fault segments providing a convincing case in which Coulomb stress modeling gives insight into the temporal and spatial manifestation of seismic activity. We extend the stress calculations to the year 2025 and provide an assessment for future seismic hazard by identifying the fault segments that are possible sites of future strong earthquakes.  相似文献   

11.
The high likelihood of a gap-filling thrust earthquake in the Alaska subduction zone within this decade is indicated by two independent methods: analysis of historic earthquake recurrence data and time-to-failure analysis applied to recent decades of instrumental data. Recent (May 1993) earthquake activity in the Shumagin Islands gap is consistent with previous projections of increases in seismic release, indicating that this segment, along with the Alaska Peninsula segment, is approaching failure. Based on this pattern of accelerating seismic release, we project the occurrence of one or moreM7.3 earthquakes in the Shumagin-Alaska Peninsula region during 1994–1996. Different segments of the Alaska-Aleutian seismic zone behave differently in the decade or two preceding great earthquakes, some showing acceleration of seismic release (type A zones), while others show deceleration (type D zones). The largest Alaska-Aleutian earthquakes—in 1957, 1964, and 1965—originated in zones that exhibit type D behavior. Type A zones currently showing accelerating release are the Shumagin, Alaska Peninsula, Delarof, and Kommandorski segments. Time-to-failure analysis suggests that the large earthquakes could occur in these latter zones within the next few years.  相似文献   

12.
By comparing seasonal rainfall data from the past 90 years with the occurrence of large (M6) earthquakes along an arid stretch of the San Andreas fault system in southern California, certain correlations have been observed. Most large earthquakes are preceded by a pattern consisting of a few years of below normal precipitation (drought) terminated by one or more consecutive seasons of heavy (above normal) rainfall. While this drought-above normal rainfall cycle can be seen at times other than prior to major earthquakes, it precedes, to varying degrees, all of the twelve M6 events. This new precursor evidence, when combined with other premonitory signals, may offer a helpful diagnostic measure that could be useful in earthquake prediction in arid regions.  相似文献   

13.
The largest earthquake (M0=4.9·1027 dyn·cm) of the 20th century in the territory of Greece occurred south of Amorgos Island, causing extensive destruction in the southern Aegean area. It occurred on an ENE–trending normal fault that is seated parallel to the Islands southern coastline. Changes in the rates of moderate–size earthquakes (M 5.0) that occurred before and after the Amorgos earthquake, within circular regions centered on its epicenter with radii of 100, 150 and 200 km, are investigated. The rate for moderate–size events just before the main shock appears to be considerably increased when compared to those of either preceding or subsequent periods. Further inspection reveals that more evident seismicity fluctuations are attributed to distances exceeding 100 km. These changes may be indicative of a broad region that is approaching a high stress state prior to an eventual large earthquake. Close to the main event, that is, within the 100–km radius, a remarkable quiescence period lasting about two decades before its occurrence was observed. Changes in seismicity are discussed in combination with static stress changes calculated by the application of the stress evolutionary model that takes into account the coseismic slip associated with the larger events (M 6.5) since the beginning of the 20th century and the tectonic loading on the major faults in the study area. These larger events, as with the intermediate magnitude seismicity taking place at distances exceeding 100 km and which encircled the quiescent area observed during the last 22 years before the Amorgos earthquake, are well correlated with stress-enhanced areas in each stage of the evolutionary model.  相似文献   

14.
The distributions of discrete frequency, N, versus interoccurrence time, t (in days), of M 7 earthquakes in the Taiwan region during the 1900–1994 period, M 6 earthquakes in the north-south seismic belt of China during the 1900–1990 period, and M 5.5 earthquakes in Southern California, U.S.A., during the 1914–1995 period are studied through two statistical models (gamma function and exponential function). Results show that both the exponential function and gamma function can describe the distributions. However, the former is more appropriate than the latter. This indicates that the three time series of earthquakes have a significant component of Poisson processes, even though the tectonic conditions, the fault distributions and the size of the three seismic regions are different.  相似文献   

15.
Strong ground motions are estimated for the Pacific Northwest assuming that large shallow earthquakes, similar to those experienced in southern Chile, southwestern Japan, and Colombia, may also occur on the Cascadia subduction zone. Fifty-six strong motion recordings for twenty-five subduction earthquakes ofM s7.0 are used to estimate the response spectra that may result from earthquakesM w<81/4. Large variations in observed ground motion levels are noted for a given site distance and earthquake magnitude. When compared with motions that have been observed in the western United States, large subduction zone earthquakes produce relatively large ground motions at surprisingly large distances. An earthquake similar to the 22 May 1960 Chilean earthquake (M w 9.5) is the largest event that is considered to be plausible for the Cascadia subduction zone. This event has a moment which is two orders of magnitude larger than the largest earthquake for which we have strong motion records. The empirical Green's function technique is used to synthesize strong ground motions for such giant earthquakes. Observed teleseismicP-waveforms from giant earthquakes are also modeled using the empirical Green's function technique in order to constrain model parameters. The teleseismic modeling in the period range of 1.0 to 50 sec strongly suggests that fewer Green's functions should be randomly summed than is required to match the long-period moments of giant earthquakes. It appears that a large portion of the moment associated with giant earthquakes occurs at very long periods that are outside the frequency band of interest for strong ground motions. Nevertheless, the occurrence of a giant earthquake in the Pacific Northwest may produce quite strong shaking over a very large region.  相似文献   

16.
IntroductionI.thasbeenpayingattentiontotheanomalousphenomenonbyseismologiststhattheearthquakeactivityoftenstrengthensinacertainspaceandduringacertaintimearoundthesourceareabeforeastrongearthquake.Mogi(1969)pointedoutthataringdistributiveareaofthestrengthenedearthquakeactivityoftenappearsaroundthesourcearea10to20yearsbeforegreatshock,whichcalledaringphenomenon(or"doughnut").Afterthat,otherscholarsreportedinsuccessionthatringdistributionofearthquakeactivityoccurredbeforeastrongertquakeorevenbefo…  相似文献   

17.
To better understand the mechanics of subduction and the process of breaking a mature seismic gap, we study seismic activity along the western New Britain subduction segment (147°E–151°E, 4°S–8°S) through earthquakes withm b 5.0 in the outer-rise, the upper area of subducting slab and at intermediate depths to 250 km, from January 1964 to December 1990. The segment last broke fully in large earthquakes of December, 28, 1945 (M s =7.9) and May 6, 1947 (M s =7.7.), and its higher seismic potential has been recognized byMcCann et al., (1979). Recently the segment broke partially in two smaller events of February, 8, 1987 (M s =7.4) and October 16, 1987 (M s =7.4), leaving still unbroken areas.We observe from focal mechanisms that the outer-rise along the whole segment was under pronounced compression from the late 60's to at least October 1987 (with exception of the tensional earthquake of December 11, 1985), signifying the mature stage of the earthquake cycle. Simultaneously the slab at intermediate depths below 40 km was under tension before the earthquake of October 16, 1987. That event, with a smooth rupture lasting 32 sec, rupture velocity of 2.0 km/sec, extent of approximately 70 km and moment of 1.2×1027 dyne-cm, did not change significantly the compressive state of stress in the outer-rise of that segment. The earthquake did not fill the gap completely and this segment is still capable of rupturing either in an earthquake which would fill the gap between the 1987 and 1971 events, or in a larger magnitude event (M s =7.7–7.9), comparable to earthquakes observed in that segment in 1906, 1945 and 1947.  相似文献   

18.
Microzonation of Bucharest: State-of-the-Art   总被引：3，自引：0，他引：3  
— The 1940 (Mw=7.7) and 1977 (Mw=7.4) Vrancea earthquakes (Romania) inflicted heavy damage and casualties in Bucharest and the statistics indicate a recurrence interval of 25 years for Mw 7.0 events. Under these circumstances, the seismic microzonation represents important information for detailed urban planning that establishes an appropriate level of preparedness to the earthquake threat. This paper reviews the main studies concerning the seismicity of the Vrancea region, the site conditions of the city, the characterization of the building stock, and the codes of practice that regulate the antiseismic design. The first-order microzonation of Bucharest was performed starting from the existing database of structural and geotechnical parameters. New insights originating from direct instrumental observation and interpretation of the local effects as well as realistic numerical modeling that update and improve the input data necessary for a detailed microzoning map of the city are also discussed.  相似文献   

19.
We review earthquake distributions associated with hydrocarbon fields in the context of pore pressure diffusion models, poroelastic stress transfer and isostasy theory. These three mechanisms trigger or induce seismic instabilities at both local scale (D5 km) and at regional scale (D20 km). The modeled changes in stress are small (1 MPa), whatever the tectonic setting. Each mechanism corresponds to different production processes. (1) Local hydraulic fracturing due to fluid injection induces seismic-slip on cracks (M L3) within the injected reservoir through decreasing the effective stress. (2) Pure fluid withdrawal causes pore pressure to decrease within the reservoir. It triggers adjustments of the geological structure to perturbations related to the reservoir response to depletion. Poroelastic mechanisms transfer this stress change from the reservoir to the surrounding levels whereM L5 seismic instabilities occur either above or below the reservoir. (3) Massive hydrocarbon recovery induces crustal readjustments due to the removal of load from the upper crust. It can induce larger earthquakes (M L6) at greater distance from the hydrocarbon fields than the two other mechanisms.Due to the mechanical properties of the shallow rock matrices involved, seismic slip triggered either by mechanism (1) or (2), is a second-order process of the main elastoplastic deformation. for a minimum of 80% of commercially productive basins, most of the local deformation is reported as aseismic, i.e., there is no evidence forM L3 earthquakes. Nevertheless, the induced stresses vary as a function of time in a manner that depends on the hydraulic diffusivity (i.e., permeability) of the reservoir and surrounding rocks. Because small earthquakes (M L3) indicate changes in stress and pore pressure, monitoring of seismicity is a means of assessingin situ reservoir behavior.The less constrained seismic response to hydrocarbon recovery is the possible connection between local fluid manipulations, triggered earthquakes and major regional earthquakes. Positive feedback mechanisms suggest that the region of seismic hazard changes is much larger than the area where hydrocarbons are extracted. These observations and models testify that fluid movement and pore pressure changes (increase or decrease) play important roles in the mechanics of earthquakes and in the triggering of natural earthquakes.  相似文献   

20.
The Bayesian extreme-value distribution of earthquake occurrences has been adopted to estimate the seismic hazard in some seismogenic zones in Greece and surrounding regions. Seismic moment, slip rate, earthquake recurrence rate and magnitude were considered as basic parameters for computing the prior estimates of the seismicity. These estimates are then updated in terms of Bayes' theorem and historical estimates of seismicity associated with each zone.High probabilities for earthquakes withM6.0 have been obtained for the northwestern part of Greece as well as for the southwestern part of the Hellenic arc.  相似文献   

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

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