首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 34 毫秒
1.
It is noticed that few geophysical studies have been carried out to decipher the crustal structure of southwestern part of the Northeast India comprising of Tripura fold belt and Bengal basin as compared to the Shillong plateau and the Brahmaputra basin. This region has a long history of seismicity that is still continuing. We have determined first-order crustal features in terms of Moho depths (H) and average VP/VS ratios (κ) using H-κ stacking technique. The inversion of receiver functions data yields near surface thick sedimentary layer in the Bengal basin, which is nearly absent in the Shillong plateau and Tripura fold belt. Our result suggests that the crust is thicker (38–45 km) in the Tripura fold belt region with higher shear-wave velocity in the lower crust than the Shillong plateau. The distribution of VP/VS ratio indicates heterogeneity throughout the whole region. While low to medium value of Poisson’s ratio (1.69–1.75) indicates the presence of felsic crust in the Shillong plateau of the extended Indian Archean crust. The medium to high values of VP/VS ratio (> 1.780) in the Bengal basin and the Tripura fold belt region represent mafic crust during the formation of the Bengal delta and the Tripura fold belt creation in the Precambrian to the Permian age. The depth of the sediments in the Bengal basin is up to 8 km on its eastern margin, which get shallower toward its northeastern and southeastern margins.  相似文献   

2.
We analyzed receiver function of teleseismic events recorded at twelve Indonesian-GEOFON (IA-GE) broadband stations using nonlinear Neighbourhood Algorithm (NA) inversion and H-k stacking methods to estimate crustal thickness, V p /V s ratios and S-wave velocity structure along Sunda-Banda arc transition zone. We observed crustal thickness of 34–37 km in Timor Island, which is consistent with the previous works. The thick crust (> 30 km) is also found beneath Sumba and Flores Islands, which might be related to the arc-continent collision causing the thickened crust. In Timor and Sumba Islands, we observed high V p /V s ratio (> 1.84) with low velocity zone that might be associated with the presence of mafic and ultramafic materials and fluid filled fracture zone. The high V p /V s ratio observed at Sumbawa and Flores volcanic Islands might be an indication of partial melt related to the upwelling of hot asthenosphere material through the subducted slab.  相似文献   

3.
The western Sichuan hydrothermal area is located at the northeastern margin of the eastern syntaxis of the Qinghai-Tibet Plateau, which is also the eastern end of the Mediterranean-Himalayan geothermal activity zone. There are 248 warm or hot springs in this area, and 11 have temperatures beyond the local boiling temperature. Most of these hot springs are distributed along the Jinshajiang, Dege-Xiangcheng, Ganzi-Litang, and Xianshuihe faults, forming a NW-SE hydrothermal belt. A geothermal analysis of this high-temperature hydrothermal area is an important basis for understanding the deep geodynamic process of the eastern syntaxis of the Qinghai-Tibet Plateau. In addition, this study offers an a priori view to utilize geothermal resources, which is important in both scientific research and application. We use gravity, magnetic, seismic, and helium isotope data to analyze the crust-mantle heat flow ratio and deep geothermal structure. The results show that the background terrestrial heat flow descends from southwest to northeast. The crustal heat ratio is not more than 60%. The high temperature hydrothermal active is related to crustal dynamics processes. Along the Batang-Litang-Kangding line, the Moho depth increases eastward, which is consistent with the changing Qc/Qm(crustal/mantle heat flow) ratio trend. The geoid in the hydrothermal zone is 4–6 km higher than the surroundings, forming a local "platform". The NW-SE striking local tensile stress zone and uplift structure in the upper and middle crust corresponds with the surface hydrothermal active zone. There is an average Curie Point Depth(CPD) of 19.5–22.5 km in Batang, Litang, and Kangding. The local shear-wave(S-wave) velocity is relatively low in the middle and lower crust. The S-wave shows a low velocity trap(Vs3.2 km s.1) at 15–30 km, which is considered a high-temperature partial melting magma, the crustal source of the hydrothermal active zone. We conclude that the hydrothermal system in this area can be divided into Batang-type and Kangding-type, both of which rely on a crustal heating cycle of atmospheric precipitation and surface water along the fracture zone. The heat is derived from the middle and lower crust: groundwater penetrates the deep faults bringing geothermal energy back to the surface and forming high-temperature springs.  相似文献   

4.
Crustal receiver functions have been calculated from 128 events for two three-component broadband seismomenters located on the south coast (FOMA) and in the central High Plateaux (ABPO) of Madagascar. For each station, crustal thickness and V p /V s ratio were estimated from H- κ plots. Self-consistent receiver functions from a smaller back-azimuthal range were then selected, stacked and inverted to determine shear wave velocity structure as a function of depth. These results were corroborated by guided forward modeling and by Monte Carlo error analysis. The crust is found to be thinner (39 ± 0.7 km) beneath the highland center of Madagascar compared to the coast (44 ± 1.6 km), which is the opposite of what would be expected for crustal isostasy, suggesting that present-day long wavelength topography is maintained, at least in part, dynamically. This inference of dynamic support is corroborated by shear wave splitting analyses at the same stations, which produce an overwhelming majority of null results (>96 %), as expected for vertical mantle flow or asthenospheric upwelling beneath the island. These findings suggest a sub-plate origin for dynamic support.  相似文献   

5.
The Q-factor estimates of the Earth’s crust and upper mantle as the functions of frequency (Q(f)) are obtained for the seismic S-waves at frequencies up to ~35 Hz. The estimates are based on the data for ~40 earthquakes recorded by the Kislovodsk seismic station since 2000. The magnitudes of these events are MW > 3.8, the sources are located in the depth interval from 1 to 165 km, and the epicentral distances range from ~100 to 300 km. The Q-factor estimates are obtained by the methods developed by Aki and Rautian et al., which employ the suppression of the effects of the source radiation spectrum and local site responses in the S-wave spectra by the coda waves measured at a fixed lapse time (time from the first arrival). The radiation pattern effects are cancelled by averaging over many events whose sources are distributed in a wide azimuthal sector centered at the receiving site. The geometrical spreading was specified in the form of a piecewise-continuous function of distance which behaves as 1/R at the distances from 1 to 50 km from the source, has a plateau at 1/50 in the interval from 50–70 km to 130–150 km, and decays as \({\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 {\sqrt R }}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{${\sqrt R }$}}\) beyond 130–150 km. For this geometrical spreading model and some of its modifications, the following Q-factor estimates are obtained: Q(f) ~ 85f0.9 at the frequencies ranging from ~1 to 20 Hz and Q(f) ~ 75f1.0 at the frequencies ranging from ~1 to 35 Hz.  相似文献   

6.
The interpretation of seismic refraction and wide angle reflection data usually involves the creation of a velocity model based on an inverse or forward modelling of the travel times of crustal and mantle phases using the ray theory approach. The modelling codes differ in terms of model parameterization, data used for modelling, regularization of the result, etc. It is helpful to know the capabilities, advantages and limitations of the code used compared to others.This work compares some popular 2D seismic modelling codes using the dataset collected along the seismic wide-angle profile DOBRE-4, where quite peculiar/uncommon reflected phases were observed in the wavefield.The ~505 km long profile was realized in southern Ukraine in 2009, using 13 shot points and 230 recording stations. Double PMP phases with a different reduced time (7.5–11 s) and a different apparent velocity, intersecting each other, are observed in the seismic wavefield. This is the most striking feature of the data. They are interpreted as reflections from strongly dipping Moho segments with an opposite dip. Two steps were used for the modelling. In the previous work by Starostenko et al. (2013), the trial-and-error forward model based on refracted and reflected phases (SEIS83 code) was published. The interesting feature is the high-amplitude (8–17 km) variability of the Moho depth in the form of downward and upward bends. This model is compared with results from other seismic inversion methods: the first arrivals tomography package FAST based on first arrivals; the JIVE3D code, which can also use later refracted arrivals and reflections; and the forward and inversion code RAYINVR using both refracted and reflected phases. Modelling with all the codes tested showed substantial variability of the Moho depth along the DOBRE-4 profile. However, SEIS83 and RAYINVR packages seem to give the most coincident results.  相似文献   

7.
Based on the Anapa (ANN) seismic station records of ~40 earthquakes (MW > 3.9) that occurred within ~300 km of the station since 2002 up to the present time, the source parameters and quality factor of the Earth’s crust (Q(f)) and upper mantle are estimated for the S-waves in the 1–8 Hz frequency band. The regional coda analysis techniques which allow separating the effects associated with seismic source (source effects) and with the propagation path of seismic waves (path effects) are employed. The Q-factor estimates are obtained in the form Q(f) = 90 × f 0.7 for the epicentral distances r < 120 km and in the form Q(f) = 90 × f1.0 for r > 120 km. The established Q(f) and source parameters are close to the estimates for Central Japan, which is probably due to the similar tectonic structure of the regions. The shapes of the source parameters are found to be independent of the magnitude of the earthquakes in the magnitude range 3.9–5.6; however, the radiation of the high-frequency components (f > 4–5 Hz) is enhanced with the depth of the source (down to h ~ 60 km). The estimates Q(f) of the quality factor determined from the records by the Sochi, Anapa, and Kislovodsk seismic stations allowed a more accurate determination of the seismic moments and magnitudes of the Caucasian earthquakes. The studies will be continued for obtaining the Q(f) estimates, geometrical spreading functions, and frequency-dependent amplification of seismic waves in the Earth’s crust in the other regions of the Northern Caucasus.  相似文献   

8.
The parameters of the earthquake that took place February 3, 2015, near the city of Sumy, Ukraine, were calculated from an analysis of records obtained by both Russian and Ukrainian seismic stations (Poltava, Skvira, Nikolaev, Dneropetrovsk, and Desna). The calculated hypocenter depth of 54 km was verified by several approaches: isolation of deep PP, SP phases from the records of remote stations and solution of the kinematic problem for the Poltava station. The focal mechanism as shear with a complex fault component was determined by the first arrivals of P-waves. The data on the azimuthal travel-time curve confirm the focal mechanism. We have calculated the earthquake parameters; they are as follows: length gap L1 = 8.08 km, L2 = 6.68 km, a destruction rate of C = 2 km/s. We have obtained the dynamic parameters of the event. The calculated fault length (L = 5.46 km) within the accuracy limits of the method coincides with the early result obtained by the azimuthal travel-time curve. On the basis of these results, we suggest that elastic energy release and formation of the dislocations in the earthquake source occurred on a smooth, prefractured fault (σr > 0). Association of the hypocenter with the tectonic node of the northern marginal fault of the Dnieper–Donets graben and northern branch of Kryvyi Rih–Kremenchuk suture confirm this. Here, we observe a considerable Moho depth, structural alteration, and high gradients of the temperature and magnetic and electric rock properties in the lower Earth’s crust and upper mantle. These circumstances are favorable for the earthquake occurring here.  相似文献   

9.
We made an attempt to assess the shear wave velocity values V S and, to a lesser extent, the V P values from ambient noise recordings in an array configuration. Five array sites were situated in the close proximity to borehole sites. Shear wave velocity profiles were modeled at these five array sites with the aid of two computational techniques, viz. spatial autocorrelation (SPAC) and H/V ellipticity. Out of these five array sites, velocity estimates could be reliably inferred at three locations. The shear wave velocities estimated by these methods are found to be quite consistent with each other. The computed V S values up to 30 m depth are in the range from 275 to 375 m/s in most of the sites, which implies prevalence of a low velocity zone at some pocket areas. The results were corroborated by evidence of site geology as well as geotechnical information.  相似文献   

10.
The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K sat and K sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V P/V S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K sat–μ, EEI, V P/V S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration.  相似文献   

11.
The electric field generation at the front of the current pulse, which originates in a coronal magnetic loop owing to the development of the Rayleigh–Taylor magnetic instability at loop footpoints, has been considered. During the τA ≈ l/V A ≈ 5?25 s time (where l is the plasma plume height entering a magnetic loop as a result of the Rayleigh–Taylor instability), a disturbance related to the magnetic field tension B ?(r,t), “escapes” the instability region with the Alfvén velocity in this case. As a result, an electric current pulse Iz(z ? V A t), at the front of which an induction magnetic field E z, which is directed along the magnetic tube axis and can therefore accelerate particles, starts propagating along a magnetic loop with a characteristic scale of Δξ ≈ l. In the case of sufficiently large currents, when B ? 2/8π > p, an electric current pulse propagates nonlinearly, and a relatively large longitudinal electric field originates E z ≈ 2I z 3 V A/c 4a2Bz 2l, which can be larger than the Dreicer field, depending on the electric current value.  相似文献   

12.
Parameters of the interplanetary magnetic field and solar wind plasma during periods of 163 isolated substorms have been studied. It is shown that the solar wind velocity V and plasma density N remain approximately constant for at least 3 h before substorm onset Т o and 1 h after Т o . On average, the velocity of the solar wind exhibits a stable trend toward anticorrelation with its density over the whole data array. However, the situation is different if the values of V and N are considered with respect to the intensity of substorms observed during that period. With the growth of substorm intensity, quantified as the maximum absolute value of AL index, an increase in both the solar wind plasma velocity and density, at which these substorms appear, is obsreved. It has been found that the magnitude of the solar wind dynamic pressure P is closely related to the magnetosphere energy load defined as averaged values of the Kan–Lee electric field EKL and Newell parameter dΦ/dt averaged for 1 h interval before Т o . The growth of the dynamic pressure is accompanied by an increase in the load energy necessary for substorm generation. This interrelation between P and values of EKL and dΦ/dt is absent in other, arbitrarily chosen periods. It is believed that the processes accompanying increasing dynamic pressure of the solar wind result in the formation of magnetosphere conditions that increasingly impede substorm generation. Thus, the larger is P, the more solar wind energy must enter the Earth’s magnetosphere during the period of the growth phase for substorm generation. This energy is later released during the period of the substorm expansion phase and creates even more intense magnetic bays.  相似文献   

13.
Stacking velocity V C2, vertical velocity ratio γ 0, effective velocity ratio γ eff, and anisotropic parameter χ eff are correlated in the PS-converted-wave (PS-wave) anisotropic prestack Kirchhoff time migration (PKTM) velocity model and are thus difficult to independently determine. We extended the simplified two-parameter (stacking velocity V C2 and anisotropic parameter k eff) moveout equation from stacking velocity analysis to PKTM velocity model updating and formed a new four-parameter (stacking velocity V C2, vertical velocity ratio γ 0, effective velocity ratio γ eff, and anisotropic parameter k eff) PS-wave anisotropic PKTM velocity model updating and process flow based on the simplified two-parameter moveout equation. In the proposed method, first, the PS-wave two-parameter stacking velocity is analyzed to obtain the anisotropic PKTM initial velocity and anisotropic parameters; then, the velocity and anisotropic parameters are corrected by analyzing the residual moveout on common imaging point gathers after prestack time migration. The vertical velocity ratio γ 0 of the prestack time migration velocity model is obtained with an appropriate method utilizing the P- and PS-wave stacked sections after level calibration. The initial effective velocity ratio γ eff is calculated using the Thomsen (1999) equation in combination with the P-wave velocity analysis; ultimately, the final velocity model of the effective velocity ratio γ eff is obtained by percentage scanning migration. This method simplifies the PS-wave parameter estimation in high-quality imaging, reduces the uncertainty of multiparameter estimations, and obtains good imaging results in practice.  相似文献   

14.
We investigated ground response for Baku (Azerbaijan) from two earthquakes of magnitude M6.3 occurred in Caspian Sea (characterized as a near event) and M7.5 in Shamakhi (characterized as a remote extreme event). S-wave velocity with the average shear wave velocity over the topmost 30 m of soil is obtained by experimental method from the V P values measured for the soils. The downtown part of Baku city is characterized by low VS30 values (< 250 m/s), related to sand, water-saturated sand, gravel-pebble, and limestone with clay. High surface PGA of 240 gal for the M7.5 event and of about 190 gal for the M6.3 event, and hence a high ground motion amplification, is observed in the shoreline area, through downtown, in the north-west, and in the east parts of Baku city with soft clays, loamy sands, gravel, sediments.  相似文献   

15.
The results of studying the process of underthrusting in the Pacific plate under Kamchatka are presented, and the phenomena associated with this process are described. The seismic component of the velocity of the plate underthrust is estimated on the basis of (1) data from the CMT (Centroid Moment Tensor) catalog and (2) the sequence of the strongest Kamchatka earthquakes. A flat site with a strike azimuth of 217° and a dip angle of 25° located at depths of 30–70 km is assumed to be the interface between the plunging plate and Kamchatka. From CMT focal mechanisms, the underthrust velocity is estimated at V = 0.7 cm/yr for the southern zone (located south of Cape Shipunski) and at V = 1.1 cm/yr for the central zone (from Cape Shipunski to Cape Kronotski). From pairs of the strongest earthquakes that have occurred successively since 1737, the underthrust velocity for the southern zone is estimated at V = 6.6–7.1 cm/yr (from two pairs) and for the central zone, at V = 6.6 cm/yr. The creep portion of the underthrust amounts to 5–15% of the total velocity (the velocity of motion of the Pacific plate is 8 cm/yr).  相似文献   

16.
Terrestrial heat flow, Q=K×ΔT/ΔZ cal/cm2 sec has been determined at 51 localities (39 on land and 12 in the sea) in and around the Japanese Islands. The average values of observed heat flow in land and sea are 1.53µ cal/cm2sec and 1.48µcal/cm2sec respectively. These value do not differ greatly from the world’s averages. The outstanding features of the heat flow distribution are as follows:a) High heat flow region (Q>2.0µcal/cm2sec) exists in the Inner Zone of the Honshu Arc. This region of high heat flow is more distinct in the northeastern Japan than in the southwestern Japan.b) The High heat flow region seems to extend, through the Fossa Magna area, down to the Izu-Mariana Arc.c) It is also probable that a similar high heat flow zone exists in the inner side of the Kurile Arc.d) These zones of high heat flow precisely coincide with the zones of the Cenozoic orogeny in the area concerned.e) Far off the coast of the northeastern Japan, the area at about 150° E may be a high heat flow region.f) Low heat flow (Q<1.0µcal/cm2sec) prevails in the Pacific coast side of the northeastern Japan and in the oceanic area directly east of it, including the area of the Japan Trench.g) The region bounded by the above mentioned high and low heat flow regions has heat flow which is more or less normal. Based on these measurements, a « steady state ” temperature distribution in the crust has been calculated for each of the above regions of high, low and intermediate heat flow, and it was found that there is a large temperature differences between the bottom of the crust of the high and low heat flow regions: the temperature at the Moho boundary in the high heat flow regions should be as high as some 800～1000°C (d=27 km), whereas that under the low heat flow region should be only about 200°C (d=23 km). The high general temperature at the Moho under the high heat flow region seems to favor a production of magma in the upper mantle. Calculated Moho temperatures disfavor the hypothesis that the Moho boundary is due to phase transition.  相似文献   

17.
We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h Am , which corresponds to the maximum increase, and layer peak height h m F, while h Am is always less than h m F. The difference between h Am and h m F linearly increases with increasing h m F. Whereas the difference is ~38 km for h m F = 280 km, it is ~54 km for h m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.  相似文献   

18.
The S wave velocity distribution in the Earth’s crust and the first two hundred kilometers of the upper mantle is inferred from data of a seismological linear network including 18 broadband stations installed in the framework of the international teleseismic experiment carried out in 2003 in the south of Siberia and in Mongolia. Models were constructed by using P-to-S received function inversion beneath each station. Vertical cross sections of S wave velocities from the surface to depths of 65 and 270 km covering the entire 1000-km profile are constructed by the linear spline interpolation of individual velocity models. The vertical sections are also represented as anomalies relative to the standard velocity model. The most intense low velocity anomalies (from ?3 to ?6%) in the crust and upper mantle are identified beneath the Sayan, Khamar-Daban, and Khangai highlands and the Djida fold zone and agree both with other geophysical data and with the distribution of Late Cenozoic volcanic fields. The results of this work suggest that the activation of Mongolian-Siberian highlands is largely connected with uplift of the asthenosphere to the base of the crust.  相似文献   

19.
Attenuation characteristics in the New Madrid Seismic Zone (NMSZ) are estimated from 157 local seismograph recordings out of 46 earthquakes of 2.6?≤?M?≤?4.1 with hypocentral distances up to 60 km and focal depths down to 25 km. Digital waveform seismograms were obtained from local earthquakes in the NMSZ recorded by the Center for Earthquake Research and Information (CERI) at the University of Memphis. Using the coda normalization method, we tried to determine Q values and geometrical spreading exponents at 13 center frequencies. The scatter of the data and trade-off between the geometrical spreading and the quality factor did not allow us to simultaneously derive both these parameters from inversion. Assuming 1/R 1.0 as the geometrical spreading function in the NMSZ, the Q P and Q S estimates increase with increasing frequency from 354 and 426 at 4 Hz to 729 and 1091 at 24 Hz, respectively. Fitting a power law equation to the Q estimates, we found the attenuation models for the P waves and S waves in the frequency range of 4 to 24 Hz as Q P?=?(115.80?±?1.36) f (0.495?±?0.129) and Q S?=?(161.34?±?1.73) f (0.613?±?0.067), respectively. We did not consider Q estimates from the coda normalization method for frequencies less than 4 Hz in the regression analysis since the decay of coda amplitude was not observed at most bandpass filtered seismograms for these frequencies. Q S/Q P?>?1, for 4?≤?f?≤?24 Hz as well as strong intrinsic attenuation, suggest that the crust beneath the NMSZ is partially fluid-saturated. Further, high scattering attenuation indicates the presence of a high level of small-scale heterogeneities inside the crust in this region.  相似文献   

20.
Solutions of P-SV equations of motion in a homogeneous transversely isotropic elastic layer contain a factor exp(±ν j z), where z is the vertical coordinate and j?=?1, 2. For computing Rayleigh wave dispersion in a multi-layered half space, ν j is computed at each layer. For a given phase velocity (c), ν j becomes complex depending on the transversely isotropic parameters. When ν j is complex, classical Rayleigh waves do not exist and generalised Rayleigh waves propagate along a path inclined to the interface. We use transversely isotropic parameters as α H , β V , ξ, ? and η and find their limits beyond which ν j becomes complex. It is seen that ν j depends on ? and η, but does not depend on ξ. The complex ν j occurs when ? is small and η is large. For a given c/β V , the region of complex ν j in a ? -η plane increases with the increase of α H /β V . Further, for a given α H /β V , the complex region of ν j increases significantly with the decrease of c/β V . This study is useful to compute dispersion parameters of Rayleigh waves in a layered medium.  相似文献   

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

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