首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Five erosion devices were compared using five intertidal estuarine sites covering a range of sediment stability from newly settled mud to very cohesive mud at the margins of a saltmarsh. The erosion devices use different methods of fluid shearing from horizontal currents/bed shear stresses to vertical water jets, and have different ‘footprint’ areas. The devices included: (1) the annular flumes (AFs—diameter 64 cm; footprint area 0.17 m2) of the Plymouth Marine Laboratory (PML); (2) PML's mini-annular flume (MAF—diameter 19 cm; area 0.026 m2); (3) the annular mini-flume (AMF—diameter 30.5 cm; area 0.032 m2) of the National Oceanography Centre Southampton (NOC); (4) NOC's Cohesive Strength Meter (CSM—diameter 3 cm; area 0.0007 m2); (5) NOC's EROMES (ER—diameter 10 cm; area 0.0079 m2). The quantification of threshold shear stress for bed erosion (τe) and sediment erosion rate was complemented by the measurement of physical, chemical and biological properties of the sediment (grain size, bulk density, water content, organic content, chlorophyll a, carbohydrates, macrofauna). The results demonstrated a significant correlation (r2=0.98) between the PML AF (laboratory measurement of undisturbed cored sediment) and PML MAF (in situ) for measurement of erosion thresholds for bed sediment. However, there were no significant correlations between AFs, the CSM and EROMES. There were no consistent correlations with physical or biological sediment properties due to the spatially unrelated sites and the marked differences in benthic assemblages. The sources of differences and the lack of correlations between erosion devices were due to several factors, including operational procedures (e.g., sediment resuspension during filling with water), definition of erosion threshold, the nature of the force applied to the bed, and method of calibration. In contrast to the CSM and EROMES, both types of AFs were able to record significant differences in the erodability of soft sediments from four sites. This indicates that the CSM and EROMES may not be very effective at measuring the differences in erosion thresholds of soft estuarine sediments.  相似文献   

2.
The reaction of the lower ionosphere to the solar proton events that occurred in 2011–2012 is studied in this paper based on the results of measurements of the propagation velocity and the E z /H τ ratio of the low-frequency electromagnetic pulses (atmospherics) in the ELF range at the high-latitude observatories Lovozero and Barentsburg. With numerical modeling methods, it is shown that horizontal local irregularities of the lower ionosphere conductivity profile could be a cause of the splashes in the E z /H τ ratio observed in the experiment during the solar proton event of March 7, 2012, which was a unique event in both the proton flux value and energy.  相似文献   

3.
The observations of the variations in the vertical component of the atmospheric electric field (E z ) at Swider midlatitude Poland observatory (geomagnetic latitude 47.8°) under the conditions of fair weather during 14 magnetic storms have been analyzed. The effect of the magnetic storm main phase in the daytime midlatitude variations in E z in the absence of local geomagnetic disturbances has been detected for the first time. Considerable (~100–300 V m?1) decreases in the electric field strength (E z ) at Swider observatory were observed in daytime simultaneously with the substorm onset in the nighttime sector of auroral latitudes (College observatory). The detected effects indicate that an intensification of the interplanetary electric field during the magnetic storm main phase, the development of magnetospheric substorms, and precipitation of energetic electrons into the nighttime auroral ionosphere can result in considerable disturbances in the midlatitude atmospheric electric field.  相似文献   

4.
A simple law of wear rate is examined for the process of gouge generation during the frictional sliding of simulated faults in rocks, by use of the Pi theorem method (dimensional analysis) and existing experimental data. The relationship between wear rate (t/d) and the applied stress can be expressed by the power-law relations $$\frac{t}{d} = C_\sigma \sigma ^{m\sigma } ,\frac{t}{d} = C_\tau \tau ^{m\tau }$$ wheret is the thickness of the gouge generated on the frictional surfaces,d is the fault displacement, σ and τ are normal stress and shear stress, respectively, andC σ,C τ,m σ andm τ are constants. These results indicate that the exponent coefficientsm σ andm τ and the coefficientsC σ andC τ depend on the material hardness of the frictional surfaces. By using the wear rates of natural faults, these power-law relationships may prove to be an acceptable palaeopiezometer of natural faults and the lithosphere.  相似文献   

5.
Semi-diurnal and fortnightly surveys were carried out to quantify the effects of wind- and navigation-induced high-energy events on bed sediments above intertidal mudflats. The mudflats are located in the upper fluvial part (Oissel mudflat) and at the mouth (Vasière Nord mudflat) of the macrotidal Seine estuary. Instantaneous flow velocities and mudflat bed elevation were measured at a high frequency and high resolution with an acoustic doppler velocimeter (ADV) and an ALTUS altimeter, respectively. Suspended particulate matter concentrations were estimated by calibrating the ADV acoustic backscattered intensity with bed sediments collected at the study sites. Turbulent bed shear stress values were estimated by the turbulent kinetic energy method, using velocity variances filtered from the wave contribution. Wave shear stress and maximum wave–current shear stress values were calculated with the wave–current interaction (WCI) model, which is based on the bed roughness length, wave orbital velocities and the wave period (TS). In the fluvial part of the estuary, boat passages occurred unevenly during the surveys and were characterized by long waves (TS>50 s) induced by the drawdown effect and by short boat-waves (TS<10 s). Boat waves generated large bottom shear stress values of 0.5 N m−2 for 2–5 min periods and, in burst of several seconds, larger bottom shear stress values up to 1 N m−2. At the mouth of the estuary, west south-west wind events generated short waves (TS<10 s) of HS values ranging from 0.1 to 0.3 m. In shallow-water environment (water depth <1.5 m), these waves produced bottom shear stress values between 1 and 2 N m−2. Wave–current shear stress values are one order of magnitude larger than the current-induced shear stress and indicate that navigation and wind are the dominant hydrodynamic forcing parameters above the two mudflats. Bed elevation and SPM concentration time series showed that these high energy events induced erosion processes of up to several centimetres. Critical erosion shear stress (τce) values were determined from the SPM concentration and bed elevation measurements. Rough τce values were found above 0.2 N m−2 for the Oissel mudflat and about 1 N m−2 for the Vasière Nord mudflat.  相似文献   

6.
Limited field and flume data suggests that both uniform and graded beds appear to progressively stabilize when subjected to inter-flood flows as characterized by the absence of active bedload transport. Previous work has shown that the degree of bed stabilization scales with duration of inter-flood flow, however, the sensitivity of this response to bed surface grain size distribution has not been explored. This article presents the first detailed comparison of the dependence of graded bed stability on inter-flood flow duration. Sixty discrete experiments, including repetitions, were undertaken using three grain size distributions of identical D50 (4.8 mm); near-uniform (σg = 1.13), unimodal (σg = 1.63) and bimodal (σg = 2.08). Each bed was conditioned for between 0 (benchmark) and 960 minutes by an antecedent shear stress below the entrainment threshold of the bed (τ*c50). The degree of bed stabilization was determined by measuring changes to critical entrainment thresholds and bedload flux characteristics. Results show that (i) increasing inter-flood duration from 0 to 960 minutes increases the average threshold shear stress of the D50 by up to 18%; (ii) bedload transport rates were reduced by up to 90% as inter-flood duration increased from 0 to 960 minutes; (iii) the rate of response to changes in inter-flood duration in both critical shear stress and bedload transport rate is non-linear and is inversely proportional to antecedent duration; (iv) there is a grade dependent response to changes in critical shear stress where the magnitude of response in uniform beds is up to twice that of the graded beds; and (v) there is a grade dependent response to changes in bedload transport rate where the bimodal bed is most responsive in terms of the magnitude of change. These advances underpin the development of more accurate predictions of both entrainment thresholds and bedload flux timing and magnitude, as well as having implications for the management of environmental flow design. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.  相似文献   

7.
Erodibility of cohesive streambeds in the loess area of the midwestern USA   总被引:2,自引:0,他引:2  
Excess stress parameters, critical shear stress (τc) and erodibility coefficient (kd), for degrading channels in the loess areas of the midwestern USA are presented based on in situ jet‐testing measurements. Critical shear stress and kd are used to define the erosion resistance of the streambed. The jet‐testing apparatus applies hydraulic stresses to the bed and the resulting scour due to the impinging jet is related to the excess stress parameters. Streams tested were primarily silt‐bedded in texture with low densities, which is typical of loess soils. Results indicate that there is a wide variation in the erosion resistance of streambeds, spanning six orders of magnitude for τc and four orders of magnitude for kd. Erosion resistance was observed to vary within a streambed, from streambed to streambed, and from region to region. An example of the diversity of materials within a river system is the Yalobusha River Basin in Mississippi. The median value of τc for the two primary bed materials, Naheola and Porters Creek Clay Formations, was 1·31 and 256 Pa, respectively. Streambeds composed of the Naheola Formation are readily eroded over the entire range of shear stresses, whereas only the deepest flows generate boundary stresses great enough to erode streambeds composed of the Porters Creek Clay Formation. Therefore, assessing material resistance and location is essential in classifying and modelling streambed erosion processes of these streams.  相似文献   

8.
This paper investigates dynamics of a spherical bubble surrounded by a viscoelastic fluid. The purpose of the study is to understand the parameters which control expansion and fragmentation of bubbly magma by decompression. In particular, we focus on which occurs first, fragmentation or expansion. Supposing that rupture of the bubble wall occurs in a critical stress condition, we calculate the change of the bubble radius and tensile stress at the bubble wall for various decompression rates. Conditions in which tensile stress is stored in the shell are represented in terms of dimensionless parameters. The results are interpreted as follows: when magma viscosity is larger than a critical value, and the decompression time is shorter than viscous expansion time, tensile stress is stored before expansion; when magma viscosity is smaller than the critical value, tensile stress is not stored, no matter how rapid the decompression. Although it is a generally accepted theory that fragmentation is effected by stress conditions and decompression time, exactly how decompression time (t1) effects the fragmentation is not yet fully understood. This study demonstrates that the stress condition is controlled by the length of the decompression time not relative to the viscoelastic relaxation time (t1 / τ), but relative to the viscous expansion time (t1 / τlrlx). As suggested by recent experimental studies, the decompression time relative to viscoelastic relaxation time (t1 / τ) is also significant to the fragmentation process itself. It indicates that the decompression time effects the fragmentation not through the stress condition. However more work must be completed to fully understand the particular relationship between the decompression time and relaxation time in terms of its influence on fragmentation.  相似文献   

9.
Measurements on acoustic pulses propagating in massive rock lead to a simple empirical relationship between the pulse rise time, τ and the time of propagation of a pulse, t:
τ=τ0+C)TQ?1dt
where τ0 is the initial rise time (at t = 0), Q is the anelastic parameter which may be expressed in terms of the fractional loss of energy per cycle of a sinusoidal wave, Q = 2π(ΔE/E)?1, and is assumed to be essentially independent of frequency, and C is a constant whose value we estimate experimentally to be 0.53 ± 0.04. Of the linear theories of seismic pulse attenuation, model 2 of Azimi et al. (1968) is favoured. Pulse shapes computed from equations of Futterman (1962) also give C = 0.5, but the pulse arrives earlier than in a non-attenuating medium with the same elasticity and density. Pulse shapes calculated using Strick's (1967, 1970, 1971) theory give values of C incompatible with our results. The observations suggest that a method of estimating the Q-structure of the earth from seismic pulse rise times may have a particular advantage over the spectral ratio method.  相似文献   

10.
Recent research has started to focus on how prolonged periods of sub‐threshold flows may be capable of imparting structural changes that contribute to increased bed stability. To date, this effect (termed ‘stress history’) has been found to be significant in acting to increase a bed's critical shear stress at entrainment threshold. However, it is supported by only limited, qualitative and often speculative information on the mechanisms of this stabilization process in grade‐specific studies. As such, this paper uses high resolution laser scanning to quantitatively ascertain the granular mechanics underpinning the relationship between stress history and entrainment threshold for beds of a range of grain size distributions. Employing a bed slope of 1/200, three grain size distributions with median grain sizes (D50) of 4·8 mm [uniform (σg = (D84/D16)0.5 = 1·13; bimodal (σg = 2·08); and, unimodal (σg = 1·63)] were exposed to antecedent stress histories of 60 and 960 minutes duration. Antecedent shear stress magnitude was set at 50% of the critical shear stress for the D50 when no stress history period was employed. Two laser displacement scans of the bed surface (approximate area 100 mm × 117 mm) were taken, one prior to the antecedent period and one after this period, so that changes to surface topography could be quantified (resolution of x = 0·10 mm, y = 0·13 mm and z = 0·24 mm). Rearrangement of bed surface structure is described using statistical analysis and two‐dimensional (2D) semi‐variograms to analyse scaling behaviour. Results reveal vertical settlement, changes to bed roughness and particle repositioning. However, the bed grain size distribution influences the relative importance of each mechanism in determining stress history induced bed stability; this is the focus of discussion in this paper. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

11.
Vertical profiles of the critical erosion threshold (τcrit) in sediment have been measured at 11 stations along the axis of the Tamar Estuary and at a single station in a tributary of the Tamar at St. John's Ford. The τcrit of surface sediment increased from 0.04 Pa in the upper, brackish estuary to 0.09 Pa in the lower estuary. In the upper estuary τcrit only increased slightly with depth whereas in the marine estuary τcrit increased rapidly from 0.09 Pa at the surface to 0.25 Pa at 15 cm below the sediment surface. The results showed that the relationship between τcrit and bulk density (ρb) obtained previously for surface sediment was also applicable to sediments from depths of 10–15 cm and probably deeper. Profiles of ρb were measured to depths of 70 cm using a corer. In the lower (marine) estuary ρb increased with depth in the sediment from 1580 kg m−3 at the surface to 1720 kg m−3 at 70 cm. In the upper estuary ρb values were lower at 1170–1200 kg m−3 and profiles were almost homogeneous indicating that consolidation was not occurring. The mid-estuary was transitional between these two situations. These results are consistent with the seasonal accumulation and loss of ‘mobile’ sediment observed previously in the upper estuary with changes in river flow, and with the apparent stability of intertidal mud in the lower marine estuary deduced from historical bathymetric survey records. The slopes of the intertidal mud banks ranged from 1–2% in the lower estuary to 20–25% in mid-estuary but, instead of continuing to increase in steepness towards the head as the estuary became narrower, the measured slopes actually decreased. It is speculated that the lack of consolidation through continual mobilisation and settlement cycles combined with an increase in silt content in the upper estuary resulted in sediment that lacked the mechanical strength to maintain steep slopes.  相似文献   

12.
Statistically significant (at the 95% significance level) changes in daily cloud cover are found to occur globally over land coincident with extreme increases in ‘fair-weather’ measurements of vertical electric field (Ez) measured at Vostok, Antarctica. Using global cloud products from the International Satellite Cloud Climatology Project (ISCCP) D1 data series, superposed epoch analyses were made of both increases and decreases in Ez. Field significance testing revealed that, both before and after extreme increases in Ez, significant absolute cloud cover changes (of 13–15%) occur in the tropics and high latitudes. While the linkages in the tropics may reflect changes in the main convective cloud generators of current flow in the global circuit, the linkages at high latitudes appear to represent responses of clouds to the current flow. This linkage offers a possible explanation of a possible solar–terrestrial climate amplification mechanism.  相似文献   

13.
Abstract

River science and management often require a design or reference discharge. The common (and sometimes unavoidable) use of such discharges may, however, obscure the fact that the magnitude and frequency of critical flows can differ due to various hydrological, geomorphological, and ecological criteria. Threshold stages and discharges were identified for six lower Brazos River, Texas gaging stations corresponding to thalweg connectivity, bed inundation, high sub-banktop flows, channel–floodplain connectivity (CFC), and overbank flooding. Critical flows were also identified for estimated thresholds for sandy bedform and medium gravel mobility, critical specific stream power for potential channel modifications, and cohesive-bank channel erosion. These thresholds have variable relationships to mean, median, and maximum flows. For four of the six stations, daily recurrence probabilities for all but flood flows are at least 1%, and as high as 11%. All stations achieve channel–floodplain connectivity at stages less than banktop. Estimated threshold flows for sediment mobility and channel erosion occur relatively frequently, with daily probabilities of 2–77%. Critical flows for bank erosion occur least often, and for sandy bedform and gravel mobility most often. Thalweg connectivity is always maintained at all sites, while bed inundation flows have a daily probability of about 80% or more. Overall, results suggest that no single flow level is dominant in hydrological or geomorphic dynamics, and that the frequency of a given threshold varies considerably even along a single river. The results support the idea that multiple flow levels and ranges are necessary to create and maintain the hydrological, geomorphological, and ecological characteristics of rivers, and that no single flow level is a reliable determinant of fluvial state.
Editor Z.W. Kundzewicz; Associate editor Q. Zhang  相似文献   

14.
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 τAl/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.  相似文献   

15.
We present the distribution of226Ra in eight vertical profiles from the eastern Pacific. The profiles are located along a meridional trend near 125°W, from 43°S to 29°N. Surface226Ra concentrations are about 7 dpm/100 kg, except for the two stations south of 30°S where the higher values are due to the Antarctic influence. Deep waters show a distinctive south-to-north increase in the226Ra content, from about 26 to 41 dpm/100 kg near the bottom. Unlike in the Atlantic and Antarctic Oceans, the effect of226Ra injection from bottom sediments is clearly discernible in the area. The presence of this primary226Ra can be traced up to at least 1–1.5 km above the ocean floor, making this part of the sea bed among the strongest source regions for the oceanic226Ra. Numerical solutions of a two-dimensional vertical advection-diffusion model applied to the deep (1.2–4 km)226Ra data give the following set of best fits: upwelling velocity(Vz) = 3.5m/yr, vertical eddy diffusivity(Kz) = 0.6cm2/s, horizontal (north-south) eddy diffusivity(Ky) = 1 × 107cm2/s, and water-column regeneration flux of226Ra(J) = 3.3 × 10?5dpmkg?1yr?1 as an upper limit. These parametric values are in general agreement with one-dimensional (vertical) model fits for the Ra-Ba system. However, consideration of226Ra balance leads us to suspect the appropriateness of describing the vertical exchange processes in the eastern Pacific with constantVz and Kz. If future modeling is attempted, it may be preferable to treat the area as a diffusion-dominant mixing regime with depth-dependent diffusivities.  相似文献   

16.
Laboratory tests using Jet Erosion Testing (JET) apparatus, impinging normally on a horizontal boundary, were conducted to determine the critical shear stress (τc) of non‐cohesive soil samples. A three‐dimensional (3D) SonTek/YSI 16 MHz Micro‐Acoustic Doppler Velocimeter (MicroADV) was used to measure turbulent kinetic energy (TKE) at a radial limit of entrainment in the wall jet zone and the measurements were used to calculate τc of the samples. The results showed that TKE increases exponentially with increasing particle size. The τc from this study were comparable (R2 = 0.8) to the theoretical τc from Shields diagram after bed roughness scale ratio (D/ks), due to the non‐uniform bed conditions, was accounted for. This study demonstrated that JET and TKE can be used to determine τc of non‐cohesive soils. The use of JET and TKE was found to be faster and easier when compared to the conventional approach of using flumes. A relationship of TKE at the onset of incipient motion (TKEc) and samples’ D50 developed in this study can be used to predict τc of non‐cohesive soils under similar non‐uniform conditions. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

17.
Z. Shi  J. M. R. Hughes 《水文研究》2002,16(16):3279-3289
The microflow environments of aquatic plants with reference to Myriophyllum and Hydrilla are simulated in a laboratory flume. A Nix Streamflow microflow meter was used to measure the mean velocity profiles of flow at different densities of plants, flow ranges and measurement positions. Each mean velocity profile consists of three hydrodynamic regimes (i.e. within‐canopy zone, above‐canopy zone and a transitional zone between them), which indicate the presence of two benthic boundary layers (internal and external ones). Out of 38 measured mean velocity profiles, most do not fit a logarithmic relationship. The following hydrodynamic parameters are used in characterizing the flow regimes: local shear velocity (u*), roughness length (zo), canopy roughness Reynolds number (Re*), bed shear stress (τo) and laminar sublayer (σ). Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

18.
The time variations in the CR geomagnetic cutoff rigidity and their relation to the interplanetary parameters and the Dst index during a strong magnetic storm of November 18–24, 2003, have been analyzed. The Tsyganenko (Ts03) model of a strongly disturbed magnetosphere [Tsyganenko, 2002a, 2002b; Tsyganenko et al., 2003] have been used to calculate effective geomagnetic thresholds with the help of the method for tracing CR particle trajectories in the magnetospheric magnetic field. The geomagnetic thresholds have been calculated using the method of global spectrographic survey (GSS), based on the data from the global network of CR stations, and the results have been compared with the effective geomagnetic cutoff rigidities. The daily anisotropy of effective geomagnetic thresholds during the Dst variation minimum have been estimated. The relation of the theoretical and experimental geomagnetic thresholds, obtained using the GSS method, to the interplanetary parameters and Dst variation is analyzed. The Dst variations, IMF B z , and solar wind density are most clearly defined in the geomagnetic thresholds during this storm. The correlation between B y and experimental geomagnetic thresholds is higher than such a correlation between this parameter and theoretical thresholds by a factor 2–3, which suggests that a real dawn-dusk asymmetry during this storm was stronger than such an asymmetry represented by the Ts03 model.  相似文献   

19.
The activation coefficient equations in the"activation criterion of pre-existing weakness"are relatively complex and not easy to apply to specific applications.The relative activity of pre-existing weaknesses is often critical in geological analysis.The Mohr circle can be used only in two-dimensional stress analysis.By applying the"activation criterion of pre-existing weakness"and combining it with numerical analysis,we establish the correspondence between the pole(n,n)of a pre-existing weakness plane and its orientation in"Mohr space".As a result,the normal stress(n)and shear stress(n)of a pre-existing weakness plane can be readily expressed in Mohr space.Furthermore,we introduce the method and procedures for predicting the activation and relative activation of pre-existing weaknesses in Mohr space.Finally,we apply the Mohr space method and compare the predictions to sandbox modeling results and 3D seismic data.The results show that Mohr space can be used in stress analysis to estimate the activation of a pre-existing weakness in any triaxial stress state.  相似文献   

20.
Erosion of cohesive soils in fluvial environments is dependent on physical, geochemical and biological properties, which govern inter‐particle attraction forces and control detachment rates from stream beds and banks. Most erosion rate models are based on the excess shear stress equation where the soil erodibility coefficient (kd) is multiplied by the difference between the boundary hydraulic shear stress (τb) and the soil critical shear stress (τc). Both kd and τc are a function of soil properties and must be obtained through in situ field or laboratory testing. Many studies have generated predictive relationships for kd and τc derived from various soil properties. These studies typically were conducted in watersheds within a single physiographic region with a common surficial geology and/or investigated a limited number of soil properties, particularly geochemical properties. With widely reported differences in relationships between τc and soil properties, this study investigated differences in predictive relationships for τc among different physiographic provinces in Tennessee, USA. Erodibility parameters were determined in the field using a mini‐jet test device. Among these provinces, statistically four unique clusters were identified from a dataset of 128 observations and these data clusters were used to develop predictive models for τc to identify dominant properties governing erosion. In these clusters, 16 significant physical and geochemical soil properties were identified for τc prediction. Among these soil properties, water content and passing #200 sieve (percentage soil less than 75 μm) were the dominant controlling parameters to predict τc in addition to clay percentage (< 2 μm), bulk density, and soil pore water chemistry. This study suggests that unique relationships exist for physiographic provinces that are likely due to soil physical‐geochemical processes associated with surficial geology that determine minerology of the cohesive soil. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

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

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