首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 27 毫秒
1.
We have imaged earthquake source zones beneath the northeast India region by seismic tomography, fractal dimension and b value mapping. 3D P-wave velocity (Vp) structure is imaged by the Local Earthquake Tomography (LET) method. High precision P-wave (3,494) and S-wave (3,064) travel times of 980 selected earthquakes, m d ≥ 2.5, are used. The events were recorded by 77 temporary/permanent seismic stations in the region during 1993–1999. By the LET method simultaneous inversion is made for precise location of the events as well as for 3D seismic imaging of the velocity structure. Fractal dimension and seismic b value has been estimated using the 980 LET relocated epicenters. A prominent northwest–southeast low Vp structure is imaged between the Shillong Plateau and Mikir hills; that reflects the Kopili fault. At the fault end, a high-Vp structure is imaged at a depth of 40 km; this is inferred to be the source zone for high seismic activity along this fault. A similar high Vp seismic source zone is imaged beneath the Shillong Plateau at 30 km depth. Both of the source zones have high fractal dimension, from 1.80 to 1.90, indicating that most of the earthquake associated fractures are approaching a 2D space. The spatial fractal dimension variation map has revealed the seismogenic structures and the crustal heterogeneities in the region. The seismic b value in northeast India is found to vary from 0.6 to 1.0. Higher b value contours are obtained along the Kopili fault (~1.0), and in the Shillong Plateau (~0.9) The correlation coefficient between the fractal dimension and b value is found to be 0.79, indicating that the correlation is positive and significant. To the south of Shillong Plateau, a low Vp structure is interpreted as thick (~20 km) sediments in the Bengal basin, with almost no seismic activity in the basin.  相似文献   

2.
We study the influence of different source characteristics (depth, distance, type and azimuth) on the site effect in Acapulco and the Valley of Mexico. Site amplification was estimated by means of spectral ratios (both horizontal-to-vertical spectral ratio and standard spectral ratio techniques were applied) from earthquake recordings at soft and hard sites. In Acapulco, 125 Mexican earthquakes covering a hypocentral range of 7–290 km and a depth range (H) of 3–61 km were analyzed in three different groups of hypocentral distances. In the Valley of Mexico, we estimate site effect at five locations using recordings from shallow-coastal interplate (200?Δ?570 km; H?35 km) and normal-faulting, intermediate-depth inslab (132?Δ?738 km; 32?H?178 km) earthquakes, as well as from teleseismic events. Our results seem to point to negligible dependence of site effects on the source location and type.  相似文献   

3.
This study analyzes the S-wave envelope broadening characteristics of 290 earthquakes recorded by 14 stations of the Spanish National Seismograph Network in the Canary Islands region. The S-wave peak delay time (t p ) and envelope duration (t q ) parameters are evaluated phenomenologically to infer the strength of velocity inhomogeneities of the medium along each seismic ray path. Crustal (0?≤?h?≤?18 km) and upper mantle (18?<?h?≤?80 km) events are analyzed separately. Results in the frequency range 1 to 12 Hz for hypocentral distances from 30 to 600 km show that both t p and t q increase according to a power of hypocentral distance and they are independent of frequency. The spatial distribution of the peak delay time reveals weak strength of heterogeneity in most of the region at shallow depths. Relatively strong inhomogeneous zones are generated under the island of Tenerife and Gran Canaria at depths of 11–22 km.  相似文献   

4.
—?The digital data acquired by 16 short-period seismic stations of the Friuli-Venezia-Giulia seismic network for 56 earthquakes of magnitude 2.3–4.7 which occurred in and near NE Italy have been used to estimate the coda attenuation Q c and seismic source parameters. The entire area under study has been divided into five smaller regions, following a criterion of homogeneity in the geological characteristics and the constrains imposed by the distribution of available events. Standard IASPEI routines for coda Q c determination have been used for the analysis of attenuation in the different regions showing a marked anomaly in the values measured across the NE border between Friuli and Austria for Q 0 value. A large variation exists in the coda attenuation Q c for different regions, indicating the presence of great heterogeneities in the crust and upper mantle of the region. The mean value of Q c (f) increases from 154–203 at 1.5?Hz to 1947–2907 at 48?Hz frequency band with large standard deviation estimates.¶Using the same earthquake data, the seismic-moment, M 0, source radius, r and stress-drop, Δσ for 54 earthquakes have been estimated from P- and S-wave spectra using the Brune's seismic source model. The earthquakes with higher stress-drop (greater than 1?Kbar) occur at depths ranging from 8 to 14?km.  相似文献   

5.
We investigate the scaling relationships among earthquake source parameters using more than 300 good quality broad band seismograms from 30 small earthquakes in the Kumaon Himalaya from the spectral analysis of P and S waves. The average ratio of P/S wave corner frequency is found to be 1.13, which is suggestive of shift in the corner frequency. The estimated seismic moment range from 1.6?×?1013–5.8?×?1015 N?m, while the stress drop varies from 0.6 to 16 bars with 80 % of the events below 10 bars. An analysis of stress drop and apparent stress drop indicates the partial stress drop mechanism in the region. The source radii are between 0.17 and 0.88 km. The total seismic energy varies from 1.79?×?108 to 7.30?×?1011 J. We also observe the variation in seismic energy for a given seismic moment. The scaling relation between the seismic moment and stress drop suggests the breakdown of constant stress drop scaling for the range of seismic moments obtained here for the region. This shows the anomalous behavior of small earthquakes in the region. The study indicates that the stress drop is the dominant scaling factor for the moments studied here.  相似文献   

6.
We study the microseismicity (M L ?<?2) in the region of Landau, SW Germany. Here, due to thick sediments (~3?km) and high cultural seismic noise, the signal-to-noise ratio is in general very low for microearthquakes. To gain new insights into the occurrence of very small seismic events, we apply a three-step detection approach and are able to identify 207 microseismic events (?1?<?M L ?<?~1) with signal-to-noise ratios smaller than 3. Recordings from a temporary broadband network are used with station distances of approximately 10?km. First, we apply a short-term to long-term average detection algorithm for data reduction. The detection algorithm is affected severely by transient noise signals. Therefore, the most promising detections, selected by coinciding triggers and high-amplitude measures, are reviewed manually. Thirteen seismic events are identified in this way. Finally, we conduct a cross-correlation analysis. As master template, we use the stacked waveforms of five manually detected seismic events with a repeating waveform. This search reveals additional 194 events with a cross-correlation coefficient exceeding 0.65 which ensures a stable identification. Our analysis shows that the repeating events occurred during the stimulation of a geothermal reservoir within a source region of only about 0.5?km3. Natural background seismicity exceeding our detection level of M L ?~?0.7 is not found in the region of Landau by our analysis.  相似文献   

7.
During the present tectonic activity in the volcanic rift zone in NE-Iceland it has become apparent that the attenuation of seismic waves is highly variable in the central region of the Krafla volcano. Earthquakes associated with the inflation of the volcano have been used to delineate two regions of high attenuation of S-waves within the caldera. These areas are located near the center of inflation have horizontal dimensions of 1–2 km and are interpreted as the expression of a magma chamber. The top of the chamber is constrained by hypocentral locations and ray paths to be at about 3 km depth. Small pockets of magma may exist at shallower levels. The bottom of the chamber is not well constrained, but appears to be above 7 km depth. Generally S-waves propagate without any anomalous aftenuation through laver 3 (vp=0.5 km sec?1) across the volcanic rift zone in NE-Iceland. The rift zone therefore does not appear to be underlain by an estensive magma chamber at crustal levels. The Krafla magma chamber is a localized feature of the Krafla central volcano.  相似文献   

8.
Spain is a low-to-moderate seismicity area with relatively low seismic hazard. However, several strong shallow earthquakes have shaken the country causing casualties and extensive damage. Regional seismicity is monitored and surveyed by means of the Spanish National Seismic Network, maintenance and control of which are entrusted to the Instituto Geográfico Nacional. This array currently comprises 120 seismic stations distributed throughout Spanish territory (mainland and islands). Basically, we are interested in checking the noise conditions, reliability, and seismic detection capability of the Spanish network by analyzing the background noise level affecting the array stations, errors in hypocentral location, and detection threshold, which provides knowledge about network performance. It also enables testing of the suitability of the velocity model used in the routine process of earthquake location. To perform this study we use a method that relies on P and S wave travel times, which are computed by simulation of seismic rays from virtual seismic sources placed at the nodes of a regular grid covering the study area. Given the characteristics of the seismicity of Spain, we drew maps for M L magnitudes 2.0, 2.5, and 3.0, at a focal depth of 10 km and a confidence level 95 %. The results relate to the number of stations involved in the hypocentral location process, how these stations are distributed spatially, and the uncertainties of focal data (errors in origin time, longitude, latitude, and depth). To assess the extent to which principal seismogenic areas are well monitored by the network, we estimated the average error in the location of a seismic source from the semiaxes of the ellipsoid of confidence by calculating the radius of the equivalent sphere. Finally, the detection threshold was determined as the magnitude of the smallest seismic event detected at least by four stations. The northwest of the peninsula, the Pyrenees, especially the westernmost segment, the Betic Cordillera, and Tenerife Island are the best-monitored zones. Origin time and focal depth are data that are far from being constrained by regional events. The two Iberian areas with moderate seismicity and the highest seismic hazard, the Pyrenees and Betic Cordillera, and the northwestern quadrant of the peninsula, are the areas wherein the focus of an earthquake is determined with an approximate error of 3 km. For M L 2.5 and M L 3.0 this error is common for almost the whole peninsula and the Canary Islands. In general, errors in epicenter latitude and longitude are small for near-surface earthquakes, increasing gradually as the depth increases, but remaining close to 5 km even at a depth of 60 km. The hypocentral depth seems to be well constrained to a depth of 40 km beneath the zones with the highest density of stations, with an error of less than 5 km. The M L magnitude detection threshold of the network is approximately 2.0 for most of Spain and still less, almost 1.0, for the western sector of the Pyrenean region and the Canary Islands.  相似文献   

9.
糯扎渡水库地区地震活动和震源参数研究   总被引：1，自引：0，他引：1  
潘颖  张晓东  付虹 《地震》2015,35(3):31-43
本文利用糯扎渡水库台网和景洪水库台网记录的地震波形反演震源谱, 计算了2011年1月至2014年3月糯扎渡水库附近143个1.0级以上地震的震源参数。 分析研究发现: ① 糯扎渡水库蓄水后, 库区内地震活动明显增强, 尤其是在水位上升3~4个月后; ② 水库开始蓄水后, 库区内地震震源深度变浅, 一段时间以后则与库区外的差别减小; ③ 地震拐角频率随地震矩增大而减小, 且蓄水以后, 库区内地震拐角频率的对数与地震矩的对数线性关系更明显; ④ 应力降、 视应力与地震矩存在正相关关系, 且相同地震矩时库区内的应力降与视应力低于蓄水前和库区外的地震应力降与视应力值; ⑤ 蓄水对库区内地震的影响可延伸至库底10 km深度处, 且在3~6 km深度内影响最大、 库区内外地震频度与视应力均值差异最显著。  相似文献   

10.
An overview of the S-wave velocity (V s) structural model of the Caribbean with a resolution of 2°?×?2° is presented. New tomographic maps of Rayleigh wave group velocity dispersion at periods ranging from 10 to 40?s were obtained as a result of the frequency time analysis of seismic signals of more than 400 ray-paths in the region. For each cell of 2°?×?2°, group velocity dispersion curves were determined and extended to 150?s by adding data from a larger scale tomographic study (Vdovin et al., Geophys. J. Int 136:324–340, 1999). Using, as independent a priori information, the available geological and geophysical data of the region, each dispersion curve has been inverted by the “hedgehog” non-linear procedure (Valyus, Determining seismic profiles from a set of observations (in Russian), Vychislitielnaya Seismologiya 4, 3–14. English translation: Computational Seismology (V.I. Keylis-Borok, ed.) 4:114–118, 1968), in order to compute a set of V s versus depth models up to 300?km of depth. Because of the non-uniqueness of the solutions for each cell, a local smoothness optimization has been applied to the whole region in order to choose a three-dimensional model of V s, satisfying this way the Occam's razor concept. Several known and some new main features of the Caribbean lithosphere and asthenosphere are shown on these models such as: the west directed subduction zone of the eastern Caribbean region with a clear mantle wedge between the Caribbean lithosphere and the subducted slab; the complex and asymmetric behavior of the crustal and lithospheric thickness in the Cayman ridge; the predominant oceanic crust in the region; the presence of continental type crust in Central America, and the South and North America plates; as well as the fact that the bottom of the upper asthenosphere gets shallower going from west to east.  相似文献   

11.
Intensity attenuation for active crustal regions   总被引：1，自引：0，他引：1  
We develop globally applicable macroseismic intensity prediction equations (IPEs) for earthquakes of moment magnitude M W 5.0?C7.9 and intensities of degree II and greater for distances less than 300?km for active crustal regions. The IPEs are developed for two distance metrics: closest distance to rupture (R rup) and hypocentral distance (R hyp). The key objective for developing the model based on hypocentral distance??in addition to more rigorous and standard measure R rup??is to provide an IPE which can be used in near real-time earthquake response systems for earthquakes anywhere in the world, where information regarding the rupture dimensions of a fault may not be known in the immediate aftermath of the event. We observe that our models, particularly the model for the R rup distance metric, generally have low median residuals with magnitude and distance. In particular, we address whether the direct use of IPEs leads to a reduction in overall uncertainties when compared with methods which use a combination of ground-motion prediction equations and ground motion to intensity conversion equations. Finally, using topographic gradient as a proxy and median model predictions, we derive intensity-based site amplification factors. These factors lead to a small reduction of residuals at shallow gradients at strong shaking levels. However, the overall effect on total median residuals is relatively small. This is in part due to the observation that the median site condition for intensity observations used to develop these IPEs is approximately near the National Earthquake Hazard Reduction Program CD site-class boundary.  相似文献   

12.
The earthquake risk on Romania is one of the highest in Europe, and seismic hazard for almost half of the territory of Romania is determined by the Vrancea seismic region, which is situated beneath the southern Carpathian Arc. The region is characterized by a high rate of occurrence of large earthquakes in a narrow focal volume at depth from 70 to 160 km. Besides the Vrancea area, several zones of shallow seismicity located within and outside the Romanian territory are considered as seismically dangerous. We present the results of probabilistic seismic hazard analysis, which implemented the “logic tree” approach, and which considered both the intermediate-depth and the shallow seismicity. Various available models of seismicity and ground-motion attenuation were used as the alternative variants. Seismic hazard in terms of macroseismic intensities, peak ground acceleration, and response spectra was evaluated for various return periods. Sensitivity study was performed to analyze the impact of variation of input parameters on the hazard results. The uncertainty on hazard estimates may be reduced by better understanding of parameters of the Vrancea source zone and the zones of crustal seismicity. Reduction of uncertainty associated with the ground-motion models is also very important issue for Romania.  相似文献   

13.
The Seismic Hazard Harmonization in Europe (SHARE) project, which began in June 2009, aims at establishing new standards for probabilistic seismic hazard assessment in the Euro-Mediterranean region. In this context, a logic tree for ground-motion prediction in Europe has been constructed. Ground-motion prediction equations (GMPEs) and weights have been determined so that the logic tree captures epistemic uncertainty in ground-motion prediction for six different tectonic regimes in Europe. Here we present the strategy that we adopted to build such a logic tree. This strategy has the particularity of combining two complementary and independent approaches: expert judgment and data testing. A set of six experts was asked to weight pre-selected GMPEs while the ability of these GMPEs to predict available data was evaluated with the method of Scherbaum et al. (Bull Seismol Soc Am 99:3234?C3247, 2009). Results of both approaches were taken into account to commonly select the smallest set of GMPEs to capture the uncertainty in ground-motion prediction in Europe. For stable continental regions, two models, both from eastern North America, have been selected for shields, and three GMPEs from active shallow crustal regions have been added for continental crust. For subduction zones, four models, all non-European, have been chosen. Finally, for active shallow crustal regions, we selected four models, each of them from a different host region but only two of them were kept for long periods. In most cases, a common agreement has been also reached for the weights. In case of divergence, a sensitivity analysis of the weights on the seismic hazard has been conducted, showing that once the GMPEs have been selected, the associated set of weights has a smaller influence on the hazard.  相似文献   

14.
Subduction earthquakes along the Pacific Coast of Central America generate considerable seismic risk in the region. The quantification of the hazard due to these events requires the development of appropriate ground-motion prediction equations, for which purpose a database of recordings from subduction events in the region is indispensable. This paper describes the compilation of a comprehensive database of strong ground-motion recordings obtained during subduction-zone events in Central America, focusing on the region from 8 to 14° N and 83 to 92° W, including Guatemala, El Salvador, Nicaragua and Costa Rica. More than 400 accelerograms recorded by the networks operating across Central America during the last decades have been added to data collected by NORSAR in two regional projects for the reduction of natural disasters. The final database consists of 554 triaxial ground-motion recordings from events of moment magnitudes between 5.0 and 7.7, including 22 interface and 58 intraslab-type events for the time period 1976–2006. Although the database presented in this study is not sufficiently complete in terms of magnitude–distance distribution to serve as a basis for the derivation of predictive equations for interface and intraslab events in Central America, it considerably expands the Central American subduction data compiled in previous studies and used in early ground-motion modelling studies for subduction events in this region. Additionally, the compiled database will allow the assessment of the existing predictive models for subduction-type events in terms of their applicability for the Central American region, which is essential for an adequate estimation of the hazard due to subduction earthquakes in this region.  相似文献   

15.
We present the regional ground-motion prediction equations for peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration (PSA), and seismic intensity (MSK scale) for the Vrancea intermediate depth earthquakes (SE-Carpathians) and territory of Romania. The prediction equations were constructed using the stochastic technique on the basis of the regional Fourier amplitude spectrum (FAS) source scaling and attenuation models and the generalised site amplification functions. Values of considered ground motion parameters are given as the functions of earthquake magnitude, depth and epicentral distance. The developed ground-motion models were tested and calibrated using the available data from the large Vrancea earthquakes. We suggest to use the presented equations for the rapid estimation of seismic effect after strong earthquakes (Shakemap generation) and seismic hazard assessment, both deterministic and probabilistic approaches.  相似文献   

16.
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.  相似文献   

17.
We derive S-wave attenuation characteristics, earthquake source parameters and site amplification functions at seismic stations used for earthquake early warning in the Irpinia–Basilicata region, using non-parametric spectral inversion of seismograms from 49 local events with M L = 1.5–3.1. We obtain relatively low Q values (Q 0 = 28 at a frequency of 1 Hz) in conjunction with a strong frequency-dependence (close to linear). The source spectra can be satisfactorily modeled using the omega-square model, with stress drops ranging between 0.01–2 MPa, and in the narrow magnitude range available for analysis, the source spectra seem to scale self-similarly. The local magnitude M L shows a linear correlation with moment magnitude M W, however with a systematic underestimation by about 0.5-magnitude units. The results obtained in this work provide important insights into the ground-motion characteristics that are required for appropriate seismic hazard assessment and are of practical relevance for a suite of applications, such as the calibration of ground-motion prediction equations or the correction for site amplification in earthquake early warning and rapid calculation of shake-maps for seismic emergency management.  相似文献   

18.
We document strong seismic scattering from around the top of the mantle Transition Zone in all available high resolution explosion seismic profiles from Siberia and North America. This seismic reflectivity from around the 410 km discontinuity indicates the presence of pronounced heterogeneity in the depth interval between 320 and 450 km in the Earth’s mantle. We model the seismic observations by heterogeneity in the form of random seismic scatterers with typical scale lengths of kilometre size (10-40 km by 2-10 km) in a 100-140 km thick depth interval. The observed heterogeneity may be explained by changes in the depths to the α-β-γ spinel transformations caused by an unexpectedly high iron content at the top of the mantle Transition Zone. The phase transformation of pyroxenes into the garnet mineral majorite probably also contributes to the reflectivity, mainly below a depth of 400 km, whereas we find it unlikely that the presence of water or partial melt is the main cause of the observed strong seismic reflectivity. Subducted oceanic slabs that equilibrated at the top of the Transition Zone may also contribute to the observed reflectivity. If this is the main cause of the reflectivity, a substantial amount of young oceanic lithosphere has been subducted under Siberia and North America during their geologic evolution. Subducted slabs may have initiated metamorphic reactions in the original mantle rocks.  相似文献   

19.
Aftershock locations, source parameters and slip distribution in the coupling zone between the overriding North American and subducted Rivera and Cocos plates were calculated for the 22 January 2003 Tecomán earthquake. Aftershock locations lie north of the El Gordo Graben with a northwest-southeast trend along the coast and superimposed on the rupture areas of the 1932 (M w?=?8.2) and 1995 (M w?=?8.0) earthquakes. The Tecomán earthquake ruptured the northwest sector of the Colima gap, however, half of the gap remains unbroken. The aftershock area has a rectangular shape of 42?±?2 by 56?±?2?km with a shallow dip of roughly 12° of the Wadati-Benioff zone. Fault geometry calculated with the Náb??lek (1984) inversion procedure is: (strike, dip, rake)?=?(277°, 27°, 78°). From the teleseimic body wave spectra and assuming a circular fault model, we estimated source duration of 20?±?2?s, a stress drop of 5.4?±?2.5?MPa and a seismic moment of 2.7?±?.7?×?1020?Nm. The spatial slip distribution on the fault plane was estimated using new additional near field strong motion data (54?km from the epicenter). We confirm their main conclusions, however we found four zones of seismic moment release clearly separated. One of them, not well defined before, is located toward the coast down dip. This observation is the result of adding new data in the inversion. We calculated a maximum slip of 3.2?m, a source duration of 30?s and a seismic moment of 1.88?×?1020?Nm.  相似文献   

20.
?—?In the last ten years (1990–1999), 21 discrete variations of continuous tilt signal have been recorded on Mount Etna, among which one episode was caused by the opening of the eruptive fracture. The remaining 20 anomalies can be classified into two categories: the first comprises 5 “instantaneous” tilt variations recorded in correspondence to the most energetic seismic events (M L ?≥?3.3) localized on the high western part of the volcano; the second consists of 15 transient anomalies ranging from some hours to 1–2 days, observed at different times at the various tilt stations, with no correlation to seismic events or other evident volcanic episodes. The aseismic variations propagate through the volcanic edifice with a velocity between 4.5–6.0?km/day. Modeling studies suggest that the deformation is generated by a tensile source located 3–6?km SW from Etna volcano summit and 5–10?km depth.  相似文献   

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

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