首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 671 毫秒
1.
In our previous works, based on numerical models, it was shown that under certain conditions a hot material can rise in portions in the tails of thermal mantle plumes. The spectrum of these pulsations can correspond to the observed spectra of catastrophic hotspot eruptions. Since most of the existing numerical models of thermal convection for the mantle of the present Earth do not reveal these pulsations, in this work, we analyze the physical cause and initiation conditions of pulsations of thermal plumes. The results of a numerical solution of the thermal convection equations for a material with varying parameters in the extended Boussinesq approximation are presented. It is shown how the structure of the convection is transformed with the increase of convection intensity. At the Rayleigh numbers Ra > 106, convection becomes unsteady, and the configuration of the ascending and descending flows changes. The new flow emerging at the mantle bottom acquires a mushroom shape with a head and a tail. After the rise of the plume’s head to the surface, the tail remains in the mantle in the form of a quasi-stationary hot steam. It turns out that at Ra ~ 5 × 107, the thermal mantle plume becomes pulsating and its tail is in fact a heated channel through which the hot material rises in successive portions. At the Rayleigh numbers Ra > 5 × 108, the tail of the thermal plume breaks and the plume becomes a regular conveyor of separate ascending portions of the hot material, which are referred to as thermals. Thus, thermal convection with pulsating plumes takes place at the transitional stage from the regime of quasi-stationary plumes to the regime of thermals.  相似文献   

2.
The coexistence of stationary mantle plumes with plate-scale flow is problematic in geodynamics. We present results from laboratory experiments aimed at understanding the effects of an imposed large-scale circulation on thermal convection at high Rayleigh number (106≤Ra≤109) in a fluid with a temperature-dependent viscosity. In a large tank, a layer of corn syrup is heated from below while being stirred by large-scale flow due to the opposing motions of a pair of conveyor belts immersed in the syrup at the top of the tank. Three regimes are observed, depending on the ratio V of the imposed horizontal flow velocity to the rise velocity of plumes ascending from the hot boundary, and on the ratio λ of the viscosity of the interior fluid to the viscosity of the hottest fluid in contact with the bottom boundary. When V≪1 and λ≥1, large-scale circulation has a negligible effect on convection and the heat flux is due to the formation and rise of randomly spaced plumes. When V>10 and λ>100, plume formation is suppressed entirely, and the heat flux is carried by a sheet-like upwelling located in the center of the tank. At intermediate V, and depending on λ, established plume conduits are advected along the bottom boundary and ascending plumes are focused towards the central upwelling. Heat transfer across the layer occurs through a combination of ascending plumes and large-scale flow. Scaling analyses show that the bottom boundary layer thickness and, in turn, the basal heat flux q depend on the Peclet number, Pe, and λ. When λ>10, q∝Pe1/2 and when λ→1, q∝(Peλ)1/3, consistent with classical scalings. When applied to the Earth, our results suggest that plate-driven mantle flow focuses ascending plumes towards upwellings in the central Pacific and Africa as well as into mid-ocean ridges. Furthermore, plumes may be captured by strong upwelling flow beneath fast-spreading ridges. This behavior may explain why hotspots are more abundant near slow-spreading ridges than fast-spreading ridges and may also explain some observed variations of mid-ocean ridge basalt (MORB) geochemistry with spreading rate. Moreover, our results suggest that a potentially significant fraction of the core heat flux is due to plumes that are drawn into upwelling flows beneath ridges and not observed as hotspots.  相似文献   

3.
本文分别在直角坐标系和柱坐标系下,研究瑞利数从104逐渐增大到107对热对流的影响,数值计算结果表明:瑞利数越大,地幔柱越窄,地幔柱上升速度也越快;源自上地幔的地幔柱半径的范围为9到210 km.根据峨眉山内带的半径推算出地幔的黏性系数约为3.8×1021Pa·s,地幔柱平均流动速度为2.5 cm/a.  相似文献   

4.
The pattern and style of mantle convection govern the thermal evolution, internal dynamics, and large-scale surface deformation of the terrestrial planets. In order to characterize the nature of heat transport and convective behaviour at Rayleigh numbers, Ra, appropriate for planetary mantles (between 104 and 108), we perform a set of laboratory experiments. Convection is driven by a temperature gradient imposed between two rigid surfaces, and there is no internal heating. As the Rayleigh number is increased, two transitions in convective behaviour occur. First we observe a change from steady to time-dependent convection at Ra≈105. A second transition occurs at higher Rayleigh numbers, Ra≈5×106, with large-scale time-dependent flow being replaced by isolated rising and sinking plumes. Corresponding to the latter transition, the exponent β in the power law relating the Nusselt number Nu to the Rayleigh number (NuRaβ) is reduced. Both rising and sinking plumes always consist of plume heads followed by tails. There is no characteristic frequency for the formation of plumes.  相似文献   

5.
A collection of numerical simulations of 2D axi-symmetric thermal convection is presented here. The aim is to investigate the shape of geoid anomalies and dynamic topography above a plume. The simulation is based on the Boussinesq approximation and infinite Prandtl number and is carried out in the spherical shell with strongly temperature- and depth-dependent Arrhenius-type viscosity. According to the Arrhenius law, plume models with purely depth-dependent rheology are unphysical and should be taken with care. The strongly coupled temperature- and depth-dependent viscosity enables us to better understand the plume's behavior inside the Earth.The topography and geoid anomalies produced from plumes are sensitive to rheology of the mantle and rheology of the plume; both have effects on shape and amplitude of the geoid anomalies. We determined different categories of the geoid which are related to various rheology. Depth-dependent viscosity models show a geoid with a negative sign above the plume, and temperature-dependent viscosity models depict a bell-shaped geoid. We identified different behaviors in the combined model with temperature- and pressure-dependent viscosity.  相似文献   

6.
The process of multiple self-nucleation and ascent of mantle plumes is studied in the numerical models of thermal convection. The plumes are observed even in the simplest isoviscous models of thermal convection that leave aside the more complex rheology of the material, thermochemical effects, phase transformations, etc., which, although controlling the features of plumes, are not necessary for their formation. The origin of plumes is mainly due to the instability of the mantle flows at highly intense (low-viscous) thermal convection. At high viscosity, convective flows form regular cells. As viscosity decreases, the ascending and descending flows become narrower and unsteady. At a further decrease in viscosity, the ascending plumes assume a mushroom-like shape and occasionally change their position in the mantle. The lifetime of each flow can attain 100 Ma. Using markers allows visualizing the evolution of the shape of the mantle plumes.  相似文献   

7.
热和组分密度异常共同驱动的动态金幔柱解析模式给出金幔柱在金幔中上升的历史.当金幔柱到达岩石层底时其头部的特征量是5个独立变量的函数:(1)金幔柱起源的深度;(2)金幔粘性系数;(3)源温度异常值;(4)源组分与热密度异常之比;(5)金幔柱的浮力通量.基于新的Mapellan数据,金星表面上有360多个冠状和类冠状构造已被发现,其中65%直径小于300km.这类小型冠状构造被认为是由具有下述特征的较小的金幔柱所形成:(1)最大直径小于300km;(2)当其头部到达岩石层底时,过剩温度足以产生部分熔融层,应高于150K;(3)被冠状构造下面的金幔柱带上来的总浮力有能力支撑冠状构造隆起的总质量.用这3个条件分析数值结果并约束金幔柱的源参数,根据本文的数值实验结果,金星上的小型冠状构造可能是起源于上金幔小于1000km深度的动态金幔柱形成的.  相似文献   

8.
地幔对流的实验研究:非立柱状幔柱和地幔涡旋   总被引:7,自引:0,他引:7       下载免费PDF全文
地幔对流的物理模拟实验结果表明 ,在地幔介质和温度非均匀分布的复杂条件下 ,热卷流 (地幔柱 )往往由立柱状转变为非立柱状 (含斜柱状、涡旋状等 )。在忽略科里奥利力的情况下 ,板块的下插和滞积下沉、岩石圈根的存在以及地幔介质粘度的非均匀分布等都可能构成不同形状的障碍 -导流体 ,导致地幔的涡旋运动。软流圈中的水平涡旋环带属于对数螺线型 ,环带旋转半径及线速度逐渐减小 ,最终在旋转中心处下沉 ,而旋转角速度大致保持恒定  相似文献   

9.
Laboratory experiments were performed to study the influence of density and viscosity layering on the formation and stability of plumes. Viscosity ratios ranged from 0.1 to 6400 for buoyancy ratios between 0.3 and 20, and Rayleigh numbers between 105 and 2.108. The presence of a chemically stratified boundary layer generates long-lived thermochemical plumes. These plumes first develop from the interface as classical thermal boundary layer instabilities. As they rise, they entrain by viscous coupling a thin film of the other layer and locally deform the interface into cusps. The interfacial topography and the entrainment act to further anchor the plumes, which persist until the chemical stratification disappears through entrainment, even for Rayleigh numbers around 108. The pattern of thermochemical plumes remains the same during an experiment, drifting only slowly through the tank. Scaled to an Earth’s mantle without plate tectonics, our results show that: (1) thermochemical plumes are expected to exist in the mantle, (2) they could easily survive hundreds of millions of years, depending on the size and magnitude of the chemical heterogeneity on which they are anchored, and (3) their drift velocity would be at most 1-2 mm/yr. They would therefore produce long-lived and relatively fixed hotspots on the lithosphere. However, the thermochemical plumes would follow any large scale motion imposed on the chemical layer. Therefore, the chemical heterogeneity acts more as a ‘floating anchor’ than as an absolute one.  相似文献   

10.
If the interpretation of the D″ layer at the base of the mantle as a thermal boundary layer, with a temperature increment in the order of 800 K, is correct, then the formation of deep-mantle plumes to vent material from it appears inevitable. We demonstrate quantitatively that the strong temperature dependence of viscosity guides the upward flow into long-lived chimneys that are ~ 20 km in diameter near the base of the mantle and decrease in width with progressive upward softening and partial melting of plume material. The speed of flow up the axis of the plume is correspondingly fast; 1.6 m y?1 at the base and 4.8 m y?1 at 670 km depth. Thermal diffusive spreading of a heated plume is compensated by a slow horizontal convergence of mantle material toward the chimney in response to the lower pressure there. This convergence, which contributes only a small increment to the flux of material up the plume, also serves to throttle the flow in the chimney. The global plume mass flux necessary to transport 1.6 × 1012 W of core heat upward through the mantle is 1.8 × 106 kg s?1. At its base, plume material is probably still significantly below its solidus or eutectic temperature, but substantial partial melting is very likely as it rises. We speculate that a small fraction of this fluid component eventually emerges at the surface in “hot spots”, with the fate of the remainder being unknown. The behaviour and properties of D″ and of plumes are closely coupled. Not only are plumes a necessary consequence of a thermal boundary layer, but their existence is impossible without that layer.  相似文献   

11.
Based on the former workers study results such as numerical simulation of fluid mechanics,seismic tomography of the whole earth and igneous rocks,the basic characteristics of mantle plumes are summarized in detail,namely the mantle plume,from the D″layer near the core-mantle bouundary(CMB)of 2900 km deep,is characterized by the spape of large head and thin narrow conduit,by the physical property of high temperature and low viscosity.The LIP(large igneous province)is the best exhibition when the mantle plume ascends to the surface.According to the basic characteristics of the mantle plumes and the LIP,as well as the temporal-spatial relationships between the mantle plume and continental breakup,the detailed research on petrology,geochemistry,temporal-spatial distribution,tectonic background of the Cenozoic-Mesozoic igneous rocks and gravity anomaly distribution in East China has been done.As a result,the Mesozoic igneous rocks in Southeast China should not be regarded as an example of typical LIP related to mantle plumes.for their related characteristics are not consitent with those of the typical LIPs related to mantle plumes.The Cenozoic igneous rocks in Northeast China have no the typical characteristics of mantle plumes and hotspots,so the Cenozoic volcanism in Northeast China might have no the direct relationships with the activity of mantle plumes.  相似文献   

12.
Trace element relationships of near-primary alkalic lavas from La Grille volcano, Grande Comore, in the Indian Ocean, as well as those of the Honolulu volcanic series, Oahu, Hawaii, show that their sources contain amphibole and/or phlogopite. Small amounts of each mineral (2% amphibole in the source of La Grille and 0.5% phlogopite plus some amphibole in the source of the Honolulu volcanics) and a range of absolute degrees of partial melting from 1 to 5% for both series are consistent with the observed trace element variation. Amphibole and phlogopite are not stable at the temperatures of convecting upper mantle or upwelling thermal plumes from the deep mantle; however, they are stable at pressure-temperature conditions of the oceanic lithospheric mantle. Therefore, the presence of amphibole and/or phlogopite in the magma source region of volcanics is strong evidence for lithospheric melting, and we conclude that the La Grille and the Honolulu series formed by melting of the oceanic lithospheric mantle.

The identification of amphibole ± phlogopite in the source region of both series implies that the metasomatism by fluids or volatile-rich melts occurred prior to melting. The presence of hydrous phases results in a lower solidus temperature of the lithospheric mantle, which can be reached by conductive heating by the thermal plumes. Isotope ratios of the La Grille and the Honolulu series display a restricted range in composition and represent compositional end-members for each island. Larger isotopic variations in shield lavas, represented by the contemporaneous Karthala volcano on Grande Comore and the older Koolau series on Oahu, reflect interaction of the upwelling thermal plumes with the lithospheric mantle rather than the heterogeneity of deep-seated mantle plume sources or entrainment of mantle material in the rising plume. Literature OsSr isotope ratio covariations constrain the process of plume-lithosphere interaction as occurring through mixing of plume melts and low-degree melts from the metasomatized oceanic lithospheric mantle.

The characterization of the lithospheric mantle signature allows the isotopic composition of the deep mantle plume components to be identified, and mixing relationships show that the Karthala and Koolau plume end-members have nearly uniform isotopic compositions. Based on independent arguments, isotopic variations on Heard and Easter islands have been shown to be a result of mixing between deep plume sources having distinct isotopic compositions with lithosphere or shallow asthenospheric mantle. To the extent that these case studies are representative of oceanic island volcanism, they indicate that interaction with oceanic lithospheric mantle plays an important role in the compositions of lavas erupted during the shield-building stage of plume magmatism, and that isotopic compositions of deep mantle plume sources are nearly uniform on the scale that they are sampled by melting.  相似文献   


13.
动态地幔柱尾管结构   总被引:2,自引:0,他引:2       下载免费PDF全文
动态地幔柱模式被广泛用于讨论地球科学中的一些重要课题,如巨大火成区的成因,冈瓦纳古陆解体的原因,板块内中小尺度动力过程的驱动因素等.但是这个基于实验研究而建立的模式中,忽略了地幔柱尾管特征及其作用.地幔柱尾管内温度和速度分布是研究地幔柱上升过程的必要条件.本文从控制尾管结构的基本方程出发,给出了一个定常轴对称地幔柱温度和速度分布的近似分析解.从而得到尾管结构的基本特征:影响尾管内温度分布的主导因素是地幔柱的热流通量,而尾管内上涌速度的大小则不仅取决于热流通量,主要是取决于地幔粘度随深度的变化方式.结果表明,对弱地幔柱,尾管的热损失可能是不可忽略的,而对强地幔柱,径向质量传递可能是不可忽略的.  相似文献   

14.
The ∼0.2 mm/yr uplift of Hawaiian islands Lanai and Molokai and Hawaiian swell topography pose important constraints on the structure and dynamics of mantle plumes. We have formulated 3-D models of mantle convection to investigate the effects of plume-plate interactions on surface vertical motions and swell topography. In our models, the controlling parameters are plume radius, excess plume temperature, and upper mantle viscosity. We have found that swell height and swell width constraints limit the radius of the Hawaiian plume to be smaller than 70 km. The additional constraint from the uplift at Lanai requires excess plume temperature to be greater than 400 K. If excess plume temperature is 400 K, models with plume radius between 50 and 70 km and upper mantle viscosity between 1020 and 3×1020 Pa s satisfy all the constraints. Our results indicate that mantle plume in the upper mantle may be significantly hotter than previously suggested. This has important implications for mantle convection and mantle melting. In addition to constraining plume dynamics, our models also provide a mechanism to produce the observed uplift at Lanai and Molokai that has never been satisfactorily explained before.  相似文献   

15.
The emplacement of kimberlites in the North American and African continents since the early Palaeozoic appears to have occurred during periods of relatively slow motion of these continents. The distribution of kimberlites in time may reflect the global pattern of convection, which forces individual plates to move faster or slower at different times. Two-dimensional numerical experiments on a convecting layer with a moving upper boundary show two different regimes: in the first, when the upper boundary velocity is high, heat is transferred by the large-scale circulation and in the second, when the upper boundary velocity is lower, heat is predominantly transferred by thermal plumes rising from the lower boundary layer. For a reasonable mantle solidus, this second regime can give rise to partial melting beneath the moving plate, far from the plate boundaries. The transition between these modes takes place over a small range of plate velocities; for a Rayleigh number of 106 it occurs around 20 mm yr?1. We suggest that the generation of kimberlite magmas may result from thermal plumes incident on the base of a slowly moving plate.  相似文献   

16.
Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit (in the frontal part of the conduit). The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v.The local coefficients of heat transfer at the boundary of the plume conduit are theoretically determined and the balance of heat fluxes through the side of the plume conduit per linear meter of the conduit height. The total heat generation rate, transmitted by the Hawaiian plume into the upper and lower mantle, is evaluated. With the use of regular patterns of heat transfer in the lower mantle, which is modeled on the horizontal layer, heated from below and cooled from above, the diameter of the plume source, the kinematic viscosity of the melt in the plume conduit, and the velocity of horizontal lower-mantle flows are evaluated and the dependences of the temperature drop, viscosity and Rayleigh number for the lower mantle on the diameter of the plume source are presented.  相似文献   

17.
This paper presents a study on the effects of phase transitions on the mantle convection of Venus in a three-dimensional(3D)spherical shell domain.Our model includes strong depth-and temperature-dependent viscosity and exothermic phase change from olivine to spinel as well as endothermic phase change from spinel to perovskite.From extensive numerical simulations of the effects of Rayleigh number(Ra),and the Clapeyron slopes and depths of phase changes,we found the following:(1)The endothermic phase change prevents mass flow through the interface.Increasing the absolute value of the Clapeyron slopes decreases radial mass flux and normalized radial mass flux at the endothermic phase boundary,and decreases the number of mantle plumes.In other words,mass flow through the phase boundary decreases.The inhibition influence of phase changes increases,as do convective wavelengths.(2)Increasing Ra also increases the convective wavelength and decreases the number of mantle plumes,but it has less influence on the mass exchange.As Ra increases,the convective vigor increases along with the radial mass flux and the mass flow through the phase boundary;however,the normalized mass flux through the phase boundary varies little with Ra,which is different from the conclusion that increasing Ra will greatly increase the inhibition of mass flow through the phase boundary based on two-dimensional(2D)modeling.(3)Increasing the depth of endothermic phase change will slightly decrease the number of mantle plumes,but has little effect on the mass flow through the phase boundary.Consistent with previous studies,our results show that the phase change from spinel to perovskite could inhibit the mass flow through the phase boundary,but they also show that the buildup of hot materials under the endothermic phase boundary in the 3D model could not be so large as to cause strong episodic overturns of mantle materials,which is quite different from previous 2D studies.Our results suggest that it is difficult for phase changes to cause significant magmatism on Venus;in other words,phase changes may not be the primary cause of catastrophic resurfacing on Venus.  相似文献   

18.
The mantle convection model with phase transitions, non-Newtonian viscosity, and internal heat sources is calculated for two-dimensional (2D) Cartesian geometry. The temperature dependence of viscosity is described by the Arrhenius law with a viscosity step of 50 at the boundary between the upper and lower mantle. The viscosity in the model ranges within 4.5 orders of magnitude. The use of the non-Newtonian rheology enabled us to model the processes of softening in the zone of bending and subduction of the oceanic plates. The yield stress in the model is assumed to be 50 MPa. Based on the obtained model, the structure of the mantle flows and the spatial fields of the stresses σxz and σxx in the Earth’s mantle are studied. The model demonstrates a stepwise migration of the subduction zones and reveals the sharp changes in the stress fields depending on the stage of the slab detachment. In contrast to the previous model (Bobrov and Baranov, 2014), the self-consistent appearance of the rigid moving lithospheric plates on the surface is observed. Here, the intense flows in the upper mantle cause the drift and bending of the top segments of the slabs and the displacement of the plumes. It is established that when the upwelling plume intersects the boundary between the lower and upper mantle, it assumes a characteristic two-level structure: in the upper mantle, the ascending jet of the mantle material gets thinner, whereas its velocity increases. This effect is caused by the jump in the viscosity at the boundary and is enhanced by the effect of the endothermic phase boundary which impedes the penetration of the plume material from the lower mantle to the upper mantle. The values and distribution of the shear stresses σxz and superlithostatic horizontal stresses σxx are calculated. In the model area of the subducting slabs the stresses are 60–80 MPa, which is by about an order of magnitude higher than in the other mantle regions. The character of the stress fields in the transition region of the phase boundaries and viscosity step by the plumes and slabs is analyzed. It is established that the viscosity step and endothermic phase boundary at a depth of 660 km induce heterogeneities in the stress fields at the upper/lower mantle boundary. With the assumed model parameters, the exothermic phase transition at 410 km barely affects the stress fields. The slab regions manifest themselves in the stress fields much stronger than the plume regions. This numerically demonstrates that it is the slabs, not the plumes that are the main drivers of the convection. The plumes partly drive the convection and are partly passively involved into the convection stirred by the sinking slabs.  相似文献   

19.
From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.  相似文献   

20.
The data on the amplitude of variations in the direction and paleointensity of the geomagnetic field and the frequency of reversals throughout the last 50 Myr near the Paleozoic/Mesozoic and Mesozoic/Cenozoic boundaries, characterized by peaks of magmatic activity of Siberian and Deccan traps, and data on the amplitude of variations in the geomagnetic field direction relative to contemporary world magnetic anomalies are generalized. The boundaries of geological eras are not fixed in recorded paleointensity, polarity, reversal frequency, and variations in the geomagnetic field direction. Against the background of the “normal” field, nearly the same tendency of an increase in the amplitude of field direction variations is observed toward epicenters of contemporary lower mantle plumes; Greenland, Deccan, and Siberian superplumes; and world magnetic anomalies. This suggests a common origin of lower mantle plumes of various formation times, world magnetic anomalies, and the rise in the amplitude of geomagnetic field variations; i.e., all these phenomena are due to a local excitation in the upper part of the liquid core. Large plumes arise in intervals of the most significant changes in the paleointensity (drops or rises), while no correlation exists between the plume generation and the reversal frequency: times of plume formation correlate with the very diverse patterns of the frequency of reversals, from their total absence to maximum frequencies, implying that world magnetic anomalies, variations in the magnetic field direction and paleointensity, and plumes, on the one hand, and field reversals, on the other, have different sources. The time interval between magmatic activity of a plume at the Earth’s surface and its origination at the core-mantle boundary (the time of the plume rise toward the surface) amounts to 20–50 Myr in all cases considered. Different rise times are apparently associated with different paths of the plume rise, “delays” in the plume upward movement, and so on. The spread in “delay” times of each plume can be attributed to uncertainties in age determinations of paleomagnetic study objects and/or the natural remanent magnetization, but it is more probable that this is a result of the formation of a series of plumes (superplumes) in approximately the same region at the core-mantle boundary in the aforementioned time interval. Such an interpretation is supported by the existence of compact clusters of higher field direction amplitudes between 300 and 200 Ma that are possible regions of formation of world magnetic anomalies and plumes.  相似文献   

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

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