首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 812 毫秒
1.
Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M L ?≤?4.4. Up to five seismograph stations were utilised to determine Q c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q 0 and α were determined for the attenuation relation Q c(f)?=?Q 0 f α . The result was Q 0?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t s???t 0) (time from origin time t 0 to the start of coda analysis window is 1.9 times the S-travel time, t s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q c was not calculated. The derived Q c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q S was Q 0?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q c result.  相似文献   

2.
—SKS phases observed at broadband stations in Germany show significant shear-wave splitting. We have analyzed SKS and SKKS phases for shear-wave splitting from 13 stations of the German Regional Seismic Network (GRSN), from 3 three-component stations of the Gräfenberg array (GRF) and from one Austrian station (SQTA). The data reveal strong differences in the splitting parameters (fast direction φ and delay time δt from a single event at various stations as well as variations at the individual stations for events with different backazimuths. The backazimuthal variations of the splitting parameters at some stations can be explained by two-layer anisotropy models with horizontal symmetry axes. The best resolved two-layer model is the GRA1 model (upper layer φ = 40°, δt = 1.15s; lower layer φ = 115°, δt = 1.95s). The upper layer can be attributed to the lithosphere. Because of the magnitude of the delay time of the upper layer, the lower layer must lie within the asthenosphere. At other stations splitting parameters are consistent with an anisotropic one-layer model for the upper mantle. Stations near the Bohemian Massif show fast directions near EW. Throughout NE Germany the directions are oriented NW/SE. The reason for this direction is probably the nearby Tornquist-Teisseyre line. The observed fast axes are subparallel to this prominent Transeuropean suture zone. At stations in southern Germany near the Alps we observed ENE/WSW directions. Below some stations we also found indications of inclined anisotropic layers.  相似文献   

3.
We present new one-dimensional SH-wave velocity models of the upper mantle beneath the Kalahari craton in southern Africa obtained from waveform inversion of regional seismograms from an Mw = 5.9 earthquake located near Lake Tanganyika recorded on broadband seismic stations deployed during the 1997–1999 Southern African Seismic Experiment. The velocity in the lithosphere beneath the Kalahari craton is similar to that of other shields, and there is little evidence for a significant low velocity zone beneath the lithosphere. The lower part of the lithosphere, from 110 to 220 km depth, is slightly slower than beneath other shields, possibly due to higher temperatures or a decrease in Mg number (Mg#). If the slower velocities are caused by a thermal anomaly, then slightly less than half of the unusually high elevation of the Kalahari craton can be explained by shallow buoyancy from a hot lithosphere. However, a decrease in the Mg# of the lower lithosphere would increase the density and counteract the buoyancy effect of the higher temperatures. We obtain a thickness of 250 ± 30 km for the mantle transition zone, which is similar to the global average, but the velocity gradient between the 410 and 660 km discontinuities is less steep than in global models, such as PREM and IASP91. We also obtain velocity jumps of between 0.16 ± 0.1 and 0.21 ± 0.1 km/s across the 410 km discontinuity. Our results suggest that there may be a thermal or chemical anomaly in the mantle transition zone, or alternatively that the shear wave velocity structure of the transition zone in global reference models needs to be refined. Overall, our seismic models provide little support for an upper mantle source of buoyancy for the unusually high elevation of the Kalahari craton, and hence the southern African portion of the African Superswell.  相似文献   

4.
We used a wavelet formulation of the classical spectral isostatic analysis to invert satellite-derived gravity and topography/bathymetry for elastic thickness (Te) over South America and its surrounding plates. To provide a homogeneous representation of the gravity field for this vast region, we corrected free-air anomalies derived from a combination of terrestrial/marine gravity data with data from the GRACE and CHAMP satellite missions (model EIGEN-CG03C) by a simple Bouguer slab using a smoothed representation of surface relief (wavelengths > 125 km). The resulting Bouguer anomaly compares well with terrestrial data acquired in the Central Andes and allows Te to be confidently estimated for values greater than 10 km. The Te map resolves regional-scale features that are well-correlated with known surface structures and shows maximum values of 100 ± 15 km over the Archean–Neoproterozoic core of the continent, decreasing to less than 30 km around continental margins. Several regions of the oceanic plates and continental margins have an elastic thickness less than 10 km. We performed a quantitative analysis by comparing the elastic thickness with the thermal structure predicted from the age of oceanic crust and igneous–metamorphic rocks. This demonstrates that oceanic plates have been weakened by thermal interaction with hotspots and locally by fracturing and hydration near the trench. We observe that only the nucleus of the continent has resisted the thermomechanical weakening induced by the rifting of Africa and South America along the passive margin and the Andean orogeny along the active margin. This latter region shows along-strike variations in Te that correlate with the geotectonic segmentation of the margin and with the pattern of crustal seismicity. Our results reveal that the rigidity structure follows the segmentation of the seismogenic zone along the subduction fault, suggesting a causal relationship that should be investigated in order to improve the understanding and predictability of great earthquakes and tsunamis.  相似文献   

5.
This paper deals with the interpretation of Bouguer gravity anomalies measured along a 250 km long Suhaitu-Etuokeqi gravity profile located at the transitional zone of the Alxa and Ordos blocks where geophysical characteristics are very complex. The analysis is carried out in terms of the ratio of elevation and Bouguer gravity anomaly, the normalized full gradient of a section of the Bouguer gravity anomaly (G h ) and the crustal density structure reveal that (1) the ratio of highs and lows of elevation and Bouguer gravity anomaly is large between Zhengyiguan fault (F4) and Helandonglu fault (F6), which can be explained due to crustal inhomogeneities related to the uplift of the Qinghai-Tibet block in the northeast; (2) the main active faults correspond to the G h contour strip or cut the local region, and generally show strong deformation characteristics, for example the Bayanwulashan mountain front fault (F1) or the southeast boundary of Alxa block is in accord with the western change belt of G h , a belt about 10 km wide that extends to about 30 km; (3) Yinchuan-Pingluo fault (F8) is the seismogenic structure of the Pingluo M earthquake, and its focal depth is about 15 km; (4) the Moho depth trend and Bouguer gravity anomaly variation indicates that the regional gravity field is strongly correlated with the Moho discontinuity.  相似文献   

6.
We present a 3D model of shear velocity of crust and upper mantle in China and surrounding regions from surface wave tomography.We combine dispersion measurements from ambient noise correlation and traditional earthquake data.The stations include the China National Seismic Network,global networks,and all the available PASSCAL stations in the region over the years.The combined data sets provide excellent data coverage of the region for surface wave measurements from 8 to 120 s,which are used to invert for 3D shear wave velocity structure of the crust and upper mantle down to about150 km.We also derive new models of the study region for crustal thickness and averaged S velocities for upper,mid,and lower crust and the uppermost mantle.The models provide a fundamental data set for understanding continental dynamics and evolution.The tomography results reveal significant features of crust and upper mantle structure,including major basins,Moho depth variation,mantle velocity contrast between eastern and western North China Craton,widespread low-velocity zone in midcrust in much of the Tibetan Plateau,and clear velocity contrasts of the mantle lithosphere between north and southern Tibet with significant E–W variations.The low velocity structure in the upper mantle under north and eastern TP correlates with surface geological boundaries.A patch of high velocity anomaly is found under the eastern part of the TP,which may indicate intact mantle lithosphere.Mantle lithosphere shows striking systematic change from the western to eastern North China Craton.The Tanlu Fault appears to be a major lithosphere boundary.  相似文献   

7.
The aim of this study is to assess the availability and quality of data from the International GNSS Service (IGS) Global Positioning System (GPS) network in Africa, especially for retrieving zenith tropospheric delay (ZTD), from which precipitable water vapour (PWV) can be derived, in view of application to the African Monsoon Multidisciplinary Analysis (AMMA) project. Three major error sources for the GPS data analysis evaluating PWV in Africa are the accuracy of the satellite orbits, the correction for the radio delay induced by the ionosphere and the vertical site displacements due to ocean loading. The first part of this study examines these error sources and the validity of GPS data for meteorological applications in Africa in dedicated analyses spanning the year 2001. These analyses were performed using the IGS precise orbits. Weak degradation of baseline precision with increasing baseline lengths suggests that the average orbital error is not limiting the GPS analysis in Africa. The impact of the ionosphere has been evaluated during a maximum of solar activity in 2001. The loss of L2 data has actually been observed. It amounts to 2% on average for 2001, with maxima of 8% during magnetic storm events. A slight decrease in formal accuracy of ZTD seems to be related to the loss of L2 data at the end of the day. This indicates that scintillation effects are present in the GPS observations but however are not a major limitation. The impact of ocean loading is found to be significant on ZTD estimates (up to ±2 mm in equivalent PWV). The use of a proper ocean loading model eliminates this effect.The second aspect of this study concerns the IGS analysis quality for the African stations. The accuracy has been assessed through position dispersion between individual solutions and the most recent version of the IGS combined solution IGb00, and residuals from the transformation of the IGS combined solution into the International Terrestrial Reference Frame 2005. The positioning performance of the IGS analysis is consistent with an accuracy in ZTD of ±6 mm (±1 mm in PWV), as requested for meteorological applications such as planned in AMMA.  相似文献   

8.
Sharp linear gradients in maps of potential field data are generally assumed to result from sharp discontinuities or boundaries between rocks having different densities or magnetic susceptibilities and are usually associated with faults or other geological contacts. The computation of the horizontal gradients of the gravity field permits us to localize the limits of such blocks and then the fault locations. The horizontal derivative maxima of the Bouguer anomaly and its upward continuation at several heights show lineaments that could reflect the layout of faults and/or contacts and their dip directions. The application of this method to the Bouguer anomaly map of Morocco (with 19,571 points, using an average crustal density ρ = 2.67 g/cm3) allowed us to perform a multiscale analysis of the gravimetric lineaments of the country. The obtained structural map is consistent with several faults already identified in previous studies, and highlights five new major subsurface faults systems with location and dip: the Saghro fault system; Bou-Arfa Midelt fault system; Sidi Slimane Mezquitem fault; Ksar El Kebir–Chefchaouen fault and the Rifan West Mediterranean fault. In addition, this study suggests a new shape and localization for the Agadir-Oujda trans-Moroccan major fault with a NE-SW direction and 900 km length, subdividing Morocco into two main domains. The results of this study contribute to the improvement of the regional structural map of the north western part of Africa, which is situated within the convergence zone between Africa and Eurasia.  相似文献   

9.
Bathymetric and gravity surveys were carried out from 1988 to 1994, in the Gulf of Naples (Southern Italy) to offshore extend the already existing Bouguer anomaly map. In order to improve the knowledge of the structural setting beneath the active Neapolitan volcanoes (Vesuvio, Campi Flegrei and Ischia), 862 stations were surveyed within the isobath of 400 m; at the same time, and about 2000 on-land gravity values were also collected. A new Bouguer anomaly map spanning the whole volcanic region was drawn from the final data set. Gravity anomalies were referred to the new absolute gravity station in Naples and computed according to 1980 Geodetic Reference System. Finally, a density value of 2200 kg/m3 was used in the computation of the Bouguer and terrain effects. We carried out the inversion of the gravity anomalies adopting a 2.5-D modelling along selected profiles crossing the investigated area. The interpretative models were constrained to data obtained from deep wells and other geophysical investigations.  相似文献   

10.
中国东北新生代板内火山广泛发育,其中诺敏河火山由于上地幔结构研究的匮乏,火山成因尚不明确.利用布设在诺敏河火山周围的40个流动台站所记录到的远震剪切波数据,测量得到82对各向异性参数和219个无效分裂结果.结果表明,研究区快慢波延迟时间变化范围为0.4～1.4s,平均0.78±0.21s;各向异性快波方向范围为N77°W—N18°E,绝大多数快波方向集中在N6.9°W±9.87°,平行于中生代晚期岩石圈伸展变形方向,推测由残留在岩石圈中的化石各向异性所引起.同时,在火山中心及周边部分台站,只观测到无效分裂而没有观测到有效分裂结果,可能是由于残存在岩石圈内的古老形变被上涌的热地幔物质所侵蚀.  相似文献   

11.
The western China lies in the convergence zone between Eurasian and Indian plates. It is an ideal place to study the lithosphere dynamics and tectonic evolutions on the continental Earth. The lithospheric strength is a key factor in controlling the lithosphere dynamics and deformations. The effective elastic thickness (Te) of the lithosphere can be used to address the lithospheric strength. Previous researchers only used one of the admittance or coherence methods to investigate the Te in the western China. Moreover, most of them ignored the internal loads of the lithosphere during the Te calculation, which can produce large biases in the Te estimations. To provide more reliable Te estimations, we used a new joint inversion method that integrated both admittance and coherence techniques to compute the Te in this study, with the WGM2012 gravity data, the ETOPO1 topographic data, and the Moho depths from the CRUST1.0 model. The internal loads are considered and investigated using the load ratio (F). Our results show that the joint inversion method can yield reliable Te and F values. Based on the analysis of Te and F distributions, we suggest (1) the northern Tibetan Plateau could be the front edge of the plate collision of Eurasian and Indian plates; (2) the southern and part of central Tibetan Plateau have a strong lithospheric mantle related to the rigid underthrusting Indian plate; (3) the southeastern Tibetan Plateau may be experiencing the delamination of lithosphere and upwelling of asthenosphere.  相似文献   

12.
The effective elastic thickness, Te, has been calculated in the collision zone between Arabia and Eurasia in Iran from the wavelet coherence. The wavelet coherence is calculated from Bouguer anomalies and topography data using the isotropic fan wavelet method, and gives Te values between 14.2 and 62.2 km. The lower value is found in the Central Iranian Blocks and the East Iranian Belt which are bounded by several large strike-slip faults with lithospheric origin. The higher value occurs in the east of the South Caspian Sea Basin. The resulting Te map shows positive and negative correlation with shear wave velocity and surface heat flow, respectively. A comparison between the seismogenic thickness (Ts) and Te in Iran suggests that Te > Ts. Results of the load ratio in Iran indicate that in most of the study area surface loads are much more prevalent than subsurface loads, except in the Central Iranian Blocks and NW of Iran. Intermediate to low Te values in Iran were inherited from multiple rifting and orogenic activities from Late Precambrian (∼650 Ma) to present day which are not only reflected in thin and warm lithosphere but also an increasing seismicity rate.  相似文献   

13.
The variations in the density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift perpendicularly to the magnetic (B) and electric (E) fields in the direction geomagnetic west-geomagnetic east have been studied using the three-dimensional nonstationary theoretical model of electron and ion densities (N e and N i ) and temperatures (T e and T i ) in the low-latitude and midlatitude ionospheric F region and plasmasphere. The method of numerical calculations of N e , N i , T e , and T i , including the advantages of the Lagrangian and Eulerian methods, is used in the model. A dipole approximation of the geomagnetic field (B), taking into account the non-coincidence of the geographic and geomagnetic poles and differences between the positions of the Earth’s and geomagnetic dipole centers, is accepted in the calculations. The calculated NmF2 and altitudes of the F2 layer maximum (hmF2) have been compared with these quantities measured at 16 low-latitude ionospheric sounding stations during the geomagnetically quiet period October 11–12, 1958. This comparison made it possible to correct the input model parameters: the NRLMSISE-00 model [O], the meridional component of the neutral wind velocity according to the HWW90 model, and the meridional component of the equatorial plasma drift due to the electric field specified by the empirical model. It has been indicated that the effect of the zonal E × B plasma drift on NmF2 can be neglected under daytime conditions and changes in NmF2 and hmF2 under the action of this drift are insignificant under nighttime conditions north of 25° and south of ?26° geomagnetic latitude. The effect of the zonal E × B plasma drift on NmF2 and hmF2 is most substantial in the nightside ionosphere approximately from ?20° to 20° geomagnetic latitude, and the neglect of this drift results in an up to 2.4-fold underestimation of NmF2. The found dependence of the effect of the zonal E × B plasma drift on NmF2 and hmF2 on geomagnetic latitude is related to the longitudinal asymmetry of B, asymmetry of the neutral wind about the geomagnetic equator, and changes in the meridional E × B plasma drift at a change in geomagnetic longitude.  相似文献   

14.
The lithosphere is interpreted as a thermal boundary layer. Approximate solutions of the boundary layer cooling problem are developed which include mantle radioactivity, partial melt in the asthenosphere, a temperature gradient in the asthenosphere, and a non-zero lithospheric thickness at the ridge crests. The cooling history of oceanic lithosphere is found to be remarkably insensitive to assumptions about the amount of radioactivity in the upper mantle and the extent of melting in the asthenosphere. Determinations of the thickness of oceanic lithosphere and the depths of oceans as a function of age are in excellent agreement with boundary layer predictions which include a heat flux from the asthenosphere. However, the determinations do not resolve how much of the total asthenospheric heat flux might be caused by a temperature gradient in the asthenosphere. Simple thermal arguments indicate that the initial lithospheric thickness, L0, at ridge crests should depend on the local half-spreading rate, V, as L0 = 3 km/V(cm/year).  相似文献   

15.
I present the results of statistical hypothesis testing of Grand’s (2002) global tomography model of three-dimensional shear velocity variations for the middle mantle underneath eastern and southern Africa. I apply an F test to evaluate the validity of a model where a tilted, slow-velocity anomaly in the deepest mantle under southern Africa, known as the African superplume, is continuous with a slow-velocity anomaly in the upper mantle under eastern Africa. This null hypothesis is tested against alternative hypotheses, in which various “obstruction volumes” in the middle mantle are constrained to zero perturbation from the one-dimensional reference velocity during the tomographic inversion. I find that there is an equal probability of accepting an alternative hypothesis with a thin “obstruction volume” at 850–1,000 km depth, whereas volumes at other depths are rejected. But the alternative hypothesis, where a connection is forced at 850–1,000 km depth, is rejected. I conclude that the African superplume rises to at least 1,150 km depth, and that the upper mantle slow-velocity anomaly continues from the surface to below the mantle transition zone. I interpret the “obstruction volume” as a weakening of the superplume in the middle mantle.  相似文献   

16.
我国海洋地学编图现状、计划与主要进展   总被引：1，自引：1，他引：0       下载免费PDF全文
我国海洋地质地球物理工作起步晚,整体调查程度较低,进一步加强我国海洋地学编图显得十分必要.我国实施的"海洋地质保障工程",将分"中国海陆"、"中国海及邻域"和"中国各海区"三个层次来进行地学编图.第一层次编制了空间重力异常图、布格重力异常图、磁力异常图、地震层析成像图、莫霍面深度图、地质图、大地构造格架图和大地构造格架演化图等8种图件.编图反映出中国海陆重力异常是"线性异常带纵横交错,块状异常区坐落其间"、重力梯级带主要为块体结合带;陆地磁力异常大致正负磁异常相间分布,海域北部为沟-弧-盆相关异常、南部磁条带异常;莫霍面总体特征为"东西分带,南北分块";地震层析成像反映出中国东部和西部岩石层厚度的差异以及上地幔软流层的分布特点;大地构造格架和演化图再现了块体体制向板块体制的转换过程.本次编图强调地球深部结构的变化对表层构造的关联,重视地球物理资料与地质构造的结合,以深-浅层结合来划分中国大地构造格架.已有成果表明,中国海陆大地构造格局可以用"块体构造学说"来描述.  相似文献   

17.
The Snake River Plain (SRP), a 90-km-wide topographic depression in southern Idaho, is a topographically anomalous feature in the western U.S. Previous seismic studies focused on the northeastern SRP to study its relationship with the Yellowstone hotspot. We present new teleseismic shear-wave splitting data from six broadband seismic stations deployed along the axis of the SRP from June 2000 to September 2001. We also analyze splitting at HLID, a permanent station of the National Seismic Network located ∼100 km north of the plain. Splitting of individual teleseismic phases is consistent at all stations within 2σ errors, and we favor the interpretation of anisotropy with a single horizontal fast axis, although a dipping-axis interpretation is statistically permitted at two of the stations. Our station fast directions, as well as shear-wave splitting data from numerous other stations throughout the Basin and Range, are best explained by a lattice preferred orientation of olivine due to horizontal shear along the base of the plate associated with the gravitational spreading of buoyant plume-like upwelling material beneath eastern Nevada into a southwestward flowing asthenosphere (with respect to a fixed hotspot reference frame). This parabolic asthenospheric flow (PAF) model for the Great Basin is attractive because it explains the observed high elevations, high mantle buoyancy, low-velocity anomaly beneath eastern Nevada, high heat flow, and depleted geochemistry of some erupted basalts. The lack of Pliocene-Recent major volcanism in eastern Nevada suggests that a significant amount of the buoyancy flux is due to compositional buoyancy. Our splitting station delay times vary in a way not predicted by the PAF model, and can be explained by: a zone of aligned magma-filled lenses and/or partially molten dikes beneath the SRP lithosphere, a depleted olivine-rich residuum underneath the sides of the eastern SRP, and/or the effect of lateral lower crustal flow from beneath the eSRP toward its adjacent flanks.  相似文献   

18.
The variations in the electron number density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift in the geomagnetic west-geomagnetic east direction perpendicularly to the electric (E) and geomagnetic (B) fields during a geomagnetically quiet period on December 7, 1989, at high solar activity have been studied based on a three-dimensional nonstationary theoretical model of electron number densities and temperatures in the ionospheric F region. Calculated and measured NmF2 values for 12 low-latitude ionospheric sounding stations have been compared. When the zonal E × B plasma drift is ignored, the NmF2 values become smaller by up to a factor of 3 under nighttime conditions in the low-latitude ionosphere. The average effect of the zonal E × B plasma drift on NmF2 in the low-latitude ionosphere is larger during winter nights than under summer nighttime conditions.  相似文献   

19.
利用73个固定台站记录的163个远震事件数据,采用多道互相关技术挑选了5524条S波到时数据,并对S波到时数据进行地壳校正,在此基础上采用天然地震层析成像方法和远震S波到时信息,获得了长江中下游成矿带上地幔的三维S波速度结构模型.研究结果表明:(1)研究区域上地幔存在着明显的低速异常,且走向与成矿带相同,可能为上涌的软流圈热物质;(2)研究区域地幔过渡带和上地幔底部存在着明显的高速异常,可能为俯冲的古太平洋板块和拆沉的岩石圈;(3)成矿带上地幔的低速异常呈现由南向北逐渐变浅的空间分布特征,该特征表明软流圈热物质由南向北上涌.综合分析认为,成矿带中生代大规模岩浆活动和成矿作用的深部过程主要与岩石圈的拆沉密切相关.  相似文献   

20.
Of 16 sites collected in the Taru grits (Permian) and Maji ya Chumvi beds (Permo-Triassic) of East Africa only 6 sites from the Maji ya Chumvi sediments gave meaningful palaeomagnetic results. After thermal cleaning the 6 sites (32 samples) give an Early Triassic pole at 67°N, 269°E with A95 = 17° in excellent agreement with other African Mesozoic poles. There are now 26 Mesozoic palaeomagnetic poles for Africa from widely diverse localities ranging in present latitude from 35°N to 30°S. The poles subdivide into Triassic (17 poles) and Cretaceous (9 poles) groups whose means are not significantly different. The palaeomagnetic pole for Africa thus remained in much the same position for 170 m.y. from Early Triassic to Late Cretaceous. The data form an especially good set for estimating the palaeoradius using Ward's method. Values of 1.08 ± 0.15 and 1.03 ± 0.19 times the present radius are deduced for the Triassic and Cretaceous respectively with a mean value of 1.08 ± 0.13 for all the Mesozoic data combined. The analysis demonstrates that hypotheses of earth expansion are very unattractive.  相似文献   

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

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