首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 203 毫秒
1.
For feasibility studies and preliminary design estimates, field measurements of shear wave velocity, V s, may not be economically adequate and empirical correlations between V s and more available penetration measurements such as cone penetration test, CPT, data turn out to be potentially valuable at least for initial evaluation of the small-strain stiffness of soils. These types of correlations between geophysical (Vs) and geotechnical (N-SPT, q c-CPT) measurements are also of utmost importance where a great precision in the calculation of the deposit response is required such as in liquefaction evaluation or earthquake ground response analyses. In this study, the stress-normalized shear wave velocity V s1 (in m/s) is defined as statistical functions of the normalized dimensionless resistance, Q tn-CPT, and the mean effective diameter, D 50 (in mm), using a data set of different uncemented soils of Holocene age accumulated at various sites in North America, Europe, and Asia. The V s1Q tn data exhibit different trends with respect to grain sizes. For soils with mean grain size (D 50) < 0.2 mm, the V s1/Q tn 0.25 ratio undergoes a significant reduction with the increase in D 50 of the soil. This trend is completely reversed with further increase in D 50 (D 50 > 0.2 mm). These results corroborate earlier results that stressed the use of different CPT-based correlations with different soil types, and those emphasized the need to impose particle-size limits on the validity of the majority of available correlations.  相似文献   

2.
The evaluation of the undrained cyclic resistance of sandy deposits is required to forecast the soil behaviour during an earthquake (liquefaction, cyclic mobility); due to the difficulties in obtaining undisturbed samples of most liquefiable soils, it is usually deduced from field test results such as cone penetration tests. This paper proposes a methodology to evaluate the undrained cyclic resistance from normalised cone resistance of two well-studied silica sands (Ticino and Toyoura), with different mineralogy, one mainly composed of feldspar, the other of quartz. The determination of the cyclic resistance of Ticino and Toyoura sands was achieved through undrained cyclic triaxial tests on reconstituted specimens. The tip resistance was deduced from CPTs performed in centrifuge with a miniaturised piezocone on homogeneous reconstituted models. Both the undrained cyclic and tip resistances were correlated with the state parameter ψ. Results of centrifuge and triaxial tests were combined through ψ to deduce the cyclic resistance ratio CRR directly from the normalised cone resistance. The shape of the curve relating CRR to the normalised cone resistance resulted unusual respect to all the recognised curves widespread in the geotechnical literature. The aim of the proposed correlations is to provide a useful instrument to improve the actual knowledge on liquefaction and to give a contribution based on the critical state soil mechanics framework to the development of refined correlations between the cyclic resistance of a sand and the results of cone penetration tests.  相似文献   

3.
Net to gross (NTG) is an important parameter in evaluating the volumes of hydrocarbon in place (VHIP) in reservoirs, and proper evaluation of this parameter will lead to significant accuracy in the estimated reserves. The use of conventional petrophysical log evaluations often does not provide sufficient resolution for net sand analysis of reservoirs, especially in laminated reservoir rocks. The uncertainties associated with log-derived net-to-gross estimates arise from the petrophysical shale volume (V sh) cutoffs used in deriving the net sand count over a reservoir interval. One way of improving the net-to-gross estimates is by using a model that calibrates log-derived net to gross to a core equivalent, which is the more accurate representation of net sand counts. In this study, a model for calibrating log-derived net to gross to a core equivalent based on genetic units and facies associations of reservoir rocks in three wells (wells 7, 36, and K008) was established. Ultraviolet core photographs show a good contrast between hydrocarbon-stained sandstones and shales, and combining it with white light-slabbed core images facilitated a manual net sand count of the core photographs on a bed-by-bed basis. Petrophysical shale volume (V sh) cutoff derived from two volumes of shale indicators was applied to generate net sand counts, which were used to get log-derived net-to-gross values. Then, the net to gross from core images and petrophysical clay volume (V sh) analysis were compared by facies associations, and this comparison yielded a reliable core-calibrated net-to-gross model, which reduces the uncertainties in net-to-gross values estimated from V sh cutoffs. The results show that for the distributary channel sands, net to gross derived using V sh cutoff results in an underestimation of net to gross by about 6–10 % when compared with the core-calibrated net to gross, while for the upper shoreface units, V sh cutoff overestimates net to gross by about 7–10 % when compared to its core-calibrated equivalent.  相似文献   

4.
The purpose of this study is to quantify the magnitudes of surface conduction and pore water conduction from the measured electrical conductivity of kaolinite, with the ultimate goal of estimating the electrical conductivity of kaolinite with a wide range of pore water conductivities (σ w = 0.013–3.356 S/m) and porosities (n = 0.368–1.0). Therefore, the theoretical background of the electrical conductivity in soils was reviewed, and electrical conductivity measurements on kaolinite were performed using both slurry and consolidation tests in this study. The results of this study demonstrate that the variations of measured electrical conductivity (σ mix) with n are debatable according to the values of σ w, because a decrease in n results in both an increase in surface conduction (K s) and a decrease in pore water conduction (K w); this causes the relative magnitude of K s compared to that of K w to vary with σ w and n. Consequently, this study develops the relation between the porosity-normalized K s/K w and 1/σ w. Additionally, the surface conductivity of the tested kaolinite is back-calculated and compared with the previous relationship between K s and zeta potential of kaolinite. The measured and estimated σ mix values are compared with the varying pore water conductivity and porosity values.  相似文献   

5.
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.  相似文献   

6.
This paper presents side-by-side comparisons of blowcount values for the Texas cone penetration (TCP) test and the standard penetration test (SPT). The comparisons yielded statistically-significant regression models for both coarse-grained soils and fine-grained soils. Consistent with expected trends and published data, the TCP–SPT relationship is nonlinear, with weak to fair correlation strength (R2 = 23–44%). For TCP blowcounts (N60, TCP) varying from 25 to 200 blows/30 cm (1 ft), corresponding SPT blowcounts (N60, SPT) are typically 30–60% lower than N60, TCP in fine-grained soils. Likewise, corresponding N60, SPT blowcounts are 10–70% lower than N60, TCP in coarse-grained soils, all other things being equal. Comparative data were obtained from published sources and from project-specific field research sites used for full-scale deep foundation load tests. The final dataset consisted of 225 test pairs obtained in similar soils and geomaterials, at equivalent depths, with all blowcounts normalized to 30 cm (12 in.) penetration (i.e., blows/30 cm or blows/ft) within the bounds of typical test precision, and corrected to 60% hammer efficiency. The generally weak correlations do not support conversion of N60, TCP to N60, SPT (or vice versa) to compute foundation capacity for final design. But, engineers can certainly get an intuitive feel about site conditions and preliminary foundation capacity by using the correlation equations to translate their knowledge of one test to the other. This study extends previous work by formally comparing and contrasting the similar yet different SPT and TCP test methods in such a way as to make the results useful to users of both tests and to the broader geotechnical engineering community.  相似文献   

7.
The South Jingyang Plateau, with a total area of 70 km2, is located in Shaanxi Province, China. Since 1976, more than 50 landslides of different types have occurred repeatedly on the edge slopes of the plateau due to the start of diversion irrigation on the plateau, resulting in great loss of lives and property. To better understand the initiation and movement mechanisms of these loess landslides, we surveyed them and carried out a detailed investigation of a large landslide in the Xihetan area. Our field survey results revealed that although most of these landslides had a long runout with high mobility, most of the landslide materials originating from the edge slopes may have been in an unsaturated state when the landslide occurred. This suggests that the materials at the toe of the edge slope as well as on the travel path along the river terrace might have played a key role in landslide movement. To examine how the materials on the travel path were involved in the landsliding, we used a multichannel surface wave technique and surveyed shear wave velocity (V s ) profiles of the landslide deposits. We also examined the internal geometry of the deposits that outcropped on the right-side slope of the landslide foot. The longitudinal profile of V s along the direction of movement showed that terrace deposits near the toe of the edge slope may have been sheared upward, indicating that at the toe, the surface of rupture might be located inside the terrace deposits. The V s contours showed an A-shaped fold within the landslide deposits in the middle part of the travel path and became greater in the most distal toe part. The V s profile across the deposits showed a U-shaped belt, in which the soil layers have smaller V s . This belt may be the boundary between the sliding landslide debris and terrace deposits. The observed internal geometry of the landslide deposits indicates that a sliding surface developed within the sandy layer underlying the gravel layer. Therefore, we inferred that after failure, the displaced landslide materials overrode and sheared the terrace deposits along its main sliding direction, resulting in the formation of thrust folds within the terrace deposits, and greater V s on the distal toe part of the landslide.  相似文献   

8.
This study quantifies the influence of various intrinsic soil properties including particle roundness, R, sphericity, S, 50% size by weight, D 50, coefficient of uniformity, C u, and the state property of relative density, D r, on the compression and recompression indices, C c and C r, of sands of various geologic origins at pre-crushing stress levels. Twenty-four sands exhibiting a wide range of particle shapes, gradations, and geologic origins were collected for the study. The particle shapes were determined using a computational geometry algorithm which allows characterization of a statistically large number of particles in specimens. One dimensional oedometer tests were performed on the soils. The new data was augmented with many previously published results. Through statistical analyses, simple functional relationships are developed for C c and C r. In both cases, the models utilized only R and D r since other intrinsic properties proved to have lesser direct influence on the compression indices. However, previous studies showed that the contributions of S and C u are felt through their effects on index packing void ratios and thus on D r. The accuracy of the models was confirmed by comparison of predicted and observed C c and C r values.  相似文献   

9.
Estimates of the masses of supermassive black holes (M bh ) in the nuclei of disk galaxies with known rotation curves are compared with estimates of the rotational velocities V m and the “indicative” masses of the galaxies M i . Although there is a correlation between M bh and V m or M i , it is appreciably weaker than the correlation with the central velocity dispersion. The values of M bh for early-type galaxies (S0-Sab), which have more massive bulges, are, on average, higher than the values for late-type galaxies with the same rotational velocities. We conclude that the black-hole masses are determined primarily by the properties of the bulge and not the rotational velocity or the mass of the galaxy.  相似文献   

10.
Before starting seismic cycle of Ahar–Varzaghan 2012 event, a partial gap in the form of a pre-seismic calm sequence (seismicity rate, r = 0.46 event/year, b = 1.4) with duration of 303 days spatially has dominated over the entire seismogenic area. From April 17, 2012, to May 31, 2012, r significantly increased to 2.16, indicating strong foreshock sequence, and b value changed to 1.9, remarkably. In the last two months before the mainshock, foreshocks have partially migrated toward the earthquake fault (with a decrease in size, b = 2.0). Significantly, high rate of seismicity and low V P /V S (1.64) in the foreshocks sequence and also very high seismicity rate (17.3) and high V P /V S (1.76) in the aftershocks sequence make substantial differences between the seismic cycle and the background seismicity. Moreover, a significant E–W migration of the microseismicity was confirmed in the study area.  相似文献   

11.
Multi-epoch observations with high spectral resolution acquired in 1998–2008 are used to study the time behavior of the spectral-line profiles and velocity fields in the atmosphere and circumstellar shell of the post-AGB star V448 Lac. Asymmetry of the profiles of the strongest absorption lines with lower-level excitation potentials χ low < 1 eV and time variations of these profiles have been detected, most prominently the profiles of the resonance lines of BaII, YII, LaII, SiII. The peculiarities of these profiles can be explained using a superposition of stellar absorption line and shell emission lines. Emission in the (0; 1) 5635 Å Swan system band of the C2 molecule has been detected in the spectrum of V448 Lac for the first time. The core of the Hα line displays radial-velocity variations with an amplitude of ΔV r ≈ 8 km/s. Radial-velocity variations displayed by weak metallic lines with lower amplitudes, ΔV r ≈ 1–2 km/s, may be due to atmospheric pulsations. Differential line shifts, ΔV r = 0–8 km/s have been detected on various dates. The position of the molecular spectrum is stationary in time, indicating a constant expansion velocity of the circumstellar shell, V exp = 15.2 km/s, as derived from the C2 and NaI lines.  相似文献   

12.
We analyze previously published geodetic data and intensity values for the M s = 8.1 Shillong (1897), M s = 7.8 Kangra (1905), and M s = 8.2 Nepal/Bihar (1934) earthquakes to investigate the rupture zones of these earthquakes as well as the amplification of ground motions throughout the Punjab, Ganges and Brahmaputra valleys. For each earthquake we subtract the observed MSK intensities from a synthetic intensity derived from an inferred planar rupture model of the earthquake, combined with an attenuation function derived from instrumentally recorded earthquakes. The resulting residuals are contoured to identify regions of anomalous intensity caused primarily by local site effects. Observations indicative of liquefaction are treated separately from other indications of shaking severity lest they inflate inferred residual shaking estimates. Despite this precaution we find that intensites are 1–3 units higher near the major rivers, as well as at the edges of the Ganges basin. We find evidence for a post-critical Moho reflection from the 1897 and 1905 earthquakes that raises intensities 1–2 units at distances of the order of 150 km from the rupture zone, and we find that the 1905 earthquake triggered a substantial subsequent earthquake at Dehra Dun, at a distance of approximately 150 km. Four or more M = 8 earthquakes are apparently overdue in the region based on seismic moment summation in the past 500 years. Results from the current study permit anticipated intensities in these future earthquakes to be refined to incorporate site effects derived from dense macroseismic data.  相似文献   

13.
Copper Adsorption by Chernozem Soils and Parent Rocks in Southern Russia   总被引:1,自引:0,他引:1  
Laboratory data in Cu2+ adsorption by chernozems and parent rocks in Rostov region show that adsorption isotherms can be approximated by the Langmuir equation, whose parameters (Kl and C) were calculated for all of the samples. The values of C show a strong negative correlation with the values of cationexchange capacity (CEC) (r =–0.88 at Р = 0.95), and Kl is correlated with the content of physical clay (particles <0.01 mm) (r = 0.78) and with clay (particles <0.001 mm) content in ordinary chernozem and southern chernozems of various particle size distribution (r = 0.80). Even stronger correlations were detected between these parameters in southern chernozems (r = 0.89 for the physical clay (PC) and r = 0.91 for the silt). However, none of the samples displays a significant correlation of C and Kl with the contents of physical clay and silt. This led us to conclude that the composition of the samples, for example, their organic matter, can affect Cu2+ adsorption by the soils and parent rocks. Acidification mechanisms of the equilibrium solutions during the Cu2+ adsorption by soils are discussed, as also are the reasons for the absence of balance between Cu2+ adsorbed by soils and exchangeable cations transferred into solution. Analysis of the fine structures of the XANES and EXAFS spectra suggests that Cu2+ can form coordinated chelate complex compounds with humic acids (HA) of soils and can substitute Al3+ at octahedral sites when interacting with clay minerals in soils.  相似文献   

14.
Radio interferometric observations of an H2O maser flare in the Orion Nebula at epoch 1982.9 have been used to determine the flare’s spatial structure. Antennas in the Crimea, Effelsberg, and Onsala were used. The emission region consists of three groups of components. The angular sizes of the components are 0.2–0.9 mas, and the widths of the emitted lines are 0.2–0.7 km/s. The velocities of the components are correlated with their relative positions, which correspond to expanding concentric rings. Assuming a 1 M protostar in a Keplerian approximation, the radius of the inner ring R is 15 AU, the velocity of its rotation Vrot is 8.98 km/s, and the radial component of the velocity Vrad is 1.79 km/s. For the outer ring, R=15.7 AU, Vrot=8.79 km/s, and Vrad=2.61 km/s.  相似文献   

15.
We have determined the main parameters of the old precataclysmic variable stars MS Peg and LM Com. The radial velocities of the components, reflection effects in the spectra, and light curves of the systems are studied based on model stellar atmospheres subject to external irradiation. Forty-seven moderate-resolution spectra for MS Peg and 57 for LM Com obtained with the 6-m telescope of the Special Astrophysical Observatory are used to derive the refined orbital periods of 0.1736660 days and 0.2586873 days, respectively; the orbital eccentricities do not exceed e=0.04. The mass (M w =0.49e) and radius (e w =0.015R) of the MS Peg primary calculated using the gravitational redshift correspond to those for a cooling carbon white dwarf with a thin hydrogen envelope. The parameters of the red dwarf (M r =0.19M, Teff=3560 K, R r =0.18R) are close to those derived from evolutionary tracks for main-sequence M stars with solar chemical composition. The radius (R r =0.22R) and temperature (Teff=3650 K) of the LM Com secondary exceed theoretical estimates for main-sequence stars with masses of M r =0.17M. The luminosity excess of the red dwarf in LM Com can be explained by a prolonged (T>5×106 yrs) relaxation of the M star to its normal state after the binary leaves the common-envelope stage. For both systems, theoretical U, B, V, and R light curves and spectra calculated using the adopted sets of parameters are generally consistent with the observations. This confirms the radiative origin of the hot spots, the unimportance of horizontal radiative transport, and the absence of large-scale velocity fields with high values (Vtrans>50 km/s) at the surfaces of the secondaries. Most of the emission lines in the spectra of these objects are formed under conditions close to thermalization, enabling modeling of their pro files in an LTE approximation. A strong λ3905 Å emission line has been identified as the 3s23p4s 1P0-3s23p2 1S SiI λ3905.52 Å line formed in the atmosphere of the hot spot. The observed intensity can be explained by non-LTE “superionization” of SiI atoms by soft UV radiation from the white dwarf. We suggest a technique for identifying binaries whose cool components are subject to UV irradiation based on observations of λ3905 Å emission in their spectra.  相似文献   

16.
Synthetic aperture radar (SAR) has operational applications in crop mapping and monitoring in many countries due to the distinct backscatter signature at various stages of crop growth. Polarimetric analysis of SAR data from different satellites was used for information extraction from different types of scatters in imaged terrain. The scattering processes were analyzed through the received scatter matrix derived from the target decomposition of SAR data. Three decomposition techniques, namely Freeman–Durden, Cloude–Pottier and Touzi decomposition of the ALOS PALSAR-1 data, were used in this study to extract land use/cover information with a specific emphasis on agriculture. The decomposed output parameters from these techniques were classified with supervised classifier of support vector machine (SVM) using region of interest (ROI) selected land use/cover classes. An accuracy assessment for the classified output was carried out using the ROI. The Ramgarh village in Jaisalmer district of Rajasthan with the predominance of agricultural land, sand dunes and settlements was chosen as the study area. Freeman–Durden decomposition resulted in the highest overall accuracy of about 85% in the land use/cover classification among the three decomposition techniques adopted in the study. It was also observed that the accuracy of land use/cover mapping derived from Cloude–Pottier and Touzi decompositions improved with the use of eigenvalues in the SVM classification. Higher accuracies in the classification of agriculture land were noted with all the three decomposition techniques. The four parameters of Cloude–Pottier (H, A, α, β) and Touzi (α s, Φ s, ψ, τ) decompositions improved the classification accuracy for all the classes due to eigenvalues. The overall classification accuracy was above 88% for both the decomposition techniques with four parameters. The soil moisture values for agriculture land and sand dunes were validated through soil moisture maps generated using Oh 1992 and 2004 models.  相似文献   

17.
The stress regime in a Rotliegend reservoir of the Northeast German Basin   总被引:2,自引:0,他引:2  
In-situ stresses have significant impact, either positive or negative, on the short and long term behaviour of fractured reservoirs. The knowledge of the stress conditions are therefore important for planning and utilization of man-made geothermal reservoirs. The geothermal field Groß Schönebeck (40 km north of Berlin/Germany) belongs to the key sites in the northeastern German Basin. We present a stress state determination for this Lower Permian (Rotliegend) reservoir by an integrated approach of 3D structural modelling, 3D fault mapping, stress ratio definition based on frictional constraints, and slip-tendency analysis. The results indicate stress ratios of the minimum horizontal stress S hmin being equal or increasing 0.55 times the amount of the vertical stress S V (S hmin ≥ 0.55S V ) and of the maximum horizontal stress S Hmax ≤ 0.78–1.00S V in stress regimes from normal to strike slip faulting. Thus, acting stresses in the 4,100-m deep reservoir are S V  = 100 MPa, S hmin = 55 MPa and S Hmax = 78?100 MPa. Values from hydraulic fracturing support these results. Various fault sets of the reservoir are characterized in terms of their potential to conduct geothermal fluids based on their slip and dilatation tendency. This combined approach can be adopted to any other geothermal site investigation.  相似文献   

18.
The deformations of the sedimentary cover at near-latitudinal geotraverses west and east of the Mid-Atlantic Ridge in the equatorial part of ocean are compared with potential fields and variations of the V p/V s attribute at a depth of ~470 km. The features of sedimentary cover deformations in abyssal basins are formulated, as well as their differences from the undisturbed bedding of sediments. The elements of chain of phenomena with common spatial manifestations and cause-and-effect relationships have been established, including heterogeneous horizontal movements, which make up macrojointing above “cold” mantle blocks at a depth of ~470 km; serpentinization of upper-mantle rocks; the formation of superposed magnetic anomalies; the release of the fluids, which acoustically bleach out the sedimentary sequence in seismic imaging; and decompaction of rocks leading to vertical motions and forced folding. The origin of the Atlantic marginal dislocation zone is explained. The coincidence of the deformation boundary in the equatorial Atlantic with the zero contour line of the V p/V s attribute is revealed. This coincidence is an indicator of the rheological state of the upper mantle.  相似文献   

19.
Subcritical crack growth plays an important role in evaluating the long-term stability of structures in rocks. By applying the constant-displacement-relax method, two groups of test specimens that one immersed in groundwater and the other in air were tested to get the parameters of subcritical crack growth in double torsion test. The relations of the stress intensity factor K I versus the subcritical crack growth velocity V were obtained under the two different environments, and the behavior of subcritical crack growth was also analyzed. The results showed: the relations of lg K I  ? lg V accorded with linear rules, which is in good agreement with Charles theory; Compared with specimens in nature state, the lg K I  ? lg V curves of saturated water specimens locate top left comer of those of air specimens. The slope of curve is smaller, and the intercept is bigger, which shows that the water–rock interaction speeds up the subcritical crack growth. And A increases 2.9 × 1018 folds but n decreases from 85.12 to 40.83 because of the water–rock interaction. Meanwhile, the fracture toughness K IC also decreases obviously from 2.55 in air to 2.26 in water due to water rock interaction. The testing results provide a basis for time-dependence of rock engineering stability.  相似文献   

20.
It is shown that the approximation of the complex, tidally distorted shape of a star as a circular disc with local line profiles and a linear limb-darkening law, which is usually applied when deriving equatorial stellar rotation velocities from line profiles, leads to overestimation of the equatorial velocity V rot sin i and underestimation of the component mass ratio q = M x /M v . A formula enabling correction of the effect of these simplifying assumptions on the shape of a star is used to re-determine the mass ratios q and the masses of the black holes M x and visual components M v in low-mass X-ray binary systems containing black holes. Taking into account the tidal–rotational distortion of the stellar shape can significantly increase the mass ratios q = M x /M v , reducing M v , while M x changes only slightly. The resulting distribution of M v attains its maximum near M v ? 0.35M , in disagreement with the results of population synthesis computations realizing standard models for Galactic X-ray novae with black holes. Possible ways to overcome this inconsistency are discussed. The derived distribution of M x also differs strongly from the mass distribution for massive stars in the Galaxy.  相似文献   

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

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