首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Debris flow density determined by grain composition   总被引:1,自引:1,他引:0  
Density is one of the most important parameters of debris flows. Because observing an active debris flow is very difficult, finding a method to estimate debris flow density is urgently needed for disaster mitigation engineering. This paper proposes an effective empirical equation in terms of grain size distribution (GSD) parameters based on observations in Jiangjia Gully, Yunnan Province, China. We found that the GSD follows P(D) = KD exp(? D/Dc), with μ and Dc representing the fine and coarse grains, respectively. In particular, μ is associated with some characteristic porosity of soil in the natural state and increases with increased porosity. Dc characterizes the grain size range of the flow and increases with the grain concentration. Studies show that flow density is related to both parameters in power law. Here, we propose an empirical equation for estimating flow density: ρ = 1.26μ -0.132 + 0.049Dc0.443, which provides not only an estimation of the density for a flow, but also describes the variation in density with the GSD of material composition; this provides important information related to the design of debris flow engineering structures.  相似文献   

2.
The equilibrium thickness of the isothermal layers of interstellar gas and volume gas densities ρ gas in the plane of the disk as a function of galactocentric distance R are computed for seven spiral galaxies (including the Milky Way) using an axisymmetrical model. In this model, the thickness of the stellar disk varies with R and remains approximately equal to the minimum thickness of a stable equilibrium disk. We found the disk thickness to increase toward the periphery in at least five of the seven galaxies. The density of the stellar disk decreases with R faster than ρ gas , so that gas dominates at the disk peripheries in terms of density. A comparison of the azimuthally averaged star formation rate SFR and the gas density shows the absence of a universal Schmidt law SFR ~ρ gas n for galaxies. However, the SFRs in various galaxies are better correlated with the volume than the gas surface density. The parameter n in the Schmidt law formally calculated using the least-squares method lies in the interval 0.8–2.4, being, on average, close to 1.5. The values of n calculated separately for the molecular gas display substantial scatter, but are, on average, close to unity. The value of n appears to increase with decreasing ρ gas , so that the fraction of gas that actively participates in star formation decreases with n.  相似文献   

3.
Shear heating by reverse faulting on a sharp straight fault plane is modelled. Increase in temperature (T i ) of faulted hangingwall and footwall blocks by frictional/shear heating for planar rough reverse faults is proportional to the coefficient of friction (μ), density and thickness of the hangingwall block (ρ). T i increases as movement progresses with time. Thermal conductivity (K i ) and thermal diffusivity (\(k_{\mathrm {i}}^{\prime }\)) of faulted blocks govern T i but they do not bear simple relation. T i is significant only near the fault plane. If the lithology is dry and faulting brings adjacent hangingwall and footwall blocks of the same lithology in contact, those blocks undergo the same rate of increase in shear heating per unit area per unit time.  相似文献   

4.
The sound velocity (V P) of liquid Fe–10 wt% Ni and Fe–10 wt% Ni–4 wt% C up to 6.6 GPa was studied using the ultrasonic pulse-echo method combined with synchrotron X-ray techniques. The obtained V P of liquid Fe–Ni is insensitive to temperature, whereas that of liquid Fe–Ni–C tends to decrease with increasing temperature. The V P values of both liquid Fe–Ni and Fe–Ni–C increase with pressure. Alloying with 10 wt% of Ni slightly reduces the V P of liquid Fe, whereas alloying with C is likely to increase the V P. However, a difference in V P between liquid Fe–Ni and Fe–Ni–C becomes to be smaller at higher temperature. By fitting the measured V P data with the Murnaghan equation of state, the adiabatic bulk modulus (K S0) and its pressure derivative (K S ) were obtained to be K S0 = 103 GPa and K S  = 5.7 for liquid Fe–Ni and K S0 = 110 GPa and K S  = 7.6 for liquid Fe–Ni–C. The calculated density of liquid Fe–Ni–C using the obtained elastic parameters was consistent with the density values measured directly using the X-ray computed tomography technique. In the relation between the density (ρ) and sound velocity (V P) at 5 GPa (the lunar core condition), it was found that the effect of alloying Fe with Ni was that ρ increased mildly and V P decreased, whereas the effect of C dissolution was to decrease ρ but increase V P. In contrast, alloying with S significantly reduces both ρ and V P. Therefore, the effects of light elements (C and S) and Ni on the ρ and V P of liquid Fe are quite different under the lunar core conditions, providing a clue to constrain the light element in the lunar core by comparing with lunar seismic data.  相似文献   

5.
The radial dependence of the pseudo phase-space density, ρ(r)/σ 3(r) is studied. We find that the pseudo phase-space density for halos consisting both of dark matter and baryons is approximately a power-law only down to 0.1% of the virial radius while it has a non-power law behavior below the quoted scale, with inner profiles changing with mass. Halos consisting just of dark matter, as the one in dark matter only simulations, are characterized by an approximately power-law behavior. The results argue against universality of the pseudo phase-space density, when the baryons effect are included, and as a consequence argue against universality of density profiles constituted by dark matter and baryons as also discussed in [1].  相似文献   

6.
Channel sand acts as a stratigraphic trap for hydrocarbon accumulation in many parts of the world. Delineation of this type of reservoir is crucial as channel sand may be scarce, and inaccurate location of the drilling wells could lose a huge currency. The Hassi Messaoud (HMD) field was subjected to multiphase tectonic events, where deep-seated structures were rejuvenated leading to intensive fault complexity. The main effective tectonic events upon the studied area are the Hercynian compression and deep erosion till the Ordovician Hamra Quartzite (HQZ) oil reservoir, followed by active Triassic rifting and filling the deeply eroded areas or the graben areas by eruptive volcanic rocks at Triassic time. Hercynian erosion and volcanic rocks distribution introduce a big uncertainty to the reservoir structural model. Amplitude versus offset (AVO) method is used as a helpful tool to differentiate channel sand from surrounding formations. Several attributes (P-impedance, S-impedance, longitudinal velocity Vp, shear velocity Vs and density ρ) are estimated from pre-stack seismic inversion. They have different sensitivity to the reservoir properties. Derived attributes such as Lamé parameters, incompressibility × density (λρ) and rigidity × density (μρ) can provide key lithology and fluid indicators (Goodway et al. 1997, Goodway CSEG Rec 26(6):39-60 2001). Petrophysically relating AVO attributes both to λρ and μρ and to each other in Lambda–Mu–Rho (LMR) cross-plot space can be a good tool for AVO interpretation (Rutherford and Williams Geophysics 54:680–688 1989 and Castagna and Swan Lead Edge 16(4):337–342 1997). After proper data conditioning, simultaneous inversion of pre-stack angle gathers is performed to get acoustic wave impedance (P-impedance), elastic wave impedances (S-impedance) and density ρ, then to calculate λρ and μρ volumes. In the studied area, λρ and μρ are used as a very important key to separate reservoir sands. The λρ and μρ curves are generated at each well location. Cross plots showed a fair separation of sand in the formation, i.e. higher μρ and lower λρ can detect sand. The output λρ and μρ volumes after simultaneous inversion follow the distribution of the sand which is consistent with the wells penetrating the target reservoir. This finding on the extension of the sand reservoir in terms of λρ and μρ. 3D cross-plot zonations are used for lithology discrimination. In this study, well logs were used to constrain lithology and to control the zonation filters by reducing the limits ambiguity. Other types of advanced attributes are calculated and tested. The obtained (μρλρ) volume acts as a good indicator for the sand distribution. It was finally used as sand presence index in the area. Also μρ has shown a good linear relationship with porosity. To note that the porosity volume is created based on the linear relationship with μρ. A product of derived porosity and the sand presence index (μρλρ) provides a good tool for reservoir characterization and lead to reservoir management, future planning of the field, and setting location for new wells.  相似文献   

7.
Fine sediment inputs can alter estuarine ecosystem structure and function. However, natural variations in the processes that regulate sediment transport make it difficult to predict their fate. In this study, sediments were sampled at different times (2011–2012) from 45 points across intertidal sandflat transects in three New Zealand estuaries (Whitford, Whangamata, and Kawhia) encompassing a wide range in mud (≤63 μm) content (0–56 %) and macrofaunal community structure. Using a core-based erosion measurement device (EROMES), we calculated three distinct measures of sediment erosion potential: erosion threshold (? c ; N m?2), erosion rate (ER; g m?2 s?1), and change in erosion rate with increasing bed shear stress (m e ; g N?1 s?1). Collectively, these measures characterized surface (? c and ER) and sub-surface (m e ) erosion. Benthic macrofauna were grouped by functional traits (size and motility) and data pooled across estuaries to determine relationships between abiotic (mud content, mean grain size) and biotic (benthic macrofauna, microbial biomass) variables and erosion measures. Results indicated that small bioturbating macrofauna (predominantly freely motile species <5 mm in size) destabilized surface sediments, explaining 23 % of the variation in ? c (p ≤ 0.01) and 59 % of the variation in ER (p ≤ 0.01). Alternatively, mud content and mean grain size cumulatively explained 61 % of the variation in m e (p ≤ 0.01), where increasing mud and grain size stabilized sub-surface sediments. These results highlight that the importance of biotic and abiotic predictors vary with erosion stage and that functional group classifications are a useful way to determine the impact of benthic macrofauna on sediment erodibility across communities with different species composition.  相似文献   

8.
The magnetic-field structure in regions of stationary, planar accretion disks around active galactic nuclei where general-relativistic effects can be neglected (from 10 to 200 gravitational radii) is considered. It is assumed that the magnetic field in the outer edges of the disk, which forms in the magnetosphere of the central black hole during the creation of the relativisitic jets, corresponds to the field of a magnetic dipole perpendicular to the plane of the disk. In this case, the azimuthal field component Bφ in the disk arises due to the presence of the radial field Bρ and the azimuthal velocity component Uφ. The value of the magnetic field at the inner radius of the disk is taken to correspond to the solution of the induction equation in a diffusion approximation. Numerical solutions of the induction equation are given for a number of cases.  相似文献   

9.
Equations of motion containing a small parameter μ are derived for stars at the peripheries of open star clusters. The parameter μ is obtained for six numerical open-cluster models. The general analytical solution of these equations of motion for μ = 0 is found. An iterative method is used to derive the frequencies of the stellar motions for first-order expansions in μ of the solutions of the equations of motion for stars at the cluster periphery. Applications of the results are discussed.  相似文献   

10.
Stable isotope data of precipitation (δ18Op and deuterium excess), drip water (δ18Od), and modern calcite precipitates (δ18Oc and δ13Cc) from Yongxing Cave, central China, are presented, with monthly sampling intervals from June 2013 to September 2016. Moderate correlations between the monthly variation of δ18Op values (from ??11.5 to ??0.7‰) and precipitation amount (r = ??0.59, n?=?34, p?<?0.01) and deuterium excess (r?=?0.39, n?=?31, p?<?0.01) imply a combined effect of changes in precipitation amount and atmospheric circulation. At five drip sites, the δ18Od values have a much smaller variability (from ??9.1 to ??7.5‰), without seasonal signals, probably a consequence of the mixing in the karst reservoir with a deep aquifer. The mean δ18Od value (??8.4‰) for all drip waters is significantly more negative than the mean δ18Op value (??6.9‰) weighted by precipitation amount, but close to the wet season (May to September) mean value (??8.3‰), suggesting that a threshold of precipitation amount must be exceeded to provide recharge. Calculation based on the equilibrium fractionation factor indicates that the δ18Oc values are not in isotopic equilibrium with their corresponding drip waters, with a range of disequilibrium effects from 0.4 to 1.4‰. The δ18Oc and δ13Cc values generally increase progressively away from the locus of precipitation on glass plates. The disequilibrium effects in the cave are likely caused by progressive calcite precipitation and CO2 degassing related to a high gradient of CO2 concentration between drip waters and cave air. Our study provides an important reference to interpret δ18Oc records from the monsoon region of China.  相似文献   

11.
The excitation of methanol in the absence of external radiation is analyzed, and LTE methods for probing interstellar gas considered. It is shown that rotation diagrams correctly estimate the gas kinetic temperature only if they are constructed using lines whose upper levels are located in the same K-ladders, such as the J0?J?1E lines at 157 GHz, the J1?J0E lines at 165 GHz, and the J2?J1E lines at 25 GHz. The gas density must be no less than 107 cm?3. Rotation diagrams constructed from lines with different K values for their upper levels (e.g., 2K?1K at 96 GHz, 3K?2K at 145 GHz, 5K?4K at 241 GHz) significantly underestimate the temperature, but enable estimation of the density. In addition, diagrams based on the 2K?1K lines can be used to estimate the methanol column density within a factor of about two to five. It is suggested that rotation diagrams should be used in the following manner. First, two rotation diagrams should be constructed, one from the lines at 96, 145, or 241 GHz, and another from the lines at 157, 165, or 25 GHz. The former diagram is used to estimate the gas density. If the density is about 107 cm?3 or higher, the latter diagram reproduces the temperature fairly well. If the density is around 106 cm?3, the temperature obtained from the latter diagram should be multiplied by a factor of 1.5–2. If the density is about 105 cm?3 or lower, then the latter diagram yields a temperature that is lower than the kinetic temperature by a factor of three or more, and should be used only as a lower limit for the kinetic temperature. The errors in the methanol column density determined from the integrated intensity of a single line can be more than an order of magnitude, even when the gas temperature is well known. However, if the J0?(J ? 1)0E lines, as well as the J1?(J ? 1)1A+ or A? lines are used, the relative error in the column density is no more than a factor of a few.  相似文献   

12.
The He/H abundances in HII regions in Blue Compact Dwarf Galaxies are redetermined using the new recombination coefficients of Benjamin et al. The electron number density n e in the He+ region, optical depth τ3889 in the HeI λ3889 Å line, and coefficient of underlying stellar absorption aHeI are determined using the self-consistent method of Olive and Skillman. The primordial helium abundance and its enrichment are found using the helium abundances obtained in this paper and heavy-element abundances from our recent paper I: Y p = 0.244 ± 0.004 and dY/dZ = 8.8 ± 4.6.  相似文献   

13.
Petrophysical evaluation and rock physics analysis are the important tools to relate the reservoir properties like porosity, permeability, pore fluids with seismic parameters. Nevertheless, the uncertainties always exist in the quantification of elastic and seismic parameters estimated through wireline logs and rock physics analysis. A workflow based on statistical relationships of rock physics and logs derived elastic and seismic parameters with porosity and the percentage error exist between them is given. The statistical linear regressions are developed for early Eocene Chorgali Formation between various petrophysically factors determined from borehole logging of well Ratana–03 drilled in tectonically disturbed zone and the seismic and elastic parameters estimated through rock physics modeling. The rock physics constraints such as seismic velocities, effective density and elastic moduli calculated from Gassmann fluid substation analysis are in harmony and close agreement to those estimated from borehole logs. The percentage errors between well logs and rock physics computed saturated bulk modulus (K sat ), effective density (ρ eff ), compressional and shear wave velocities (V P and V S) are 1.31%, 4.23 %, 5.25% and 4.01% respectively. The permeability of reservoir intervals show fairly strong linear relationship with the porosity, indicating that the reservoir interval of the Chorgali Formation is permeable and porous thus having large potential of hydrocarbon accumulation and production.  相似文献   

14.
A high-spatial resolution study design was used to investigate the relationship between land use practices, stream physicochemistry, hydroclimate, and stream Escherichia (E) coli concentrations in a mixed-land-use watershed in the Appalachian region. Stream samples were collected daily from six monitoring sites and analyzed for total E. coli counts using an enzyme metabolism indicator method. Statistical comparison of E. coli concentration time series showed significant (p?<?0.05) differences between study sites. Although highest average E. coli concentrations were observed at two agricultural sites (534 and 582 colony-forming counts (CFU) per 100 mL, respectively), highest total loadings were observed within the receiving stream, with values increasing downstream (2?×?1012 and 4.2?×?1012 study total CFU for bracketed upstream and downstream sites, respectively). No single physical variable displayed a significant correlation (p?<?0.05) with observed E. coli concentration at every site. However, sites displayed different patterns of significant correlations (p?<?0.05) between E. coli concentration and both physicochemical (e.g. pH, dissolved oxygen saturation) and hydroclimate variables (e.g. streamflow and precipitation). Percent agricultural land cover was the only land use category that showed significant (p?<?0.04) correlation with study average E. coli concentrations, thereby emphasizing the importance of land use practices to stream pathogen regimes. Results validate the analytical method and provide high-resolution, detailed, quantitative characterizations of stream E. coli regimes, thereby supplying land and water resource managers with science-based information to advance management decisions and improve public health.  相似文献   

15.
The electrical conductivity of aqueous fluids containing 0.01, 0.1, and 1 M NaCl was measured in an externally heated diamond cell to 600 °C and 1 GPa. These measurements therefore more than double the pressure range of previous data and extend it to higher NaCl concentrations relevant for crustal and mantle fluids. Electrical conductivity was generally found to increase with pressure and fluid salinity. The conductivity increase observed upon variation of NaCl concentration from 0.1 to 1 M was smaller than from 0.01 to 0.1 M, which reflects the reduced degree of dissociation at high NaCl concentration. Measured conductivities can be reproduced (R 2 = 0.96) by a numerical model with log \(\sigma\) = ?1.7060– 93.78/T + 0.8075 log c + 3.0781 log \(\rho\) + log \(\varLambda\) 0(T, \(\rho\)), where \(\sigma\) is the conductivity in S m?1, T is temperature in K, c is NaCl concentration in wt%, \(\rho\) is the density of pure water (in g/cm3) at given pressure and temperature, and \(\varLambda\) 0 (T, \(\rho\)) is the molar conductivity of NaCl in water at infinite dilution (in S cm2 mol?1), \(\varLambda\) 0 = 1573–1212 \(\rho\) + 537 062/T–208 122 721/T 2. This model allows accurate predictions of the conductivity of saline fluids throughout most of the crust and upper mantle; it should not be used at temperatures below 100 °C. In general, the data show that already a very small fraction of NaCl-bearing aqueous fluid in the deep crust is sufficient to enhance bulk conductivities to values that would be expected for a high degree of partial melting. Accordingly, aqueous fluids may be distinguished from hydrous melts by comparing magnetotelluric and seismic data. H2O–NaCl fluids may enhance electrical conductivities in the deep crust with little disturbance of v p or v p/v s ratios. However, at the high temperatures in the mantle wedge above subduction zones, the conductivity of hydrous basaltic melts and saline aqueous fluids is rather similar, so that distinguishing these two phases from conductivity data alone is difficult. Observed conductivities in forearc regions, where temperatures are too low to allow melting, may be accounted for by not more than 1 wt% of an aqueous fluid with 5 wt% NaCl, if this fluid forms a continuous film or fills interconnected tubes.  相似文献   

16.
A technique for IR spectroscopic determination of the total nitrogen content N S in the form of A-and B 1-defects is suggested. It provides for the computer processing and decomposition of IR spectra into constituent bands, calculation of the total absorption band area S N and individual areas S A and S B1 and their normalization with respect to the total area of the diamond intrinsic absorption S 0, with the normalization coefficients K S , K A , and K B1 being calculated. Based on the analysis of the IR spectra of 60 octahedral diamond crystals from the Mir and Yubileinaya pipes (Sakha-Yakutiya), the empirical functions N S = 911.85 K S 0.9919 ppm (R 2 = 0.9859), N A = 1185.6 K A 1.1511 ppm (R 2 = 0.8703), and N B1 = 911.85 K S 0.9919 ? 1185.6 K A 1.1511 ppm have been defined.  相似文献   

17.
The isotopic systematics of noble gases (He and Ar) were studied in Neoarchean and Paleoproterozoic lower crustal xenoliths from the Belomorian mobile belt. The xenolith suite is dominated by garnet granulites (Grt + Cpx + Pl ± Opx ± Qtz ± Kfs ± Phl ± Hbl) and two-pyroxene or garnet pyroxenites (Cpx + Pl ± Grt ± Opx ± Hbl ± Qtz). The xenoliths and the host Devonian ultramafic lamprophyres forming diatremes and explosion dikes contain fluid with similar He and Ar isotopic compositions. It was found that the fluid was trapped by the rocks and xenolith minerals approximately simultaneously with the formation of the lamprophyres. This conclusion is based on the identical K-Ar ages of the majority of xenoliths and ultramafic lamprophyres. When the xenoliths were transported toward the surface by high-temperature ultrabasic melts, the noble gases occurring in them were partly (Ar) or completely (He) lost. The melts were contaminated by meteoric waters during their emplacement in the upper crust, which resulted in that the Ar isotopic composition of lamprophyres approached the composition of atmospheric Ar. The fluid phase that was liberated during melt crystallization severely affected the xenoliths, diminishing the difference between the isotope compositions of He and Ar in the xenoliths and ultramafic lamprophyres. The He isotope composition includes an admixture of mantle 3He, which is suggested by the high measured 3He concentrations, exceeding the calculated values, and high 3He/36Ar ratios in the xenoliths and their host lamprophyres. The fraction of the mantle component in the fluid trapped by the xenoliths was estimated as ~20%.  相似文献   

18.
There is a need for research that advances understanding of flow alterations in contemporary watersheds where natural and anthropogenic interactions can confound mitigation efforts. Event-based flow frequency, timing, magnitude, and rate of change were quantified at five-site nested gauging sites in a representative mixed-land-use watershed of the central USA. Statistically independent storms were paired by site (n = 111 × 5 sites) to test for significant differences in event-based rainfall and flow response variables (n = 17) between gauging sites. Increased frequency of small peak flow events (i.e., 64 more events less than 4.0 m3 s?1) was observed at the rural–urban interface of the watershed. Differences in flow response were apparent during drier periods when small rainfall events resulted in increased flow response at urban sites in the lower reaches. Relationships between rainfall and peak flow were stronger with decreased pasture/crop land use and increased urban land use by approximately 20%. Event-based total rainfall explained 40–68% of the variance in peak flow (p < 0.001). Coefficients of determination (r2) were negatively correlated with pasture/crop land use (r2 = 0.92; p = 0.007; n = 5) and positively correlated with urban land use (r2 = 0.90; p = 0.008; n = 5). Significant differences in flow metrics were observed between rural and urban sites (p < 0.05; n = 111) that were not explained by differences in rainfall variables and drainage area. An urban influence on flow timing was observed using median time lag to peak centroid and time of maximum precipitation to peak flow. Results highlight the need to establish manageable flow targets in rapidly urbanizing mixed-land-use watersheds.  相似文献   

19.
The results of hydrodynamical calculations of radially pulsating helium stars with masses 0.5MM≤0.9M, bolometric luminosities 600L≤5×103L, and effective temperatures 1.5×104 K≤Teff≤3.5×104 K are presented. The pulsation instability of these stars is due to the effects of ionization of iron-group elements in layers with temperatures T~2×105 K. The calculations were carried out using opacities for the relative mass abundances of hydrogen and heavy elements X=0 and Z=0.01, 0.015, and 0.02. Approximate formulas for the pulsation constant Q over the entire range of pulsation instability of the hot helium stars in terms of the mass M, radius R, effective temperature Teff, and heavy-element abundance Z are derived. The instability of BX Cir to radial pulsations with the observed period Π=0.1066 d occurs only for a mass M≥0.55M, effective temperature Teff≥23000 K, and heavy-element abundance Z≥0.015. The allowed mass of BX Cir is in the range 0.55MM≤0.8M, which corresponds to luminosities 800LM≤1400L and mean radii 1.7R?R?2.1R.  相似文献   

20.
Thermal inertia (I) is an important parameter in the Earth’s thermal study. There is no doubt that correct and up-to-date knowledge of thermal inertia particularly as it is affected by the land use/cover will provide good and useful information to agriculturists and environmental scientists. In this work, thermal inertia of Abeokuta City of Ogun State, southwestern Nigeria, was determined. Map of the study area was gridded using 2-min resolution which gave 5 × 5 sampling points from where core samples were collected. Bulk density (ρ s ), thermal conductivity (λ s ), and heat capacity (C s ) of the samples were determined. The pattern of variation of the city view of I (in the order of ×103 Jm?2 s?1/2 K?1) showed that the main urban built-up part of the study area, Abeokuta South, had the highest I mean, 0.76160 with standard deviation of 0.032547 and standard error of 0.01455541. The trend also varied along each sampling latitudinal line. It ranges between 0.509 and 0.756 on latitude 7° 14′, 0.557 and 0.768 on latitude 7° 12′, 0.642 and 0.782 on latitude 7° 10′, 0.7 and 0.794 on latitude 7° 08′, and between 0.642 and 0.728 on latitude 7° 06′. Hence, we have estimated thermal inertia of Abeokuta using the thermophysical properties of the study area. Thermal inertia tends to gradually decrease with distance from the city center (Abeokuta South). This result is promising in possible future consideration of urban ground heat energy conversion to other forms of energy.  相似文献   

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

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