首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Ionospheric scintillations are fluctuations in the phase and/or amplitude of trans-ionospheric radio signals caused by electron density irregularities in the ionosphere. A better understanding of the scintillation pattern is important to make a better assessment of GPS receiver performance, for instance. Additionally, scintillation can be used as a tool for remote sensing of ionospheric irregularities. Therefore, the study of ionospheric scintillation has both scientific as well as technological implications. In the past few years, there has been a significant advance in the methods for analysis of scintillation and in our understanding of the impact of scintillation on GPS receiver performance. In this work, we revisit some of the existing methods of analysis of scintillation, propose an alternative approach, and apply these techniques in a comprehensive study of the characteristics of amplitude scintillation. This comprehensive study made use of 32?days of high-rate (50?Hz) measurements made by a GPS-based scintillation monitor located in S?o José dos Campos, Brazil (23.2°S, 45.9°W, ?17.5° dip latitude) near the Equatorial Anomaly during high solar flux conditions. The variability of the decorrelation time (τ0) of scintillation patterns is presented as a function of scintillation severity index (S 4). We found that the values of τ0 tend to decrease with the increase of S 4, confirming the results of previous studies. In addition, we found that, at least for the measurements made during this campaign, averaged values of τ0 (for fixed S 4 index values) did not vary much as a function of local time. Our results also indicate a significant impact of τ0 in the GPS carrier loop performance for S 4?≥?0.7. An alternative way to compute the probability of cycle slip that takes into account the fading duration time is also presented. The results of this approach show a 38% probability of cycle slips during strong scintillation scenarios (S 4 close to 1 and τ0 near 0.2?s). Finally, we present results of an analysis of the individual amplitude fades observed in our set of measurements. This analysis suggests that users operating GPS receivers with C/N 0 thresholds around 30?dB-Hz and above can be affected significantly by low-latitude scintillation.  相似文献   

2.
Observations of non-Fickian transport in sandbox experiments [Levy M, Berkowitz B. Measurement and analysis of non-Fickian dispersion in heterogeneous porous media. J Contam Hydrol 2003;64:203–26] were analyzed previously using a power law tail ψ(t) ∼ t−1−β with 0 < β < 2 for the spectrum of transition times comprising a tracer plume migration. For each sandbox medium a choice of β resulted in an excellent fit to the breakthrough curve (BTC) data, and the value of β decreased slowly with increasing flow velocity. Here, the data are reanalyzed with the full spectrum of ψ(t) gleaned from analytical calculations [Cortis A, Chen Y, Scher H, Berkowitz B. Quantitative characterization of pore-scale disorder effects on transport in “homogeneous” granular media. Phys Rev E 2004;10(70):041108. doi: 10.1103/PhysRevE.70.041108], numerical simulations [Bijeljic B, Blunt MJ. Pore-scale modeling and continuous time random walk analysis of dispersion in porous media. Water Resour Res 2006;42:W01202. doi: 10.1029/2005WR004578] and permeability fields [Di Donato G, Obi E-O, Blunt MJ. Anomalous transport in heterogeneous media demonstrated by streamline-based simulation. Geophys Res Lett 2003;30:1608–12s. doi: 10.1029/2003GL017196]. We represent the main features of the full spectrum of transition times with a truncated power law (TPL), ψ(t) ∼ (t1 + t)−1−βexp(−t/t2), where t1 and t2 are the limits of the power law spectrum. An excellent fit to the entire BTC data set, including the changes in flow velocity, for each sandbox medium is obtained with a single set of values of t1, β, t2. The influence of the cutoff time t2 is apparent even in the regime t < t2. Significantly, we demonstrate that the previous apparent velocity dependence of β is a result of choosing a pure power law tail for ψ(t). The key is the change in the log–log slope of the TPL form of ψ(t) with a shifting observational time window caused by the change in the mean velocity. Hence, the use of the full spectrum of ψ(t) is not only necessary for the transition to Fickian behavior, but also to account for the dynamics of these laboratory observations of non-Fickian transport.  相似文献   

3.
Effects of temporally correlated infiltration on water flow in an unsaturated–saturated system were investigated. Both white noise and exponentially correlated infiltration processes were considered. The moment equations of the pressure head (ψ) were solved numerically to obtain the variance and autocorrelation functions of ψ at 14 observation points. Monte Carlo simulations were conducted to verify the numerical results and to estimate the power spectrum of ψ (S ψψ ). It was found that as the water flows through the system, the variance of the ψ (\( \sigma_{\psi }^{2} \)) were damped by the system: the deeper in the system, the smaller the \( \sigma_{\psi }^{2} \), and the larger the correlation timescale of the infiltration process (λ I ), the larger the \( \sigma_{\psi }^{2} \). The unsaturated–saturated system gradually filters out the short-term fluctuations of ψ and the damping effect is most significant in the upper part of the system. The fluctuations of ψ is non-stationary at early time and becomes stationary as time progresses: the larger the value of λ I , the longer the non-stationary period. The correlation timescale of the ψ (λ ψ ) increases with depth and approaches a constant value at depth: the larger the value of λ I , the larger the value of λ ψ . The results of the estimated S ψψ is consistent with those of the variance and autocorrelation function.  相似文献   

4.
By modelling the observed distribution of210Pb and210Po in surface waters of the Pacific, residence times relative to particulate removal are determined. For the center of the North Pacific gyre these are τPo = 0.6years andτPb = 1.7years. The surface ocean τPb is determined by particulate transport rather than plankton settling. The fact that it is about two orders of magnitude smaller than τPb for the deep ocean implies a sharp change in the adsorptive quality of particles during descent through the water column.  相似文献   

5.
Repeat times of strong shallow mainshocks have been determined by the use of instrumental and historical data for 68 seismogenic sources in the Aegean and surrounding area (34°N–43°N, 18°E–30°E). For 49 of these sources at least two interevent times (three mainshocks) are available for each source. By using the repeat times for these 49 sources the following relation has been determined: $$\log T_t = 0.36M_{\min } + 0.35M_p + a$$ whereT t is the repeat time, measured in years,M p the surface wave magnitude of the preceding mainshock,M min the magnitude of the smallest earthquake considered and “a” parameter which varies from source to source. A multilinear correlation coefficient equal to 0.89 was determined for this relation. By using the same repeat times for the 49 seismogenic sources, the following relation has been determined between the magnitude,M f , of the following mainshock andM min andM p . $$M_f = 0.95M_{\min } - 0.49M_p + m$$ wherem is a constant which varies from source to source. A multilinear correlation coefficient equal to 0.80 was found for this relation. The model expressed by these two relations is represented by a scheme of a time variation of stress under constant tectonic loading. In this scheme, the maximum stress values during the different seismic cycles fluctuate around a value, τ1, in a relatively narrow stress interval, expressing the high correlation coefficient of the relation between LogT andM p . On the contrary, the minimum stress values fluctuate around a value, τ2, in a much broader stress interval. However, each of these minimum stress values becomes lower or higher than τ2 if the previous one is higher or lower than τ2, respectively, expressing the negative correlation betweenM f andM p .  相似文献   

6.
State parameter, ψ, has been widely used to combine the influence of void ratio, e, and stress level, p′, on the soils behavior. Stress ratio, Rs, and modified state parameter, ψm, have also been proposed for the same purpose. This paper aims to evaluate and compare the different state indices in combining the effect of fines content, density and stress level for five different types of sands, by processing a large number of previously published experimental data. The use of the recently established concept of equivalent interparticle void ratio, e, in definition of the state indices is also evaluated. The results indicate that the influence of fines presence, in addition to the e and p′, on the behavior is favorably reflected by the state indices. Unique correlations were derived between the cyclic or monotonic strength and each of the state indices, independent of the fines content. The correlations, for all the different types of soils, fell into limited types of common formulations. ψ and Rs worked generally better than ψm, whether defined in terms of e or e. The extension of straight part of critical state line was found to be an appropriate reference line for calculating Rs used in conjunction with e or e.  相似文献   

7.
8.
Cadmium,zinc, copper,and barium in foraminifera tests   总被引:1,自引:0,他引:1  
The concentrations of cadmium, zinc, copper and barium have been determined on 2-mg samples of single-species foraminifera populations. Cleaning techniques were tested using North Atlantic core tops, and followed by a detailed downcore study for the last 30,000 years in South Atlantic core V22-174. Raw foram tests cleaned by crushing followed by ultrasonic removal of fine-grained material, and dissolved in a pH 5.5 acetate buffer, contain appreciable amounts of trace elements associated with ferromanganese and alumino-silicate contaminants. A reductive/complexing cleaning treatment reduces ferromanganese contamination by 1–2 orders of magnitude. Acetate buffers at pH 5.5 complex iron and raise the solubility of ferromanganese oxides; these buffers are unsuitable for separating carbonate lattice-bound trace elements from the fraction associated with ferromanganese phases. Improved mechanical and ultrasonic reductive cleaning combined with a mild dissolution in distilled water under 1 atm.PCO2 reduces contaminant levels another order of magnitude. The Cd and Zn concentrations (order 10?8 mole Cd/mole Ca and 10?5 mole Zn/mole Ca) of species with low surface area show an increase with decreasing isotopic temperatures. This increase is consistent with the increasing concentrations of these metals from low values in surface waters to higher values at depth. The variance of Cd and Zn over the last 30,000 years in the central South Atlantic is consistent with the probable variability of the dissolved trace elements at the calcification levels of the species analyzed. Cu and Ba are irreproducible and probably sensitive to residual contaminant phases. The trace element content of the tests differs from levels observed in a recent coprecipitation study. Foraminifera may be a significant vector in zinc cycling in the ocean.  相似文献   

9.
The particulate concentrations of 17 trace metals, Al, Sc, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Se, Ag, Sb, Au, Hg, Pb and Th have been measured in the marine atmosphere (58 samples) and in the deep waters (35 samples) of the Tropical North Atlantic. For oceanic suspended matter, our results are similar to those in samples from the Atlantic and the Pacific Oceans collected during the GEOSECS Program. Based on these results, we have made a flux balance for the mixed layer between input via the atmosphere and removal through small and large particles. These data show that the primary flux of suspended aluminosilicates in the Tropical North Atlantic is attributable to the atmospheric input. Elements Sc, Th, Fe, V, Mn, Co and Cr show high correlation with Al in the marine atmosphere. Of these elements, Fe, Mn, V, Co and Cr are influenced by additional processes such as biological, in the marine environment. For elements Ni, Cu, Zn, Se, Ag, Sb, Au, Hg and Pb, we observe high enrichments (relative to average crustal material) in the marine atmosphere which may be due, at least partially, to the influence of anthropogenic sources. These metals also show similar enrichments in deep ocean suspended matter. Model calculations indicate that the atmospheric flux may not control the deep ocean particulate chemistry of Ni, Cu, Zn, Ag, Sb, Au and Hg. Hence it is likely that, for these elements, the enrichment in the ocean is due to processes within the marine regime, for example their involvement in the biological cycle of the ocean. For Se and Pb, the atmospheric source looks to be the dominant contribution to their particulate concentration in seawater. In the deep North Atlantic, particulate Pb appears to be mostly of anthropogenic origin, which is not the case for Se.  相似文献   

10.
This study utilises the equivalent granular state parameter, ψ, as a key parameter for studying static and cyclic instability and their linkage. ψ can be considered as a generalisation of the state parameter as first proposed by Been and Jefferies so that the influence of fines content in addition to stress and density state can be captured. Test results presented in this study conclusively showed that ψ at the start of undrained shearing and ηIS, the stress ratio at onset of static instability, can be described by a single relationship irrespective of fines content for both compression and extension shearing. This single relationship is referred as instability curve. However, the instability curve in extension shearing is different from that of compression. In this paper, the capacity of the instability curve in predicting triggering of cyclic instability was evaluated experimentally. An extensive series of undrained one-way (compression) and non-symmetric two-way cyclic triaxial tests, in addition to monotonic triaxial tests in both compression and extension were conducted for this evaluation. Furthermore, a published database for Hokksund sand with fines was also used. Test results show that cyclic instability was triggered shortly after the cyclic effective stress path crossed the estimated ηIS-zone(s) as obtained from instability curve(s) irrespective of whether instability occurs in the compression or extension side.  相似文献   

11.
In the western part of the Gardar Igneous Province of southern Greenland, lamprophyre dykes intruded at ca. 1276-1254 m.y. RbSr biotite ages yield a palaeomagnetic pole at 206.5°E,3°N (nine sites, dψ = 5.1°, dχ = 10.1°) Slightly younger dolerite dykes with RbSr biotite ages in the range 1278-1263 m.y. give a pole at 201.5°E,8.5°N (24 sites, dψ = 4.7°, dχ = 9.4°), and the syeno-gabbro ring dyke of the Kûngnât complex (RbSr isochron age 1245 ± 17 m.y.) cutting both of these dykes swarms, gives a pole at 198.5°E, 3.5°N (four sites, dψ = 2.3°,dχ = 4.4°). All these rock units have the same polarity and the poles are identical to those from Mackenzie and related igneous rocks of North America (1280-1220 m.y.) after closure of the Davis Strait; they confirm that this part of the Gardar Province is a lateral extension of the Mackenzie igneous episode within the Laurentian craton.In the Tugtutôq region of the eastern part of the Gardar Province 47 NNE-trending dykes of various petrologic types, and intruded between 1175 ± 9 and 1168 ± 37 m.y. (RbSr isochron ages) yield a palaeomagnetic pole at 223.9° E, 36.4°N (dψ = 4.1°, dχ = 6.1°). Fifteen other dykes in this swarm were intruded during a transitional phase of the magnetic field which, however, does not appear to have achieved a complete reversal over a period of several millions of years. The majority of dykes studied are highly stable to AF and thermal demagnetisation and contain single high blocking temperature components with single Curie points in the range 380–560°C.Palaeomagnetic poles from the Gardar Province between ca. 1330 and 1160 m.y. in age define the earlier part of the Great Logan apparent polar-wander loop; they correlate closely with contemporaneous North American results and confirm the coherence of the Laurentian craton in Upper Proterozoic times.  相似文献   

12.
Marine limestones from southern Germany were analysed by the non-destructive neutron activation/Ge(Li) spectrometric method for La, Ce, Sm, Eu, Tb, Yb, Lu, Sr, Ta, U, Sc, Th, Cr, Hf, Cs, Fe and Rb. On the basis of carbonate and non-carbonate portions of the rock, the elements could be divided into (1) those associated with clastic material (Sc, Cr, Hf, etc.), (2) those which were partly associated with authigenic material (REE), and (3) those with pronounced diagenetic history (Sr, U). The calculated effective REE distribution coefficient (Di,eff) between CaCO3 and seawater varied from about 1400 for La to about 460 for Lu suggesting fractionation of the lighter and the heavier REE's in favour of the former.  相似文献   

13.
Coral-reef terraces on tectonically elevated coasts provide important data on variations of sea level. In this paper a systematic approach is presented for the construction of synthetic reef terraces. Idealized reef terraces are generated assuming a constant rate of tectonic uplift ν and a harmonic variation of sea level with amplitudehwo and period τ. It is further assumed that coral terraces grow only during periods of relative sea-level rise. The behavior is governed by the parameter u = τυ/2π hwo. If u > 1 no reef terraces form. If 0.217 < u < 1 a discontinuous series of terraces form on the uplifting basement. If 0 < u < 0.217 the basement is covered with coral and a series of coral terraces forms on this coral. Uplifted coral-reef terraces on the Huon Peninsula, Papua New Guinea, compare reasonably well with terraces generated by our model with τ = 20,000 years andhwo = 21.4 m. We also compare the observed terraces with model terraces using the paleo sea level inferred from18O/16O data. The observed terraces are qualitatively similar to those predicted using the oxygen isotope data, however, the model terraces require larger relative sea-level variations with a near 20,000 year period than those inferred directly from the18O/16O data.  相似文献   

14.
— We have evaluated how the parameters prescribing the slip-dependent constitutive law are affected by temperature and effective normal stress, by conducting the triaxial fracture experiments on Tsukuba-granite samples in seismogenic environments, which correspond to a depth range to 15 km. The normalized critical slip displacement D c almost remains constant below 300oC (insensitive to both temperature and effective normal stress σ n eff); D c increases with increasing temperature above 300 °C, and the rate of D c increase with temperature tends to be largest at higher σ n eff. The breakdown stress drop Δτ b for the granite at constant σ n eff is roughly 80 MPa below 300 °C, and does not depend on σ n eff. Above 300 °C, Δτ b decreases gradually with increasing temperature, and the rate of Δτ b reduction with temperature increases at higher σ n eff. The peak shear strength τ p increases nearly linearly with increasing σ n eff below 300 °C. However, τ p becomes lower above 300 °C, deviating from the linear relation extrapolated from below 300 °C. This is consistent with the onset of crystal plastic deformation mechanisms of Tsukuba granite.  相似文献   

15.
The objective of this work is to analyse the storage–flux hysteretic behaviour of a simplified model for subsurface flow processes. The subsurface flow dynamics is analysed by means of a model based on the kinematic wave assumptions and by using a width weighting/depth averaging scheme which allows to map the three-dimensional soil mantle into a one-dimensional profile. Continuity and a kinematic form of Darcy’s law lead to a hillslope-storage kinematic wave equation for subsurface flow, solvable with the method of characteristics. Adopting a second order polynomial function to describe the bedrock slope and an exponential function to describe the variation of the width of the hillslope with hillslope distance, we derive general solutions to the hillslope-storage kinematic wave equations, applicable to a wide range of hillslopes. These solutions provide a physical basis for deriving two geometric parameters α and ψ which define the hydrological similarity between hillslopes with respect to their characteristic response and hysteresis. The hysteresis η, quantified by the area of the hysteretic dimensionless loop, has been therefore computed for a range of values of parameters α and ψ. Slopes exhibit generally clockwise hysteretic loop in the flux-storage plot, with higher groundwater mean volume for given discharge on rising limb than at same discharge on falling limb. It has been found that hysteresis increases with decreasing α and ψ, i.e. with increasing convergence (for the shape) and concavity (for the profile), and vice versa. For relatively large values of α and ψ the hysteresis may take a complex pattern, with combination of clockwise to anticlockwise loop cycles. Application of the theory to three hillslopes in the Eastern Italian Alps provides an opportunity to examine how natural topographies are represented by the two hillslope hydrological similarity parameters.  相似文献   

16.
In order to investigate the pollution levels, sources and ecological risks of arsenic (As) and heavy metals (Cr, Ni, Cu, Zn, Pb and Cd) in inshore sediments of the Yellow River estuary, the surface sediment in areas of inshore coastal waters were sampled in October 2014 as the flow-sediment regulation project (FSRP) was implemented for 13 years. Results showed that the concentrations of As and heavy metals in inshore sediments of the Yellow River estuary were in the order of Zn?>?Cr?>?Cu?>?Ni?>?Pb?>?As?>?Cd. Higher levels of As, Cr, Ni, Cu, Zn and Pb generally occurred in fine-grained sediments of the Yellow River estuary and the southeast region, which was consistent with the spatial distribution of clay. In contrast, higher concentrations of Cd were generally observed in northwest area of the Yellow River estuary and near the Qingshuigou estuary, which showed similarly spatial distribution with that of sand. The sediment quality guidelines (SQGS) and geoaccumulation indices (Igeo) indicated that the inshore sediments were polluted by Cu, Cd, As, Pb and Zn, and, among them, Cd pollution was more serious. Ecological risk indices (E r i ) demonstrated low risks for Cr, Ni, Cu, Zn, Pb and As, and high potential toxicity by Cd. The integrated ecological risk index implied that 6.8% of stations presented moderate risk, 4.5% of stations exhibited disastrous risk, and 88.7% of stations demonstrated considerable risk. Principal component analysis indicated that Ni, Cu, Zn, Pb and As might originate from common pollution sources, while Cr and Cd might share another similar sources. With the continuous implementation of FSRP, As and heavy metal levels in inshore sediments of the Yellow River estuary could be classified as stage I (2002–2010) and stage II (2010–2014). In the stage I, As, Cr, Ni, Cu, Zn and Pb levels fluctuated but decreased significantly, whereas Cd concentrations showed little variation. In the stage II, As and heavy metal levels significantly increased although some little fluctuations occurred. The continuous accumulation of As and heavy metals (especially for Cd) in inshore sediments of the Yellow River estuary would occur again as the FSRP was implemented for 9 years (since 2010). The ecotoxicological risk of Cd, As, Ni and Cu in inshore sediments might be more serious since the accumulation of the four elements would be continuously occurred in future years. Next step, there will be long-term potential consequences for marine organism if effective measures are not taken to control the loadings of metal pollutants into estuary.  相似文献   

17.
To avoid soil erosion, lters often are installed in the downstream cross section of dams. However, the probability of piping failure and the lter effect both are related to the soil pore network. Previous erosion probability models for a base soil- lter system do not consider the effect of local ow. Therefore, in this study, an improved erosion probability model is established, in which the deviation between the main ow direction of the soil and the local ow direction in the lter was considered based on a previous model. The improved model was validated by numerical simulation. The erosion probability was found to affect the reasonable evaluation of lter effects and the optimal selection of lter thickness. The controlling con- striction size in a granular lter was increased to Dc* Dc5 with a lter thickness of 59 Dm.  相似文献   

18.
Favourable oceanographic and environmental conditions allow the formation of a fine-grained deposit (Tagus prodelta) located at the mouth of the Tagus River. This fine-grained deposit results from the sink and accumulation of terrestrial and marine-derived materials. Three short sediment cores collected in the Tagus prodelta were investigated through the variability in grain-size, major and trace elements, Corg, Ntot, δ13C and 210Pb dating to characterise the historical development of trace metal contamination. Historical trends indicated significant anthropogenic enrichments for Hg, Pb, Zn, Cu, Sb and Sn since the 1930s. Hg presents the highest level of anthropogenic enrichment (ca. EFHg=20) in cores 3576 (PO287-26-1B) and 3579 (PO287-27-1B) collected closer to the shore. Despite these elements were derived from distinct industrial sources, all of them presented similar temporal trends, which points to the importance of estuarine mixing processes originated by the tidal regime and wind before transfer to the adjoining coastal areas. The cores were characterised by smoothness of down-core variations and river flood events occurred in the last decades were not registered. Furthermore, the improvement of effluents treatment together with the closing of some industrial point sources in the last two decades is not evidenced in recent sediment composition of the Tagus prodelta. Besides sediment reworking in the upper sediment layers, sediment mixing inside the estuary may explain minor abrupt alterations with the depth.  相似文献   

19.
Sb, As, Cu, and Zn toxicity and contamination have become a growing concern in recent years. Phytoremediation, a plant based and cost effective technology, may be an effective approach in the cleanup of water contaminated by these metals. In this study, the aquatic macrophyte Eleocharis acicularis was used in laboratory and field experiments to assess its capability to accumulate Sb, As, Cu, and Zn, and thereby investigate its potential application in phytoremediation. The results showed that E. acicularis adapted well to water contaminated by these metals. The removal rates of Sb, As, Cu, and Zn in the laboratory experiment were 3.04, 2.75, 0.417, and 1.49 μg/L/day, respectively. The highest concentrations of these metals accumulated in E. acicularis after 10 days of the laboratory experiment were 6.29, 6.44, 20.5, and 73.5 mg/kg dry weight, respectively. Only 8% of As, 12% of Sb, 87% of Cu and 93% of Zn removed from the water were used by E. acicularis. The highest concentrations of Sb, As, Cu, and Zn accumulated in E. acicularis after 10 wk of the field experiment were 76.0, 22.4, 33.9, and 266 mg/kg dry weight, respectively. The results indicate that E. acicularis has the ability to accumulate Sb, As, Cu, and Zn from contaminated water.  相似文献   

20.
Preliminary data on major elements, Cs, Ba, Rb, Pb, Sr, REE, Y, Th, U, Zr, Ht, Sn, Nb, W, Mo, Cr, V, Sc, Ni, Co and Cu contents for eight samples coming from the Upper Cretaceous volcanic belt of the Pontic Chain (Northern Turkey) are reported. SiO, versus K2O relationship shows that the analyzed samples belong to the calc-alkaline and shoshonite series. The calc-alkaline rocks appear to represent two distinct magma types one close in composition to typical island are calc-alkaline magmas and one with high incompatible elements concentration and tractionated heavy REE patterns which suggest a genesis by partial melting at high pressure with a garnet bearing residue. Shoshonitic rocks show Na2O/K2O close to one, high incompatible elements concentration, and TiO2%. Al2O3%, Ni and Co contents, Ni/Co and V/Ni ratios and REE patterns similar to typical island are andesites which suggest for these rocks similar genetical processes as the island are calc-alkaline magmas.  相似文献   

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

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