首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Thermodynamic properties of PbO-SiO2 melts, obtained from published data and calculated from freezing point depressions, reflect the gradual polymerization of silicate anions in the melt as the SiO2PbO ratio is increased. The free energy of mixing curve at 1000°C has a minimum at 40 mole % SiO2 and is convex-upward between 72 and 98 mole % SiO2. The latter is an indication of metastable liquid immiscibility. The free energy minimum is correlated with the maximum in the distribution of nonbridging oxygens in the melt. In SiO2-poor melts, the activities of PbO and SiO2 (pure liquid standard states) show sharp negative deviations from ideality. The PbO activity reflects the paucity of free oxygen species in the melt whereas the SiO2 activity reflects the depolymerized state of the silicate anions. In more SiO2-rich melts, the activity of SiO2 shows a positive deviation from ideality which is qualitatively correlated to a polymerization parameter. The heat of mixing term has a minimum of ?2000 cal at 35 mole % SiO2 and a maximum of +200 cal at 90 mole % SiO2. The minimum is associated with the exothermic heat effect obtained during the reaction (O0) + (O2?) = 2(O?), whereas the maximum corresponds to the endothermic heat effect obtained when coordination polyhedra of oxygens form around the Pb cation. The entropy of mixing curve has the same form but is systematically smaller than a theoretical curve calculated on the assumption of random mixing of oxygen species. The discrepancy is due to the entropy loss obtained by the clustering of oxygen species to form complex silicate species.  相似文献   

2.
The argon solubility of 38 liquids in the system Na2O-CaO-MgO-Al2O3-SiO2 (NCMAS) has been determined at 1873 K and 1 bar, the argon concentration of presaturated glasses being measured using a static mass spectrometer. For compositions in the subsystem diopside (CaMgSi2O6), nepheline (NaAlSiO4), albite (NaAlSi3O8), anorthite (CaAl2Si2O8), argon solubility is generally a linear function of the relative proportion of each end member, solubility being lowest in diopside melt (1.53 10−5 cm3 STP · g−1 · bar−1) and highest in albite melt (2.88 10−4 cm3 STP · g−1 · bar−1). For the tectosilicate joins studied (SiO2-Na2Al2O4, SiO2-CaAl2O4, SiO2-MgAl2O4) solubility decreases with decreasing silica content in all cases, being highest for Na-bearing liquids and lowest for Mg-bearing liquids at constant molar silica content. Where comparison is possible our results are in good agreement with data from the literature. When our data are considered in isolation we find that argon solubility shows an excellent correlation with calculated ionic porosity. The covariation of argon solubility and liquid density is also reasonable, that with molar volume less convincing and that with polymerization state (as defined by the ratio of the number of nonbridging oxygens and tetrahedral network forming cations; NBO/T) nonexistent. However, when our data are combined with those from the literature no well constrained correlation between argon solubility and ionic porosity is apparent. Based upon this observation and consideration of the temperature dependence of noble gas solubility it is concluded that ionic porosity is not a universally applicable parameter which may be used to predict noble gas solubility as a function of composition, temperature and pressure. Two new models for calculating argon solubility are proposed, both employing the notion of partial molar argon solubilities. The first uses oxide components, for which partial molar argon solubility is directly proportional to partial molar ionic porosity calculated at 1873 K, irrespective of the temperature of experimental equilibration. The second model, which offers the best fit to the available data, employs tetrahedral units rather than oxides as the proposed melt components. This latter model successfully accounts for reported argon solubilities in simple Al-free systems, in simple Al-bearing systems and in natural liquids. This is interpreted to infer that argon is incorporated in large sites in the liquid structure (such as the space within rings of n-tetrahedra) although further work is required to understand the quantitative links between melt structure and noble gas solubility.  相似文献   

3.
4.
Phase equilibria data in the systems SiO2-P2O5, P2O5-MxOy, and P2O5-MxOy-SiO2 are employed in conjunction with Chromatographic and spectral data to investigate the role of P2O5 in silicate melts. Such data indicate that the behavior of P2O5 is complex. P2O5 depolymerizes pure SiO2 melts by entering the network as a four-fold coordinated cation, but polymerizes melts in which an additional metal cation other than silicon is present. The effect of this polymerization is apparent in the widening of the granite-ferrobasalt two-liquid solvus. In this complex system P2O5 acts to increase phase separation by further enrichment of the high charge density cations Ti, Fe, Mg, Mn, Ca, in the ferrobasaltic liquid. P2O5 also produces an increase in the ferrobasalt-granite REE liquid distribution coefficients. These distribution coefficients are close to 4 in P2O5-free melts, but close to 15 in P2O5-bearing melts.The dual behavior of P2O5 is explained in a model which requires complexing of phosphate anions (analogous to silicate anions) and metal cations in the melt. This interaction destroys Si-O-M-O-Si bonds polymerizing the melt. The higher concentration of Si-O-M-O-Si bond complexes in immiscible ferrobasaltic liquids relative to their conjugate immiscible granite liquids explains the partitioning of P2O5 into the ferrobasaltic liquid.  相似文献   

5.
The solubility behavior of K2O, Na2O, Al2O3, and SiO2 in silicate-saturated aqueous fluid and coexisting H2O-saturated silicate melts in the systems K2O-Al2O3-SiO2-H2O and Na2O-Al2O3-SiO2-H2O has been examined in the 1- to 2-GPa pressure range at 1100°C. Glasses of Na- and K-tetrasilicate compositions with 0, 3, and 6 mol% Al2O3 were used as starting materials. In both systems, the oxides dissolve incongruently in aqueous fluid and silicate melt. When recalculated to an anhydrous basis, the aqueous fluids are enriched in alkalis and depleted in silica and alumina relative to their proportions in the starting materials. The extent of incongruency is more pronounced in the Na2O-Al2O3-SiO2-H2O system than in the K2O-Al2O3-SiO2-H2O system.The partition coefficients of the oxides, Doxidefluid/melt, are linear and positive functions of the oxide concentration in the fluid for each composition. There is a slight dependence of the partition coefficients on bulk composition. No effect of pressure could be discerned. For alkali metals, the fluid/melt partition coefficients range from 0.06 to 0.8. For Al2O3 this range is 0.01 to 0.2, and for SiO2, it is 0.01 to 0.32. For all compositions, DK2Ofluid/melt∼DNa2Ofluid/melt>DSiO2fluid/melt>DAl2O3fluid/melt for the same oxide concentration in the fluid. DK2Ofluid/melt, DNa2Ofluid/melt, and DSiO2fluid/melt correlate negatively with the Al2O3 content of the systems. This correlation is consistent with a solubility model of alkalis that involve associated KOH°, NaOH°, silicate, and aluminate complexes.  相似文献   

6.
The mixing properties of aluminosilicate melts in the pseudobinary system NaAlSi3O8-KAlSi 3O8 have been determined by measuring the compositions of their saturated vapours by hightemperature Knudsen cell mass spectrometry. The melts mix close to ideally over most of the composition range with small positive deviations from ideality for K-rich compositions. These may be related to incipient partial ordering of melt constituents into leucite-like and SiO2-like structures above the feldspar liquidus.  相似文献   

7.
An inversion of SiO2 and MgO volatility occurs during high-temperature melt evaporation in the CaO–MgO–Al2O3–SiO2 (CMAS) system. This results in that SiO2, which is usually more volatile than MgO, becomes less volatile during the evaporation of melts enriched in the refractory oxides CaO and Al2O3. The volatility inversion is adequately explained within the theory of acid–base interaction of silicate melt components developed by D.S. Korzhinskii. The compositions of high-Al2O3 and SiO2-poor glasses (known as HASP glasses) from the lunar regolith show a systematic decrease in MgO/SiO2 with increasing CaO content, which is a direct consequence of the influence of acid–base effects.  相似文献   

8.
Nanometric solid inclusions in diamond incorporated in garnet and zircon from felsic gneiss of the Kokchetav massif, Kazakhstan, have been examined utilizing electron microscopy and focused ion beam techniques. Host garnet and zircon contain numerous pockets of multiple inclusions, which consist of 1–3 diamond crystals intergrown with quartz, phengite, phlogopite, albite, K‐feldspar, rutile, apatite, titanite, biotite, chlorite and graphite in various combinations. Recalculation of the average chemical composition of the entrapped fluid represented by multiple inclusion pockets indicates that such fluid contained a low wt% of SiO2, suggesting a relatively low‐temperature fluid rather than a melt. Transmission electron microscopy revealed that the diamond contains abundant nanocrystalline inclusions of oxides, rare carbonates and silicates. Within the 15 diamond crystals studied, abundant inclusions were found of SiO2, TiO2, FexOy, Cr2O3, ZrSiO4, and single grains of ThxOy, BaSO4, MgCO3, FeCr2O4 and a stoichiometric Fe‐rich pyroxene. The diversity of trace elements within inclusions of essentially the same stoichiometry suggests that the Kokchetav diamond crystallized from a fluid containing variable amounts of Si, Fe, Ti, Cr, Zr, Ba, Mg and Th and other minor components such as K, Na, P, S, Pb, Zn, Nb, Al, Ca, Cl. Most of the components in crystals included in diamond appear to have their origin in the subducted metasediments, but some of them probably originate from the mantle. It is concluded that Kokchetav diamond most likely crystallized from a COH‐rich multicomponent supercritical fluid at a relatively low temperature (hence the apparently low content of rock‐forming elements), and that the diversity of major and minor components suggests interactions between subducted metasediments and mantle components.  相似文献   

9.
The enthalpies of solution of synthetic Mg2SiO4-Fe2SiO4 olivine solid solutions have been measured in Pb2B2O5 melt at 970 K. The heat of solution of forsterite was found to be 15.62 ± 0.3 kcal mol?1 and that of fayalite 9.39 ± 0.14 kcal mol?1. Solid solutions between these end-members exhibit small positive deviations from mixing ideality, asymmetric towards the Fe end-member. In terms of the sub-regular solution model, excess enthalpies of intermediate olivine are adequately represented by the equation Hxs = 2(1000 + 1000XFe) XFeXMgThe enthalpies of solution at 970 K are consistent with high temperature phase equilibrium measurements of activity-composition relationships in the olivine series. Excess entropy terms are not needed to relate the phase equilibrium data to the calorimetric data presented here.The enthalpy of solution of FeSiO3 ferrosilite at 970 K was found to be 4.36 ± 0.10 kcal mol?1. This value, when taken together with calorimetric measurements on fayalite and quartz, is consistent with phase equilibrium investigations of the reaction: 2FeSiO3 = Fe2SiO4 + SiO2 Ferrosilite Fayalite QuartzThese provide a check on the internal consistency of the calorimetric data presented here.  相似文献   

10.
Principal components analysis is used to study the chemical compositions of pyroxenes of five Apollo 12 specimens. Important correlations recognized in the variation of oxide weight per cent are: MGO, Al2O3, SiO2| CaO, TiO2, FeO MgO, Al2O3, SiO2| FeO MgO, SiO2, FeO | Al2O3, CaO, TiO2 where the oxides on one side of the bar are correlated positively with each other and negatively with the oxides on the other side. Several other similarly distinct relationships with significantly less variance could be noted. These correlations indicating substitutional relationships can be interpreted as representative of stable and metastable trends of crystallization by using crystal-chemical and thermodynamic information. The per cent variance of pyroxene groups with characteristic trends in each specimen can be evaluated and interpreted in terms of history of crystallization. Distribution of Fe and Mg in certain pairs of olivine and pyroxene, which are found in contact in the rock and which may have crystallized simultaneously, is useful in recognizing the tendency towards chemical equilibrium in FeMg distribution during a limited interval in the liquidus or subsolidus stages.  相似文献   

11.
Solution enthalpies of synthetic olivine solid solutions in the system Mg2SiO4-Fe2SiO4 have been measured in molten 2PbO·B2O3 at 979 K. The enthalpy data show that olivine solid solutions have a positive enthalpy of mixing and the deviation from ideality is approximated as symmetric with respect to composition, in contrast to the previous study. Applying the symmetric regular solution model to the present enthalpy data, the interaction parameter of ethalpy (WH) is estimated to be 5.3±1.7 kJ/mol (one cation site basis). Using this Wh and the published data on excess free energy of mixing, the nonideal parameter of entropy (Ws) of olivine solid solutions is estimated as 0.6±1.5 J/mol·K.  相似文献   

12.
Ab initio STO-3G molecular orbital theory has been used to calculate energy-optimized Si-O bond lengths and angles for molecular orthosilicic and pyrosilicic acids. The resulting bond length for orthosilicic acid and the nonbridging bonds for pyrosilicic acid compare well with Si-OH bonds observed for a number of hydrated silicate minerals. Minimum energy Si-O bond lengths to the bridging oxygen of the pyrosilicic molecule show a close correspondence with bridging bond length data observed for the silica polymorphs and for gas phase and molecular crystal siloxanes when plotted against the SiOSi angle. In addition, the calculations show that the mean Si-O bond length of a silicate tetrahedron increases slightly as the SiOSi angle narrows. The close correspondence between the Si-O bond length and angle variations calculated for pyrosilicic acid and those observed for the silica polymorphs and siloxanes substantiates the suggestion that local bonding forces in solids are not very different from those in molecules and clusters consisting of the same atoms with the same coordination numbers. An extended basis calculation for H4SiO4 implies that there are about 0.6 electrons in the 3d-orbitals on Si. An analysis of bond overlap populations obtained from STO-3G* calculations for H6Si2O7 indicates that Si-O bond length and SiOSi angle correlations may be ascribed to changes in the hybridization state of the bridging oxygen and (dp) π-bonding involving all five of the 3d AO's of Si and the lone-pair AO's of the oxygen. Theoretical density difference maps calculated for H6Si2O7 show a build-up of charge density between Si and O, with the peak-height charge densities of the nonbridging bonds exceeding those of the bridging bonds by about 0.05 e Å?3. In addition, atomic charges (+1.3 and ?0.65) calculated for Si and O in a SiO2 moiety of the low quartz structure conform reasonably well with the electroneutrality postulate and with experimental charges obtained from monopole and radial refinements of diffraction data recorded for low quartz and coesite.  相似文献   

13.
The solubility mechanism of fluorine in quenched SiO2-NaF and SiO2-AlF3 melts has been determined with Raman spectroscopy. In the fluorine abundance range of F/(F+Si) from 0.15 to 0.5, a portion of the fluorine is exchanged with bridging oxygen in the silicate network to form Si-F bonds. In individual SiO4-tetrahedra, one oxygen per silicon is replaced in this manner to form fluorine-bearing silicate complexes in the melt. The proportion of these complexes is nearly linearly correlated with bulk melt F/(F+Si) in the system SiO2-AlF3, but its abundance increases at a lower rate and nonlinearly with increasing F/(F+Si) in the system SiO2-NaF. The process results in the formation ofnonbridging oxygen (NBO), resulting in stabilization of Si2O 5 2? units as well as metal (Na+ or Al3+) fluoride complexes in the melts. Sodium fluoride complexes are significantly more stable than those of aluminum fluoride.  相似文献   

14.
Incremental amounts of Na2O and K2O added to immiscible melts in the MgO-CaO-TiO2-Al2O3 SiO2 system cause a decrease in critical temperature, phase separation and change in the pattern of Al2O3 partitioning. Al2O3, which is concentrated in the low SiO2 immiscible melts in the alkali-free system, is increasingly partitioned into the high-SiO2 immiscible melt as the alkali/aluminium ratio is increased. However, K2O is more effective than Na2O in stabilizing Al2O2 in the SiO2-rich melt. The coordination changes occurring in the aluminosilicate melts upon the addition of the alkali oxides are described by CaAl2O4+2SiOK=2KAlO2+SiOCaOSi where K (or Na) displaces Ca as the charge-balancing cation for the networkforming AlO4 tetrahedra. The increased stability of the AlO4 species in the highly polymerized SiO2-rich melt and the consequent shrinkage of the miscibility gap is ascribed to positive configurational entropy and negative enthalpy changes associated with the formation of K, Na-AlO4 species. Element partition systematics indicate that (Na, K)AlO2 species favor the more polymerized, CaAl2O4 and TiO2 species, the less polymerized silicate structure in the melt.  相似文献   

15.
Twelve samples of Nigerian laterites were obtained from Ilorin, a rapidly growing urban center, and capital of Kwara State, Nigeria. Three varieties of laterites (clay, gravel and crust) were identified and subjected to mineralogical, chemical and geotechnical analyses which included: identification of clay and non-clay minerals by X-ray diffraction (XRD) techniques; chemical composition by X-ray fluorescence spectrometer analysis; pH of soil in water; moisture contents and specific gravity determinations, grain size analysis; compaction test by Harvard Compaction Apparatus and unconfined compressive strength determination.The laterite soil samples are composed of kaolinite and illite clay minerals with some quartz and feldspar. They were found to be rich in SiO2 (45%) Fe2O3, (16%) and Al203 (10%).These soils yielded maximum strength when compacted on the dry side of their optimum moisture content (OMC).The soils are not expected to perform very well as concrete aggregates since they contain high amounts of SiO2 and Fe2O3. These oxides are known to have deleterious effects on construction materials, particularly concrete aggregates.  相似文献   

16.
Revealing the atomic structure and disorder in oxide glasses, including sodium silicates and aluminosilicates, with varying degrees of polymerization, is a challenging problem in high-temperature geochemistry as well as glass science. Here, we report 17O MAS and 3QMAS NMR spectra for binary sodium silicate and ternary sodium aluminosilicate glasses with varying degrees of polymerization (Na2O/SiO2 ratio and Na2O/Al2O3 ratio), revealing in detail the extent of disorder (network connectivity and topological disorder) and variations of NMR parameters with the glass composition. In binary sodium silicate glasses [Na2O-k(SiO2)], the fraction of non-bridging oxygens (NBOs, Na-O-Si) increases with the Na2O/SiO2 ratio (k), as predicted from the composition. The 17O isotropic chemical shifts (17O δiso) for both bridging oxygen (BO) and NBO increase by about 10-15 ppm with the SiO2 content (for k = 1-3). The quadrupolar coupling products of BOs and NBOs also increase with the SiO2 content. These trends suggest that both NBOs and BOs strongly interact with Na; therefore, the Na distributions around BOs and NBOs are likely to be relatively homogenous for the glass compositions studied here, placing some qualitative limits on the extent of segregation of alkali channels from silica-enriched regions as suggested by modified random-network models. The peak width (in the isotropic dimension) and thus bond angle and length distributions of Si-O-Si and Na-O-Si increase with the SiO2 content, indicating an increase in the topological disorder with the degree of polymerization. In the ternary aluminosilicate glasses [Na2O]x[Al2O3]1−xSiO2, the NBO fraction decreases while the Al-O-Si and Al-O-Al fractions apparently increase with increasing Al2O3 content. The variation of oxygen cluster populations suggests that deviation from “Al avoidance” is more apparent near the charge-balanced join (Na/Al = 1). The Si-O-Si fraction, which is closely related to the activity coefficient of silica, would decrease with increasing Al2O3 content at a constant mole fraction of SiO2. Therefore, the activity of silica may decrease from depolymerized binary silicates to fully polymerized sodium aluminosilicate glasses at a constant mole fraction of SiO2.  相似文献   

17.
When melts of Ca–Al inclusions in chondrites, which are dominated by the oxides SiO2, MgO, CaO, and Al2O3, evaporate at high temperatures, the SiO2 and MgO fugacities are inverted: SiO2, which is more volatile than MgO, becomes less volatile when melts rich in refractory CaO and Al2O3 evaporate. This fugacity inversion can be realistically explained within the framework of D.S. Korzhinskii’s theory of acid–base interaction between components in silicate melts. According to this theory, an increase in CaO concentration in the melt increases its basicity, and this, in turn, increases the activity (and hence, also fugacity) of MgO and decreases those of SiO2. In the real compositions of the Ca–Al inclusions in chondrites, the MgO/SiO2 ratio systematically decreases with an increase in the CaO concentration under the effect of acid–base interaction.  相似文献   

18.
Isobaric (200 MPa) experiments have been performed to investigate the effects of H2O alone or in combination with P, S, F or Cl on liquid-phase separation in melts in the systems Fe2SiO4–Fe3O4–KAlSi2O6–SiO2, Fe3O4–KAlSi2O6–SiO2 and Fe3O4–Fe2O3–KAlSi2O6–SiO2 with or without plagioclase (An50). Experiments were heated in a rapid-quench internally heated pressure vessel at 1,075, 1,150 or 1,200 °C for 2 h. Experimental fO2 was maintained at QFM, NNO or MH oxygen buffers. H2O alone or in combination with P, S or F increases the temperature and composition range of two-liquid fields at fO2 = NNO and MH buffers. P, S, F and Cl partition preferentially into the Fe-rich immiscible liquid. Two-liquid partition coefficients for Fe, Si, P and S correlate well with the degree of polymerization of the SiO2-rich liquid and plot on similar but distinct power-law curves compared with equivalent anhydrous or basaltic melts. The addition of 2 wt% S to the system Fe3O4–Fe2O3–KAlSi2O6–SiO2 stabilizes three immiscible melts with Fe-, FeS- and Si-rich compositions. H2O-induced suppression of liquidus temperatures in the experimental systems, considered with the effects of pressure on the temperature and composition ranges of two-liquid fields in silicate melts, suggests that liquid-phase separation may be stable in some H2O-rich silicate magmas at pressures in excess of 200 MPa.  相似文献   

19.
A new statistical model is proposed for the molecular mass distributions (MMD) of polymerized anions in silicate melts. The model is based on the known distribution of Q n species in the MeO-Me2O-SiO2 system. In this model, chain and ring complexes are regarded as a random series of Q n structons with various concentrations of bridging bonds (1 ≤ n ≤ 4, Q 0 corresponds to SiO 4 4? ). This approach makes it possible to estimate the probability of formation of various ensembles of polymer species corresponding to the general formula (Si i O3i+1?j )2(i+1?j)?, where i is the size of the ion, and j is the cyclization number of intrachain bonds. The statistical model is utilized in the STRUCTON computer model, which makes use of the Monte Carlo method and is intended for the calculation of the composition and proportions of polyanions at a specified degree of polymerization of silicate melts (STRUCTON, version 1.2; 2007). Using this program, we simulated 1200 MMD for polyanions in the range of 0.52 ≤ p ≤98, where p is the fraction of nonbridging bonds in the silicon-oxygen matrix. The average number of types of anions in this range was determined to increase from three (SiO 4 4? , Si2O 7 6? , and Si3O 10 8? ) to 153, and their average size increases from 1 to 7.2. A special option of the STRUCTON program combines MMD reconstructions in silicate melts with the formalism of the Toop-Samis model, which enables the calculation of the mole fraction of the O2? ion relative to all anions in melts of specified composition. It is demonstrated that, with regard for the distribution and average size of anion complexes, the concentration of the O2? ion in the MeO-SiO2 system is characterized by two extrema: a minimum at 40–45 mol % SiO2, which corresponds to the initial stages of the gelenization of the polycondensated silicate matrix, and a maximum, which is predicted for the range of 60–80 mol % SiO2.  相似文献   

20.
New members of the crichtonite mineral series are described in which K, Ba, Ca and REE are in significant concentrations (5 wt% oxides) filling the A formula position in AM21O38. These phases are chromium (16 wt% Cr2O3) titanates (58 wt% TiO2) enriched in ZrO2 (5 wt%) and constitute a mineral repository for refractory and large ion lithophile elements in the upper mantle. The mineral senes coexists with Mg-Cr-ilmenite, Nb-Cr-rutile, and Ca-Cr (NbZr) armalcolite that have equally unusual chemistries. Kimberlitic crichtonites are depleted in the intermediate lanthanides but highly enriched in LREE and HREE with chondrite normalized abundances of 103 to 105. Crichtonite, armalcolite, and Nb-Cr-rutile occupy a compositional range in TiO2 contents bridging the gap between ilmenite and rutile, two minerals having a widespread distribution in kimberlites and mantle-derived nodule suites.In common with other associations, and based on similarities in mineral chemistry, it is concluded that these minerals formed at P = 20–30 kb, 900–1100°C by reaction of peridotite with metasomatizing fluids. Kimberlitic crichtonite may be expressed as spinel + Cr-ferropseudobrookite, and armalcolite is equivalent to Cr-geikielite + rutile in the system (FeMg)-TiO2-Cr2O3. This system contains a number of Cr-Ti compounds not found as minerals but it is proposed that the ubiquitous occurrence of ilmenite intergrowths in kimberlitic rutile results from decomposition of high pressure αPbO2-type crystallographic shear structures. The new minerals have exotic chemistries and the high K-affinities broaden the scope for the origin of alkalic rocks, the generation of highly potassic magmas in the upper mantle, and suggest that alkali metasomatism may be pervasive.  相似文献   

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

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