首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 437 毫秒
1.
Hydraulic properties of deeply weathered basement rocks and variably weathered sedimentary materials were measured by pumping and slug-test methods. Results from over 200 bores in 13 catchments, and eight pumping-test sites across the eastern and central wheatbelt of Western Australia were analysed. Measurements were made in each of the major lithological units, and emphasis placed on a ubiquitous basal saprolite aquifer. Comparisons were made between alternative drilling and analytical procedures to determine the most appropriate methods of investigation.

Aquifers with an average hydraulic conductivity of 0.55 m day−1 occur in variably weathered Cainozoic sediments and poorly weathered saprolite grits (0.57 m day−1). These aquifers are separated by an aquitard (0.065 m day−1) comprising the mottled and pallid zones of the deeply weathered profile. Locally higher values of hydraulic conductivity occur in the saprolite aquifer, although after prolonged periods of pumping the values decrease until they are similar to those obtained from the slug-test methods. Hydraulic conductivities measured in bores drilled with rotary auger rigs were approximately an order of magnitude lower than those measured in the same material with bores drilled by the rotary air-blast method.

Wheatbelt aquifers range from predominantly unconfined (Cainozoic sediments), to confined (saprolite grit aquifer). The poorly weathered saprolite grit aquifer has moderate to high transmissivities (4–50 m2 day−1) and is capable of producing from less than 5 to over 230 kl day−1 of ground water, which is often of a quality suitable for livestock. Yields are influenced by the variability in the permeability of isovolumetrically weathered materials from which the aquifer is derived.

The overlying aquitard has a low transmissivity (< 1 m2 day−1), especially when deeply weathered, indurated and silicified. The transmissivity of the variably weathered sedimentary materials ranges from less than 0.5 m2 day−1 to over 10 m2 day−1, depending on the texture of the materials and their position within the landscape. Higher transmissivity zones may occur as discrete layers of coarser textured materials. The salinity of the saprolite and sedimentary aquifers ranges from less than 2000 mgl−1 to greater than 250000 mgl−1 (total dissolved solids; TDS), depending on position within the landscape. Secondary soil salinization develops when groundwater discharge occurs from either saprolite or sedimentary aquifers.  相似文献   


2.
Microencapsulation of degradative organisms enhances microorganism survivability (Stormo and Crawford 1994). The use of encapsulated cell microbeads for in situ biodegradation depends not only on microorganism survival but also on microbead transport characteristics. Two forced-gradient, recirculating-loop tracer experiments were conducted to evaluate the feasibility of encapsulated cell transport and bioremediation on the basis of polystyrene microsphere transport results. The tracer tests were conducted in a shallow, confined, unconsolidated, heterogeneous, sedimentary aquifer using bromide ion and 2 μm, 5 μn, and 15 μm microsphere tracers. Significant differences were observed in the transport of bromide solute and polystyrene microspheres. Microspheres reached peak concentrations in monitoring wells before bromide, which was thought to reflect the influence of aquifer heterogeneity. Greater decreases in microsphere C/Co ratios were observed with distance from the injection wells than in bromide C/Co ratios, which was attributed to particle filtration and/or settling. Several methods might be considered for introducing encapsulated cell microbeads into a subsurface environment, including direct injection into a contaminated aquifer zone, injection through a recirculating ground water flow system, or emplacement in a subsurface microbial curtain in advance of a plume. However, the in situ use of encapsulated cells in an aquifer is probably limited to aquifers containing sufficiently large pore spaces, allowing passage of at least some encapsulated cells. The use of encapsulated cells may also be limited by differences in solute and microbead transport patterns and flowpath clogging by larger encapsulated cell microbeads.  相似文献   

3.
The laterites in Burundi, which are formed by weathering of ultrabasic rocks, show a complete profile with the following horizons: canga, the ferruginous crust capping, ferralite, consisting essentially of iron hydroxides, and saprolite, which contains a large quantity of hydrosilicate minerals. Nickel bearing minerals occur in the saprolite and the lower portion of ferralite. Resistivity well-logging and resistivity sounding indicated that the electrical properties of rocks depend upon their composition: Canga and ferralite showed high resistivities of 6,500 Ωm and 800 Ωm, respectively. The resistivity of saprolite was found to be much lower, between 10 Ωm and 20 Ωm. The laterite is underlain by resistive peridotite. The chargeability of saprolite was found to be lower than that of the upper horizons and the bedrock. Electrolytic conductivity of laterite, which depends on the geometry of the deposit, was found to be low, because the laterite contains moisture and ground water, which are highly resistive. The relatively high conductivity of saprolite is caused by nickeliferous hydrosilicates, which exhibit the electrical properties of clay minerals, with an apparent maximum conductivity of 0.25 S/m. The conductivity of saprolite corresponds to a concentration between 30% and 50% of conductive silicate minerals distributed in the pore space of deposit. A nickel enrichment of up to 6% was estimated from the resistivity of the saprolite. Prospecting for laterites by electrical sounding showed that the development of laterite horizons in a nickel deposit correlates with the surface morphology of weathered ultrabasic massif. Thus the method can be used in preliminary exploration of such deposits.  相似文献   

4.
Applicability of spectral analysis to determine hydraulic diffusivity   总被引：1，自引：1，他引：0  
This study is to evaluate the applicability of estimating the one-dimensional horizontal hydraulic diffusivity of an unconfined aquifer with time-dependent fluctuation of lateral head and vertical recharge boundaries using observed water level spectra. Different models of boundary condition are imposed to evaluate the statistical significance between the calculated hydraulic diffusivity (ξ′) with the given hydraulic diffusivity (ξ). The auto-spectra of the water level in observation wells tapping the same aquifer are closely related to those at the disturbed boundaries. For an aquifer with a constant hydraulic diffusivity, the water level fluctuation in the monitoring wells is linearly related to the water level spectra observed at the boundaries. The spectral density function of aquifer hydraulic head varies inversely with specific yield (S y) and directly with recharge. Given small variation in water level spectra at the disturbed boundaries, the water level fluctuation in the aquifer is affected by the recharge condition and the aquifer spectral density function is sensitive to S y. Using an iterative technique to estimate ξ from 1400 sets of given parameters, 99% of the ξ′/ξ values deviated within only one order of magnitude with the model length (L) being equal to 1 km and 10 km. For L equal to 100 m, approximately 82% of the ξ′/ξ population falls within two orders of magnitude. Therefore, spectral analysis of aquifer hydraulic head response can be used to estimate the hydraulic diffusivity of an unconfined aquifer which is affected by periodic variations in recharge and head at boundaries.  相似文献   

5.
Ground water from springs and public supply wells was investigated for hydrochemistry and environmental isotopes of 3H, 18O and D in Jeju volcanic island, Korea. The wells are completed in a basaltic aquifer and the upper part of hydrovolcanic sedimentary formation. Nitrate contamination is conspicuous in the coastal area where most of the samples have nitrate concentrations well above 1 mg NO3 N/l. Agricultural land use seems to have a strong influence on the distribution of nitrate in ground water. Comparison of stable isotopic compositions of precipitation and ground water show that ground water mostly originates from rainy season precipitation without significant secondary modification and that local recharge is dominant. 3H concentration of ground water ranged from nearly zero to 5 TU and is poorly correlated with vertical location of well screens. The occurrence of the 3H‐free, old ground water is due to the presence of low permeability layers near the boundary of the basaltic aquifer and the hydrovolcanic sedimentary formation, which significantly limits ground water flow from the upper basaltic aquifer. The old ground water exhibited background‐level nitrate concentrations despite high nitrate loadings, whereas young ground water had considerably higher nitrate concentrations. This correlation of 3H and nitrate concentration may be ascribed to the history of fertilizer use that has increased dramatically since the early 1960s in the island. This suggests that 3H can be used as a qualitative indicator for aquifer vulnerability to nitrate contamination. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

6.
Abstract

The assessment of groundwater vulnerability to pollution has proved to be an effective tool for water resource management, especially in arid and semi-arid regions like Mahdia and Ksour Essaf. The main objective of this study is to assess the aquifer vulnerability by applying the DRASTIC method as well as using sensitivity analysis to evaluate the effect of each DRASTIC parameter on the final vulnerability map. An additional objective is to demonstrate the role of the GIS techniques in the vulnerability assessment. The DRASTIC method assigns a high vulnerability to the coast of the Mahdia-Ksour Essaf. The lowest values are observed in the southern part of the study area. A sensitivity analysis applied in this study suggests that net recharge, aquifer media and depth of groundwater are the key factors determining vulnerability. The model is validated with groundwater quality data and the results have shown strong relationships between modified DRASTIC Vulnerability Index and nitrate and chloride concentrations.

Citation Saidi, S., Bouri, S. & Ben Dhia, H. (2011) Sensitivity analysis in groundwater vulnerability assessment based on GIS in the Mahdia-Ksour Essaf aquifer, Tunisia: a validation study. Hydrol. Sci. J. 56(2), 288–304.  相似文献   

7.
Landscapes evolve in response to external forces, such as tectonics and climate, that influence surface processes of erosion and weathering. Internal feedbacks between erosion and weathering also play an integral role in regulating the landscapes response. Our understanding of these internal and external feedbacks is limited to a handful of field‐based studies, only a few of which have explicitly examined saprolite weathering. Here, we report rates of erosion and weathering in saprolite and soil to quantify how climate influences denudation, by focusing on an elevation transect in the western Sierra Nevada Mountains, California. We use an adapted mass balance approach and couple soil‐production rates from the cosmogenic radionuclide (CRN) 10Be with zirconium concentrations in rock, saprolite and soil. Our approach includes deep saprolite weathering and suggests that previous studies may have underestimated denudation rates across similar landscapes. Along the studied climate gradient, chemical weathering rates peak at middle elevations (1200–2000 m), averaging 112·3 ± 9·7 t km–2 y–1 compared to high and low elevation sites (46·8 ± 5·2 t km?2 y?1). Measured weathering rates follow similar patterns with climate as those of predicted silica fluxes, modeled using an Arrhenius temperature relationship and a linear relationship between flux and precipitation. Furthermore, chemical weathering and erosion are tightly correlated across our sites, and physical erosion rates increase with both saprolite weathering rates and intensity. Unexpectedly, saprolite and soil weathering intensities are inversely related, such that more weathered saprolites are overlain by weakly weathered soils. These data quantify exciting links between climate, weathering and erosion, and together suggest that climate controls chemical weathering via temperature and moisture control on chemical reaction rates. Our results also suggest that saprolite weathering reduces bedrock coherence, leading to faster rates of soil transport that, in turn, decrease material residence times in the soil column and limit soil weathering. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

8.
The ultrasonic profiling method of measuring the compressional and shear wave velocities in cylindrical rock samples is extended to measurements in some weathered and fresh granite blocks collected from the Hyderabad (India) region. This possibility of the method provides a means of investigating the elastic properties of the less compact rocks, of which the near-surface formations are particularly important. In this article the important parts of the ultrasonic profiling instrument developed are described and the relevant aspects of the seismic wave fields and identification of the individual waves in the wavetrain responses to longitudinal excitation are considered. Compressional, shear and surface (Rayleigh) wave velocities in some fresh and weathered granites are detailed. The compressional velocities range from 4.8 km/s to 5.5 km/s in fresh granites and lie between 1.1 km/s and 2.5 km/s in weathered granites. Young's modulus and Poisson's ratios computed from the measured velocities are also presented. An empirical relation of the form log E= 4.27 + 2.11 log Vp between Young's modulus E and compressional velocities Vp in the fresh granites studied is deduced. The versatility of the approach is thus demonstrated.  相似文献   

9.
Abstract

The dependence of the recession of the ground water levels and the ground water discharge upon the initial state of the aquifer is examined for deep unconfined aquifers. It is shown that only in the early stages of the recession does the initial state exert a limited influence on the recession. An estimate of the upper limit of the time t 0 for which for t > t 0 the recession becomes effectively independent of the initial state of the aquifer, valid for physically realistic initial states can be gained from inequalities (11) and (12a) and equation (16). t 0 depends essentially on the parameters of the aquifer and it is estimated that for useful aquifers t 0 can not be expected to exceed one month in relatively adverse cases. This explains why empirical recessions often are found to be consistent, of an exponential form.  相似文献   

10.
Heat as a tracer in fractured porous aquifers is more sensitive to fracture-matrix processes than a solute tracer. Temperature evolution as a function of time can be used to differentiate fracture and matrix characteristics. Experimental hot (50 °C) and cold (10 °C) water injections were performed in a weathered and fractured granite aquifer where the natural background temperature is 30 °C. The tailing of the hot and cold breakthrough curves, observed under different hydraulic conditions, was characterized in a log–log plot of time vs. normalized temperature difference, also converted to a residence time distribution (normalized). Dimensionless tail slopes close to 1.5 were observed for hot and cold breakthrough curves, compared to solute tracer tests showing slopes between 2 and 3. This stronger thermal diffusive behavior is explained by heat conduction. Using a process-based numerical model, the impact of heat conduction toward and from the porous rock matrix on groundwater heat transport was explored. Fracture aperture was adjusted depending on the actual hydraulic conditions. Water density and viscosity were considered temperature dependent. The model simulated the increase or reduction of the energy level in the fracture-matrix system and satisfactorily reproduced breakthrough curves tail slopes. This study shows the feasibility and utility of cold water tracer tests in hot fractured aquifers to boost and characterize the thermal matrix diffusion from the matrix toward the flowing groundwater in the fractures. This can be used as complementary information to solute tracer tests that are largely influenced by strong advection in the fractures.  相似文献   

11.
M. C. Tom Kuo 《Ground water》2022,60(4):510-517
Few published data are available for two-phase flow in fractures from field studies. All measurements of relative permeability reported in the literature were done in laboratory-scale. The in situ water saturations are normally not known for multiphase flow in natural fractures; therefore, the direct measurements of relative permeability are difficult in field-scale. With the help of a case study before and after the 2008 Mw 5.4 Antung earthquake, groundwater radon was used as a tracer to determine the gas and water saturations in a small naturally fractured aquifer. Well tests were also conducted to estimate aquifer transmissivity before and after the 2008 Antung earthquake. Anomalous declines in both groundwater radon concentration and transmissivity were observed precursory to the 2008 Antung earthquake. Both declines are two precursory phenomena having a common effect of gas bubbles. Using the data from well tests and radon tracer, one data point of water relative permeability can be obtained for in situ fractures. This data point reveals strong phase interference between water and gas bubbles for multiphase flow in natural fractures. Both the data of well tests and radon tracer are essential to gain an improved understanding of mass transfer behavior of groundwater-dissolved gases between water and gas phases.  相似文献   

12.
In this study, we used an archive of borehole logs from the British Geological Survey to collect information on the spatial structure of weathering that extends from the surface to competent bedrock across the Triassic Sherwood Sandstone Group outcrop (750 km2), in the East Midlands, UK. The borehole logs were used to estimate the thickness of the soil (n = 280) and soil and saprolite (S&S) to competent rock (n = 500). The weathering profile of the sandstone consisted of soil (median thickness ～ 1·5 m) overlying a transition zone of compacted and weakly cemented weathered sandstone saprolite over bedrock. Topographic analysis using a NEXTMAP 5 m × 5 m digital elevation model (DEM) revealed no significant relationships between slope properties (relief, flow length, flow accumulation or slope angle) and soil or S&S thickness. A weak, but statistically significant correlation was found between the thickness of the soil and S&S (rs = 0·25, p < 0·001, n = 192). The variation in soil thickness may be related to changes in current and historic and land‐use, variation in sandstone properties and the influence of glacial/peri‐glacial processes. The thickness of the saprolite was more variable towards the southern part of the study area, where it increased to a maximum 40 m. We hypothesize and provide evidence that the greater weathering thickness is related to the occurrence of increased faulting in this part of the study region, allowing increased access to meteoric waters. A possible source of increased water supply is meltwater from Quaternary ice sheets; the overburden of ice may have increased sub‐glacial pore water pressure, with the fractures and faults acting as a drainage system for the removal of dissolved weathering products. British Geological Survey © NERC 2010  相似文献   

13.
The position of the auroral luminosity equatorward boundary during the interaction between the Earth’s magnetosphere and isolated solar wind streams from different solar sources has been statistically studied based on the ground and satellite observations of auroras. These studies continue the series of the works performed in order to develop the technique for predicting auroras based on the characteristics of the interplanetary medium and auroral disturbances. The dependences of the minimal position of the auroral luminosity equatorward boundary (Φ′) on the values of the azimuthal component of the interplanetary electric field (E y ) and AL indices of magnetic activity, averaged over 6 and 24 h, are presented. The distribution limits for each type of isolated solar wind streams on the Φ′-E y and Φ′-AL planes have been determined.  相似文献   

14.
Traditional aquifer vulnerability techniques primarily rely on spatial property data for a region and are limited by their ability to directly or indirectly assess flow and transport processes occurring from the surface to depth within an aquifer system. The main objective of this study was to investigate groundwater vulnerability in terms of aquifer interconnectivity and flow dynamics. A combination of stable isotopes, groundwater age‐dating (radiocarbon), and geomorphic/geogenic spatial analyses was applied to a regional, highly developed coastal aquifer to explain the presence of nitrate at depth. The average δ13C value (?17.3 ± 2‰ VPDB, n = 27) is characteristic of groundwater originating from locally infiltrated precipitation through extensively cultivated soils. The average δ18O and δD values (?4.0 ± 0.1‰ VSMOW, n = 27; δD: ?19.3 ± 1‰ VSMOW, n = 27, respectively) are similar to precipitation water derived from maritime sources feeding the region's surface water and groundwater. Stable and radioactive isotopes reveal significant mixing between shallow and deep aquifers due to high velocities, hydraulic connection, and input of local recharge water to depths. Groundwater overdevelopment enhances deeper and faster modern water downward flux, amplifying aquifer vulnerability. Therefore, aquifer vulnerability is a variable, dependent on the type and degree of stress conditions experienced by a groundwater system as well as the geospatial properties at the near surface.  相似文献   

15.
To enhance the understanding of solute dynamics within the stream‐to‐riparian continuum during flood event‐driven water fluctuation (i.e., flood wave), a variable saturated groundwater flow and solute transport model were developed and calibrated against in situ measurements of the Inbuk stream, Korea, where seasonal flooding prevails. The solute dynamics were further investigated for flood waves (varying by amplitude [A], duration [T], roundness [r], and skewness [tp]) that were parameterised by real‐time stream stage fluctuations. We found that the solute transferred faster and farther in the riparian zone, especially within the phreatic zone, above which in the variable saturated zone the concentration required a significantly longer time, particularly at higher altitudes, to return to the initial state. By comparison, solute transferred shallowly in the streambed where the solute plume exhibited an exponential growth trend from the centre to the bank. The dynamic changes of solute flux and mass along the stream–aquifer interface and stream concentration were linked to the shape of flood wave. As the flood wave became higher (A↗), wider (T↗), rounder (r↘), and less skewed (tp↗), the maximum solute storage in aquifer increased. Maximum stream concentration (Cstr?max) not only presented a positive linear relationship with A or tp but also showed a negative logarithmic trend with increasing T or r. The sensitivity of Cstr_max to A was approximately two times that of tp, and between these values, the r was slightly more sensitive than T. Cstr?max linearly increased as hydraulic conductivity increased and logarithmically increased as longitudinal dispersivity increased. The former relationship was more sensitive than the latter.  相似文献   

16.
Changes in Entrapped Gas Content and Hydraulic Conductivity with Pressure   总被引：1，自引：0，他引：1  
Water table fluctuations continuously introduce entrapped air bubbles into the otherwise saturated capillary fringe and groundwater zone, which reduces the effective (quasi‐saturated) hydraulic conductivity, Kquasi, thus impacting groundwater flow, aquifer recharge and solute and contaminant transport. These entrapped gases will be susceptible to compression or expansion with changes in water pressure, as would be expected with water table (and barometric pressure) fluctuations. Here we undertake laboratory experiments using sand‐packed columns to quantify the effect of water table changes of up to 250 cm on the entrapped gas content and the quasi‐saturated hydraulic conductivity, and discuss our ability to account for these mechanisms in ground water models. Initial entrapped air contents ranged between 0.080 and 0.158, with a corresponding Kquasi ranging between 2 and 6 times lower compared to the Ks value. The application of 250 cm of water pressure caused an 18% to 26% reduction in the entrapped air content, resulting in an increase in Kquasi by 1.16 to 1.57 times compared to its initial (0 cm water pressure) value. The change in entrapped air content measured at pressure step intervals of 50 cm, was essentially linear, and could be modeled according to the ideal gas law. Meanwhile, the changes in Kquasi with compression–expansion of the bubbles because of pressure changes could be adequately captured with several current hydraulic conductivity models.  相似文献   

17.
The conversion of bedrock to regolith marks the inception of critical zone processes, but the factors that regulate it remain poorly understood. Although the thickness and degree of weathering of regolith are widely thought to be important regulators of the development of regolith and its water‐storage potential, the functional relationships between regolith properties and the processes that generate it remain poorly documented. This is due in part to the fact that regolith is difficult to characterize by direct observations over the broad scales needed for process‐based understanding of the critical zone. Here we use seismic refraction and resistivity imaging techniques to estimate variations in regolith thickness and porosity across a forested slope and swampy meadow in the Southern Sierra Critical Zone Observatory (SSCZO). Inferred seismic velocities and electrical resistivities image a weathering zone ranging in thickness from 10 to 35 m (average = 23 m) along one intensively studied transect. The inferred weathering zone consists of roughly equal thicknesses of saprolite (P‐velocity < 2 km s?1) and moderately weathered bedrock (P‐velocity = 2–4 km s?1). A minimum‐porosity model assuming dry pore space shows porosities as high as 50% near the surface, decreasing to near zero at the base of weathered rock. Physical properties of saprolite samples from hand augering and push cores are consistent with our rock physics model when variations in pore saturation are taken into account. Our results indicate that saprolite is a crucial reservoir of water, potentially storing an average of 3 m3 m?2 of water along a forested slope in the headwaters of the SSCZO. When coupled with published erosion rates from cosmogenic nuclides, our geophysical estimates of weathering zone thickness imply regolith residence times on the order of 105 years. Thus, soils at the surface today may integrate weathering over glacial–interglacial fluctuations in climate. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

18.
A field tracer test was carried out in a light nonaqueous phase liquid (LNAPL) source zone using a well pattern consisting of one injection well surrounded by four extraction wells (5‐spot well pattern). Multilevel sampling was carried out in two observation wells located inside the test cell characterized by heterogeneous lithology. Tracer breakthrough curves showed relatively uniform flow within soil layers. A numerical flow and solute transport model was calibrated on hydraulic heads and tracer breakthrough curves. The model was used to estimate an average accessible porosity of 0.115 for the swept zone and an average longitudinal dispersivity of 0.55 m. The model was further used to optimize the relative effects of viscous forces versus capillary forces under realistic imposed hydraulic gradients and to establish optimal surfactant solution properties. Maximum capillary number (NCa) values between injection and extraction wells were obtained for an injection flow rate of 16 L/min, a total extraction flow rate of 20 L/min, and a surfactant solution with a viscosity of 0.005 Pa?s. The unconfined nature of the aquifer limited further flow rate or viscosity increases that would have led to unrealistic hydraulic gradients. An NCa range of 3.8 × 10?4 to 7.6 × 10?3 was obtained depending on the magnitude of the simulated LNAPL‐water interfacial tension reduction. Finally, surfactant and chase water slug sizing was optimized with a radial form of the simplified Ogata‐Banks analytical solution (Ogata and Banks 1961) so that injected concentrations could be maintained in the entire 5‐spot cell.  相似文献   

19.
Estimation of specific yield (Sy) of an aquifer is of great importance in water resource management. In this study, an experimental drainage method was developed to determine specific yield of an alluvial aquifer of the Platte River valley, Nebraska, USA. Sets of sediment cores with continuous interval depth were collected to plastic tubes using the direct push technique and then taken to the laboratory. During the Sy experiment, those sediment cores were re-saturated by placing them vertically in a large water tank. Sy was determined by the water drained from the sediments by force of gravity in a bracket. Our results show that the values of Sy varied largely with depth at each site and the variability for Sy with interval depth between the test sites is also observed. This spatial heterogeneity in Sy might result from the variation of grain size, grain shape, sorting and compaction of sediments in different cores with interval depth. The Sy for all sediment cores ranged from 0.01 to 0.18 and the mean value was 0.08±0.04. Our drainage method can functionally preserve the sedimentary structures in their original state and it is easier to experiment at a lower expense.  相似文献   

20.
  相似文献   

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

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