首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 38 毫秒
1.
We revised an equation for estimating palaeostress magnitude using the microboudin technique by incorporating the influence of time on the fracture strength of minerals. The equation was used to estimate triaxial palaeostresses from a rare sample of metachert from Turkey that contains microboudinaged, columnar tourmaline grains in a wide range of orientations within the foliation plane. The estimated principal palaeostresses are σ1 = 605 MPa, σ2 = 598 MPa, and σ3 = 597 MPa. As the microboudinage is considered to have occurred immediately before the rock encountered the brittle-plastic transition during exhumation, these stress values correspond to conditions at approximately 18 km depth and 300 °C within a Cretaceous orogenic belt.  相似文献   

2.
Talc is one of the weakest minerals that is associated with fault zones. Triaxial friction experiments conducted on water-saturated talc gouge at room temperature yield values of the coefficient of friction, μ (shear stress, τ/effective normal stress, σ′N) in the range 0.16–0.23, and μ increases with increasing σ′N. Talc gouge heated to temperatures of 100°–400 °C is consistently weaker than at room temperature, and μ < 0.1 at slow strain rates in some heated experiments. Talc also is characterized by inherently stable, velocity-strengthening behavior (strength increases with increasing shear rate) at all conditions tested. The low strength of talc is a consequence of its layered crystal structure and, in particular, its very weak interlayer bond. Its hydrophobic character may be responsible for the relatively small increase in μ with increasing σ′N at room temperature compared to other sheet silicates.Talc has a temperature–pressure range of stability that extends from surficial to eclogite-facies conditions, making it of potential significance in a variety of faulting environments. Talc has been identified in exhumed subduction zone thrusts, in fault gouge collected from oceanic transform and detachment faults associated with rift systems, and recently in serpentinite from the central creeping section of the San Andreas fault. Typically, talc crystallized in the active fault zones as a result of the reaction of ultramafic rocks with silica-saturated hydrothermal fluids. This mode of formation of talc is a prime example of a fault-zone weakening process. Because of its velocity-strengthening behavior, talc may play a role in stabilizing slip at depth in subduction zones and in the creeping faults of central and northern California that are associated with ophiolitic rocks.  相似文献   

3.
The co-seismic deformations produced during the September 27, 2003 Chuya earthquake (Ms = 7.5) that affected the Gorny Altai, Russia, are described and discussed along a 30 km long segment. The co-seismic deformations have manifested themselves both in unconsolidated sediments as R- and R′-shears, extension fractures and contraction structures, and in bedrock as the reactivation of preexisting schistosity zones and individual fractures, as well as development of new ruptures and coarse crushing zones. It has been established that the pattern of earthquake ruptures represents a typical fault zone trending NW–SE with a width reaching 4–5 km and a dextral strike–slip kinematics. The initial stress field that produced the whole structural pattern of co-seismic deformations during the Chuya earthquake, is associated with a transcurrent regime with a NNW–SSE, almost N–S, trending of compressional stress axis (σ1), and a ENE–WSW, almost E–W, trending of tensional stress axis (σ3). The state of stress in the newly-formed fault zone is relatively uniform. The local stress variations are expressed in insignificant deviation of σ1 from N–S to NW–SE or NE–SW, in short-term fluctuations of relative stress values in keeping their spatial orientations, or in a local increase of the plunge angle of the σ1. The geometry of the fault zone associated with the Chuya earthquake has been compared with the mechanical model of fracturing in large continental fault zones with dextral strike–slip kinematics. It is apparent that the observed fracture pattern corresponds to the late disjunctive stage of faulting when the master fault is not fully developed but its segments are already clearly defined. It has been shown that fracturing in widely different rocks follows the common laws of the deformation of solid bodies, even close to the Earth surface, and with high rates of movements.  相似文献   

4.
A note on fault reactivation   总被引：2，自引：0，他引：2  
Reactivation of existing faults whose normal lies in the σ1′σ3′ plane of a stress field with effective principal compressive stresses σ1′ >σ2′ >σ3′ is considered for the simplest frictional failure criterion, τ = μσn′ = μ(σn − P), where τ and σn are respectively the shear and normal stresses to the existing fault, P is the fluid pressure and μ is the static friction. For a plane oriented at θ to σ1′, the stress ratio for reactivation is (σ1′/σ3′) = (1 + μ cot θ)/(1 − μ tan θ). This ratio has a minimum positive value at the optimum angle for reactivation given by (1/μ) but reaches infinity when θ = 2θ*, beyond which σ3′ < 0 is a necessary condition for reactivation. An important consequence is that for typical rock friction coefficients, it is unlikely that normal faults will be reactivated as high-angle reverse faults or thrusts as low-angle normal faults, unless the effective least principal stress is tensile.  相似文献   

5.
A statistical analysis for the Rf/φ method is developed in the paper. The technique provides estimates of Rs, Ri and θ according to Rf and φ. Moreover, we can restore the pattern of initial orientations of marker objects by means of their frequency histogram, provided the strain is homogeneous. It is unnecessary to postulate random initial orientations if the principal directions of strain ellipsoid are known.  相似文献   

6.
This review of geological, seismological, geochronological and paleobotanical data is made to compare historic and geologic rates and styles of deformation of the Sierra Nevada and western Basin and Range Provinces. The main uplift of this region began about 17 m.y. ago, with slow uplift of the central Sierra Nevada summit region at rates estimated at about 0.012 mm/yr and of western Basin and Range Province at about 0.01 mm/yr. Many Mesozoic faults of the Foothills fault system were reactivated with normal slip in mid-Tertiary time and have continued to be active with slow slip rates. Sparse data indicate acceleration of rates of uplift and faulting during the Late Cenozoic. The Basin and Range faulting appears to have extended westward during this period with a reduction in width of the Sierra Nevada.The eastern boundary zone of the Sierra Nevada has an irregular en-echelon pattern of normal and right-oblique faults. The area between the Sierra Nevada and the Walker Lane is a complex zone of irregular patterns of hörst and graben blocks and conjugate normal-to right- and left-slip faults of NW and NE trend, respectively. The Walker Lane has at least five main strands near Walker Lake, with total right-slip separation estimated at 48 km. The NE-trending left-slip faults are much shorter than the Walker Lane fault zone and have maximum separations of no more than a few kilometers. Examples include the 1948 and 1966 fault zone northeast of Truckee, California, the Olinghouse fault (Part III) and possibly the almost 200-km-long Carson Lineament.Historic geologic evidence of faulting, seismologic evidence for focal mechanisms, geodetic measurements and strain measurements confirm continued regional uplift and tilting of the Sierra Nevada, with minor internal local faulting and deformation, smaller uplift of the western Basin and Range Province, conjugate focal mechanisms for faults of diverse orientations and types, and a NS to NE—SW compression axis (σ1) and an EW to NW—SE extension axis (σ3).  相似文献   

7.
This paper describes a new procedure for assessing the ratio between in situ stresses in rock masses by means of K (K = σH / σv, being σH and σv principal stress) and tectonics for purposes of engineering geology and rock mechanics. The method combines the use of the logic decision tree and the empirical relationship between the Tectonic Stress Index, TSI, and a series of K in situ values obtained from an extensive database. The decision tree considers geological and geophysical factors affecting stress magnitudes both on the regional and local scale. The TSI index is defined by geological and geomechanical parameters. The method proposed provides an assessment of the magnitude of horizontal stresses of tectonic origin. Results for several regions of Europe are presented and the possible applications of the procedure are discussed.  相似文献   

8.
Neotectonic evolution of the Central Betic Cordilleras (Southern Spain)   总被引：1，自引：0，他引：1  
Paleostress orientations were calculated from fault-slip data of 36 sites located along a traverse through the Central Betic Cordilleras (southern Spain). Heterogeneous fault sets, which are frequent in the area, have been divided into homogeneous subsets by cross-cutting relationships observed in the field and by a paleostress stratigraphy approach applied on each individual fault population. The state of stress was sorted according to main tectonic events and a new chronology is presented of the Miocene to Recent deformation in the central part of the Betic Cordilleras. The deviatoric stress tensors fall into four distinct groups that are regionally consistent and correlate with three Late Oligocene–Aquitanian to Recent major tectonic events in the Betic Cordilleras. The new chronology of the neotectonic evolution includes, from oldest to youngest, the following main tectonic phases:
(1) Late Oligocene–Aquitanian to Early Tortonian: σ1 subhorizontal N–S, partly E–W directed, σ3 subvertical; compressional structures (thrusting of nappes, large-scale folding) and strike-slip faulting in the Alborán Domain and the External Zone of the Betic Cordilleras;
(2) Early Tortonian to Pliocene–Pleistocene: σ1 subvertical, σ3 subhorizontal NW–SE, partly N–S directed or E–W-directed (radial extension); large-scale normal faulting in the Central Betic Cordilleras and in the oldest Neogene formations of the Granada Basin related to the gravitational collapse of the Betic Cordilleras and the exhumation of the intensely metamorphosed rock series of the Internal Zones, at the same time formation of the Alborán Basin and intramontane basins such as the Granada Basin;
(3) Pleistocene to Recent: (3a) σ1 subvertical, σ3 subhorizontal NE–SW with prominent normal faulting, but coevally; (3b) σ1 subhorizontal NW directed, σ3 NE–SW subhorizontal with strike-slip faulting. Extensional structures and strike-slip faulting are related to the ongoing convergence of the Eurasian and African Plates and coeval uplift of the Betic Cordilleras. Reactivation of pre-existing fractures and faults was frequently observed. Phase 3 is interpreted as periodic strike-slip and normal faulting events due to a permutation of the principal stress axes, mainly σ1 and σ2.
Keywords: Neotectonics; Paleostress; Fault-slip data; Deformation history; Betic Cordilleras  相似文献   

9.
Summary. Most currently used techniques for analysing the stability of near surface structures, such as rock slopes, are based on the application of the effective Coulomb shear strength parameters cohesion c′, and the angle of friction φ′ on some known or anticipated shear surface subjected to an effective normal stress σ′n. The most widely used of these techniques are the variants of the method of slices and related upper bound techniques. If the Hoek-Brown criterion is to be used to model the strength of near surface fractured rocks, it is necessary to determine equivalent Coulomb shear strength parameters for the specified level of effective normal stress. Calculation of the equivalent Coulomb parameters for the Hoek-Brown criterion for cases when a ≠ 0.5 is not a straightforward matter. A simple procedure for calculating instantaneous values of c′i and φ′i has been developed based on spreadsheet calculations and the application of a numerical optimisation routine. This procedure can also be applied to calculating the Hoek-Brown envelope plotted in shear stress/normal stress space. A simple closed form solution for c′i and tan φ′i has also been developed for the special case when a = 1. A three-dimensional version of the Hoek-Brown criterion has been developed by combining it with the Drucker-Prager criterion. This new yield criterion has been implemented by numerical solution of the governing equations. A simplification of this three-dimensional yield criterion has been developed by introducing an intermediate principal stress weighting factor. Comparison with published results demonstrates that this simplified criterion has the capacity to model the results of true triaxial tests for a range of different rock types over a wide range of stress levels. The new three-dimensional yield criterion has the advantage that its input parameters can be determined from routine uniaxial compression tests and mineralogical examination.  相似文献   

10.
Field data from an orogenic foreland and an orogenic belt (the Mesozoic rocks of southern England and the Umbria‐Marche Apennines of Italy respectively) indicate the following. Firstly, stress evolution during the tectonic cycle, between maximum compressive stress (σ1) being vertical during extension and least compressive stress (σ3) being vertical during contraction, can involve phases when the intermediate compressive stress (σ2) is vertical, promoting strike‐slip deformation. Secondly, variations in the relative magnitudes of the stress axes are caused by variations in overburden and tectonic forces. Thirdly, overpressure can develop because of compaction during burial, and, as overburden is reduced during uplift and erosion, the vertical stress (σV) reduces but fluid pressure (Pf) remains approximately constant. Brittle deformation, including transient strike‐slip faults, reverse‐reactivated normal faults and normal‐reactivated thrusts, is preferentially developed in overpressured areas because high Pf promotes faulting.  相似文献   

11.
In the last decade, even in areas that had been considered tectonically stable, a great amount of Cenozoic, including the Quaternary period, structural data have been collected throughout Brazil. The main goal of this study is to describe the Cenozoic structures and tectonic evolution of an area that is located at the border of the Paraná Basin in the state of São Paulo.The research methods consisted of the analysis of: (1) brittle structure data, mainly conjugate fractures and fault slip data; (2) lineaments traced on air photos and TM Landsat and radar images; and (3) a second-order base surface map.The study area, during the Cenozoic, has been affected by five strike–slip tectonic events, which generated mainly strike–slip faults, and secondarily normal and reverse ones. The events were named, from the oldest to the youngest, E1-NE, E2-EW, E3-NW, E4-NS, and E5-NNE; and the maximum principal stresses σ1 strike approximately NE–SW, E–W, NW–SE, N–S, and NNE–SSW, respectively. Event E2-EW seems to have been contemporaneous with the deposition of the Rio Claro Formation, the most important Cenozoic deposit of probable Neogenic age, and also to have controlled the distribution of its deposits. Event E3-NW was the strongest one in the area, as is pointed out by structural data, and the maximum principal stress σ1 of event E5-NNE is partially concordant with the orientation of σH-max of well break-out data in the Paraná Basin, suggesting a Neotectonic activity for this event. Finally, discontinuities parallel and correlated to the directions of strike–slip faults of the Cenozoic events seem to have actively controlled the sculpturing of the relief in the study area.  相似文献   

12.
Reactions of CO2 with carbonate and silicate minerals in continental sediments and upper part of the crystalline crust produce HCO3− in river and ground waters. H2SO4 formed by the oxidation of pyrite and reacting with carbonates may produce CO2 or HCO3−. The ratio, ψ, of atmospheric or soil CO2 consumed in weathering to HCO3− produced depends on the mix of CO2 and H2SO4, and the proportions of the carbonates and silicates in the source rock. An average sediment has a CO2 uptake potential of ψ = 0.61. The potential increases by inclusion of the crystalline crust in the weathering source rock. A mineral dissolution model for an average river gives ψ = 0.68 to 0.72 that is within the range of ψ = 0.63 to 0.75, reported by other investigators using other methods. These results translate into the CO2 weathering flux of 20 to 24 × 1012mol/yr.  相似文献   

13.
Fault slip analysis of Quaternary faults in southeastern Korea   总被引：1，自引：0，他引：1  
The Quaternary stress field has been reconstructed for southeast Korea using sets of fault data. The subhorizontal direction of the maximum principal stress (σ1) trended ENE and the direction of the minimum principal stress (σ3) was nearly vertical. The stress ratio (Φ = (σ2 − σ3) / (σ1 − σ3)) value was 0.65. Two possible interpretations for the stress field can be made in the framework of eastern Asian tectonics; (1) The σHmax trajectory for southeast Korea fits well with the fan-shaped radial pattern of maximum principal stress induced by the India–Eurasia collision. Thus, we suggest that the main source for this recent stress field in southeast Korea is related to the remote India–Eurasia continental collision. (2) The stress field in Korea shows a pattern similar to that in southwestern Japan. The origin for the E–W trending σHmax in Japan is known to be related to the mantle upwelling in the East China Sea. Thus, it is possible that Quaternary stress field in Korea has evolved synchronously with that in Japan. We suggest further studies (GPS and in situ stress measurement) to test these hypotheses.  相似文献   

14.
We estimated the stress fields of the aftershocks of the 2000 western Tottori earthquake (Mw 6.6) and the northern Hyogo swarm (max Mw 5.2) by a stress tensor inversion of moment tensor solutions reported from the National Research Institute for Earth Science and Disaster Prevention (Japan). The maximum principal stress direction of the western Tottori sequence was estimated as N107°E with a strike–slip regime. In the northern Hyogo swarm, the orientations of the principal stress directions could not be well constrained by the observed data, but after examining the detailed characteristics of the solution, we obtained a most probable solution of N113°E for the σ1 direction. These solutions are consistent with the maximum horizontal directions roughly estimated from the strike directions of large earthquakes occurring geographically between these two seismic activities. We measured the angle between each fault–slip direction and maximum principal stress direction to investigate the frictional properties of earthquakes. The distribution of the angles was forward modeled to estimate the coefficient of friction and the stress ratio, assuming uniformly distributed fault orientations. For the western Tottori sequence, a homogeneous stress field with a coefficient of friction less than 0.4 was estimated. A high stress level was also suggested because very little change occurred in the stress field during the mainshock. For the northern Hyogo sequence, the coefficient of friction was estimated to be between 0.5 and 1.0.  相似文献   

15.
A new method to determine stress directions using the preferential orientation of plagioclase mechanical twins has been applied to high-temperature mylonitic rocks from the Além Paraíba shear zone, Ribeira fold belt, southeastern Brazil. We have measured the lattice-preferred orientation of plagioclase grains and calculated the orientation of the stress axes possible for the observed twin orientations. The maximum compressive stress direction (σ1), determined for all studied samples, is a function of the mechanical twin orientations of a number of distinct plagioclase populations. The σ1 direction is generally subperpendicular to the (010) plane. The statistical treatment for most of the plagioclase grains examined for each sample shows that σ1 is almost perpendicular to the foliation plane, suggesting a significant coaxial component in the deformation process of these rocks.  相似文献   

16.
A rare metachert pebble containing amphibole grains with microboudin structures in a wide range of orientations provides an opportunity to perform stress analysis in two orthogonal orientations on the foliation surface. The sample was analysed by the microboudin method to infer the triaxial stress state during microboudinage. Stress parameters proportional to the far-field differential stress were determined for sodic amphibole grains in the two orientations. The ratio of the stresses in the two orthogonal orientations (σ1−σ2)/(σ1−σ3) was calculated to be 0.64, indicating that σ2 lies closer to the midpoint between σ1 and σ3 than to σ3.  相似文献   

17.
To understand the generation mechanism of the Bam earthquake (Mw 6.6), we studied three-dimensional VP, VS and Poisson's ratio (σ) structures in the Bam area by using the seismic tomography method. We inverted accurate arrival times of 19490 P waves and 19015 S waves from 2396 aftershocks recorded by a temporal high-sensitivity seismic network. The 3-D velocity structure of the seismogenic region was well resolved to a depth of 14 km with significant velocity variations of up to 5%. The general pattern of aftershock distribution was relocated by using the 3-D structure to delineate a source fault for a length of approximately 20 km along a line 4.5 km west of the known geological Bam fault; this source fault dips steeply westward and strikes a nearly north–south line. The main shallow cluster of aftershocks south of the city of Bam is distributed just under the minor surface ruptures in the desert. The 3-D velocity structure shows a thick layer of high VS and low σ (minimum: 0.20) at a depth range of 2–6 km. The deeper layer, with a thickness of about 2 km, appears to have a low VS and high σ (maximum: 0.28) from 6 km depth beneath Bam to a depth of 9 km south of the city. The inferred increase of Poisson's ratio from 2 to 10 km in depth may be associated with a change from rigid and SiO2-rich rock to more mafic rock, including the probable existence of fluids. The main seismic gap of aftershock distribution at the depth range of 2 to 7 km coincides well with the large slip zone in the shallow thick layer of high VS and low σ. The large slip propagating mainly in the shallow rigid layer may be one of the main reasons why the Bam area suffered heavy damage.  相似文献   

18.
Klaus-G. Hinzen   《Tectonophysics》2003,377(3-4):325-356
Fault plane solutions (FPS) from 110 earthquakes in the northern Rhine area with local magnitudes, ranging from 1.0 to 6.1, and occurring between 1976 and 2002 are determined. FPS are retrieved from P-wave first motions using a grid search approach allowing a detailed exploration of the parameter space. The influence of the 1D velocity model on take-off angles and resulting FPS is examined. All events were relocated with a recently developed minimum 1D model of the velocity structure [J. Geophys. Res. (2003)]. Rose diagrams of the orientation of P, T and B axes show a clear preference of trends of P and T axes at N292°E and N27°E, respectively. The majority of B axes trend in northerly directions. Plunges of P and T axes are mostly around 45° while most B axes are subhorizontal. The main direction of the maximum horizontal stress directly inferred from the fault plane solutions is N118°E.To calculate the orientations of the principal stress axes and the shape of the stress tensor, the inversion method of Gephard and Forsyth [J. Geophys. Res. 89 (1984) 9305] was applied to the whole data set and to several subsets of data. The subsets were formed by grouping events from various geological and tectonic areas and by grouping events into different depth ranges. The subset areas include the Lower Rhine Embayment, the Rhenish Massif, the middle Rhine area, the Neuwied Basin and the area known as the Stavelot–Venn Massif. Inversion of the entire data set shows some ambiguity between a strike-slip and extensional stress regime, with a vertical axis for the medium principal stress and a trend of N305°E and N35°E for the σ1 and σ3 axis, respectively, as the best fitting tensor. Earthquakes from the Lower Rhine Embayment and, to some degree, from the middle Rhine area indicate an extensional stress regime. In the Lower Rhine Embayment, plunge and trend of the σ1 axis are 76° and N162°E and for the σ3 axis 7° and N42°E. The best fitting solution for the area of the Stavelot–Venn Massif is a strike-slip regime with subhorizontal σ1 and σ3 axes with a trend of N316°E and N225°E, respectively. Stress orientations found here agree overall with the results from earlier studies based on smaller data sets. The directions of the maximum and minimum horizontal stresses inverted from focal mechanisms agree well with the stress field predicted by the European Stress Map. This confirms earlier interpretations that the stress field of the Rhine Graben system is controlled by plate driving forces acting on the plate boundaries. However, amplitudes of the stresses change on a local scale and with depth. Estimates of the absolute magnitude of principal stresses favor a normal faulting regime in the shallow crust (above 12-km depth) and a strike-slip regime in the lower crust.  相似文献   

19.
The magnetic susceptibility anisotropy of 275 specimens comprising 38 sites from the Cambrian slate belt in North Wales was measured to determine the magnetic fabric of the slates. The susceptibility ellipsoid is oblate for all sites, and the maximum/intermediate susceptibility plane always coincides with the cleavage plane of the slates which has a Caledonian strike and is nearly vertical. The maximum axes align sub-vertically and the intermediate axes sub-horizontally, trending NE-SW. The minimum susceptibility axes are normal to this foliation plane and coincide with the poles to the slaty cleavage. The orientations of the principal susceptibility axes are found to be in excellent agreement with the orientations of the principal strain directions, determined by X-ray goniometry on one of the samples from almost all of the sites. Correlation of the magnetic susceptibility anisotropy with predicted March strains (March, 1932) shows that the principal magnitudes of susceptibility can be related to those of the strain by: (for i = 1, 2, 3. The orthogonal principal axes), where χf and χ0 are the final and initial susceptibilities along a given axis i and lf and li are final and initial axial dimensions in the same direction i of a principal strain axis. The exponent a for the North Wales slates was found to be 0.145 ± 0.005. Knowledge of such a relationship may permit rapid approximate determinations of a petrofabric in similar rocks from their magnetic fabrics. However, the exponent a will probably have to be recalibrated for each rock type.  相似文献   

20.
The Húsavík–Flatey Fault (HFF) is an oblique dextral transform fault, part of the Tjörnes Fracture Zone (TFZ), that connects the North Volcanic Zone of Iceland and the Kolbeinsey Ridge. We carry out stress inversion to reconstruct the paleostress fields and present-day stress fields along the Húsavík–Flatey Fault, analysing 2700 brittle tectonic data measured on the field and about 700 earthquake focal mechanisms calculated by the Icelandic Meteorological Office. This allows us to discuss the Latest Cenozoic finite deformations (from the tectonic data) as well as the present-day deformations (from the earthquake mechanisms). In both these cases, different tectonic groups are reconstructed and each of them includes several distinct stress states characterised by normal or strike-slip faulting. The stress states of a same tectonic group are related through stress permutations (σ1−σ2 and σ2−σ3 permutations as well as σ1−σ3 reversals). They do not reflect separate tectonic episodes. The tectonic groups derived from the geological data and the earthquake data have striking similarity and are considered to be related. The obliquity of the Húsavík–Flatey Fault implies geometric accommodation in the transform zone, resulting mainly from a dextral transtension along an ENE–WSW trend. This overall mechanism is subject to slip partitioning into two stress states: a Húsavík–Flatey Fault-perpendicular, NE–SW trending extension and a Húsavík–Flatey Fault-parallel, NW–SE trending extension. These three regimes occur in various local tectonic successions and not as a regional definite succession of tectonic events. The largest magnitude earthquakes reveal a regional stress field tightly related to the transform motion, whereas the lowest magnitude earthquakes depend on the local stress fields. The field data also reveal an early extension trending similar to the spreading vector. The focal mechanism data do not reflect this extension, which occurred earlier in the evolution of the HFF and is interpreted as a stage of structural development dominated by the rifting process.  相似文献   

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

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