首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
蒙古及周边地区重力异常和地壳不均匀体分布   总被引:4,自引:1,他引:3       下载免费PDF全文
基于全球EGM2008自由空气重力异常模型,本文计算了蒙古及周边地区的布格重力异常和AiryHeiskanen均衡重力异常.在此基础上,本文采用Crust 1.0地壳模型为参考,通过重力正演方法,对蒙古及周边地区不同深度地壳密度结构模型的重力异常进行了计算,并对得到的正演布格重力异常与实际重力异常进行了对比和分析.研究结果表明:蒙古西部杭爱山地区与阿尔泰山地区的构造变形差异性明显,现今均衡重力异常中杭爱山周边没有明显的均衡异常高值区,而阿尔泰山地区西南方向存在均衡重力异常高值分布,分析与新构造运动密切相关;Crust 1.0模型给出的壳幔横向密度不均匀体分布对于计算Moho面起伏引起的重力异常作用明显;Crust 1.0给出的地壳内界面变形可以反映深大活动断裂的深部构造变形.研究结果对于认识蒙古东西部构造特征差异,以及现今西部活动断裂的地球物理场特征具有参考意义,也可以为进一步应用Crust 1.0模型为参考开展三维密度结构反演提供一定帮助.  相似文献   

2.
In one geophysical survey over a salt dome in the Gulf of Mexico, data taken from a single magnetic and seismic shipborne traverse demonstrated significantly the importance of magnetic measurements. The reduced magnetic and seismic profiles were integrated with gravity measurements from an earlier survey demonstrating the value of obtaining simultaneous geophysical data. The traverse was approximately 25 miles in length, 135 miles southeast of Galveston, Texas. The observed gradient averaged approximately 2 gammas per 1000 feet and crosses the southwest corner of a large shallow salt dome called the Way Dome. The recorded magnetic profile over the salt dome shows a negative magnetic anomaly extending for approximately 12 miles with a maximum amplitude of approximately 20 gammas. This negatively anomaly is due to the negative susceptibility contrast between the salt and cap rock that truncates the surrounding parent sediments. The shape of the negative anomaly correlates well with the continous seismic record and clearly delineates the topographical and structural features of the salt dome. The recorded analog magnetic profile was replotted, with both the regional gradient and time variations removed, which makes the observed anomaly much sharper than is noticeable on the original recorded strip chart record. A theoretical magnetic anomaly was calculated using a computer program for comparison to the observed anomaly. The results showed good correlation.  相似文献   

3.
 Lascar Volcano (5592 m; 23°22'S, 67°44'W) entered a new period of vigorous activity in 1984, culminating in a major explosive eruption in April 1993. Activity since 1984 has been characterised by cyclic behaviour with recognition of four cycles up to the end of 1993. In each cycle a lava dome is extruded in the active crater, accompanied by vigorous degassing through high-temperature, high-velocity fumaroles distributed on and around the dome. The fumaroles are the source of a sustained steam plume above the volcano. The dome then subsides back into the conduit. During the subsidence phase the velocity and gas output of the fumaroles decrease, and the cycle is completed by violent explosive activity. Subsidence of both the dome and the crater floor is accommodated by movement on concentric, cylindrical or inward-dipping conical fractures. The observations are consistent with a model in which gas loss from the dome is progressively inhibited during a cycle and gas pressure increases within and below the lava dome, triggering a large explosive eruption. Factors that can lead to a decrease in gas loss include a decrease in magma permeability by foam collapse, reduction in permeability due to precipitation of hydrothermal minerals in the pores and fractures within the dome and in country rock surrounding the conduit, and closure of open fractures during subsidence of the dome and crater floor. Dome subsidence may be a consequence of reduction in magma porosity (foam collapse) as degassing occurs and pressurisation develops as the permeability of the dome and conduit system decreases. Superimposed upon this activity are small explosive events of shallow origin. These we interpret as subsidence events on the concentric fractures leading to short-term pressure increases just below the crater floor. Received: 12 December 1996 / Accepted: 6 May 1997  相似文献   

4.
Rb-Sr whole-rock analyses yield a Cambro-Ordovician (495 ± 11 m.y.) sedimentation age for the supposed Precambrian Greenland Group and a late Precambrian age, 680 ± 21 m.y., for parts of the Constant Gneiss, the first confirmation of Precambrian rocks in New Zealand. A Precambrian age for the Greenland Group is thus unlikely and the large area of Upper Cambrian-Lower Ordovician rocks now established can be considered as a lateral equivalent of the fossiliferous Lower Palaeozoic succession of northwest Nelson to the east. The Greenland Group, especially in the Paparoa Range has been affected subsequently by a thermal metamorphic overprint about 360 m.y. ago during the Tuhuan Orogeny. Although the Constant Gneiss must form the local basement to the Greenland Group in north Westland, the former does not appear to be the source of the sediments and the true provenance must lie elsewhere.  相似文献   

5.
The aim of the study was — besides the dating of metamorphic events — to evaluate the effects of multi-stage crystal growth, episodic and continuous Pb loss, and U gain on the discordant age patterns found for zircon populations of the polymetamorphic Baltimore Gneiss, the Precambrian basement in the Maryland Piedmont. Eight gneiss and migmatite samples were collected at two localities in the Phoenix and Towson dome, respectively. Their zircon populations were separated into twenty-three fractions of different size and optical appearance. A low-contamination method (T.E. Krogh, 1973) was used for the U-Pb analyses.Microscopy and electron-microprobe studies revealed internal heterogeneities of the zircon crystals: at least half of the grains of each population reflect more than one stage of crystal growth, with the last stage consisting of U-poor overgrowths (U: below 400 ppm, mostly below 200 ppm). Evidence for episodic U gain and overgrown material other than zircon has not been found. On a concordia diagram the “ages” obtained by upward extrapolations (1080 and 1180 m.y.) and downward extrapolations (421 and 455 m.y.) of the best-fit lines to the data points are in fair agreement with the geochronologic data found by other investigators and the probable times of metamorphic periods of Precambrian (Grenville) and early Paleozoic (Taconic) orogenies.Models of Pb loss by continuous diffusion cannot adequately explain the discordant age patterns: these are essentially the result of superposition of episodic Pb loss and zircon overgrowth during the Taconic (and Acadian?) metamorphisms. The zircon overgrowth appears to be present in all fractions, but its influence on the U-Pb systematics is generally not perceptible because it is overridden by the effect of episodic Pb loss. For the fractions showing the most discordant ages, the contribution of Pb loss to the discordancy was found to be at least 85 %.From the microscopic picture and the isotopic data, it appears that the bulk of the zircon substance crystallized during one or several high-grade metamorphisms accompanied by migmatization and granitization of the rocks in the course of the Grenville orogeny. Under consideration of zircon ages of Baltimore Gneiss rocks of Pennsylvania, the results point to a complex Grenville metamorphic history in the Maryland and Pennsylvania Piedmont, that lasted from at least 1200 m.y. until about 980 m.y. The granulite-facies metamorphism in the West Chester Prong, Pennsylvania, may be 50–200 m.y. younger than the metamorphic events in the gneiss domes of the Baltimore area. Although it seems that real differences exist with respect to the Precambrian ages of major zircon-forming events between the Phoenix and the Towson dome, the apparent difference of about 100 m.y. should be interpreted with caution, because it is impossible, so far, to evaluate quantitatively the influence of possibly much older inherited zircon components.  相似文献   

6.
A 109.91 m ice core was recovered from Dome A (or Dome Argus), the highest ice feature in Antarctica, during the 2004/05 austral summer by the 21st Chinese National Antarctic Research Expedition (CHINARE-21). Both methane profile along the core and firn densification model calculation suggest that the close-off depth is at about 102.0 m with an ice age about 4200 a. Stable isotopes (δ18O and δD) of the chips samples produced during each run of ice core drilling at Dome A, together with those of the other co...  相似文献   

7.
Hosted within the Pacquet Harbour Group (PHG) on the Baie Verte Peninsula of north-central Newfoundland, the Rambler rhyolite is a 487 Ma unit of felsic tuffs, flows and subvolcanic intrusive rocks. The PHG has been affected by multiple phases of deformation with the youngest D4 deformation event producing broad northeast plunging upright cross folds in the Rambler rhyolite. Fold culminations on the upper bounding surface of the rhyolite host Cu +/− Au volcanogenic massive sulfide deposits (e.g. Rambler and Ming mines). Geophysical inversions of recently acquired high resolution gravity and magnetic data have been implemented to determine the extent of the fold axis (dome) at depth. To direct the outcome of the inversion process towards a more geologically reasonable solution this study outlines a procedure which permits the inclusion of known geological and geophysical constraints into the input (reference) model for inversion using the MAG3D and GRAV3D algorithms provided by the University of British Columbia Geophysical Inversion Facility. Reference model constraints included surficial geological contacts as defined by aeromagnetic data, and subsurface distribution of physical property variations from a series of drill-hole logs. The output (computed) model images the surface of the rhyolite dome as dipping roughly 40° to the northeast as a series of voxels with density values ranging from 2.71 to 2.75 g/cm3. While previously published ore deposit models parallel this structure in the near surface, results from these inversions suggest deeper exploration may be favorable. Magnetic inversion modeling has not provided any insight into dome morphology however it outlines the distribution of gabbroic dykes surrounding the dome.  相似文献   

8.
Santiaguito volcano has shown a continuous slow extrusion of dacite lava since 1922. In the 50 years of activity there have been four periods of abnormally high extrusion rates, interspersed by periods of little magma production. The type of activity shown by the volcano has been varied and crudely cyclic. Dome extrusion periods are accompanied by pyroclastic activity and followed by lava flows. There are now 16 time stratigraphic units delineated on the dome. Activity since 1967 has been especially closely observed. Dome extrusion at the west end of the complex has been accompanied by pyroclastic cruptions and plug dome extrusion at the east end. The eurrent extrusion rate has remained essentially constant since 1967 at about 5×106 m3/yr, far below Santiaguito’s 1922–71 average of 14×106 m3/yr. The active vent at the east end of the volcano (Caliente vent) has been the principal vent of the volcano since the creation of the explosion crater in 1902. After its initial period of dome extrusion (1922–25), the Caliente vent has chiefly produced pyroclastic eruptions as well as at least 95% of the dome’s lumarolic activity, while lateral vents have continued to give rise to lavas. Lava flows at Santiaguito have effective viscosity values of about 106 poises, while dome lavas are significantly more viscous. The differences in viscosity are in part related to volatile content of the lava when it reaches the surface. During dome extrusion, lavas lose their volatiles through pyroclastic activity before they reach the surface. Lava flows at Santiaguito occur when lava reaches the surface with higher volatile content. Obstruction of either the central (pyroclastic) vent or the lateral (dome extrusion) vent or both vents has an important influence on succeeding activity. In June 1972, at the time of this writing, the outbreak of new lava flows at both the Caliente and lateral El Brujo vents has just occurred, resulting from obstruction of pyroclastic activity by a large plug dome at the Caliente vent.  相似文献   

9.
Numerical models of the extrusion and gravity current of a viscous incompressible fluid are studied to determine the shapes of salt structures formed on the Earth’s surface and the velocities of rock salt extrusion and lateral spreading. Two main types of salt extrusion are examined. In the case of active extrusion, salt rises to the surface at a velocity of about 30–35 cm/yr, forming a salt dome about 550 m in height. In the case of passive extrusion, the velocity of salt extrusion from under a newly developing sedimentary minibasin is significantly lower. In the course of its evolution, the salt dome becomes topographically lower and transforms into a broad plateau. The extrusion velocity of salt controls the shape, size, and velocity of its gravity current. The shapes of salt domes modeled in this study agree well with observations. The gravity current velocities in the models vary from 3 m/yr to 60 cm/yr, depending on the proximity to the current orifice. Numerical modeling of salt extrusion and gravity current in various geodynamic settings can be used for a detailed analysis of the geological and geophysical evolution of structures containing salt layers and the related oil and gas fields.  相似文献   

10.
We describe a low-cost application of digital photogrammetry using commercially available photogrammetric software and oblique photographs taken with an off-the-shelf digital camera to create sequential digital elevation models (DEMs) of a lava dome that grew during the 2004–2008 eruption of Mount St. Helens (MSH) volcano. Renewed activity at MSH provided an opportunity to devise and test this method, because it could be validated against other observations of this well-monitored volcano. The datasets consist of oblique aerial photographs (snapshots) taken from a helicopter using a digital single-lens reflex camera. Twelve sets of overlapping digital images of the dome taken during 2004–2007 were used to produce DEMs and to calculate lava dome volumes and extrusion rates. Analyses of the digital images were carried out using photogrammetric software to produce three-dimensional coordinates of points identified in multiple photos. The evolving morphology of the dome was modeled by comparing successive DEMs. Results were validated by comparison to volume measurements derived from traditional vertical photogrammetric surveys by the US Geological Survey Cascades Volcano Observatory. Our technique was significantly less expensive and required less time than traditional vertical photogrammetric techniques; yet, it consistently yielded volume estimates within 5% of the traditional method. This technique provides an inexpensive, rapid assessment tool for tracking lava dome growth or other topographic changes at restless volcanoes.  相似文献   

11.
The 14 ka Puketarata eruption of Maroa caldera in Taupo Volcanic Zone was a dome-related event in which the bulk of the 0.25 km3 of eruption products were emplaced as phreatomagmatic fall and surge deposits. A rhyolitic dike encountered shallow groundwater during emplacement along a NE-trending normal fault, leading to shallow-seated explosions characterised by low to moderate water/magma ratios. The eruption products consist of two lava domes, a proximal tuff ring, three phreatic collapse craters, and a widespread fall deposit. The pyroclastic deposits contain dominantly dense juvenile clasts and few foreign lithics, and relate to very shallow-level disruption of the growing dome and its feeder dike with relatively little involvement of country rock. The distal fall deposit, representing 88% of the eruption products is, despite its uniform appearance and apparently subplinian dispersal, a composite feature equivalent to numerous discrete proximal phreatomagmatic lapilli fall layers, each deposited from a short-lived eruption column. The Puketarata products are subdivided into four units related to successive phases of:(A) shallow lava intrusion and initial dome growth; (B) rapid growth and destruction of dome lobes; (C) slower, sustained dome growth and restriction of explosive disruption to the dome margins; and (D) post-dome withdrawal of magma and crater-collapse. Phase D was phreatic, phases A and C had moderate water: magma ratios, and phase B a low water: magma ratio. Dome extrusion was most rapid during phase B, but so was destruction, and hence dome growth was largely accomplished during phase C. The Puketarata eruption illustrates how vent geometry and the presence of groundwater may control the style of silicic volcanism. Early activity was dominated by these external influences and sustained dome growth only followed after effective exclusion of external water from newly emplaced magma.  相似文献   

12.
Pyroclastic-laden explosive eruptions from Santiaguito Volcano (Guatemala) are vented from the 200-m diameter Caliente Dome summit and result in a superposition of spatially extensive and temporally sustained (tens of seconds to minutes) acoustic sources. A network of infrasonic microphones distributed on various sides of the volcano record distinct waveforms, which are poorly correlated across the network and suggestive of acoustic interference from multiple sources. Presuming the infrasound wavefield is a linear superposition of spatially and temporally distinct sub-events, we introduce a semblance mapping technique to recover the time history of the spatially evolving sources during successive time windows. Coincident high-resolution video footage corroborates that both rapid dome uplift and individual explosive pulses are likely sources of high semblance infrasound that are identifiable during short (2 s) time windows. This study suggests that complex and network-variable infrasound waveforms are produced whenever a volcanic vent source dimension is large compared to the wavelength of the sound being produced. Non-compact infrasound radiators are probably commonplace at silicic volcanic systems, where venting often occurs across a dome surface.  相似文献   

13.
Magmatic intrusions can trigger thin-skinned compression of the adjacent sedimentary cover by three processes: (a) gravity gliding away from the topographic dome resulting from the ascending magma; (b) fluid push from the rear resulting from forceful intrusion then lateral spreading of the magma; and (c) increased loading by volcanic accumulation. The applicability of the first two mechanisms, gravity gliding and fluid push, was tested using dynamically scaled experiments. Model results help to elucidate the kinematics and structural evolution of thrusts and folds formed by such processes and determine which geological parameters control the deformation style. The results show that the presence of a weak layer within the sedimentary overburden is essential to form thrust and fold belts around the intrusion. Experiments suggest that although gravity gliding can dominate the early stages of deformation, most of the deformation is caused by lateral spreading of the fluid magma pushing the adjacent sedimentary rocks. Models also suggest that true laccoliths can form only if the sedimentary section comprises a basal low-viscosity stratum. Comparison with natural examples allows the evaluation of the applicability of experimental results.  相似文献   

14.
West of Boston, Mass., Castle and others (1976) recognized an up to 5km wide, possibly folded, NE-SW trending Burlington Mylonite Zone. We have extended mapping south into Natick and Framington quadrangles, and supplemented it by fixing local directions of tectonic motion, which are more variable than reported by Goldstein (1989). In Natick the mylonite zone is partly migmatized and converted into blastomylonites, forming the lithodemic Rice Gneiss and is intersected by the Dedham Granite dated ca 630 Ma. The granite also invades deformed, folded, and commonly mylonitized Westboro Quartzite. Thus mylonitization, folding, and formation of migmatitic blastomylonites are all earlier than ca 630 Ma, and can collectively be attributed to the main phase of the Avalonian orogeny that in Africa is referred to as the Pan-African I. The sense of movements in the Rice Gneiss is generally sinistral strike-slip with a NE-SW trend of foliation. Other local mylonites have more variable directions of motion.A narrower E-W zone of mylonitization has been recognized by Grimes (M.S. thesis 1993, Boston College) and named the Nobscot Shear Zone. It affects the Milford Granite, also about 630 Ma in age, while similar but narrow shear zones affect other local granites including the Dedham. These zones, dipping steeply north and including the Nobscot, are less intensely mylonitized and are not associated with migmatites. Their age is not known, but since they affect only Precambrian rocks, they are assumed to be late Proterozoic. We attribute these zones to the second stage of the Avalonian or the Pan-African II.The older rocks west of Boston are widely affected by numerous brittle faults. These are all of unknown age, but probably Phanerozoic. The most significant brittle fault in the Burlington area is the mid to late Paleozoic Bloody Bluff Fault. We do not associate large scale mylonitization with that fault, because the mylonites are commonly cut by undeformed or little deformed Siluro-Devonian gabbro-diorites.  相似文献   

15.
Usu volcano (Hokkaido, Japan) is a dacitic volcano, known for its high production rate of lava domes and crypto-domes. It is thus a good target to study processes of volcanic dome evolution (upheaval and/or relaxation). We carried out repeated GPS and microgravity surveys on the three most recent domes of Mt. Usu (1910: Meiji Shinzan; 1943–1945: Showa-Shinzan and 1977–1982: Usu-Shinzan). The repeat period was 1 to 2 months and extended from October 1996 to June 1997. We also compare new data with results from former studies. More than 20 years after the start of Usu-Shinzan dome growth, there is still subsidence at a maximum rate of about 7 to 8 cm/year. The reasons for this subsidence are discussed. Repeated gravity surveys revealed an increase of gravity on the domes (about 60±10 microgal/year for Usu-Shinzan, about 15 microgal at Showa-Shinzan and 10 to 20 microgal for Meiji-shinzan); this gravity increase exceeds that expected due to subsidence. We discuss and interpret the excess gravity change in terms of a density increase in the edifice, caused by a combination of processes (contraction of the edifice, water level change, devesiculisation, cooling and magma intrusion). Quantification of these processes at Usu volcano may help to understand the processes of evolution at domes on other volcanoes such as Merapi (Indonesia), Unzen (Japan) or Montserrat (West Indies).  相似文献   

16.
Taal Volcano is peculiar in its violent explosivity in spite of its low altitude. The surrounding topographies suggest that the origin of Taal Volcano is either a caldera or a graben structure. To confirm the caldera origin one must search for a vast quantity of caldera ejecta balancing with the depression. As a first step, a gravity survey was carried out on and around Taal Volcano, and high gravity anomalies were lound on Volcano Island. The distribution of the gravities may suggest a graben structure overlying a denser layer of igneous material.  相似文献   

17.
 On King George Island during latest Oligocene/earliest Miocene time, submarine eruptions resulted in the emplacement of a small (ca. 500 m estimated original diameter) basalt lava dome at Low Head. The dome contains a central mass of columnar rock enveloped by fractured basalt and basalt breccia. The breccia is crystalline and is a joint-block deposit (lithic orthobreccia) interpreted as an unusually thick dome carapace breccia cogenetic with the columnar rock. It was formed in situ by a combination of intense dilation, fracturing and shattering caused by natural hydrofracturing during initial dome effusion and subsequent endogenous emplacement of further basalt melt, now preserved as the columnar rock. Muddy matrix with dispersed hyaloclastite and microfossils fills fractures and diffuse patches in part of the fractured basalt and breccia lithofacies. The sparse glass-rich clasts formed by cooling-contraction granulation during interaction between chilled basalt crust and surrounding water. Together with muddy sediment, they were injected into the dome by hydrofracturing, local steam fluidisation and likely explosive bulk interaction. The basalt lava was highly crystallised and degassed prior to extrusion. Together with a low effusion temperature and rapid convective heat loss in a submarine setting, these properties significantly affected the magma rheology (increased the viscosity and shear strength) and influenced the final dome-like form of the extrusion. Conversely, high heat retention was favoured by the degassed state of the magma (minimal undercooling), a thick breccia carapace and viscous shear heating, which helped to sustain magmatic (eruption) temperatures and enhanced the mobility of the flow. Received: 1 August 1996 / Accepted: 15 September 1997  相似文献   

18.
Detailed gravity measurements integrated with geological data were computed to constrain the mechanisms that were active during the emplacement of the Triassic evaporite-bearing folds of Jebel Cheid from the salt-dome zone in the Atlassic region. The gravity analysis consists in mapping the contrasting gravity responses: complete Bouguer anomaly, residual anomaly and derivative maps; the main results display a positive amplitude gravity anomaly as the response of Triassic evaporite bodies and important NE–SW-trending features at the boundaries between the Triassic outcrops and their enveloping strata. In contrast with gravity calculations of a salt dome structure usually resulting in negative gravity anomaly models, the Jebel Cheid clearly expresses a positive gravity anomaly; furthermore, this result is supported by synthetic gravity interpretation.  相似文献   

19.
In oil and mineral exploration, gravity gradient tensor data include higher-frequency signals than gravity data, which can be used to delineate small-scale anomalies. However, full-tensor gradiometry (FTG) data are contaminated by high-frequency random noise. The separation of noise from high-frequency signals is one of the most challenging tasks in processing of gravity gradient tensor data. We first derive the Cartesian equations of gravity gradient tensors under the constraint of the Laplace equation and the expression for the gravitational potential, and then we use the Cartesian equations to fit the measured gradient tensor data by using optimal linear inversion and remove the noise from the measured data. Based on model tests, we confirm that not only this method removes the high-frequency random noise but also enhances the weak anomaly signals masked by the noise. Compared with traditional low-pass filtering methods, this method avoids removing noise by sacrificing resolution. Finally, we apply our method to real gravity gradient tensor data acquired by Bell Geospace for the Vinton Dome at the Texas-Louisiana border.  相似文献   

20.
The distribution, stratigraphic relationships and fragmental components of the May 8 and 20, 1902, pyroclastic flows from Mt. Pelée, Martinique, together with eyewitness accounts, suggest the following explanation for those eruptions. The eruptions were vertically directed magmatic (perhaps initiated phreatically), and contained abundant juvenile lithics from congealed magma of the dome and neck. This resulted in a two-part eruption column having (1) a dense, lithic-charged part which collapsed into the crater and flowed out of a pre-existing notch in its side, giving rise to pyrochlastic flows, and (2) a magmatically derived column containing gases, juvenile vitric material and crystals which largely by-passed the neck and dome and escaped into the atmosphere. All of the energy of the flows was apparently focused through the notch. They emerged fully turbulent and flowed down Rivière Blanche. Gravity segregation of large and abundant fragments soon resulted in a dense, high-concentration, poorly fluidized block-and-ash flow confined to the valley, leaving above a fully turbulent, high-energy ash-cloud surge. As the ash-cloud surge moved down the mountain, it continued to expand outward. The process of gravity segregation continued as the ash-cloud surge expanded, resulting in secondary block-and-ash underflows. Toward St. Pierre, the secondary block-and-ash flows developed on a gently sloping upland surface 100 m or more above the valley of Rivière Blanche. The turbulent, fragment-depleted surges above the secondary block-and-ash flows maintained sufficient energy to devastate the landscape outward to about 3000 m, including St. Pierre. The surges refracted around obstacles and in one place, moved up a small valley in a direction opposite to the main flows.  相似文献   

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

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