首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
《Planetary and Space Science》1999,47(10-11):1183-1200
Interior models of Jupiter and Saturn are calculated and compared in the framework of the three-layer assumption, which rely on the perception that both planets consist of three globally homogeneous regions: a dense core, a metallic hydrogen envelope, and a molecular hydrogen envelope. Within this framework, constraints on the core mass and abundance of heavy elements (i.e. elements other than hydrogen and helium) are given by accounting for uncertainties on the measured gravitational moments, surface temperature, surface helium abundance, and on the inferred protosolar helium abundance, equations of state, temperature profile and solid/differential interior rotation. Results obtained solely from static models matching the measured gravitational fields indicate that the mass of Jupiter’s dense core is less than 14 M (Earth masses), but that models with no core are possible given the current uncertainties on the hydrogen–helium equation of state. Similarly, Saturn’s core mass is less than 22 M but no lower limit can be inferred. The total mass of heavy elements (including that in the core) is constrained to lie between 11 and 42 M in Jupiter, and between 19 and 31 M in Saturn. The enrichment in heavy elements of their molecular envelopes is 1–6.5, and 0.5–12 times the solar value, respectively. Additional constraints from evolution models accounting for the progressive differentiation of helium (Hubbard WB, Guillot T, Marley MS, Burrows A, Lunine JI, Saumon D, 1999. Comparative evolution of Jupiter and Saturn. Planet. Space Sci. 47, 1175–1182) are used to obtain tighter, albeit less robust, constraints. The resulting core masses are then expected to be in the range 0–10 M, and 6–17 M for Jupiter and Saturn, respectively. Furthermore, it is shown that Saturn’s atmospheric helium mass mixing ratio, as derived from Voyager, Y=0.06±0.05, is probably too low. Static and evolution models favor a value of Y=0.11−0.25. Using, Y=0.16±0.05, Saturn’s molecular region is found to be enriched in heavy elements by 3.5 to 10 times the solar value, in relatively good agreement with the measured methane abundance. Finally, in all cases, the gravitational moment J6 of models matching all the constraints are found to lie between 0.35 and 0.38×10−4 for Jupiter, and between 0.90 and 0.98×10−4 for Saturn, assuming solid rotation. For comparison, the uncertainties on the measured J6 are about 10 times larger. More accurate measurements of J6 (as expected from the Cassini orbiter for Saturn) will therefore permit to test the validity of interior models calculations and the magnitude of differential rotation in the planetary interior.  相似文献   

2.
Ravit Helled  Gerald Schubert 《Icarus》2009,199(2):368-377
We present ‘empirical’ models (pressure vs. density) of Saturn's interior constrained by the gravitational coefficients J2, J4, and J6 for different assumed rotation rates of the planet. The empirical pressure-density profile is interpreted in terms of a hydrogen and helium physical equation of state to deduce the hydrogen to helium ratio in Saturn and to constrain the depth dependence of helium and heavy element abundances. The planet's internal structure (pressure vs. density) and composition are found to be insensitive to the assumed rotation rate for periods between 10h:32m:35s and 10h:41m:35s. We find that helium is depleted in the upper envelope, while in the high pressure region (P?1 Mbar) either the helium abundance or the concentration of heavier elements is significantly enhanced. Taking the ratio of hydrogen to helium in Saturn to be solar, we find that the maximum mass of heavy elements in Saturn's interior ranges from ∼6 to 20 M. The empirical models of Saturn's interior yield a moment of inertia factor varying from 0.22271 to 0.22599 for rotation periods between 10h:32m:35s and 10h:41m:35s, respectively. A long-term precession rate of about 0.754 yr−1 is found to be consistent with the derived moment of inertia values and assumed rotation rates over the entire range of investigated rotation rates. This suggests that the long-term precession period of Saturn is somewhat shorter than the generally assumed value of 1.77×106 years inferred from modeling and observations.  相似文献   

3.
Using the helium abundance measured by Galileo in the atmosphere of Jupiter and interior models reproducing the observed external gravitational field, we derive new constraints on the composition and structure of the planet. We conclude that, except for helium which must be more abundant in the metallic interior than in the molecular envelope, Jupiter could be homogeneous (no core) or could have a central dense core up to 12M. The mass fraction of heavy elements is less than 7.5 times the solar value in the metallic envelope and between 1 and 7.2 times solar in the molecular envelope. The total amount of elements other than hydrogen and helium in the planet is between 11 and 45M.  相似文献   

4.
We have calculated evolutionary and static models of Jupiter and Saturn with homogeneous solar composition mantles and dense cores of material consisting of solar abundances of SiO2, MgO, Fe, and Ni. Evolutionary sequences for Jupiter were calculated with cores of mass 2, 4, 6, and 8% of the Jovian mass. Evolutionary sequences for Saturn were calculated with cores of mass 16, 18, 20, and 22% of total mass. Two envelope mixtures, representative of the solar abundances were used: X (mass fraction of hydrogen) = 0.74, Y (mass fraction of helium) = 0.24 and X = 0.77 and Y = 0.21. For Jupiter, the observations of the temperature at 1 bar pressure (T1bar), radius and internal luminosity were best fit by evolutionary models with a core mass of ~6.5% and chemical composition of X = 0.77, Y = 0.21. The calculated cooling time for Jupiter is approximately 4.9 × 109 years, which is consistent, within our error bars, with the known age of the solar system. For Saturn, the observations of the radius, internal luminosity and T1BAR can be best fit by evolutionary models with a core mass of ~21% and chemical composition of X = 0.77, Y = 0.21. The cooling time calculated for Saturn is approximately 2.6 × 109 years, almost a factor 2 less than the present age of the solar system. Static models of Jupiter and Saturn were calculated for the above chemical compositions in order to investigate the sensitivity of the calculated gravitational moments, J2 and J4, to the mass of the dense core, T1BAR and hydrogen/helium ratio. We find for Jupiter that a model having a core mass of approximately 7% gives values of J2, J4, and T1BAR that are within observational limits, for the mixture X = 0.77, Y = 0.21. The static Jupiter models are completely consistent with the evolutionary results. For Saturn, the quantities J2, J4, and J6 determined from the static models with the most probable T1BAR of 140°K, using modeling procedures which result in consistent models for Jupiter, are considerably below the observed values.  相似文献   

5.
We present simple two-layer models of Uranus with rocky core and polytropic envelope satisfying exactly the observed mass, radius and the gravitational moments. The models show that the value of the fourth order zonal harmonic isJ 4 –38×10–6, whileJ 6 10–6. More elaborate threelayer models fail to satisfy the observational constraints of the ice/rock ratio and/or of the rotation period. We conclude that three-layer models with uniform chemical composition in each layer may be too restrictive. More realistic models should account for variable chemical composition within each layer.  相似文献   

6.
Expected characteristics of RR Lyrae stars as a function of the evolutive parameters are reported. Results from both evolutionary and pulsational investigations are collected in a suitable form, to show the general constraints to any interpretative analysis of the observations. It is shown that the spread in luminosity among the RR Lyrae stars results a function of the original chemical composition. On this basis a set of independent indications is found, suggesting that the globular cluster ω Cen is more He-rich than M 3; agreement with the whole observational frame is attained ifY ωCen~0.35,Z ωCen~5×10?4 andY M3~0.25,Z M3~10?3. No mass loss is needed to account for the RR Lyrae stars observed in ω Cen. The results are discussed, and it is shown that M 13-type clusters can be just characterized by a larger value ofZ in comparison with ω Cen. It is suggested that variations in the original helium content of the order of ΔY~0.1 and a correlationZ=Z(t) can account for some well-observed galactic globular clusters, without allowing for mass loss in the redder HB stars belonging to each cluster.  相似文献   

7.
A theory for the origin of the solar system, which is based on ideas of supersonic turbulent convection and indicates the possibility that the original Laplacian hypothesis may by valid, is presented. We suggest that the first stage of the Sun's formation consisted of the condensation of CNO ices (i.e. H2O, NH3, CH4,...) and later H2, including He as impurity atoms, at interstellar densities to from a cloud of solid grains. These grains then migrate under gravity to their common centre of mass giving up almost two orders of magnitude of angular momentum through resistive interaction with residual gases which are tied, via the ions, to the interstellar magnetic field. Grains rich in CNO rapidly dominate the centre of the cloud at this stage, both giving up almost all of their angular momentum and forming a central chemical inhomogeneity which may account for the present low solar neutrino flux (Prentice, 1976). The rest of the grain cloud, when sufficiently compressed to sweep up the residual gases and go into free fall, is not threatened by rotational disruption until its mean size has shrunk to about the orbit of Neptune. When the central opacity rises sufficiently to halt the free collapse at central density near 10?13 g cm?3, corresponding to a mean cloud radius of 104 R , we find that there is insufficient gravitational energy, for the vaporized cloud to acquire a complete hydrostatic equilibrium, even if a supersonic turbulent stress arising from the motions of convective elements becomes important, as Schatzman (1967) has proposed. Instead we suggest that the inner 3–4% of the cloud mass collapses freely all the way to stellar size to release sufficient energy to stabilize the rest of the infalling cloud. Our model of the early solar nebula thus consists of a small dense quasi-stellar core surrounded by a vast tenuous but opaque turbulent convective envelope. Following an earlier paper (Prentice, 1973) we show how the supersonic turbulent stress \((\rho _t v_t ^2 ) = \beta \rho GM(r)/r\) , where β is called the turbulence parameter, ρ is the gas density andM(r) the mass interior to radiusr causes the envelope to become very centrally condensed (i.e. drastically lowers its moment-of-inertia coefficientf) and leads to a very steep density inversion at its photosurface, as well as causing the interior to rotate like a solid body. As the nebula contracts conserving its angular momentum the ratio θ of centrifugal force to gravitational force at the equator steadily increases. In order to maintain pressure equilibrium at its photosurface, material is extruded outwards from the deep interior of the envelope to form a dense belt of non-turbulent gases at the equator which are free of turbulent viscosity. If the turbulence is sufficiently strong, we find that when θ→1 at equatorial radiusR e=R0, corresponding to the orbit of Neptune, the addition of any further mass to the equator causes the envelope to discontinuously withdraw to a new radiusR e>R0, leaving behind the circular belt of gas at the Kepler orbitR 0. The protosun continues to contract inwards, again rotationally stabilizing itself by extruding fresh material to the equator, and eventually abandoning a second gaseous ring at radiusR 1, and so on. If the collapse occurs homologously the sequence of orbital radiiR n of the system of gaseous Laplacian rings satisfy the geometric progression $$R_n /R_{n + 1} = [1 + m/Mf]^2 = constant, n = 0, 1,2, \ldots ,$$ analogous to the Titius-Bode Law of planetary distances, wherem denotes the mass of the disposed ring andM the remaining mass of the envelope. Choosing a ratio of surface to central temperature for the envelope equal to about 10?3 and adjusting the turbulence parameter β~~0.1 so thatR n/Rn+1 matches the observed mean ratio of 1.73, we typically findf=0.01 and that the rings of gas each have about the same mass, namely 1000M of the solar material. Detailed calculations which take into account non-homologous behaviour resulting from the changing mass fraction of dissociated H2 in the nebula during the collapse do not appreciably disturb this result. This model of the contracting protosun enables us to account for the observed physical structure and mass distribution of the planetary system, as well as the chemistry. In a later Paper II we shall examine in detail the condensation of the planets from the system of gaseous rings.  相似文献   

8.
This paper presents a new family of interior solutions of Einstein–Maxwell field equations in general relativity for a static spherically symmetric distribution of a charged perfect fluid with a particular form of charge distribution. This solution gives us wide range of parameter, K, for which the solution is well behaved hence, suitable for modeling of superdense star. For this solution the gravitational mass of a star is maximized with all degree of suitability by assuming the surface density equal to normal nuclear density, ρ nm=2.5×1017 kg?m?3. By this model we obtain the mass of the Crab pulsar, M Crab, 1.36M and radius 13.21 km, constraining the moment of inertia >?1.61×1038 kg?m2 for the conservative estimate of Crab nebula mass 2M . And M Crab=1.96M with radius R Crab=14.38 km constraining the moment of inertia >?3.04×1038 kg?m2 for the newest estimate of Crab nebula mass, 4.6M . These results are quite well in agreement with the possible values of mass and radius of Crab pulsar. Besides this, our model yields moments of inertia for PSR J0737-3039A and PSR J0737-3039B, I A =1.4285×1038 kg?m2 and I B =1.3647×1038 kg?m2 respectively. It has been observed that under well behaved conditions this class of solutions gives us the overall maximum gravitational mass of super dense object, M G(max)=4.7487M with radius $R_{M_{\max}}=15.24~\mathrm{km}$ , surface redshift 0.9878, charge 7.47×1020 C, and central density 4.31ρ nm.  相似文献   

9.
Models of Delta Scuti stars are tested against radial pulsations in the four lowest modes, using a linear, nonadiabatic approximation. It is shown that the theoretical blue edge locations of the instability strip on the H-R diagram depend strongly (ΔlogT e?0.03) on both opacity interpolation and variations of the envelope helium abundanceY. Several other factors have small effect (Δ logT e?0.01): variations of the heavy element abundanceZ, the use of various opacity tables, the treatment of convection, and the number of envelope mass zones. A comparison with both results obtained by various authors and observations is performed. The locations of the observed blue edge and our theoretical third-harmonic ones are consistent. A nonlinear dependence of the theoretical blue edge locations on envelope helium abundanceY is derived. A star withY?0.2 may pulsate within the instability strip, if it is not near the blue edge.  相似文献   

10.
The evolution of a rotating star with a mass of 16M at the hydrogen burning phase is considered together with the hydrodynamic processes of angular momentum transport in its interior. Shear turbulence is shown to limit the amplitude of the latitudinal variations in mean molecular weight on a surface of constant pressure in a layer with variable chemical composition. The resulting nonuniformity in the mean molecular weight distribution and the turbulent energy transport along the surface of constant pressure reduce the absolute value of the meridional circulation velocity. Nevertheless, meridional circulation remains the main mechanism of angular momentum transport in the radial direction in a layer with variable chemical composition. The intensity of the processes of angular momentum transport by meridional circulation and shear turbulence is determined by the angular momentum of the star. At a fairly high angular momentum, more specifically, at J = 3.69 × 1052 g cm2 s?1, the star during the second half of the hydrogen-burning phase in its convective core has characteristics typical of classical early Be stars.  相似文献   

11.
More than 80 giant planets are known by mass and radius. Their interior structure in terms of core mass, number of layers, and composition however is still poorly known. An overview is presented about the core mass M core and envelope mass of metals M Z in Jupiter as predicted by various equations of state. It is argued that the uncertainty about the true H/He EOS in a pressure regime where the gravitational moments J 2 and J 4 are most sensitive, i.e. between 0.5 and 4 Mbar, is in part responsible for the broad range \(M_{\mathit{core}}=0{-}18\:M_{\oplus }\), \(M_{Z}=0{-}38\:M_{\oplus }\), and \(M_{\mathit{core}}+M_{Z}=14{-}38\:M_{\oplus }\) currently offered for Jupiter. We then compare the Jupiter models obtained when we only match J 2 with the range of solutions for the exoplanet \(\mathrm{GJ}\:436\mathrm{b}\), when we match an assumed tidal Love number k 2 value.  相似文献   

12.
The paper presents a new class of parametric interior solutions of Einstein–Maxwell field equations in general relativity for a static spherically symmetric distribution of a charged perfect fluid with a particular form of electric field intensity. This solution gives us wide range of parameter, K (0.69≤K≤7.1), for which the solution is well behaved hence, suitable for modeling of superdense star. For this solution the gravitational mass of a superdense object is maximized with all degree of suitability by assuming the surface density of the star equal to the normal nuclear density ρ nm=2.5×1017kg?m?3. By this model we obtain the mass of the Crab pulsar M Crab=1.401M and the radius, R Crab=12.98 km constraining the moment of inertia I NS,38>1.61 for the conservative estimate of Crab nebula mass 2M and M Crab=2.0156M with radius, R Crab=14.07 km constraining the moment of inertia I NS,38>3.04 for the newest estimate of Crab nebula mass 4.6M which are quite well in agreement with the possible values of mass and radius of Crab pulsar. Besides this, our model yields the moments of inertia for PSR J0737-3039A and PSR J0737-3039B are I A,38=1.4624 and I B,38=1.2689 respectively. It has been observed that under well behaved conditions this class of parametric solution gives us the maximum gravitational mass of causal superdense object 2.8020M with radius 14.49 km, surface redshift z R =0.4319, charge Q=4.67×1020 C, and central density ρ c =2.68ρ nm.  相似文献   

13.
The various measurements of the linear matter density perturbation amplitude obtained from the observations of the cosmic microwave background (CMB) anisotropy, weak gravitational lensing, galaxy cluster mass function, matter power spectrum, and redshift space distortions are compared. The Planck data on the CMB temperature anisotropy spectrum at high multipoles, ? > 1000 (where the effect of gravitational lensing is most significant), are shown to give a measurement of the matter density perturbation amplitude that contradicts all other measurements of this quantity from both Planck CMB anisotropy data and other data at a significance level of about 3.7σ. Thus, at present these data should not be combined together for the calculations of constraints on cosmological parameters. Except for the Planck data on the CMB temperature anisotropy spectrum at high multipoles, all the remaining measurements of the density perturbation amplitude agree well between themselves and give the following constraints: σ8 = 0.792± 0.006 on the linear matter density perturbation amplitude, Ωm = 0.287± 0.007 on the matter density parameter, and H0 = 69.4 ± 0.6 km s?1 Mpc?1 on the Hubble constant. Various constraints on the sum of neutrino masses and the number of neutrino flavors can be obtained by additionally taking into account the data on baryon acoustic oscillations and (or) direct Hubble constant measurements in the local Universe.  相似文献   

14.
In the previous years, p-mode oscillations (pressure oscillations stochastically excited by convection) have been detected in several solar-like stars thanks to the ground-based spectroscopic and space spectroscopic and photometric observations. We study the importance of seismic constraints on stellar modeling and the impact of their accuracy on reducing the uncertainties of global stellar parameters (i.e. mass, age, etc.). We use the Singular Value Decomposition (SVD) method to analyze the sensitivity of stellar models to seismic constraints. In this context, we construct a grid of evolutionary sequences for solar-like stars with varying age and mass. Around each model of this grid, we evaluate the partial derivatives with respect to a large set of free parameters: mass ?, age τ, mixing-length parameter α, initial helium abundance Y 0, and initial metallicity Z/X 0. Masses between 0.9 and 1.55 M and central hydrogen abundances from Xc=0.7 to 0.05 have been considered in this study.  相似文献   

15.
We show that hydrostatically equilibrium models for the thin photospheres of helium stars based on new opacities κR (OPAL and OP) can be constructed only for masses M<5M . The parameter Г=κL/4πGMc, defined as the ratio of light pressure to gravity, exceeds a critical value of 1.0 for larger masses, which must result in mass outflow under light pressure. This mass limit matches the observed lower limit for the masses of Wolf-Rayet stars (M WR>5M )), which is an additional argument that the Wolf-Rayet stellar cores are actually helium stars. By solving the equation of radiative transfer in extended atmospheres, we construct a semiempirical model for a WN5 star (M WN5=10M )) with a helium core and an expanding envelope, whose physical and geometric parameters are known mainly from light-curve solution for the eclipsing binary V444 Cyg (WN5+06): outflow rate $\dot M \approx 1.0 \times 10^{ - 5} M_ \odot yr^{ - 1} $ , terminal velocity V ≈2000 km s?1, and expanding-envelope optical depth τenv≈25. The temperature at the outer boundary of the photosphere of a helium star surrounded by such an envelope is approximately 130 kK higher than that in the absence of an envelope, being T ph≈240 kK. Because of the high temperatures, the absorption coefficients at the corresponding photospheric levels are smaller than those in models with no envelope; therefore, the photosphere turns out to be in hydrostatic equilibrium and stable against light pressure (Гmax≈0.9). As a way out of this conflicting situation (an expanding envelope together with a hydrostatically equilibrium photosphere), we propose a model of discrete mass outflow, which is also supported by the observed cloudy structure of the envelopes in this type of stars. To quantitatively estimate parameters of the nonuniform outflow model requires detailed gas-dynamical calculations.  相似文献   

16.
《Planetary and Space Science》1999,47(10-11):1211-1224
The spectrum of free oscillations of Jupiter is calculated for a set of models, each of them fitting all available observational data. Diagnostic capabilities of the spectrum are studied. They could be used, as soon as relevant observations are performed, for both the identification of the observed modes and the improvement of the models themselves. The calculations were made for five-layer models. They differ in the core mass (2–10 M) and in the molecular-metallic phase transition pressure of hydrogen (1.5–3 Mbar). The spectrum of Jupiter consists of gravitational modes related to density jumps in the planetary interiors and of acoustic modes. The periods of the acoustic modes are calculated for degree up to l=30 and overtone number up to n=20. The investigated models have a characteristic frequency of ≈152–155 μHz. Two outer gravitational modes related to density jumps in the molecular envelope and at the interface with the metallic envelope have nonzero displacements at the planetary surface. These modes have good diagnostic properties. The values of the kinetic energy averaged over the period of oscillation are calculated for a 1-m amplitude of the displacement at the planetary surface. The influence of all effects of rotation on the spectrum is discussed.  相似文献   

17.
Spectroscopic observations of the recurrent Nova CI Aql in the wavelength range 4000–11000 Å are presented. Its evolution is traced from maximum light to the disappearance of nebular lines. 〈E B-V 〉 = 0.91 ± 0.11, as inferred from the Balmer decrement. The mean expansion velocity of the envelope measured near maximum light is 2800 km s?1. The helium abundance in the Nova envelope has been found to be enhanced, ?He/H?=0.22. CI Aql is similar in spectral evolution, in change of the envelope expansion velocity, and in helium abundance to other recurrent novae.  相似文献   

18.
S. Inaba  G.W. Wetherill 《Icarus》2003,166(1):46-62
We have calculated formation of gas giant planets based on the standard core accretion model including effects of fragmentation and planetary envelope. The accretion process is found to proceed as follows. As a result of runaway growth of planetesimals with initial radii of ∼10 km, planetary embryos with a mass of ∼1027 g (∼ Mars mass) are found to form in ∼105 years at Jupiter's position (5.2 AU), assuming a large enough value of the surface density of solid material (25 g/cm2) in the accretion disk at that distance. Strong gravitational perturbations between the runaway planetary embryos and the remaining planetesimals cause the random velocities of the planetesimals to become large enough for collisions between small planetesimals to lead to their catastrophic disruption. This produces a large number of fragments. At the same time, the planetary embryos have envelopes, that reduce energies of fragments by gas drag and capture them. The large radius of the envelope increases the collision rate between them, resulting in rapid growth of the planetary embryos. By the combined effects of fragmentation and planetary envelope, the largest planetary embryo with 21M forms at 5.2 AU in 3.8×106 years. The planetary embryo is massive enough to start a rapid gas accretion and forms a gas giant planet.  相似文献   

19.
The Kuiper Belt zone is unique insofar as the major heat sources of objects a few tens of kilometers in size—solar radiation on the one hand and radioactive decay on the other—have comparable power. This leads to unique evolutionary patterns, with heat waves propagating inward from the irradiated surface and outward from the radioactively heated interior. A major radioactive source that is considered in this study is 26Al. The long-term evolution of several models with characteristics typical of Kuiper Belt objects is followed by means of a 1-D numerical code that solves the heat and mass balance equations on a spherically symmetric grid. The free parameters considered are radius (10-500 km), heliocentric distance (30-120 AU), and initial 26Al content (0-5×10−8 by mass). The initial composition assumed is a porous mixture of ices (H2O, CO, and CO2) and dust. Gases released in the interior are allowed to escape to the surface. It is shown that, depending on parameters, the interior may reach quite high temperatures (up to 180 K). The models suggest that Kuiper Belt objects are likely to lose the ices of very volatile species during early evolution; ices of less volatile species are retained in a surface layer, about 1 km thick. The models indicate that the amorphous ice crystallizes in the interior, and hence some objects may also lose part of the volatiles trapped in amorphous ice. Generally, the outer layers are far less affected than the inner part, resulting in a stratified composition and altered porosity distribution. These changes in structure and composition should have significant consequences for the short-period comets, which are believed to be descendants of Kuiper Belt objects.  相似文献   

20.
Based on our spectrophotometric observations, we have studied the envelope of the HeN Nova Mon 2012. The abundances of some chemical elements in the envelope and its mass have been estimated. Our results show that the helium, nitrogen, oxygen, and neon abundances in the Nova envelope exceed the solar ones by a factor of 1.5, 33, 9, and 95, respectively. The envelope mass has been found to be 2.3 × 10?4 M .  相似文献   

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

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