首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Major, minor and trace element abundances were determined in seven Angra dos Reis samples including whole rocks, fassaite (clinopyroxene), olivine and whitlockite separates via sequential instrumental neutron activation analysis. The chondritic normalized rare earth element (REE) abundance pattern for the Angra dos Reis clinopyroxene separates shows a concave downward shape with a small negative Eu anomaly. The strong fractionation between the light and the heavy REE in olivine separates could be attributed to the presence of islands of kirschsteinite in the olivines. The large-ion lithophile trace elements were highly enriched in the whitlockite separate as expected (e.g. La ≈ 370 ppm). The lower Hf and Sc abundances in whitlockite compared to that in the equilibrium “magma” could be the result of favorable partitioning of Hf and Sc in baddeleyite, which may have crystallized prior to or with whitlockite in the interstitial liquid. Comparison of whole rock with mineral separate data shows the presence of ~3% olivine, ~2.6% spinel and small amounts of metallic Ni-Fe and troilite in the whole rock.The trace element abundances in the derivative magma from which the Angra dos Reis clinopyroxene crystallized were estimated from the clinopyroxene data and the clinopyroxene mineral-liquid partition coefficients. From the derivative magma, the trace element abundances in the possible parent magmas were calculated by assuming that these parent magmas have undergone different degrees of clinopyroxene fractional crystallization to yield the Angra dos Reis derivative magma. Using the trace element abundances in these possible parent magmas, a two-stage crystal-liquid fractionation model with source material containing olivine, orthopyroxene and clinopyroxene is presented for the genesis of Angra dos Reis. Possible combinations of the degree of equilibrium non-modal partial melting, the source mineral composition and the initial element abundances required to generate possible Angra dos Reis parent magmas are calculated by the multilinear regression analysis method. Favorable solutions for this two-stage crystal-liquid fractionation model could be that Angra dos Reis crystallized at ~70% fractional crystallization of clinopyroxene from magmas generated by reasonable degrees of equilibrium partial melting (~7–10%) of deep-seated primitive source materials (olivine ~54–30%, orthopyroxene ~33–53%, and clinopyroxene ~13–17%) with trace element (Ba, Sr, REE and Sc) abundances ~3.5–4.7 × chondrites. These calculated REE abundances in the Angra dos Reis parent body are very similar to those suggested for the primordial moon (~3–5 × chondrites).Possible genetic relationships between Angra dos Reis and other achondrites, especially cumulate eucrites and nakhlites, are studied. Apparently, the unique Angra dos Reis could not be related to those achondrites by crystal-liquid fractionation of the same parent body.  相似文献   

2.
利用CHAMP卫星矢量和标量地磁测量数据,通过三维Taylor多项式模型建立了2010.0年中国及邻近地区在300km高度附近的X、Y、Z和F分量的三维磁场模型.为了比较验证所建模型,分析了其与Taylor多项式模型(二维模型)不同截断阶数所对应的均方偏差(RMSE)、残差及分布等.结果表明,三维模型每一阶的RMSE和残差绝对值的平均值均要比二维模型的小约45%.由于采用了系数完全展开的建模方式,三维模型的系数数量约为二维模型的2倍,三维模型较低的截断阶数可以反映更多的地磁信息.本研究中的5阶三维Taylor模型基本可达到8阶Taylor模型的精度.两种模型绘制的地磁场及残差分布有较好的一致性.  相似文献   

3.
三维复杂山谷地形SV波垂直输入地震反应分析   总被引:3,自引:0,他引:3       下载免费PDF全文
本文基于显式有限元法研究了地震波垂直入射时三维复杂山谷地形对地震地面运动的影响,在数值分析中应用了三维化二维的解法和黏弹性人工边界的处理方法,实现了地震波垂直输入下三维复杂场地地震动数值模拟,并验证了该方法的合理性.以四川桃坪地区一山谷地形作为研究对象,基于地表高程数据分别建立了二维和三维场地模型,对比研究表明:在复杂地形情况下考虑二、三维模型时具有明显差异,三维模型能更真实地反映地形变化对地震动的影响,复杂地形条件下有必要考虑三维实际场地模型.本文对边界自由场的处理方法也可用于处理三维复杂场地地震动斜入射问题,为三维复杂地形场地地震效应研究提供参考.  相似文献   

4.
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性.  相似文献   

5.
D. Yu  S. N. Lane 《水文研究》2006,20(7):1541-1565
High‐resolution data obtained from airborne remote sensing is increasing opportunities for representation of small‐scale structural elements (e.g. walls, buildings) in complex floodplain systems using two‐dimensional (2D) models of flood inundation. At the same time, 2D inundation models have been developed and shown to provide good predictions of flood inundation extent, with respect to both full solution of the depth‐averaged Navier–Stokes equations and simplified diffusion‐wave models. However, these models have yet to be applied extensively to urban areas. This paper applies a 2D raster‐based diffusion‐wave model to determine patterns of fluvial flood inundation in urban areas using high‐resolution topographic data and explores the effects of spatial resolution upon estimated inundation extent and flow routing process. Model response shows that even relatively small changes in model resolution have considerable effects on the predicted inundation extent and the timing of flood inundation. Timing sensitivity would be expected, given the relatively poor representation of inertial processes in a diffusion‐wave model. Sensitivity to inundation extent is more surprising, but is associated with: (1) the smoothing effect of mesh coarsening upon input topographical data; (2) poorer representation of both cell blockage and surface routing processes as the mesh is coarsened, where the flow routing is especially complex; and (3) the effects of (1) and (2) upon water levels and velocities, which in turn determine which parts of the floodplain the flow can actually travel to. It is shown that the combined effects of wetting and roughness parameters can compensate in part for a coarser mesh resolution. However, the coarser the resolution, the poorer the ability to control the inundation process, as these parameters not only affect the speed, but also the direction of wetting. Thus, high‐resolution data will need to be coupled to a more sophisticated representation of the inundation process in order to obtain effective predictions of flood inundation extent. This is explored in a companion paper. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

6.
Trace element abundances in Ivory Coast normal microtektites and Australasian bottle-green microtektites confirm that microtektites are genetically related to tektites in the associated strewn field. Although major and compatible trace element abundances imply that bottle-green microtektites are members of a fractional crystallization sequence, the similar rare earth element distributions in Australasian normal and bottle-green microtektites and tektites cannot be explained by a simple fractionation model. The similar REE abundances in tektites and microtektites of widely different major element composition also preclude simple models calling for sedimentary rock precursors.  相似文献   

7.
V. Tayefi  S. N. Lane  R. J. Hardy  D. Yu 《水文研究》2007,21(23):3190-3202
A much understudied aspect of flood inundation is examined, i.e. upland environments with topographically complex floodplains. Although the presence of high‐resolution topographic data (e.g. lidar) has improved the quality of river flood inundation predictions, the optimum dimensionality of hydraulic models for this purpose has yet to be fully evaluated for situations of both topographic and topological (i.e. the connectivity of floodplain features) complexity. In this paper, we present the comparison of three treatments of upland flood inundation using: (a) a one‐dimensional (1D) model (HEC‐RAS v. 3·1·2) with the domain defined as series of extended cross‐sections; (b) the same 1D model, but with the floodplain defined by a series of storage cells, hydraulically connected to the main river channel and other storage cells on the floodplain according to floodplain topological characteristics; (c) a two‐dimensional (2D) diffusion wave treatment, again with explicit representation of floodplain structural features. The necessary topographic and topological data were derived using lidar and Ordnance Survey Landline data. The three models were tested on a 6 km upland reach of the River Wharfe, UK. The models were assessed by comparison with measured inundation extent. The results showed that both the extended cross‐section and the storage cell 1D modes were conceptually problematic. They also resulted in poorer model predictions, requiring incorrect parameterization of the main river to floodplain flux in order to approach anything like the level of agreement observed when the 2D diffusion wave treatment was assessed. We conclude that a coupled 1D–2D treatment is likely to provide the best modelling approach, with currently available technology, for complex floodplain configurations. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

8.
This paper addresses the time history finite element analysis of rock-structure interaction. Modeled is not only the lateral energy dissipation, but also the interaction between the far field and the numerical model itself. This is accomplished by a preliminary analysis of the far field as a shear beam (for lateral excitation), and then velocities and displacements are transferred to the model as nodal forces through damping and stiffness matrices respectively. Details of the finite element implementation are given, along with an extensive series of simulations comparing this method, with the one of Lysmer for both 2D and 3D models. The model is derived from the principle of virtual work, and its implementation does not require any modification of existing finite element codes, only clever pre and postprocessing of results are needed.  相似文献   

9.
New rare earth element (REE) data for Archaean basalts and spinifex-textured peridotites (STP) show a range of La/Sm ratios (chondrite-normalized) from 0.36 to 3.5, with the bulk of the data in the range 0.7–1.3. This supports the hypothesis, based on Sr isotope initial ratios, that the Archaean mantle was chemically heterogeneous. We suggest that the bulk mantle source for Archaean basaltic magmas was close to an undepleted earth material. An average chemical composition of the Archaean mantle is estimated using chemical regularities observed in Archaean STP and high-magnesian basalts. TiO2 and MgO data show an inverse correlation which intersects the MgO axis at about 50% MgO (Fo92). TiO2 abundance in the mantle source is measured on this plot by assigning anMgO= 38% for the mantle. Concentrations of other elements are also estimated and these data are then used to obtain a composition for the bulk earth. We suggest an earth model with about 1.35 times ordinary chondrite abundances of refractory lithophile elements and about 0.2 times carbonaceous type 1 chondrite abundances of moderately volatile elements (such as Na, Rb, K, Mn). P shows severe depletion in the model earth relative to carbonaceous chondrites, a feature either due to volatilization or core formation (preferred). Our data support the hypothesis of Ringwood that the source material for the earth is a carbonaceous chondrite-like material.The generation of mid-ocean ridge basalts (MORB) is examined in the light of the model earth composition and Al2O3/TiO2, CaO/TiO2 ratios. It is suggested that for primitive basalts, these values can be used to predict the residual phases in their source. Comparison of chemical characteristics of inferred sources for 2.7-b.y. Archaean basalts and modern “normal” MORB indicates that the MORB source is severely depleted in highly incompatible elements such as Cs, Ba, Rb, U, Th, K, La and Nb, but has comparable abundances of less incompatible elements such as Ti, Zr, Y, Yb. The cause of the depletion in the MORB source is examined in terms of crust formation and extraction of silica-undersaturated melts. The latter seems to be a more likely explanation, since the degree of enrichment of highly incompatible elements in the crust only accounts for up to 40% of their abundances in the bulk earth and cannot match the depletion pattern in normal MORB. A large volume of material, less depleted than the source for normal MORB must therefore exist in the mantle and can serve as the source for the ocean island basalts and “normal” MORB.Three different mantle evolution models are examined and each suggests that the mantle is stratified with respect to abundances of incompatible trace elements. We suggest that no satisfactory model is available to fully explain the spectrum of geochemical and geophysical data. In particular the Pb and Sr isotope data on oceanic basalts, the depletion patterns of MORB and the balance between lithophile abundances in the crust and mantle, are important geochemical constraints to mantle models. Further modelling of the mantle evolution will be dependent on firmer information on the role of subduction, mantle convection pattern, and basalt production through geologic time together with a better understanding of the nature of Archaean crustal genesis.  相似文献   

10.
考虑关断时间的回线源激发TEM三维时域有限差分正演   总被引:14,自引:9,他引:5       下载免费PDF全文
从麦克斯韦旋度方程出发可以直接导出瞬变电磁场扩散方程,然而扩散方程不含电场对时间的一阶导数,不能构成显式的时域有限差分方程,借鉴du Fort-Frankel有限差分离散方法引入虚拟位移电流项构建显式时域有限差分方程.对Wang和Hohmann的经典时域算法进行了两点改进:第一,通过将矩形回线源电流密度加入麦克斯韦方程组的安培环路定理方程,实现回线源瞬变电磁激发源加入;第二,在计算中考虑关断时间.第一点改进使时域有限差分方程考虑了一次场的计算,并且源的计算不再依赖均匀半空间模型响应作为初始条件,使算法能够适应表层电阻率不均匀时的三维复杂模型.由于实际观测中不可能出现阶跃电流的关断形式,第二点改进可以方便设置发射电流下降沿.采用改进的三维时域有限差分正演算法对均匀半空间模型、四类三层模型、均匀半空间中含有低阻块体模型进行了计算并分别与解析解、线性数字滤波解、积分方程解和Wang的三维时域有限差分解进行了对比验证.以H模型为例,采用建立的三维时域有限差分正演算法计算了不同关断时间的斜阶跃脉冲回线源瞬变电磁中心点感应电动势衰减曲线.以实际地质资料为基础,构建包含两层采空区的三维复杂模型,以1 μs的极短关断时间进行了复杂模型定回线源瞬变电磁响应计算,并计算了该复杂模型的视电阻率曲线.  相似文献   

11.
The high observed abundances of Na and Cu in chondrules indicate that the amount of loss during chondrule formation was minor and possibly negligible, consistent with the view that loss was controlled by diffusion kinetics rather than equilibrium volatility, and that the surface of the chondrule quickly cooled to temperatures at which diffusional transport was negligible. Ordinary chondrite/CI abundance ratios appear to be randomly distributed in the range 0.9-0.1. Very few values are observed in the 0.36–0.70 range, but this is not statistically significant, nor is it predicted by the two-component (chondrule-matrix) model.If CI chondrite abundances are representative of mean solar-system material, the very low chondrule content in CM chondrites (<5% of high-temperature materials) indicates that the observed volatile distribution resulted from incomplete accretion of volatile carriers (perhaps a fine aerosol). At the ordinary chondrite formation location the fraction of an element sited in unaccreted carriers increased with decreasing condensation temperature. At the CM location a similar trend is observed for elements less volatile than S, but the unaccreted fraction of more volatile elements was nearly constant.  相似文献   

12.
Abstract

Data unavailability is the main reason for limited applications of hydrodynamic models for predicting inundation in the developing world. This paper aims to generate moderately high-resolution hybrid terrain data by merging height information from low-cost Indian Remote Sensing satellite (IRS) Cartosat-1 stereo satellite images, freely-available Shuttle Radar Topograph Mission (SRTM) digital elevation model (DEM) data, and limited surveyed channel cross-sections. The study reach is characterized by anabranching channels that are associated with channel bifurcation, loops and river islands. We compared the performance of a simple 1D–2D coupled LISFLOOD-FP model and a complex fully 2D finite element TELEMAC-2D model with the hybrid terrain data. The results show that TELEMAC-2D produced significantly improved simulated inundation with the hybrid terrain data, as compared to the SRTM DEM. LISFLOOD-FP was found unsuitable to work with the hybrid DEM in a complicated fluvial environment, as it failed to efficiently divert water in the branches from the main channel.
Editor D. Koutsoyiannis; Associate editor A. Viglione  相似文献   

13.
常规的三维时间域航空电磁模拟通常采用隐式步长方法进行时间离散,需要几次矩阵分解和上百次右端源项回带,计算效率较低.为了提高正演计算效率,本文提出使用有理Krylov方法求解时间域电场扩散方程.首先使用非结构四面体网格进行空间离散,采用Nédélec矢量基函数近似四面体单元内的电场;然后基于有限元离散给出矩阵指数和矢量乘积表示的电场显式解;最后采用有理Arnoldi算法构造Krylov子空间内的正交基函数并进一步求解矩阵指数与矢量的乘积,直接得到任意时刻的电场解向量,避免步长离散过程.此外,本文还提出一种指数加权偏移参数优化方法,使得有理Arnoldi近似在瞬变衰减晚期具备更高的精度,从而降低Krylov子空间阶数并提高计算效率.通过和层状模型解析解的对比验证了有理Krylov方法的精度.针对三维异常体模型使用全局网格和局部网格剖分并和其他数值方法比较,进一步说明了有理Krylov方法的有效性.  相似文献   

14.
A correct understanding of the hydrodynamics and morphodynamics of tidal basins is of fundamental importance for the fate of the Venice Lagoon, Italy. If on one hand, the development of sophisticated numerical models is called for in order to reproduce the complexity of the mechanisms governing the morphodynamic evolution of many natural environments, including lagoons, on the other hand, a clear knowledge of the reliability and limits of the results provided by these models is crucial in order to establish the condition under which they can be safely applied. To this aim, researchers involved in numerical modeling in the framework of the recent Corila research programmes, agreed to perform an accurate comparison of results provided by three different numerical models, applying them to the test case offered by the experimental investigations performed under controlled conditions by Tambroni et al. (2005a). Here, we consider the following numerical models: (i) a 2D finite element hydrodynamic model coupled with a 2D finite volume morphodynamic model (5 and 3); (ii) a 2D finite element morphodynamic model (Ferrarin et al., 2008); (iii) a 2D depth-averaged model for the inlet region, coupled with a 1D model for the channel (Tambroni et al., 2005b). A first set of simulations concerns the fixed bed case and shows that all the models provide similar results: in particular, they are able to predict the observed free surface oscillations satisfactorily, while comparison with the measured velocity field is less satisfactory. Moreover, as far as the flow field at the inlet is concerned, the models describe accurately the potential flow into the channel during the flood phase, while they are not able to adequately reproduce the occurrence of the fine structure of the shear layers shed by the inlet edges during the ebb phase. This limit is related to the shallow water character of the models. As for the morphodynamics, the long term equilibrium configurations of the bottom of the channel and of the near inlet region show qualitative agreement with the experimental observations, although in this case the differences between the results provided by the distinct numerical approaches are more marked.  相似文献   

15.
Granodiorite and diorite plutons of the Rosetown complex, N.Y., which are associated with the nearby Cortlandt complex, have chemical and textural characteristics indicating that large-scale liquid immiscibility played a major role in their petrogenesis. Rare earth element, zirconium, niobium and phosphorus abundances are much greater in the diorite precluding the possibility that the Rosetown diorite and granodiorite are related by fractional crystallization. The trace element data also eliminate the possibility that the granodiorite represents: (1) a partial melt of crustal rocks including basalt; (2) a granitic cumulate; or (3) a residue from an aqueous fluid derived either from a silicate melt or crustal rocks. Liquid immiscibility appears to be viable model for the origin of the Rosetown granodiorite and iron-rich diorite. This model is supported by the following: (1) the major element compositions occur in a two-liquid field on a Greig diagram; (2) both bodies have similar Sr isotope compositions; (3) common phases in the two rock types have overlapping compositions; (4) the major and trace element data of the diorite and granodiorite are similar to the experimentally determined partition data of immiscible liquid pairs; and (5) possible ocelli of iron-rich diorite are found in the granodiorite.  相似文献   

16.
基于非结构网格的电阻率三维带地形反演   总被引:6,自引:3,他引:3       下载免费PDF全文
吴小平  刘洋  王威 《地球物理学报》2015,58(8):2706-2717
地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探.  相似文献   

17.
Rare earth element (REE) abundances determined by activation analysis in rocks, plagioclase and mafic separates from the Fiskenaesset Complex are presented together with data on major and trace elements in the minerals. The REE data for the rocks and plagioclases are distinct from those of many other anorthositic complexes and the abundances are some of the lowest recorded for plagioclase from terrestrial anorthosites. The bulk and trace element compositions of the Fiskenaesset plagioclases show a number of similarities to those of lunar plagioclases. The plagioclases show a positive Eu anomaly of about 10 and a depletion in the heavy REE relative to the light ones. The mafic separates are enriched in the heavy REE relative to the light ones, and show no Eu anomaly except in one sample with a positive anomaly not attributable to plagioclase contamination. It is estimated, from experimental partition coefficient data, that the REE pattern in the magma at an early stage of fractionation was La (17×) to Lu (0.7× chondrites) with a possible positive Eu anomaly. This highly fractionated REE pattern may be attributed to partial melting of a garnet-bearing source.  相似文献   

18.
Petrogenesis of basalts from the FAMOUS area: Mid-Atlantic Ridge   总被引:1,自引:0,他引:1  
Fresh basalt glasses most of which have Mg/(Mg + Fe2+) of 0.66–0.72 from outcrops within 3 km of one other in the rift valley at the Project FAMOUS locality have been analyzed for major, minor and trace elements in order to determine their petrogenesis.Transition metal abundances of the FAMOUS samples are similar to a wide variety of continental and oceanic basalts with high MgO and Ni, all of which show remarkably little variation, with the exception of Cu, Zn and Ti, on a chondrite-normalized plot. Modelling of these data suggests that the mantle beneath both continents and oceans is systematically fractionated relative to chondrites. This fractionation provides a constraint for models of earth formation and subsequent evolution.The abundances of the rare earth and the incompatible elements, Ba, La, Th, U, and Nb, vary by more than a factor of three and the La/Yb and La/Sm ratios vary by factors of 3.1 and 1.6, respectively, in samples with similar, high Mg/(Mg + Fe2+). There is no correlation between the degree of light-REE enrichment and the heavy-REE abundance. Furthermore, the trace element variations do not appear correlated with respect to location in the rift valley or to time of eruption. These trace element features demonstrate that successive eruptions in one small area of the rift valley can show wide variations in trace element chemistry over a short span of time; they preclude the derivation of these basalt glasses from a single magma chamber.Despite the heterogeneities in REE and the variable trace element abundances, a homogeneous mantle source is suggested by the similarities among the samples in the incompatible element ratios of La/Ce, Ba/Th, Zr/Nb and K/Ba and the small range in87Sr/86Sr isotope ratios observed in other samples from the FAMOUS region (White and Bryan, 1977). Thus, trace element heterogeneities appear to be generated by processes in the mantle during melting. However, processes such as batch partial melting, fractional fusion, fractional crystallization, zone refining, or mixing of magmas or sources acting alone are incapable of explaining the lack of correlation between the light and heavy REE.It is suggested that the observed variations are a consequence of dynamic partial melting of a homogeneous mantle source region. This process includes varying degrees of partial melting of an uprising mantle source with continuous but incomplete removal of melt as melting proceeds, varying extents of batch partial melting, and zone refining. Dynamic melting can produce different melts from a homogeneous source which have different degrees of light-REE enrichment and crossing REE patterns. The variable trace element abundances which may be produced through dynamic melting may be the cause of the apparent decoupling of major and trace elements (Bryan et al., 1976) which previously has been suggested for the FAMOUS region (Bryan and Moore, 1977).  相似文献   

19.
1 INTRODUCTION The Parana River in South America ranks as the sixth largest in the world in terms of average discharge (22 x 103 m3/s according to Schumm and Winkley, 1994), and its Argentinean reach stretches from Puerto Iguazu on the north to the Rio de la Plata estuary in the south (Fig. 1). Its estimated average suspended load, bed load, and total load are 23 x 106, 2.2 x 106, and 25.2 x 106 metric tons per year, respectively (Paoli and Schreider, 2001.) The economic developmen…  相似文献   

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

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