首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 84 毫秒
1.
Thermal state, rheology and seismicity in the pannonian basin, Hungary   总被引：1，自引：0，他引：1  
On the basis of data on crustal structure and terrestrial heat flow, a 3-D geothermal model for the lithosphere in the Pannonian basin, Hungary, has been calculated. This model, together with information on crustal composition, laboratory data on rock friction, and certain assumptions about fluid conditions and strain-rate levels within the lithosphere, has been used to construct a rheological model of the area.The results obtained show a layered rheological structure where an aseismic part of the crust is “sandwiched” between an upper and a lower seismogenic crustal layers. According to the proposed rheological model, seismic activity in the upper crust may be expected down to depths of 10–12 km, which is confirmed well by the observed depth distribution of seismicity. The model also predicts a lower crustal seismogenic layer down to 20–22 km. Because of infrequent occurrences of deep earthquakes and/or a generally small number of reliable hypocenter depth determinations in the study area, this seismogenic zone is less constrained by observations.The depth of the different rheologic horizons within the crust is governed mainly by thermal conditions. The lower boundary of both seismogenic layers appears isothermal. Brittle-ductile transition in the upper crust coincides with the ˜200 °C isotherm, while in the lower crust it coincides with the ˜ 375 °C isotherm. The lowermost crust and the upper mantle beneath Hungary show ductile behavior, thus the possibility of siesmic activity at these horizons can be excluded.  相似文献   

2.
Aso Volcano experienced a huge pyroclastic eruption 90 thousand years ago, and formed a large caldera (18 km × 25 km). In order to test the hypothesis of a magma body in the mid and lower crust that has been suggested geophysically and geochemically, we investigated seismic velocity discontinuities and velocity structure beneath Aso Caldera using receiver functions and a genetic algorithm inversion. We confirm the existence of the Moho at depths between 30 km and 35 km and a large velocity anomaly should exist in the deep portion of the crust beneath Aso Caldera, from imaging of receiver functions observed only at stations outside the caldera. As a result of a more detailed examination with GA inversion, a low velocity layer is detected at depths between 10 km and 24 km beneath the western part of the caldera. S-wave velocity of the layer is estimated to be 2.0–2.4 km/s. We estimate that the low velocity layer contains at most 15% melt or 30% aqueous fluid. The layer exists near the Conrad and at the same depths as the swarm of the low frequency earthquakes and a compressional and dilatational deformation source which are expected to be caused by fluid movement beneath the middle-eastern part of the caldera. Fluid contained in the layer might be related with huge pyroclastic eruptions of Aso Volcano.  相似文献   

3.
Using aeromagnetic data acquired in the area from the Cerro Prieto geothermal field, we estimated the depth to the Curie point isotherm, interpreted as the base of the magnetic sources, following statistical spectral-based techniques. According to our results the Curie point isotherm is located at a depths ranging from 14 to 17 km. Our result is somewhat deeper than that obtained previously based only in 2-D and 3-D forward modeling of previous low-quality data. However, our results are supported by independent information comprising geothermal gradients, seismicity distribution in the crust, and gravity determined crustal thickness. Our results imply a high thermal gradient (ranging between 33 and 38 °C/km) and high heat flow (of about 100 mW/m2) for the study area. The thermal regime for the area is inferred to be similar to that from the Salton trough.  相似文献   

4.
In this study, we aim to map the Curie point depth surface for the northern Red Sea rift region and its surroundings based on the spectral analysis of aeromagnetic data. Spectral analysis technique was used to estimate the boundaries (top and bottom) of the magnetized crust. The Curie point depth (CPD) estimates of the Red Sea rift from 112 overlapping blocks vary from 5 to 20 km. The depths obtained for the bottom of the magnetized crust are assumed to correspond to Curie point depths where the magnetic layer loses its magnetization. Intermediate to deep Curie point depth anomalies (10–16 km) were observed in southern and central Sinai and the Gulf of Suez (intermediate heat flow) due to the uplifted basement rocks. The shallowest CPD of 5 km (associated with very high heat flow, ~235 mW m?2) is located at/around the axial trough of the Red Sea rift region especially at Brothers Island and Conrad Deep due to its association with both the concentration of rifting to the axial depression and the magmatic activity, whereas, beneath the Gulf of Aqaba, three Curie point depth anomalies belonging to three major basins vary from 10 km in the north to about 14 km in the south (with a mean heat flow of about 85 mW m?2). Moreover, low CPD anomalies (high heat flow) were also observed beneath some localities in the northern part of the Gulf of Suez at Hammam Fraun, at Esna city along River Nile, at west Ras Gharib in the eastern desert and at Safaga along the western shore line of the Red Sea rift. These resulted from deviatoric tensional stresses developing in the lithosphere which contribute to its further extension and may be due to the opening of the Gulf of Suez and/or the Red Sea rift. Furthermore, low CPD (with high heat flow anomaly) was observed in the eastern border of the study area, beneath northern Arabia, due to the quasi-vertical low-velocity anomaly which extends into the lower mantle and may be related to volcanism in northern Arabia. Dense microearthquakes seem to occur in areas where the lateral gradients of the CPD are steep (e.g. entrance of the Gulf of Suez and Brothers Island in the Red Sea). These areas may correspond to the boundaries between high and low thermal regions of the crust. Thus, the variations in the microseismic activity may be closely related to thermal structures of the crust. Indeed, shallow cutoff depths of seismicity can also be found in some geothermal areas (e.g. western area of Safaga city along the Red Sea coastal region and at Esna city along the River Nile). These facts indicate that the changes in the thickness of the seismogenic layer strongly depend on temperature. Generally, the shallow Curie point depth indicates that some regions in our study area are promising regions for further geothermal exploration particularly in some localities along the River Nile, Red Sea and Gulf of Suez coastal regions.  相似文献   

5.
Ground and aeromagnetic data are combined to characterize the onshore and offshore magnetic properties of the central Philippines, whose tectonic setting is complicated by opposing subduction zones, large-scale strike-slip faulting and arc–continent collision. The striking difference between the magnetic signatures of the islands with established continental affinity and those of the islands belonging to the island arc terrane is observed. Negative magnetic anomalies are registered over the continental terrane, while positive magnetic anomalies are observed over the Philippine Mobile Belt. Several linear features in the magnetic anomaly map coincide with the trace of the Philippine Fault and its splays. Power spectral analysis of the magnetic data reveals that the Curie depth across the central Philippines varies. The deepest point of the magnetic crust is beneath Mindoro Island at 32 km. The Curie surface shallows toward the east: the Curie surface is 21 km deep between the islands of Sibuyan and Masbate, and 18 km deep at the junction of Buruanga Peninsula and Panay Island. The shallowest Curie surface (18 km) coincides with the boundary of the arc–continent collision, signifying the obduction of mantle rocks over the continental basement. Comparison of the calculated Curie depth with recent crustal thickness models reveals the same eastwards thinning trend and range of depths. The coincidence of the magnetic boundary and the density boundary may support the existence of a compositional boundary that reflects the crust–mantle interface.  相似文献   

6.
Two-dimensional crustal velocity models are derived from passive seismic observations for the Archean Karelian bedrock of north-eastern Finland. In addition, an updated Moho depth map is constructed by integrating the results of this study with previous data sets. The structural models image a typical three-layer Archean crust, with thickness varying between 40 and 52 km. P wave velocities within the 12–20 km thick upper crust range from 6.1 to 6.4 km/s. The relatively high velocities are related to layered mafic intrusive and volcanic rocks. The middle crust is a fairly homogeneous layer associated with velocities of 6.5–6.8 km/s. The boundary between middle and lower crust is located at depths between 28 and 38 km. The thickness of the lower crust increases from 5–15 km in the Archean part to 15–22 km in the Archean–Proterozoic transition zone. In the lower crust and uppermost mantle, P wave velocities vary between 6.9–7.3 km/s and 7.9–8.2 km/s. The average Vp/Vs ratio increases from 1.71 in the upper crust to 1.76 in the lower crust.The crust attains its maximum thickness in the south-east, where the Archean crust is both over- and underthrust by the Proterozoic crust. A crustal depression bulging out from that zone to the N–NE towards Kuusamo is linked to a collision between major Archean blocks. Further north, crustal thickening under the Salla and Kittilä greenstone belts is tentatively associated with a NW–SE-oriented collision zone or major shear zone. Elevated Moho beneath the Pudasjärvi block is primarily explained with rift-related extension and crustal thinning at ∼2.4–2.1 Ga.The new crustal velocity models and synthetic waveform modelling are used to outline the thickness of the seismogenic layer beneath the temporary Kuusamo seismic network. Lack of seismic activity within the mafic high-velocity body in the uppermost 8 km of crust and relative abundance of mid-crustal, i.e., 14–30 km deep earthquakes are characteristic features of the Kuusamo seismicity. The upper limit of seismicity is attributed to the excess of strong mafic material in the uppermost crust. Comparison with the rheological profiles of the lithosphere, calculated at nearby locations, indicates that the base of the seismogenic layer correlates best with the onset of brittle to ductile transition at about 30 km depth.We found no evidence on microearthquake activity in the lower crust beneath the Archean Karelian craton. However, a data set of relatively well-constrained events extracted from the regional earthquake catalogue implies a deeper cut-off depth for earthquakes in the Norrbotten tectonic province of northern Sweden.  相似文献   

7.
Curie Point Depth Estimates and Correlation with Subduction in Mexico   总被引：2，自引：0，他引：2  
We investigate the regional thermal structure of the crust in Mexico using Curie Point Depth (CPD) estimates. The top and bottom of the magnetized crust were calculated using the power-density spectra of the total magnetic field from the freely available ??Magnetic Anomaly Map of North America??. We applied this method to estimate the regional crustal thermal structure in overlapping square windows of 2°?×?2°. The CPD estimates range between 10 and 40?km and show several regions of relatively shallow and deep magnetic sources, with a general inverse correlation with measured heat flow. A deep CPD region (20?C30?km) is located in the fore-arc area where the subducting Cocos plate has a flat-slab geometry. This deep region is bound to the NW and SE by shallow CPD areas beneath the states of Michoacan (CPD?=?12?C16?km) and Oaxaca (CPD?=?~16?km), respectively. There is a good spatial correlation between this deep CPD area and two main fracture zones located on the incoming Cocos plate (Orozco and O??Gorman fracture zones), suggesting that subduction plays an important role in setting apart different CPD provinces along the Mexican coast. Another deep CPD (16?C32?km) area corresponds to the region where the Rivera plate subducts beneath Jalisco block. The Trans-Mexican Volcanic Belt is characterized by a decrease in Curie depths from west (16?C20?km) to east (10?C12?km). Finally, several deep CPD areas are situated in the back-arc region where old Mesozoic terrains are present. Our results suggest that the main control on the crust??s regional thermal structure in the fore-arc and volcanic arc regions is due to the subduction of the Cocos and Rivera plates beneath Mexico.  相似文献   

8.
The Cape Verde Islands are located on a mid-plate topographic swell and are thought to have formed above a deep mantle plume. Wide-angle seismic data have been used to determine the crustal and uppermost mantle structure along a ~ 440 km long transect of the archipelago. Modelling shows that ‘normal’ oceanic crust, ~ 7 km in thickness, exists between the islands and is gently flexed due to volcano loading. There is no direct evidence for high density bodies in the lower crust or for an anomalously low density upper mantle. The observed flexure and free-air gravity anomaly can be explained by volcano loading of a plate with an effective elastic thickness of 30 km and a load and infill density of 2600 kg m− 3. The origin of the Cape Verde swell is poorly understood. An elastic thickness of 30 km is expected for the ~ 125 Ma old oceanic lithosphere beneath the islands, suggesting that the observed height of the swell and the elevated heat flow cannot be attributed to thermal reheating of the lithosphere. The lack of evidence for high densities and velocities in the lower crust and low densities and velocities in the upper mantle, suggests that neither a crustal underplate or a depleted swell root are the cause of the shallower than expected bathymetry and that, instead, the swell is supported by dynamic uplift associated with the underlying plume.  相似文献   

9.
芦山与汶川地震之间存在约40 km的地震空区.震源区和地震空区的深部构造背景的研究对深入了解中强地震的深部孕育环境及地震空区的地震活动性具有重要科学意义.利用本小组布设的15个临时观测地震台以及21个芦山科考台站和21个四川省地震局固定台站记录的远震数据,用H-K叠加方法得到各个台站的地壳厚度和平均泊松比,并构建了接收函数共转换点(CCP)偏移叠加图像以及反演得到台站下方的S波速度模型.我们的结果揭示了震源区和地震空区地壳结构特征差异:(1)汶川震源区的地壳平均泊松比为~0.28;芦山震源区为~0.29;而地震空区处于泊松比变化剧烈的区域;(2)汶川地震与芦山地震的震源区以西下方的Moho面呈现深度上的突变(这与前人的研究成果基本一致),分别从~44 km突变到~59 km,~40 km突变到~50 km,而地震空区地壳平均厚度呈现渐变性变化;(3)地震空区Moho面下凹且具有低速的上地壳.综合一维S波速度结构和H-k以及CCP的初步结果,这可能显示汶川地震的发震断裂在深部方向上向西倾斜并形成切割整个地壳的大型断裂;芦山地震则可能是由于上、下地壳解耦引起的;而地震空区处于两种地震形成机制控制区域的过渡带中.  相似文献   

10.
Surface heat flows are calculated from elastic lithosphere thicknesses for the heavy cratered highlands of Mars, in terms of the fraction of the surface heat flow derived from crustal heat sources. Previous heat flow estimations for Mars used linear thermal gradients, which is equivalent to ignoring the existence of heat sources within the crust. We compute surface heat flows following a methodology that relates effective thickness and curvature of an elastic plate with the strength envelope of the lithosphere, and assuming crustal heat sources homogeneously distributed in a radioactive element-rich layer 20 or 60 km thick. The obtained results show that the surface heat flow increases with the proportion of heat sources within the crust, and with the decrease of both radioactive element-rich layer thickness and surface temperature. Also, the results permit us to calculate representative temperatures for the crust base, rock strength for the upper mantle, and lower and upper limits to the crustal magnetization depth and intensity, respectively. For Terra Cimmeria, an effective elastic thickness of 12 km implies between 30% and 80% of heat sources located within the crust. In this case the uppermost mantle would be weak at the time of loading, and temperatures in the lower crust cold enough to favor unrelaxed crustal thickness variations and to permit deep Curie depths in the highlands, as suggested by the observational evidence.  相似文献   

11.
The basal depth of the outer layer with internal magnetic sources was calculated from magnetic data available within a roughly 500 km wide and 1200 km long area, running from central Germany to southern Italy. The dataset, deriving from different aeromagnetic surveys, is reduced to the reference altitude of 3000 m a.s.l. and a reference year of 1980.0. The adopted method, which transforms the spatial data into the frequency domain, provides a relationship between the two-dimensional spectrum of the magnetic anomalies and the top and centroid depths of the magnetic sources. The magnetic layer bottom depth (MLBD) thus obtained is 29-33 km deep in the stable areas (central Europe Variscan units, Corsica-Sardinia Variscan block) and corresponds to the Moho, having an average temperature of 560 °C. From the Alps to the Apennines, MLBD ranges between 22 and 28 km and is clearly shallower than the Moho. In these units, the wide variation of MLBD appears to be compatible with the presence of shallow magnetised bodies, consisting of lower crustal rocks (Ivrea-Verbano zone), ophiolitic units (Penninic zone and Voltri Massif) and intrasedimentary basic volcanic bodies (Po Basin). An average value of 25 km can be attributed to MLBD, which corresponds to a temperature of 550 °C. In the peri-Tyrrhenian zone and the Ligurian Sea, MLBD is below the Moho, which ranges from 17 to 20 km depth, and it has a temperature matching approximately to the Curie temperature of magnetite (580 °C).  相似文献   

12.
Introduction Jiashi-Artux area in southwest Xinjiang is one of the most active earthquake provinces at pre-sent in Chinese mainland. In the last century, about 3/4 strong earthquakes in Chinese mainland hit this area, and especially from January 21 to April 16 in 1997, 7 earthquakes with the magnituderanging from 6.0 to 6.9 occurred in a very small area of 9 km×18 km near Jiashi (ZHU et al, 1998). It has never taken place before in Chinese mainland that a series of strong earthquakes shoo…  相似文献   

13.
Study on S wave velocity structure beneath part stations in Shanxi Province   总被引：2，自引：0，他引：2  
Introduction Shanxi uplift area is bounded by Ordos block to the west, Taihang Mountain front fault zone to the east. On the north side, it is also adjoined with Yinshan fold belt and Hetao-Baotou fault depression basin accompanied with its southern active ruptures, as well as on the south side, connected to Weihe fault depression basin accompanied with large ruptures on the northern margin of Qinling fold belt. The interior of the Shanxi uplift area, with complicated tectonic structures, …  相似文献   

14.
The crust and upper mantle structure beneath southeastern China   总被引：1，自引：0，他引：1  
We analyzed teleseismic waveforms recorded by 44 stations in the Fujian and Taiwan provinces of China and obtained 5344 high quality receiver functions. The crustal thickness (H) and average crustal VP/VS ratio (k) beneath every station were estimated using the H–k stacking method. Crustal thicknesses near the Fujian Province range from 28.3 to 32.8 km with an average of 31.1 km, and the corresponding VP/VS ratios vary from 1.70 to 1.84 with a mean of 1.76. From inland to offshore of the Fujian Province, the crustal thicknesses decrease and Poisson's ratios increase. These may indicate decreasing SiO2 and increasing calc-alkaline contents in the crust. The discontinuity structures such as the Moho, subducting slab, the 410- and 660-km discontinuities (hereafter we call them the 410 and the 660) are also studied using common converted point (CCP) stacking of receiver functions. Along two NW–SE lines of central and northern Taiwan, the CCP stacking results show a western dipping structure at depths above 50 km, suggesting that the Philippine Sea plate is probably subducting beneath the Eurasian continent plate near the central and northern Taiwan. The CCP stacking results show sharp and flat 410- and 660-km discontinuities, and the transition zone thickness (TZT) is the same as that of ambient mantle beneath Fujian and Taiwan Strait, but thickens in the east of Taiwan. These results suggest that (1) the subducting Eurasian continent plate is confined to the depths above 410 km beneath Fujian and Taiwan Strait; and (2) the South China Sea slab may reach the transition zone beneath the east of Taiwan.  相似文献   

15.
The seismogenic zone of subduction thrust faults   总被引：13，自引：0，他引：13  
Abstract Subduction thrust faults generate earthquakes over a limited depth range. They are aseismic in their seaward updip portions and landward downdip of a critical point. The seaward shallow aseismic zone, commonly beneath accreted sediments, may be a consequence of unconsolidated sediments, especially stable-sliding smectite clays. Such clays are dehydrated and the fault may become seismogenic where the temperature reaches 100--150°C, that is, at a 5--15 km depth. Two factors may determine the downdip seismogenic limit. For subduction of young hot oceanic lithosphere beneath large accretionary sedimentary prisms and beneath continental crust, the transition to aseismic stable sliding is temperature controlled. The maximum temperature for seismic behavior in crustal rocks is ～ 350°C, regardless of the presence of water. In addition, great earthquake ruptures initiated at less than this temperature may propagate with decreasing slip to where the temperature is ～ 450°C. For subduction beneath thin island arc crust and beneath continental crust in some areas, the forearc mantle is reached by the thrust shallower than the 350°C temperature. The forearc upper mantle probably is aseismic because of stable-sliding serpentinite hydrated by water from the underthrusting oceanic crust and sediments. For many subduction zones the downdip seismogenic width defined by these limits is much less than previously assumed. Within the narrowly defined seismic zone, most of the convergence may occur in earthquakes. Numerical thermal models have been employed to estimate temperatures on the subduction thrust planes of four continental subduction zones. For Cascadia and Southwest Japan where very young and hot plates are subducting, the downdip seismogenic limit on the subduction thrust is thermally controlled and is shallow. For Alaska and most of Chile, the forearc mantle is reached before the critical temperature, and mantle serpentinite provides the limit. In all four regions, the seismogenic zones so defined agree with estimates of the extent of great earthquake rupture, and with the downdip extent of the interseismic locked zone.  相似文献   

16.
南海西沙地块岛屿地震观测和海陆联测初步结果   总被引：6，自引：1，他引：5       下载免费PDF全文
为了研究南海西沙地块下方的地壳结构,我们在岛屿区架设了流动地震台站进行天然地震观测和海底地震仪的人工地震探测.本文利用远震接收函数方法和射线追踪方法对琛航台的远震数据和海陆联测数据进行了初步处理和分析,建立了琛航岛下方的一维横波速度模型以及横穿琛航岛的二维地壳结构模型.琛航岛顶部存在2 km厚的新生代低速沉积层,下地壳...  相似文献   

17.
玉树震区及邻近地区壳幔速度结构特征   总被引：3，自引：0，他引：3       下载免费PDF全文
对获得的玉树震区的人工地震测深资料,利用正演拟合地震测深方法进行处理与计算,获得了玉树震区的地壳速度结构和构造的基本特征。结果表明,不同区域的地壳速度结构在纵向和横向上均具有明显的非均匀性。 地壳呈层状结构,研究区结晶基底界面起伏较大,玉树附近下方较厚达8km左右,在400.0km桩号向北至温泉方向基底结构逐渐减薄为2.5km。基底界面的凹陷、隆起与不同的地质构造单元具有较好的对应关系;地壳厚度由囊谦、玉树附近的72km向南、北逐渐减薄为62km。在200.0～400.0km桩号之间二维速度等值线及壳内界面起伏变化较大,在玉树附近下方莫霍面有一弧形凹陷。  相似文献   

18.
中国中部的地壳及上地幔构造   总被引：1，自引：0，他引：1       下载免费PDF全文
本文利用中国数字地震台网的资料,对昆仑山—柴达木盆地—秦岭—大别山纬向构造带的地壳及上地幔构造进行了研究。沿这一构造带沉积层平均厚度为8km,剪切波速度β=2.62km/s,壳内平均速度较低,为3.5km/s左右,在下地壳还存在低速层。地壳平均厚度为48km,自东而西相差不很大。上地慢低速层在70—90km深处开始出现。这种构造特征有别于地台和高原,具有块体间的过渡性质,表明这是中国地壳块体中的一个独立构造单元。  相似文献   

19.
云南西部地壳深部结构特征   总被引：10，自引：3，他引：7       下载免费PDF全文
在云南西部,穿过红河、小江断裂带完成了一条长360 km、呈北东向的深地震宽角反射/折射剖面.通过对该测线的观测资料进行一维、二维模拟解释,得到了沿剖面的二维地壳速度模型.研究结果显示,沿测线Moho界面埋深横线变化大,其西南侧Moho埋深约35 km,东北侧Moho最大埋深可达43 km.沿剖面从西南到北东方向,地壳平均P波速度从5.9 km/s逐渐增加到6.13 km/s,但显著低于全球大陆平均值.结合以往的接收函数和面波联合反演结果,我们推算沿测线从西南到东北,其下方地壳泊松比介于0.23~0.25之间.剖面西南侧上地壳具有异常低的P波速度和泊松比,暗示其下方上地壳以α-相长英质组分为主;而剖面东北上地壳相对较高的P波速度和泊松比则暗示其物质组成以花岗岩-花岗闪长岩为主.研究区下地壳的P波速度和泊松比分别介于6.25~6.75 km/s和0.24~0.26 km/s之间,暗示其上部组成以花岗岩相的片麻岩为主,而下部组成则以角闪石类岩石为主.红河断裂两侧地壳速度显著不同,从浅到深其速度差异逐渐变弱,但红河断裂两侧地壳厚度变化较大.而小江断裂下方两侧地壳速度和地壳厚度变化并没有红河断裂那么明显.  相似文献   

20.
From an analysis of many seismic profiles across the stable continental regions of North America and northern Europe, the crustal and upper mantle velocity structure is determined. Analysis procedures include ray theory calculations and synthetic seismograms computed using reflectivity techniques. TheP wave velocity structure beneath the Canadian Shield is virtually identical to that beneath the Baltic Shield to a depth of at least 800 km. Two major layers with a total thickness of about 42 km characterize the crust of these shield regions. Features of the upper mantle of these region include velocity discontinuities at depths of about 74 km, 330 km, 430 km and 700 km. A 13 km thickP wave low velocity channel beginning at a depth of about 94 km is also present.A number of problems associated with record section interpretation are identified and a generalized approach to seismic profile analysis using many record sections is described. TheS wave velocity structure beneath the Canadian Shield is derived from constrained surface wave data. The thickness of the lithosphere beneath the Canadian and Baltic Shields is determined to be 95–100 km. The continental plate thickness may be the same as the lithospheric thickness, although available data do not exclude the possibility of the continental plate being thicker than the lithosphere.  相似文献   

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

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