首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
This paper focuses on using high-frequency GAP-SENSORs (GSs), accelerometers, and load cells in a laminar shear box (LSB) filled with loose Toyoura sand to understand the effects of impact loads and cyclic shaking at 1-G on soil properties. The shear wave velocity at small strain (Vs) was calculated directly from first arrival reference using displacement time-history of two GSs under impact loading. Moreover, from first peak using the reduced deformation amplitude technique, damping ratio was calculated. In addition, shaking table tests were performed under harmonic loading with amplitude of acceleration inside the model ground varying from 0.02 g to 1 g. The frequencies of excitation varied from 1 Hz to 10 Hz. GSs and inside accelerometers were used to directly measure the outside lateral deformation and shear stress at different elevations of LSB, respectively. Results show that the shear modulus (G) and the damping ratio (D) behavior of model sand are generally consistent with the behavior presented by similar tests using only accelerometers. In addition, damping ratio increases as frequency loading increases. Characteristic changes in two shear stress components in shaking loading conditions were also investigated using high precision inside load cells.  相似文献   

2.
On 6 April 2009 a Mw=6.1 earthquake produced severe destruction and damage over the historic center of L’Aquila City (central Italy), in which the accelerometer stations AQK and AQU recorded a large amount of near-fault ground motion data. This paper analyzes the recorded ground motions and compares the observed peak accelerations and the horizontal to vertical response spectral ratios with those revealed from numerical simulations. The finite element method is considered herein to perform dynamic modeling on the soil profile underlying the seismic station AQU. The subsurface model, which is based on the reviewed surveys that were carried out in previous studies, consists of 200–400 m of Quaternary sediments overlying a Meso-Cenozoic carbonate bedrock. The Martin-Finn-Seed's pore-water pressure model is used in the simulations. The horizontal to vertical response spectral ratio that is observed during the weak seismic events shows three predominant frequencies at about 14 Hz, 3 Hz and 0.6 Hz, which may be related to the computed seismic motion amplification occurring at the shallow colluvium, at the top and base of the fluvial-lacustrine sequence, respectively. During the 2009 L’Aquila main shock the predominant frequency of 14 Hz shifts to lower values probably due to a peculiar wave-field incidence angle. The predominant frequency of 3 Hz shifts to lower values when the earthquake magnitude increases, which may be associated to the progressive softening of soil due to the excess pore-water pressure generation that reaches a maximum value of about 350 kPa in the top of fluvial-lacustrine sequence. The computed vertical peak acceleration underestimates the experimental value and the horizontal to vertical peak acceleration ratio that is observed at station AQU decreases when the earthquake magnitude increases, which reveals amplification of the vertical component of ground motion probably due to near-source effects.  相似文献   

3.
High-frequency (≥ 2 Hz) Multi-channel Analysis of Love Waves (MALW) provides a practical way to determine velocity of horizontally polarized shear (SH) waves for a layered earth model up to 30 m below the ground surface in many geological settings. The information used in the MALW method is phase of Love waves. Information on amplitude of Love waves is not utilized in the MALW method. In this paper we present a method that uses information on amplitude of high-frequency Love waves to estimate quality factors (Qs) of near-surface materials. Unlike Rayleigh waves, attenuation coefficients (amplitude) of Love waves are independent of quality factors for P waves and are function of quality factors of Love waves. In theory, a fewer parameters make the inversion of attenuation coefficients of Love waves more stable and reduce the degree of nonuniqueness. We discussed sensitivity of an inversion system based on a linear relationship between attenuation coefficients and dissipation factors (1/Qs). The sensitivity analysis suggested that damping and constraints to an inversion system are necessary to obtain a smooth and meaningful quality factor model when no other information is available. We used synthetic and real-world data to demonstrate feasibility of inversion of attenuation coefficients of high-frequency Love-wave data acquired with the MALW method for quality factors with a linear, damped and constrained system.  相似文献   

4.
In western India during the Bhuj earthquake (Mw 7.6) on January 26, 2001, the Anjar City at ~30 km southwest of Bhuj experienced three types of damage scenario: severely damaged, less damaged and non-damaged. Similar damage patterns were also observed for the 1819 (Mw 7.8) and the 1956 (Mw 6.0) earthquakes. Microtremor array measurements were conducted in and around the Anjar city to examine the strength of soil structures and damage pattern. Significant differences are observed in frequencies and amplitudes in horizontal-to-vertical spectral ratio (HVSR) using microtremor measurements. The severely- damaged site shows two peak amplitudes: 2.8 at 1.2 Hz; and 4.0 at 8.0 Hz. The less-damaged site also shows two amplitudes: 2.5 and 2.1 at 1.4 Hz; and 2.0 Hz, respectively. The non-damaged site, on the other hand, shows that the HVSR curves become almost flatter. Similar results for three types of damage scenario based on analyses of earthquake records are also observed for the study area. The microtremor array measurements has revealed shear wave velocity Vs≥400 m/s at 18 m depth in the non-damaged, at 40 m in the less-damaged and at 60 m depth in the severely-damaged sites. The site amplitudes and the Vs values show a good correlation with the soil characteristics and damage pattern, suggesting that strength of soil layers at varying depths is a dictating factor for the estimate of the earthquake risk evaluation of the area under study.  相似文献   

5.
Serpentinization of the mantle wedge is an important process that influences the seismic and mechanical properties in subduction zones. Seismic detection of serpentines relies on the knowledge of elastic properties of serpentinites, which thus far has not been possible in the absence of single-crystal elastic properties of antigorite. The elastic constants of antigorite, the dominant serpentine at high-pressure in subduction zones, were measured using Brillouin spectroscopy under ambient conditions. In addition, antigorite lattice preferred orientations (LPO) were determined using an electron back-scattering diffraction (EBSD) technique. Isotropic aggregate velocities are significantly lower than those of peridotites to allow seismic detection of serpentinites from tomography. The isotropic VP/VS ratio is 1.76 in the Voigt–Reuss–Hill average, not very different from that of 1.73 in peridotite, but may vary between 1.70 and 1.86 between the Voigt and Reuss bonds. Antigorite and deformed serpentinites have a very high seismic anisotropy and remarkably low velocities along particular directions. VP varies between 8.9 km s? 1 and 5.6 km s? 1 (46% anisotropy), and 8.3 km s? 1 and 5.8 km s? 1 (37%), and VS between 5.1 km s? 1 and 2.5 km s? 1 (66%), and 4.7 km s? 1 and 2.9 km s? 1 (50%) for the single-crystal and aggregate, respectively. The VP/VS ratio and shear wave splitting also vary with orientation between 1.2 and 3.4, and 1.3 and 2.8 for the single-crystal and aggregate, respectively. Thus deformed serpentinites can present seismic velocities similar to peridotites for wave propagation parallel to the foliation or lower than crustal rocks for wave propagation perpendicular to the foliation. These properties can be used to detect serpentinite, quantify the amount of serpentinization, and to discuss relationships between seismic anisotropy and deformation in the mantle wedge. Regions of high VP/VS ratios and extremely low velocities in the mantle wedge of subduction zones (down to about 6 and 3 km.s?1 for VP and VS, respectively) are difficult to explain without strong preferred orientation of serpentine. Local variations of anisotropy may result from kilometer-scale folding of serpentinites. Shear wave splittings up to 1–1.5 s can be explained with moderately thick (10–20 km) serpentinite bodies.  相似文献   

6.
The single backscattering model was used to estimate total attenuation of coda waves (Qc) of local earthquakes recorded on eight seismological stations in the complex area of the western continental Croatia. We estimated Q0 and n, parameters of the frequency dependent coda-Q using the relation Qc = Q0fn. Lapse time dependence of these parameters was studied using a constant 30 s long time window that was slid along the coda of seismograms. Obtained Qc were distributed into classes according to their lapse time, tL. For tL = 20–50 s we estimated Q0 = 45–184 and n = 0.49–0.94, and for tL = 60–100 s we obtained Q0 = 119–316 and n = 0.37–0.82. There is a tendency of decrease of parameter n with increasing Q0, and vice versa. The rates of change of both Q0 and n seem to decrease for lapse times larger than 50–80 s, indicating an alteration in rock properties controlling coda attenuation at depths of about 100–160 km. A very good correlation was found between the frequency dependence parameter n and the Moho depths for lapse times of 50, 60 and 70 s.  相似文献   

7.
By using six 4.5 Hz geophones, surface wave tests were performed on four different sites by dropping freely a 65 kg mass from a height of 5 m. The receivers were kept far away from the source to eliminate the arrival of body waves. Three different sources to nearest receiver distances (S), namely, 46 m, 56 m and 66 m, were chosen. Dispersion curves were drawn for all the sites. The maximum wavelength (λmax), the maximum depth (dmax) up to which exploration can be made and the frequency content of the signals depends on the site stiffness and the value of S. A stiffer site yields greater values of λmax and dmax. For stiffer sites, an increase in S leads to an increase in λmax. The predominant time durations of the signals increase from stiffer to softer sites. An inverse analysis was also performed based on the stiffness matrix approach in conjunction with the maximum vertical flexibility coefficient of ground surface to establish the governing mode of excitation. For the Site 2, the results from the surface wave tests were found to compare reasonably well with that determined on the basis of cross boreholes seismic tests.  相似文献   

8.
  相似文献   

9.
In a typical seismic dam safety evaluation, standard penetration, cone penetration, Becker penetration, or shear wave velocity (Vs) tests are often first conducted near the toe of an earth dam to infer if any liquefiable soil exists in the foundation of the dam footprint. In current practice, a level-ground condition is commonly assumed when normalizing penetration resistance and Vs, and may be assumed (particularly in preliminary assessments) in applying the cyclic stress method (with or without the Kα correction) to evaluate liquefaction. However, the presence of an earth dam, or any other large embankment or structure, significantly alters the normal and shear stresses in the foundation. This paper identifies and quantifies potential errors in ignoring altered stresses near heavy structures, and presents a methodology to incorporate these effects within the framework of the simplified procedure. Specifically, the effects of these altered stresses (in comparison to the level-ground assumption with and without Kα correction) on the: (1) normalization of field measurements such as penetration resistance and Vs; (2) cyclic stress ratio (CSR); (3) cyclic resistance ratio (CRR); and (4) factor of safety against liquefaction triggering (FSliq), are evaluated by considering static and dynamic analyses of a generic earthen embankment (60 m high) resting on a saturated, cohesionless foundation (30 m deep). Our analyses indicated that ignoring the presence of induced static shear stresses can result in potentially unconservative errors in overburden correction factors of 30% to 60% at shallow depth (although this error is greatly muted at depths exceeding about 15 m), while errors in CSR potentially can range from about 20% too conservative to 40% unconservative. Potential errors in CRR can approach 50% unconservative at shallow depths, but again, this error is muted at depths exceeding about 15 m. Combining these factors, potentially unconservative errors in computing FSliq could exceed 100% at shallow depths (less than 15 m to 20 m) while at greater depth (exceeding 20 m) errors approach 20% on the conservative side.  相似文献   

10.
The spectral attenuation of solar irradiation was measured during summer in two types of coastal waters in southern Chile, a north Patagonian fjord (Seno Reloncaví) and open coast (Valdivia). In order to relate the light availability with the light requirements of upper subtidal seaweeds, the saturating irradiance for photosynthesis (Ek) from P–I curves was measured. In addition the UV risk was assessed. Based on the z1% of PAR, the lower limit of the euphotic zone in the studied systems averaged 21 m (Kd 0.24 m?1) in Seno Reloncaví and 18 m (Kd 0.27 m?1) in the coast of Valdivia. Photosynthesis of the studied seaweeds was saturated at markedly lower irradiances than found in their natural depths at the time of the study. Solar radiation penetrating into these depths at both locations largely supports the light requirements for the photosynthesis of subtidal species: 50–160 μmol m?2 s?1 for seaweeds from Seno Reloncaví (7 m tidal range) and 20–115 μmol m?2 s?1 for Valdivia assemblages (2 m tidal range). Optimal light conditions to saturate photosynthesis (Ek) were present at 10–16 m water depth. The attenuation of solar irradiation did not vary significantly between the fjord and coastal sites of this study. However, the underwater light climates to which seaweeds are exposed in these sites vary significantly because of the stronger influence of tidal range affecting the fjord system as compared with the open coastal site. The patterns of UV-B penetration in these water bodies suggest that seaweeds living in upper littoral zones such as the intertidal and shallow subtidal (<3 m) may be at risk.  相似文献   

11.
Inversion of local earthquake travel times and joint inversion of receiver functions and Rayleigh wave group velocity measurements were used to derive a simple model for the velocity crustal structure beneath the southern edge of the Central Alborz (Iran), including the seismically active area around the megacity of Tehran. The P and S travel times from 115 well-located earthquakes recorded by a dense local seismic network, operated from June to November 2006, were inverted to determine a 1D velocity model of the upper crust. The limited range of earthquake depths (between 2 km and 26 km) prevents us determining any velocity interfaces deeper than 25 km. The velocity of the lower crust and the depth of the Moho were found by joint inversion of receiver functions and Rayleigh wave group velocity data. The resulting P-wave velocity model comprises an upper crust with 3 km and 4 km thick sedimentary layers with P wave velocities (Vp) of ～5.4 and ～5.8 km s?1, respectively, above 9 km and 8 km thick layers of upper crystalline crust (Vp ～6.1 and ～6.25 km s?1 respectively). The lower crystalline crust is ～34 km thick (Vp ～ 6.40 km s?1). The total crustal thickness beneath this part of the Central Alborz is 58 ± 2 km.  相似文献   

12.
《Journal of Geodynamics》2010,49(3-5):331-339
The Free Core Nutation (FCN) is investigated with the help of its resonance effect on the tidal amplitudes in Superconducting Gravimeter (SG) records of the GGP network. The FCN resonance parameters are combined in a resonance equation involving the Earth's interior parameters. The sensitivity of the FCN parameters to the diurnal tidal waves demonstrates that the quality factor of the FCN is strongly dependent on the accuracy of the imaginary part estimates of the gravimetric factors close to the resonance. The weak amplitude of Ψ1 tidal wave on the Earth, which is the closest in frequency to the FCN, in addition to errors in ocean loading correction, explains the poor determination of the quality factor Q from surface gravimetric data. The inversion of tidal gravimetric factors leads to estimates of the period, Q and resonance strength of the FCN. We show that, by inverting log(Q) instead of Q, the results using the least-squares method optimized using the Levenberg–Marquardt algorithm are in agreement with the Bayesian probabilistic results and agree with the results obtained from VLBI nutation data. Finally, a combined inversion of 7 GGP European SG data is performed giving T = 428 ± 3 days and 7762 < Q < 31,989 (90% C.I.). An experimental estimate of the internal pressure Love number is also proposed.  相似文献   

13.
Ground motion intensity measures such as the peak ground acceleration (PGA) and the pseudo-spectral acceleration (PSA) at two sites due to the same seismic event are correlated. The spatial correlation needs to be considered when modeling ground-motion fields for seismic loss assessments, since it can have a significant influence on the statistical moments and probability distribution of aggregated seismic loss of a building portfolio.Empirical models of spatial correlation of ground motion intensity measures exist only for a few seismic regions in the world such as Japan, Taiwan and California, since for this purpose a dense observation network of earthquake ground motion is required. The Istanbul Earthquake Rapid Response and Early Warning System (IERREWS) provides one such dense array with station spacing of typically 2 km in the urban area of Istanbul. Based on the records of eight small to moderate (Mw3.5–Mw5.1) events, which occurred since 2003 in the Marmara region, we establish a model of intra-event spatial correlation for PGA and PSA up to the natural period of 1.0 s.The results indicate that the correlation coefficients of PGA and short-period PSA decay rapidly with increasing interstation distance, resulting in correlation lengths of approximately 3–4 km, while correlation lengths at longer natural periods (above 0.5 s) exceed 6 km. Finally, we implement the correlation model in a Monte Carlo simulation to evaluate economic loss in Istanbul's district Zeytinburnu due to a Mw7.2 scenario earthquake.  相似文献   

14.
The late Triassic to early Tertiary Coast Mountains Batholith (CMB) of British Columbia provides an ideal locale to study the processes whereby accreted terranes and subduction-related melts interact to form stable continental crust of intermediate to felsic composition and complementary ultramafic residuals. Seismic measurements, combined with calculated elastic properties of various CMB rock compositions, provide a window into the deep-crustal lithologies that are key to understanding the processes of continental growth and evolution. We use a combination of seismic observations and petrologic modeling to construct hypothetical crustal sections at representative locations across the CMB, then test the viability of these sections via forward modeling with synthetic seismic data. The compositions that make up our petrologic forward models are based on calculations using the free energy minimization program Perple_X to predict mineral assemblages at depth for the bulk compositions of exposed plutonic rocks collected in the study area. Seismic data were collected along two transects in west-central British Columbia: a southern line that crossed the CMB near the town of Bella Coola (near 52° N), and a northern line centered on the towns of Terrace and Kitimat (near 54° N). Along both transects, seismic receiver functions reveal high Vp/Vs ratios near the Insular/Intermontane terrane boundary and crustal thickness increasing from 26 ± 3 km to 34 ± 3 km (at the 1 sigma certainty level) from west to east across the Coast Shear Zone (CSZ). On the southern line, we observe an anomalous region of complex receiver functions and diminished Moho signals beneath the central portion of the CMB. Our petrologic and seismic profiles show that observed seismic data from much of the CMB can be well-matched in terms of crustal thickness and structure, average Vp/Vs, and amplitude of the Moho converted phase, without including ultramafic residual material in the lower crust.  相似文献   

15.
Modeling of multimode surface wave group velocity dispersion data sampling the eastern and the western Ganga basins, reveals a three layer crust with an average Vs of 3.7 km s?1, draped by ～2.5 km foreland sediments. The Moho is at a depth of 43 ± 2 km and 41 ± 2 km beneath the eastern and the western Ganga basins respectively. Crustal Vp/Vs shows a felsic upper and middle crust beneath the eastern Ganga basin (1.70) compared to a more mafic western Ganga basin crust (1.77). Due to higher radiogenic heat production in felsic than mafic rocks, a lateral thermal heterogeneity will be present in the foreland basin crust. This heterogeneity had been previously observed in the north Indian Shield immediately south of the foreland basin and must also continue northward below the Himalaya. The high heat producing felsic crust, underthrust below the Himalayas could be an important cause for melting of midcrustal rocks and emplacement of leucogranites. This is a plausible explanation for abundance of leucogranites in the east-central Himalaya compared to the west. The uppermost mantle Vs is also significantly lower beneath the eastern Ganga basin (4.30 km s?1) compared to the west (4.44 km s?1).  相似文献   

16.
Many authors have proposed that the study of seismicity rates is an appropriate technique for evaluating how close a seismic gap may be to rupture. We designed an algorithm for identification of patterns of significant seismic quiescence by using the definition of seismic quiescence proposed by Schreider (1990). This algorithm shows the area of quiescence where an earthquake of great magnitude may probably occur. We have applied our algorithm to the earthquake catalog on the Mexican Pacific coast located between 14 and 21 degrees of North latitude and 94 and 106 degrees West longitude; with depths less than or equal to 60 km and magnitude greater than or equal to 4.3, which occurred from January, 1965 until December, 2014. We have found significant patterns of seismic quietude before the earthquakes of Oaxaca (November 1978, Mw = 7.8), Petatlán (March 1979, Mw = 7.6), Michoacán (September 1985, Mw = 8.0, and Mw = 7.6) and Colima (October 1995, Mw = 8.0). Fortunately, in this century earthquakes of great magnitude have not occurred in Mexico. However, we have identified well-defined seismic quiescences in the Guerrero seismic-gap, which are apparently correlated with the occurrence of silent earthquakes in 2002, 2006 and 2010 recently discovered by GPS technology.  相似文献   

17.
Outcrops of the Cretaceous Upper sandstone formation some 375 km to the East of Addis Ababa on the motor Highway to Harar was paleomagnetically investigated. About seventy core samples were collected at various stratigraphic levels from 250–300 meters thick sedimentary formation. After standard sample preparations in the laboratory the resulting specimens were subjected to routine paleomagnetic demagnetization protocol. In the first steps of demagnetizations process the recent and viscous magnetizations were removed by heating until a temperature of level of 300 °C. Further demagnetization of the samples resulted in the isolation of the final magnetization with stable line segments that is directed towards the origin, which is interpreted as Characteristic Remanent Magnetization (ChRM). Rock – magnetic experiments have identified goethite (αFeOOH), hematite (αFe2O3), detritial hematite, and magnetite as the magnetic mineral phases carrying the remanence. The ChRM identified resulted in an average value of (Ds = 0.5°, Is = ?0.7°, α95 = 4.3°, N = 34) for the red sandstones while an average value of (Ds = 335.8°, Is = ?31.8°, α95 = 4.7°, N = 14) for the limestone intercalations. The former ChRM in the red sandstone is determined to be secondary while the latter ChRM is known to be primary. Comparison of these directional results and their pole equivalents with the African plate Apparent Polar Wander Path curve established by Besse and Courtillot (2003) give ages of between 115–130 Million years for the limestone intercalation and ages of 30 million years for red sandstone unit. These are interpreted respectively as estimates of the age of deposition and a later remagnetization respectively.  相似文献   

18.
The paper focuses on the study of the dynamic response of a wind turbine installed close to seismic fault in an earthquake prone region of Greece. The investigation of the seismic behavior of the turbine is performed with response spectrum analysis using the elastic acceleration spectrum provided by the Greek National Aseismic Code for the project area, increased by 25% due to proximity to seismic fault, as recommended by the Code. Also, dynamic analysis is performed using the spectrum obtained from the assessment of seismic hazard at the site of the project (Local), for two typical cases of earthquakes of magnitudes Ms=5.8 and Ms=7.3 with epicentral distances of 1 km and 11 km from the test position, respectively. Finally, a time-history analysis is carried out using an artificial accelerogram obtained from the local spectrum envelope, for the site of the project, considering 1% damping ratio. The direct comparison of the response spectrum analysis results with the corresponding ones obtained from the dynamic and time-history analyses indicates that the implementation of the elastic acceleration spectrum is insufficient, although increased, and therefore it could be deduced that the current recommendation of the Greek National Aseismic Code to consider an increment of 25% on the elastic acceleration spectrum in regions close to seismic faults, requires review and further improvement.  相似文献   

19.
The present work deals with 1D and 2D ground response analysis and liquefaction analysis of alluvial soil deposits from Kanpur region along Indo-Gangetic plains. Standard penetration tests and seismic down hole tests have been conducted at four locations namely IITK, Nankari village, Mandhana and Bithoor at 1.5 m interval up to a depth of 30 m below the ground surface to find the variation of penetration blows and the shear wave velocity along the depth. From the selected sites undisturbed as well as representative soil samples have been collected for detailed soil classification. The soil profiles from four sites have been considered for 1D and 2D ground response analysis by applying the free field motions of three Himalayan earthquakes namely Chamba earthquake (Mw—5.1), Chamoli earthquake (Mw—6.4) and Uttarkashi earthquake (Mw—6.5). An average value of Peak Ground Acceleration (PGA) obtained from 1D and 2D analysis is considered for liquefaction analysis and post-liquefaction settlement. The excess pore water pressure ratio is greater than 0.8 at a depth of 24 m from ground surface for IITK, Nankari village, Bithoor sites. More than 50% of post liquefaction settlement is contributed by layers from 21–30 m for all sites. In general, the soil deposits in Kanpur region have silty sand and sand deposits and are prone to liquefaction hazards due to drastic decrease of cyclic resistance ratio (CRR) at four chosen sites in Kanpur.  相似文献   

20.
The attenuation properties of the lithosphere in the Bam region, East-Central Iran, have been investigated. For this purpose, 42 local earthquakes having focal depths less than 25 km have been used. The quality factor of coda waves (Qc) has been estimated using the single back-scattering model. The quality factors Qp, Qd (P and direct S-waves) have been estimated using the extended-coda normalization method. Qi and Qs (the intrinsic and scattering attenuation parameters) have been estimated for the region. The values of Qp, Qd, Qc, Qi and Qs show a dependence on frequency in the range of 1.5–24 Hz for the Bam region. The average frequency-dependent relationships estimated for the region are Qp=(36±6)f(1.03±0.06), Qd=(59±8)f(1.00±0.03), Qc=(79±5)f(1.01±0.04), Qs=(131±4)f(1.01±0.04) and Qi=(104±6)f(1.01±0.05). A comparison between Qi and Qs shows that intrinsic absorption is predominant over scattering.The variation of Q has also been estimated at different lapse times to observe heterogeneities variation with depth. The variation of Q with frequency and lapse time shows that the lithosphere becomes more homogeneous with depth.The estimated Qo values at different stations suggest a low value of Q indicating a heterogeneous and attenuative crust beneath the entire region.  相似文献   

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

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