首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
On the basis of S wave information from Tai’an-Xinzhou DSS profile and with reference to the results from P-wave interpretation, the 2-D structures, including S-wave velocity V s, ratio γ between V p and V s; and Poisson’ s ratio σ, are calculated; the structural configuration of the profile is presented and the relevant inferences are drawn from the above results. Upwarping mantle districts (V s≈4.30 km/s) and sloping mantle districts (V s≈4.50 km/s) of the profile with velocity difference about −4% at the top of upper mantle are divided according to the differences of V s, γ and σ in different media and structures, also with reference to the information of their neighbouring regions; the existence of Niujiaqiao-Dongwang high-angle ultra-crustal fault zone is reaffirmed; the properties of low and high velocity blocks (zones) including the crust-mantle transitionalzone and the boudary indicators of North China rift valley are discussed. A comprehensive study is conducted on the relation of the interpretation results with earthquakes. It is concluded that the mantle upwarps, thermal material upwells through the high-angle fault, the primary hypocenter was located at the crust-mantle juncture 30.0∼33.0 km deep, and additional stress excited the M S=6.8 and M S=7.2 earthquakes at specific locations around 9.0 km below Niujiaqiao-Dongwang, the earthquakes took place around the high-angle ultra-crustal fault and centered in the brittle media and rock strata with low γ and low σ values. This subject is part of the 85-907-02 key project during the “8th Five-Year Plan” from the State Science and Technology Commission.  相似文献   

2.
Receiver function study in northern Sumatra and the Malaysian peninsula   总被引：1，自引：0，他引：1  
In this receiver function study, we investigate the structure of the crust beneath six seismic broadband stations close to the Sunda Arc formed by subduction of the Indo-Australian under the Sunda plate. We apply three different methods to analyse receiver functions at single stations. A recently developed algorithm determines absolute shear-wave velocities from observed frequency-dependent apparent incidence angles of P waves. Using waveform inversion of receiver functions and a modified Zhu and Kanamori algorithm, properties of discontinuities such as depth, velocity contrast, and sharpness are determined. The combination of the methods leads to robust results. The approach is validated by synthetic tests. Stations located on Malaysia show high-shear-wave velocities (V S) near the surface in the range of 3.4–3.6 km s − 1 attributed to crystalline rocks and 3.6–4.0 km s − 1 in the lower crust. Upper and lower crust are clearly separated, the Moho is found at normal depths of 30–34 km where it forms a sharp discontinuity at station KUM or a gradient at stations IPM and KOM. For stations close to the subduction zone (BSI, GSI and PSI) complexity within the crust is high. Near the surface low V S of 2.6–2.9 km s − 1 indicate sediment layers. High V S of 4.2 km s − 1 are found at depth greater than 6 and 2 km at BSI and PSI, respectively. There, the Moho is located at 37 and 40 km depth. At station GSI, situated closest to the trench, the subducting slab is imaged as a north-east dipping structure separated from the sediment layer by a 10 km wide gradient in V S between 10 and 20 km depth. Within the subducting slab V S ≈ 4.7 km s − 1. At station BSI, the subducting slab is found at depth between 90 and 110 km dipping 20° ± 8° in approximately N 60° E. A velocity increase in similar depth is indicated at station PSI, however no evidence for a dipping layer is found.  相似文献   

3.
S-wave velocity structure beneath the Ailaoshan-Red River fault was obtained from receiver functions by using teleseismic body wave records of broadband digital seismic stations. The average crustal thickness, Vp/Vs ratio and Poisson’s ratio were also estimated. The results indicate that the interface of crust and mantle beneath the Ailaoshan-Red River fault is not a sharp velocity discontinuity but a characteristic transition zone. The velocity increases relatively fast at the depth of Moho and then increases slowly in the uppermost mantle. The average crustal thickness across the fault is 36―37 km on the southwest side and 40―42 km on the northeast side, indicating that the fault cuts the crust. The relatively high Poisson’s ratio (0.26―0.28) of the crust implies a high content of mafic materials in the lower crust. Moreover, the lower crust with low velocity could be an ideal position for decoupling between the crust and upper mantle.  相似文献   

4.
There is a high possibility of reoccurrence of the Tonankai and Nankai earthquakes along the Nankai Trough in Japan. It is very important to predict the long-period ground motions from the next Tonankai and Nankai earthquakes with moment magnitudes of 8.1 and 8.4, respectively, to mitigate their disastrous effects. In this study, long-period (>2.5 s) ground motions were predicted using an earthquake scenario proposed by the Headquarters for Earthquake Research Promotion in Japan. The calculations were performed using a fourth-order finite difference method with a variable spacing staggered-grid in the frequency range 0.05–0.4 Hz. The attenuation characteristics (Q) in the finite difference simulations were assumed to be proportional to frequency (f) and S-wave velocity (V s) represented by Q = f · V s / 2. Such optimum attenuation characteristic for the sedimentary layers in the Osaka basin was obtained empirically by comparing the observed motions during the actual M5.5 event with the modeling results. We used the velocity structure model of the Osaka basin consisting of three sedimentary layers on bedrock. The characteristics of the predicted long-period ground motions from the next Tonankai and Nankai earthquakes depend significantly on the complex thickness distribution of the sediments inside the basin. The duration of the predicted long-period ground motions in the city of Osaka is more than 4 min, and the largest peak ground velocities (PGVs) exceed 80 cm/s. The predominant period is 5 to 6 s. These preliminary results indicate the possibility of earthquake damage because of future subduction earthquakes in large-scale constructions such as tall buildings, long-span bridges, and oil storage tanks in the Osaka area.  相似文献   

5.
Flow velocity is one of the most important hydrodynamic variables for both channelized (rill and gullies) and interrill erosive phenomena. The dye tracer technique to measure surface flow velocity Vs is based on the measurement of the travel time of a tracer needed to cover a known distance. The measured Vs must be corrected to obtain the mean flow velocity V using a factor αv = V/Vs which is generally empirically deduced. The Vs measurement can be influenced by the method applied to time the travel of the dye-tracer and αv can vary in different flow conditions. Experiments were performed by a fixed bed small flume simulating a rill channel for two roughness conditions (sieved soil, gravel). The comparison between a chronometer-based (CB) and video-based (VB) technique to measure Vs was carried out. For each slope-discharge combination, 20 measurements of Vs, characterized by a sample mean Vm, were carried out. For both techniques, the frequency distributions of Vs/Vm resulted independent of slope and discharge. For a given technique, all measurements resulted normally distributed, with a mean equal to one, and featured by a low variability. Therefore, Vm was considered representative of surface flow velocity. Regardless of roughness, the Vm values obtained by the two techniques were very close and characterized by a good measurement precision. The developed analysis on αv highlighted that it is not correlated with Reynolds number for turbulent flow regime. Moreover, αv is correlated neither with the Froude number nor with channel slope. However, the analysis of the empirical frequency distributions of the correction factor demonstrated a slope effect. For each technique (CB, VB)-roughness (soil, gravel) combination, a constant correction factor was statistically representative even if resulted in less accurate V estimations compared to those yielded by the slope-specific correction factor.  相似文献   

6.
Strombolian-type volcanic activity is characterized by a series of gas bubbles bursting at the top of a magma column and leading to the ejection of lava clots and gas emission at the surface. The quantitative analysis of physical parameters (e.g., velocity, size, and mass fluxes) controlling the emission dynamics of these volcanic products is very important for the understanding of eruption source mechanisms but remains difficult to obtain in a systematic fashion. Ground-based Doppler radar is found to be a very effective tool for measuring ejecta velocities at a high acquisition rate and close to the emission source. We present here a series of measurements carried out at Mt. Etna’s Southeast crater, using an L-band volcanological Doppler radar, during the 4 July 2001 Strombolian eruptions. Doppler radar data are supplemented by the analysis of video snapshots recorded simultaneously. We provide here a set of physical parameters systematically retrieved from 247 Strombolian explosions spanning 15 min and occurring during the paroxysm of the eruption from 21:30 to 21:45 UT. The time-average values give a maximum particle velocity of Vmaxp = 94.7±24 \textm/s V_{{\max }}^p = {94}.{7}\pm {24} {\text{m/s}} , a bulk lava jet velocity of V\textPW - rad = 37.6±1.9 \textm/s {V_{{{\text{PW - rad}}}}} = {37}.{6}\pm {1}.{9} {\text{m/s}} , and an initial gas velocity at the source vent of V0g = 118.4±36 \textm/s V_0^g = {118}.{4}\pm {36} {\text{m/s}} . The time-averaged particle diameter is found to be about D\textPW - rad = 4.2±2.1 \textcm {D_{{{\text{PW - rad}}}}} = {4}.{2}\pm {2}.{1} {\text{cm}} . The volume and mass gas fluxes are estimated from time-averaged source gas velocities over the sequence duration at Qvg = 3 - 11 ×103\textm3\text/s Q_v^g = {3} - {11} \times {1}{0^{{3}}}{{\text{m}}^{{3}}}{\text{/s}} and Qmg = 0.5 - 2 ×103\textkg/s Q_m^g = 0.{5} - {2} \times {1}{0^{{3}}}{\text{kg/s}} , respectively.  相似文献   

7.
Array measurements of microtremors at 16 sites in the city of Thessaloniki were performed to estimate the Vs velocity of soil formations for site effect analysis. The spatial autocorrelation method was used to determine phase velocity dispersion curves in the frequency range from 0.8–1.5 to 6–7 Hz. A Rayleigh wave inversion technique (stochastic method) was subsequently applied to determine the Vs profiles at all the examined sites. The determination of Vs profiles reached a depth of 320 m. Comparisons with Vs values from cross-hole tests at the same sites proved the reliability of the SPAC method. The accuracy of the Vs profiles, the ability to reach large penetration depths in densely populated urban areas and its low cost compared to conventional geophysical prospecting, make Mictrotremor Exploration Method very attractive and useful for microzonation and site effects studies. An example of its application for the site characterization in Thessaloniki is presented herein.  相似文献   

8.
By processing the data of vertical ionospheric sounding in Almaty for 2000–2009, we obtained the distributions of the heights of the maximum (h m F) and bottom (h bot F) of the F2-layer, incremental changes in its semi-thickness (δh), the characteristic time of losses (τ), and the vertical displacement velocity of the node of the thermospheric wind (V) during the transitional time of the day during nighttime increases in the electron concentration at the layer maximum. The comparison of the measured V and modeled V m velocities showed a certain discrepancy. The influence of the altitude gradient of the meridional thermospheric wind velocity on the behaviors of h m F, h bot F, δh, and τ during nighttime increases in the electron concentration is studied.  相似文献   

9.
景亮  吴琪  方怡 《震灾防御技术》2020,15(4):696-707
为探讨测试方法、试验条件以及级配特征对砂-粉混合料剪切波速的影响，对具有不同细粒含量FC，相对密度Dr以及初始有效围压$\sigma_{3 \mathrm{c}}^{\prime}$的砂-粉混合料进行弯曲元和共振柱试验。结果表明:当Dr =35%或50%时，剪切波速Vs随FC的增大先减小后增大；当Dr =60%时，Vs随FC的增大而减小；弯曲元试验测得的Vs明显大于共振柱试验测得的Vs，随着FC的增大，弯曲元试验与共振试验得到的Vs差值逐渐减小，而当FC＞20%时，两种试验得到的Vs基本相同。在考虑Vs弥散性之后，不同FC的混合料弯曲元与共振柱试验得到的Vs结果具有较好的一致性。基于Hardin模型建立的砂-粉混合料Vs预测方法具有较好的预测效果。  相似文献   

10.
Fragmentation of magma during Plinian volcanic eruptions   总被引：2，自引：0，他引：2  
 The ratio of the volume of vesicles (gas) to that of glass (liquid) in pumice clasts (V G /V L ) reflects the degassing and dynamic history experienced by a magma during an explosive eruption. V G /V L in pumices from a large number of Plinian eruption deposits is shown here to vary by two orders of magnitude, even between pumices at a given level in a deposit. These variations in V G /V L do not correlate with crystallinity or initial water content of the magmas or their eruptive intensities, despite large ranges in these variables. Gas volume ratios of pumices do, however, vary systematically with magma viscosity estimated at the point of fragmentation, and we infer that pumices do not quench at the level of fragmentation but undergo some post-fragmentary evolution. On the timescale of Plinian eruptions, pumices with viscosities <109 Pa s can expand after fragmentation, as long as their bubbles retain gas, at a rate inversely proportional to their viscosity. Once the bubbles connect to form a permeable network and lose their gas, expansion halts and pumices with viscosities <105 Pa s can collapse under the action of surface tension. Textural evidence from bubble sizes and shapes in pumices indicates that both expansion and collapse have taken place. The magnitudes of expansion and collapse, therefore, depend critically on the timing of bubble connectivity relative to the final moment of quenching. We propose that bubbles in different pumices become connected at different times throughout the time span between fragmentation and quenching. After accounting for these effects, we derive new information on the fragmentation process from two characteristics of pumices. The most important is a relatively constant minimum value of V G /V L of ∼1.78 (64 vol.% vesicularity) in all samples with viscosities >105 Pa s. This value is independent of magma composition and thus reflects a property of the eruptive mechanism. The other characteristic is that highly expanded pumices (>85 vol.% vesicularities) are common, which argues against overpressure in bubbles as a mechanism for fragmenting magma. We suggest that magma fragments when it reaches a vesicularity of ∼64 vol.%, but only if sheared sufficiently strongly. The intensity of shear varies as a function of velocity in the conduit, which is related to overpressure in the chamber, so that changes in overpressure with time are important in controlling the common progression from explosive to effusive activity at volcanoes. Received: 19 April 1995 / Accepted: 3 April 1996  相似文献   

11.
We estimated three-dimensional P- (Vp) and S-wave velocity (Vs) and Vp/Vs structures in and around the Onikobe volcanic area, northeastern Japan, by local travel time tomography. We used travel time data from source and receiver pairs located within and outside the study area, which plays an important role in obtaining the optimum ray coverage and in elucidating the deeper structure more accurately. Detailed information on deeper structures is essential for imaging the complete volcanic system from the magmatic source zone through areas of shallow hydrothermal circulation. More than 50 000 travel time data for the P-waves and 35 000 for the S-waves were used to image the velocity structure. Our results show the following dominant features: (1) two conduits in the upper crust with low Vp and low Vs indicative of H2O-rich fluid pathways: one lying beneath Naruko volcano, the other beneath the focal area of the 1962 Northern Miyagi earthquake (M6.5); (2) an underlying broad region in the lower crust with low Vp, low Vs and high Vp/Vs, suggestive of a zone of partial melt, from which the fluids in (1) are derived; and (3) low Vp/Vs areas near the surface of the Sanzugawa and Onikobe calderas, suggesting a diffuse vapor-saturated cap.  相似文献   

12.
We studied the applicability of two types of existing three-dimensional (3-D) basin velocity structure models of the Osaka basin, western Japan for long-period ground motion simulations. We synthesized long-period (3–20 s) ground motions in the Osaka basin during a M6.5 earthquake that occurred near the hypothetical Tonankai earthquake source area, approximately 200 km from Osaka. The simulations were performed using a 3-D finite-difference method with nonuniform staggered grids using the two basin velocity structure models. To study the ground motion characteristics inside the basin, we evaluated the wave field inside the basin using the transfer functions derived from the synthetics at the basin and a reference rock site outside the basin. The synthetic waveforms at the basin site were obtained by a convolution of the calculated transfer function and the observed waveform at the reference rock site. First, we estimated the appropriate Q values for the sediment layers. Assuming that the Q value depends on the S wave velocity V S and period T, it was set to Q = (1/3V S)(T 0/T) where V S is in m/s and the reference period T 0 is 3.0 s. Second, we compared the synthetics and the observations using waveforms and pseudovelocity response spectra, together with a comparison of the velocity structures of the two basin models. We also introduced a goodness-of-fit factor to the pseudovelocity response spectra as an objective index. The synthetics of both the models reproduced the observations reasonably well at most of the stations in the central part the basin. At some stations, however, especially where the bedrock depth varies sharply, there were noticeable discrepancies in the simulation results of the models, and the synthetics did not accurately reproduce the observation. Our results indicate that the superiority of one model over the other cannot be determined and that an improvement in the basin velocity structure models based on simulation studies is required, especially along the basin edges. We also conclude that our transfer function method can be used to examine the applicability of the basin velocity structure models for long-period ground motion simulations.  相似文献   

13.
Scaling properties of variable electric fields in the topside ionosphere have been investigated on scales s from ∼30 m to 2 km by FAST electric field observations with sample rate of 512 s−1, in sixteen events of the broadband ELF turbulence. It is shown that down to scales of a few hundred meters, the power of turbulent electric fluctuations is a power law, ∼s α. Scaling index α derived from the slope of logarithmic diagrams (LD) constructed by the discrete wavelet transform of data can be estimated as α = 2.2 ± 0.3, which is close to α estimate earlier reported for scales 1–30 km by electric field observations of the Dynamics Explorer 2 satellite. The behavior of α index is analyzed near the scale of the order of electron inertial length λe = c/ω0 (ω0 being the electron plasma frequency). At altitudes considered (700–2500 km), λe makes 100–900 m. We demonstrate that at scales ≤λe, a decrease of LD slope and deviation from the power law are typically observed. As pointed out in the discussion, this feature cannot be identified as a transition to the diffusion range, where dissipation of the turbulence occurs.  相似文献   

14.
We studied the long-period ground motions in the Osaka sedimentary basin, Japan, which contains a 1- to 3-km thickness of sediments and is the site of many buildings or construction structures with long-natural period. We simulated the broadband ground motions likely to be produced by the hypothetical Nankai earthquake: the earthquake expected to give rise to the most severe long-period ground motion within the basin. For the simulation, we constructed multiscale heterogeneous source models based on the Central Disaster Management Council of Japan (CDMC) source model and adopted a hybrid computation method in which long-period motion and short-period motion are computed using a 3-D finite difference method and the stochastic Green’s function method, respectively. In computing long-period motions, we used a 3-D structure model of the crust and the Osaka sedimentary basin. The ground motions are estimated to have peak velocities of 50–90 cm/s, prolonged durations exceeding 300 s, and long predominant periods of 5–10 s in the area with great thickness of sediments. The predominant periods are in agreement with an approximate evaluation by 4 H/V s where H and V s are the thickness of the sediment and the average S wave velocity, respectively.  相似文献   

15.
The Effect of Velocity Inversions on H/V   总被引：2，自引：0，他引：2  
We analyzed the phenomenology of microtremor H/V curves under inversions in the shear-wave velocity (Vs) profile in the subsoil. Under no Vs inversion the spectral signature of the H/V peaks is found to be ‘eye-shaped’ with the horizontal components higher than the vertical. Conversely, under negative velocity gradients, numerous of differences emerge. I) A H/V ratio below 1 is observed for a wide range of frequencies, due to the decrease of the horizontal components below the vertical one. II) In the presence of persistent H/V < 1, small bumps in the H/V ratio given by local minima in the vertical spectral component may represent the relics of the peaks indicating resonances and stratigraphic discontinuities. As a consequence, in the presence of velocity inversions the H/V > 2 SESAME (2004) criterion fails but a stratigraphic interpretation may still be possible. III) The H/V curves should always be interpreted together with the single component spectra. IV) Microtremor H/V measurements for stratigraphic/microzonation purposes on stiff artificial soils, (asphalt, concrete, cement, pavements) should always be avoided since the latter often produce velocity inversions. This may have consequences in the intermediate to high frequency domain ( > 1 Hz) also in the application of reference site methods, like Hsite/Hbedrock, to microtremor. Theoretical modeling confirms these experimental findings.  相似文献   

16.
Summary In a series of triaxial experiments we have measuredV p ,V s and volumetric strain simultaneously in dilating dry and saturated rocks. For the first time these data permit quantitative comparison of seismic velocities or their ratio and dilatant volumetric strain. In air-dry samplesV p /V s decreases by a few per cent at strains of 10–3; in saturated materials with high pore pressure,V p /V s increases by a comparable amount. Decreases in seismic velocity ratio are difficult to generate in initially saturated rocks even with low pore pressures and at strain rates of 10–4/sec. A liquid-vapor transition will not produce a significant drop inV p /V s . If dilatancy and fluid flow are responsible for seismic travel time anomalies prior to earthquakes, our results suggest that such anomalies will occur only in regions where pore fluid source to sink dimensions are of the order of 10 km or more, or in regions where the rocks are not saturated to begin with.  相似文献   

17.
The evaluation and design of stone column improvement ground for liquefaction mitigation is a challenging issue for the state of practice. In this paper, a shear wave velocity-based approach is proposed based on the well-defined correlations of liquefaction resistance (CRR)-shear wave velocity (V s)-void ratio (e) of sandy soils, and the values of parameters in this approach are recommended for preliminary design purpose when site specific values are not available. The detailed procedures of pre- and post-improvement liquefaction evaluations and stone column design are given. According to this approach, the required level of ground improvement will be met once the target V s of soil is raised high enough (i.e., no less than the critical velocity) to resist the given earthquake loading according to the CRR-V s relationship, and then this requirement is transferred to the control of target void ratio (i.e., the critical e) according to the V s-e relationship. As this approach relies on the densification of the surrounding soil instead of the whole improved ground and is conservative by nature, specific considerations of the densification mechanism and effect are given, and the effects of drainage and reinforcement of stone columns are also discussed. A case study of a thermal power plant in Indonesia is introduced, where the effectiveness of stone column improved ground was evaluated by the proposed V s-based method and compared with the SPT-based evaluation. This improved ground performed well and experienced no liquefaction during subsequent strong earthquakes.  相似文献   

18.
In this paper we study the variation of Vp/Vs and Poisson's ratio (δ) in the Yellowstone National Park region, using earthquakes which were well recorded by a local seismic network. We find that the average Vp/Vs value within the geothermally active Yellowstone caldera is about 7% lower than in the area outside the caldera. Within the caldera itself there may be a further 2–7% reduction of Vp/Vs in the hydrothermally active Norris Geyser Basin, the Upper and Lower Geyser Basins, and the Yellowstone Lake and Mud Volcano regions. After considering various possible causes for Vp/Vs changes, such as geologic and structural differences, thermal effects, partial melting, and hydrothermal activity, we conclude that the most plausible explanation for the observed Vp/Vs reduction is the presence of hot-water at temperatures and pore-pressures near the water steam transition in the caldera geothermal reservoirs.  相似文献   

19.
A non-hydrostatic model in cross-sectional form with an idealized sill is used to examine the influence of sill depth (h s) and aspect ratio upon internal motion. The model is forced with a barotropic tide and internal waves and mixing occurs at the sill. Calculations using a wide sill and quantifying the response using power spectra show that for a given tidal forcing namely Froude number F r as the sill depth (h s) increases the lee wave response and vertical mixing decrease. This is because of a reduction in across sill velocity U s due to increased depth. Calculations show that the sill Froude number F s based on sill depth and across sill velocity is one parameter that controls the response at the sill. At low F s (namely F s ≪ 1) in the wide sill case, there is little lee wave production, and the response is in terms of internal tides. At high F s, calculations with a narrow sill show that for a given F s value, the lee wave response and internal mixing increase with increasing aspect ratio. Calculations using a narrow sill with constant U s show that for small values of h s, a near surface mixed layer can occur on the downstream side of the sill. For large values of h s, a thick well-mixed bottom boundary layer occurs due to turbulence produced by the lee waves at the seabed. For intermediate values of h s, “internal mixing” dominates the solution and controls across thermocline mixing.  相似文献   

20.
The characteristic of seismic coda wave attenuation in Yunnan area in 7 frequency-bands range from 1 Hz to 20 Hz was estimated by using the local earthquake's waveform data recorded from 22 Yunnan digital seismic stations.Coda attenuation Q-c1 of each station was firstly calculated by single scattering method. Then, mean free path Le and seismic albedo Bo of each station were calculated, and scattering attenuation Q-1s and intrinsic attenuation Q-1i were separated from total attenuation Q-1t by multiple lapse time window analysis based on the multiple scattering model in uniform random isotropic scattering medium. The attenuating characteristics in Yunnan show that most value of Le are in 10～30 km, with maximal within 2～6 Hz;Bo are about 0.5 at 1～2 Hz, but less than 0.5at other frequency-bands, which means Q-1i is comparable with Q-1s at 1～2 Hz, and after 1～2 Hz, Q-1i is greater than Q-1s and dominates the attenuation process. Q-1c is close to Q-1i at other frequency bands except 1～2 Hz.Results show that Q-1 especially Qs-1 varies spatially, Q-1 in eastern Yunnan zone is a bit higher than in northwestern Yunnan zone;northwestern Yunnan zone higher than southwestern Yunnan zone. Comparing with other results in global, Qs-1 in Yunnan is lower than the global average value among these results, Q-1i is higher than the global average value, and Q-1t lies the middle among these results.  相似文献   

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

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