首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 922 毫秒
1.
A new in‐situ remediation concept termed a Horizontal Reactive Media Treatment Well (HRX Well®) is presented that utilizes horizontal wells filled with reactive media to passively treat contaminated groundwater in‐situ. The approach involves the use of large‐diameter directionally drilled horizontal wells filled with granular reactive media generally installed parallel to the direction of groundwater flow. The design leverages natural “flow‐focusing” behavior induced by the high in‐well hydraulic conductivity of the reactive media relative to the aquifer hydraulic conductivity to passively capture and treat proportionally large volumes of groundwater within the well. Clean groundwater then exits the horizontal well along its downgradient sections. Many different types of solid granular reactive media are already available (e.g., zero valent iron, activated carbon, ion exchange resins, zeolite, apatite, chitin); therefore, this concept could be used to address a wide range of contaminants. Three‐dimensional flow and transport simulations were completed to assess the general hydraulic performance, capture zones, residence times, effects of aquifer heterogeneity, and treatment effectiveness of the concept. The results demonstrate that capture and treatment widths of up to tens of feet can be achieved for many aquifer settings, and that reductions in downgradient concentrations and contaminant mass flux are nearly immediate. For a representative example, the predicted treatment zone width for the HRX Well is approximately 27 to 44 feet, and contaminant concentrations immediately downgradient of the HRX Well decreased an order of magnitude within 10 days. A series of laboratory‐scale physical tests (i.e., tank tests) were completed that further demonstrate the concept and confirm model prediction performance. For example, the breakthrough time, peak concentration and total mass recovery of methylene blue (reactive tracer) was about 2, 35, and 20 times (respectively) less than chloride (conservative tracer) at the outlet of the tank‐scale HRX Well.  相似文献   

2.
Groundwater interacts with surface water features nearly in all types of landscapes. Understanding these interactions has practical consequences on the quantity and quality of water in either system, because the depletion or contamination of one of the systems will eventually affect the other one. Many studies have shown that the use of heat as natural tracer in conjunction with water level measurements is an effective method for estimating water flow (fluxes) between groundwater and surface water. A number of studies have explored the effects of spatial and temporal variability of groundwater–surface water flux exchanges using temperature and water level measurements; however, the effect of temporal resolution of water level and temperature data on estimating flux remains unexplored. Therefore, this study investigated the effect of temporal resolution of input data on temporal variation of groundwater–surface water flux exchanges. To this end, we calibrated a variably saturated two‐dimensional groundwater flow and heat transport model (VS2DH) at hourly and daily time scales using temperatures measured at multiple depths below the riverbed of the Zenne River, located at a well‐known Belgian brownfield site. Results of the study showed that the computed water flux through the streambed ranged between ?32 mm/day and +25 mm/day using the hourly model and from ?10 mm/day to ?37 mm/day using the daily model. The hourly model resulted in detecting reversal of flow direction inducing short‐term surface water flow into the streambed. However, such events were not captured if daily temperature and water level measurements were used as input. These findings have important implications for understanding contaminant mass flux and their attenuation in the mixing zone of groundwater and surface water. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

3.
《Advances in water resources》2005,28(10):1102-1111
We present a range of analytical solutions to the combined transient water and solute transport for horizontal flow. We adopt the concept of a scale and time dependent dispersivity used for contaminant transport in aquifers and apply it to transient, unsaturated horizontal flow to develop similarity solutions for both constant solute concentration and solute flux boundary conditions. Through the use of a specific form of the water profile as used by Brutsaert [Water Resour Res 1968:4;785], the solute profiles can be reduced to a simple quadrature. We also derive a solution for the instantaneous injection of water and solute into a horizontal media for an arbitrary dispersivity. It is found that the solute concentration remains constant in both space and time as the water redistributes, suggesting that the solute does not disperse relative to the water.  相似文献   

4.
Measurement and interpretation of mass fluxes in favor of concentrations is gaining more and more interest, especially within the framework of the characterization and management of large-scale volatile organic carbon (VOC) groundwater contamination (source zones and plumes). Traditional methods of estimating contaminant fluxes and discharges involve individual measurements/calculations of the Darcy water flux and the contaminant concentrations. However, taken into account the spatially and temporally varying hydrologic conditions in complex, heterogeneous aquifers, higher uncertainty arises from such indirect estimation of contaminant fluxes. Therefore, the potential use of passive sampling devices for the direct measurement of groundwater-related VOC mass fluxes is examined. A review of current passive samplers for the measurement of organic contaminants in water yielded the selection of 18 samplers that were screened for a number of criteria. These criteria are related to the possible application of the sampler for the measurement of VOC mass fluxes in groundwater. This screening study indicates that direct measurement of VOC mass fluxes in groundwater is possible with very few passive samplers. Currently, the passive flux meter (PFM) is the only passive sampler which has proven to effectively measure mass fluxes in near source groundwater. A passive sampler for mass flux measurement in plume zones with regard to long-term monitoring (several months to a year) still needs to be developed or optimized. A passive sampler for long-term monitoring of contaminant mass fluxes in groundwater would be of considerable value in the development of risk-based assessment and management of soil and groundwater pollutions.  相似文献   

5.
Erick Carlier 《水文研究》2008,22(17):3500-3506
An analytical transport‐model was developed to simulate the propagation of a contaminant in one‐ and two‐dimensional transient flow in groundwater. It is proved that the distribution of concentration at a given time and for a given discharge is identical to that obtained for a different discharge if the volumetric flux of water is the same in the two cases. The results of simulations have been compared with results obtained using the MT3DMS numerical model. There is good agreement when the calculated concentrations are flux‐weighted concentrations. On the other hand, there is a notable divergence when the resident mode is considered. Resident mode concentrations express the mass per unit volume whereas flux mode concentrations express the ratio of mass flux to fluid flux. The solutions presented in this paper can thus be a useful alternative to code MT3DMS when the objective is to simulate concentrations in transient flow according to a resident mode. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

6.
The internal tracer method for estimating contaminant degradation rates separates the attenuation effects not associated with degradation by using a codisposed recalcitrant internal tracer to normalize the degrading contaminant concentration. The remaining attenuation between the internal tracer and degrading contaminant is attributed to degradation and the degradation rate half-life is estimated from the first-order decay equation. An analytical solution of the advection- dispersion equation was used to evaluate flow-and-transport conditions that could result in incorrect estimates of contaminant degradation rate constants. Flow-and-transport characteristics that result in overestimating degradation rates were of particular interest because the internal tracer method often used to demonstrate natural attenuation can achieve remedial objectives. The analytical solution was also used to estimate the magnitude of error associated with using the internal tracer method at an example site and to explain different degradation rates estimated using tracers with different decay rate constants.  相似文献   

7.
Field tracer tests indicate that advective porosity, the quantity relating advective velocity to Darcy flux, may exhibit directional dependence. Hydraulic anisotropy explains some but not all of the reported directional results. The present paper shows mathematically that directional variations in advective porosity may arise simply from incomplete mixing of an inert tracer between directional flow channels within a sampling (or support) volume ω of soil or rock that may be hydraulically isotropic or anisotropic. In the traditional fully homogenized case, our theory yields trivially a scalar advective porosity equal to the interconnected porosity ϕ, thus explaining neither the observed directional effects nor the widely reported experimental finding that advective porosity is generally smaller than ϕ. We consider incomplete mixing under conditions in which the characteristic time tD of longitudinal diffusion along channels across ω is much shorter than the characteristic time tH required for homogenization through transverse diffusion between channels. This may happen where flow takes place preferentially through relatively conductive channels and/or fractures of variable orientation separated by material that forms a partial barrier to diffusive transport. Our solution is valid for arbitrary channel Peclet numbers on a correspondingly wide range of time scales tD ⩽ t ≪ tH. It shows that the tracer center of mass is advected at a macroscopic velocity which is generally not collinear with the macroscopic Darcy flux and exceeds it in magnitude. These two vectors are related through a second-rank symmetric advective dispersivity tensor Φ. If the permeability k of ω is a symmetric positive-definite tensor, so is Φ. However, the principal directions and values of these two tensors are generally not the same; whereas those of k are a fixed property of the medium and the length-scale of ω, those of Φ depend additionally on the direction and magnitude of the applied hydraulic gradient. When the latter is large, diffusion has negligible effect on Φ and one may consider tracer mass to be distributed between channels in proportion to the magnitude of their Darcy flux. This is made intuitive through a simple example of an idealized fracture network. Our analytical formalism reveals the properties of Φ but is too schematic to allow predicting the latter accurately on the basis of realistic details about the void structure of ω and tracer mass distribution within it. Yet knowing the tensorial properties of Φ is sufficient to allow determining it indirectly on the basis of ω-scale hydraulic and tracer data, including concentrations that represent homogenized samples extracted from (or sensed externally across) an ω-scale plume.  相似文献   

8.
Part of a small drainage basin on the Sevilleta National Wildlife Refuge (about 25 km north of Socorro, NM) was intensively instrumented with soil monitoring equipment to estimate natural ground-water recharge. Soil-moisture data were analysed with special attention to characterizing the influence of topography on the direction of vadose water flow paths in fine to medium aeolian sand. Moisture content data were obtained by the neutron scattering technique, and hydraulic head data were obtained using tensiometers. In addition, tracer experiments were conducted on a sandy hillslope to delineate the flow paths of vadose water. The results indicate that there is a strong lateral component to unsaturated flow on a hillslope, even in the absence of apparent sublayers of much lower permeability. Darcian calculations estimate the long-term, steady, deep flux beneath a concave location to be about 4 per cent of an assumed mean annual precipitation of 20 cm. The deep soil water flux downward varied by several orders of magnitude during the 17 month period of record.  相似文献   

9.
多环芳烃在城市湖泊气-水界面上的交换   总被引：4，自引：1，他引：3  
李军  张干  祁士华  刘国卿 《湖泊科学》2004,16(3):238-244
利用大流量采样器和半渗透膜装置对广州麓湖大气中气态多环芳烃和水体中溶解态多环芳烃进行了连续一年的观测,在此数据基础上,依据国内外通用的双膜理论计算了多环芳烃在广州麓湖水面上的交换通量.结果显示萘、苊、二氢苊的交换通量方向是由湖水挥发进入大气,其它主要化合物都是从大气进入水体;在单个化合物的交换通量中,以菲的通量值为最大,平均达13137.1ng/(m2·d),最高达到21211.22ng/(m2·d),其次为萘、荧蒽、蒽、芘、芴、二氢苊和苊等.各化合物交换通量的季节变化特点大部分都是在夏季达到最大值,在冬季交换通量降到最小值,但唯一例外的是芴,其交换通量却是在夏季降到最低值.  相似文献   

10.
Due to differences in hydraulic conductivity and effects of well construction geometry, groundwater lateral flow through a monitoring well typically differs from groundwater flow in the surrounding aquifer. These differences must be well understood in order to apply passive measuring techniques, such as passive flux meters (PFMs) used for the measurement of groundwater and contaminant mass fluxes. To understand these differences, lab flow tank experiments were performed to evaluate the influences of the well screen, the surrounding filter pack and the presence of a PFM on the natural groundwater flux through a monitoring well. The results were compared with analytical calculations of flow field distortion based on the potential theory of Drost et al. (1968). Measured well flow field distortion factors were found to be lower than calculated flow field distortion factors, while measured PFM flow field distortion factors were comparable to the calculated ones. However, this latter is not the case for all conditions. The slotted geometry of the well screen seems to make a correct analytical calculation challenging for conditions where flow field deviation occurs, because the potential theory assumes a uniform flow field. Finally, plots of the functional relationships of the distortion of the flow field with the hydraulic conductivities of the filter screen, surrounding filter pack and corresponding radii make it possible to design well construction to optimally function during PFM applications.  相似文献   

11.
Innovative remediation studies were conducted between 1994 and 2004 at sites contaminated by nonaqueous phase liquids (NAPLs) at Hill and Dover AFB, and included technologies that mobilize, solubilize, and volatilize NAPL: air sparging (AS), surfactant flushing, cosolvent flooding, and flushing with a complexing-sugar solution. The experiments proved that aggressive remedial efforts tailored to the contaminant can remove more than 90% of the NAPL-phase contaminant mass. Site-characterization methods were tested as part of these field efforts, including partitioning tracer tests, biotracer tests, and mass-flux measurements. A significant reduction in the groundwater contaminant mass flux was achieved despite incomplete removal of the source. The effectiveness of soil, groundwater, and tracer based characterization methods may be site and technology specific. Employing multiple methods can improve characterization. The studies elucidated the importance of small-scale heterogeneities on remediation effectiveness, and fomented research on enhanced-delivery methods. Most contaminant removal occurs in hydraulically accessible zones, and complete removal is limited by contaminant mass stored in inaccessible zones. These studies illustrated the importance of understanding the fluid dynamics and interfacial behavior of injected fluids on remediation design and implementation. The importance of understanding the dynamics of NAPL-mixture dissolution and removal was highlighted. The results from these studies helped researchers better understand what processes and scales are most important to include in mathematical models used for design and data analysis. Finally, the work at these sites emphasized the importance and feasibility of recycling and reusing chemical agents, and enabled the implementation and success of follow-on full-scale efforts.  相似文献   

12.
For the purposes of risk mitigation and remediation design, Darcy's Law calculations, that are typically used to estimate flow direction and magnitude, can be usefully supplemented with more direct measurements. The point velocity probe (PVP), which measures groundwater velocity by conducting a mini‐tracer test across a cylindrical probe surface, was developed to address this need. Laboratory and field experiments have shown that these probes can provide accurate velocity magnitude and direction measurements at the centimeter scale. In an effort to streamline the production of PVPs, three‐dimensional (3D) printing was investigated. The 3D printer produces a plastic probe body with designated detector and injection port locations. The subsequent installation of injection lines and detector wires requires, for a single‐port probe, no longer than 1–2 h of assembly time, compared to approximately 10 h to assemble previous designs (mainly determined by injection port installation). Probes can be printed in batches with the number of units made depending on the printer and probe size. Laboratory tests of the original PVP models yielded velocity magnitudes and directions with average errors (from expected values) of ±9% to ±15% and ±8°, respectively. A printed PVP was tested in a nested storage tank system (NeST) packed with sand that mimics flow conditions through a sandy aquifer. The probe was tested at multiple pumping rates with the injector port at angles of (α =) 30°, 45°, and 75° from the expected linear flow direction. The printed PVPs provided an average magnitude percent error of ±13.5% and an average direction error of ±4°. This shows that in laboratory tests, printed PVPs performed as well as the original probes previously reported, making them viable units for field applications as well as being quicker to assemble than the earlier designs.  相似文献   

13.
The measured drainage fluxes through a layered volcanic vadose zone exhibited high spatial variability as a consequence of heterogeneous flow conditions. The drainage flux variability was quantified using automated equilibrium tension lysimeters, installed in close‐proximity and resulted in high variability in the Br masses recovered from a conservative tracer experiment. The primary cause of the heterogeneous flow was attributed to textural changes occurring at the interface between volcanic layers, resulting in development of funnel‐flow patterns, and further enhanced by the existence of hydrophobic conditions. The Br recoveries in individual automated equilibrium tension lysimeters were used to determine the corresponding variable sizes of the surface areas contributing drainage to the lysimeters. The tracer experiment confirmed the existence of unsaturated lateral transport occurring at the sloping interface of the coarse Taupo Ignimbrite material with the silty Palaeosol layer at approximately 4.2 m depth. This study demonstrates that measurements of both flux and solute concentrations at multiple locations are essential when heterogeneous flow is suspected to be present, to be able to determine reliable estimates of contaminant leaching through the vadose zone at the plot scale. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
This article provides a partial answer to the question “What is the relation between excess hydraulic head and volume flux of water in a conduit within a porous matrix?”, focusing on the case that the forcing is steady. The conduit is modelled as a horizontal circular cylinder, imbedded within a porous matrix of rectangular cross section, having constant head prescribed on the sidewalls and being confined top and bottom. Laminar flow in the matrix is assumed to obey Darcy's law, while turbulent flow in the conduit is quantified using the Darcy–Weisbach equation. Analysis of the latter equation shows that the length scale of variations in the direction of the conduit is large compared with the scale of lateral and vertical variations. This permits separation of the full three-dimensional non-linear problem into a two-dimensional linear problem for head within the matrix and a one-dimensional non-linear problem for head within the conduit. Analytic solutions are obtained for the distribution of head in the matrix and in a conduit of either infinite or finite length. In both cases, the volume flux of water is proportional to the excess head to the 2/3 power, the conduit radius to the 5/3 power, the matrix permeability to the 1/3 power and gravity to the 1/3 power. The scale of variation of head along the conduit is proportional to the excess head to the ?1/3 power, the conduit radius to the 5/3 power, the matrix permeability to the ?2/3 power and gravity to the 1/3 power.  相似文献   

15.
Resource extraction and transportation activities in subarctic Canada can result in the unintentional release of contaminants into the surrounding peatlands. In the event of a release, a thorough understanding of solute transport within the saturated zone is necessary to predict plume fate and the potential impacts on peatland ecosystems. To better characterize contaminant transport in these systems, approximately 13,000 L/day of sodium chloride tracer (200 mg/L) was released into a bog in the James Bay Lowland. The tracer was pumped into a fully penetrating well (1.5 m) between July 5 and August 18, 2015. Horizontal and vertical plume development was measured via in situ specific conductance and water table depth from an adaptive monitoring network. Over the spill period, the bulk of the plume travelled a lateral distance of 100 m in the direction of the slight regional groundwater and topographical slope. The plume shape was irregular and followed the hollows, indicating preferential flow paths due to the site microtopography. Saturated transport of the tracer occurred primarily at ~25 cm below ground surface (bgs), and at a discontinuous high hydraulic conductivity layer ~125 cm bgs due to a complex and heterogeneous vertical hydraulic conductivity profile. Plume measurement was confounded by a large amount of precipitation (233 mm over the study period) that temporarily diluted the tracer in the highly conductive upper peat layer. Longitudinal solute advection can be approximated using local water table information (i.e., depth and gradient); microtopography; and meteorological conditions. Vertical distribution of solute within the peat profile is far more complex due to the heterogeneous subsurface; characterization would be aided by a detailed understanding of the site‐specific peat profile; the degree of decomposition; and the type of contaminant (e.g., reactive/nonreactive). The results of this research highlight the difficulty of tracking a contaminant spill in bogs and provide a benchmark for the characterization of the short‐term fate of a plume in these complex systems.  相似文献   

16.
Measurement uncertainty is a key hindrance to the quantification of water fluxes at all scales of investigation. Predictions of soil‐water flux rely on accurate or representative measurements of hydraulic gradients and field‐state hydraulic conductivity. We quantified the potential magnitude of errors associated with the parameters and variables used directly and indirectly within the Darcy – Buckingham soil‐water‐flux equation. These potential errors were applied to a field hydrometric data set collected from a forested hillslope in central Singapore, and their effect on flow pathway predictions was assessed. Potential errors in the hydraulic gradient calculations were small, approximately one order of magnitude less than the absolute magnitude of the hydraulic gradients. However, errors associated with field‐state hydraulic conductivity derivation were very large. Borehole (Guelph permeameter) and core‐based (Talsma ring permeameter) techniques were used to measure field‐saturated hydraulic conductivity. Measurements using these two approaches differed by up to 3\9 orders of magnitude, with the difference becoming increasingly marked within the B horizon. The sensitivity of the shape of the predicted unsaturated hydraulic conductivity curve to ±5% moisture content error on the moisture release curve was also assessed. Applied moisture release curve error resulted in hydraulic conductivity predictions of less than ±0\2 orders of magnitude deviation from the apparent conductivity. The flow pathways derived from the borehole saturated hydraulic conductivity approach suggested a dominant near‐surface flow pathway, whereas pathways calculated from the core‐based measurements indicated vertical percolation to depth. Direct tracer evidence supported the latter flow pathway, although tracer velocities were approximately two orders of magnitude smaller than the Darcy predictions. We conclude that saturated hydraulic conductivity is the critical hillslope hydrological parameter, and there is an urgent need to address the issues regarding its measurement further. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

17.
Understanding runoff generation processes is important for flood prediction, water management, erosion control, water quality, contaminant transport and the evaluation of impacts of land use change. However, little process research has been carried out in southern Chile. In particular the young volcanic ash soils, which are typical for this area, are not well understood in their hydrologic behaviour. To establish a ‘reference study’ which can then be used for comparison with other (disturbed) sites, this study focuses on the investigation of runoff generation processes in an undisturbed, forested catchment in the Chilean Andes. The paper reports on an investigation of these processes with different tracer methods at different spatial scales. Hydrograph separation with environmental isotopes and geochemical constituents was used on the catchment scale. Thermal energy was used as a tracer to investigate groundwater–surface water interactions at the local stream reach scale and dye tracers were used to study infiltration and percolation characteristics at the plot scale. It was found that pre‐event water dominates the storm hydrograph. In the lower reaches, however, water usually exfiltrates from the stream into the adjacent aquifer. The dye tracer experiments showed that while preferential vertical flow dominates under forest, water infiltrates as a straight horizontal front in the bare volcanic ashes (no vegetation) on the catchment rim. Subsurface flow patterns in the forest differ significantly from summer to winter. All three approaches used in this study suggest an important shift in dominant processes from dry to wet season. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

18.
Large rivers have been previously shown to be vertically heterogeneous in terms of suspended particulate matter (SPM) concentration, as a result of sorting of suspended solids. Therefore, the spatial distribution of suspended sediments within the river section has to be known to assess the riverine sedimentary flux. Numerous studies have focused on the vertical distribution of SPM in a river channel from a theoretical or experimental perspective, but only a few were conducted so far on very large rivers. Moreover, a technique for the prediction of depth‐integrated suspended sediment fluxes in very large rivers based on sediment transport dynamics has not yet been proposed. We sampled river water along depth following several vertical profiles, at four locations on the Amazon River and its main tributaries and at two distinct water stages. Depending on the vertical profile, a one‐ to fivefold increase in SPM concentration is observed from river channel surface to bottom, which has a significant impact on the ‘depth‐averaged’ SPM concentration. For each cross section, a so‐called Rouse profile quantitatively accounts for the trend of SPM concentration increase with depth, and a representative Rouse number can be measured for each cross section. However, the prediction of this Rouse number would require the knowledge of the settling velocity of particles, which is dependent on the state of aggregation affecting particles within the river. We demonstrate that in the Amazon River, particle aggregation significantly influences the Rouse number and renders its determination impossible from grain‐size distribution data obtained in the lab. However, in each cross section, the Rouse profile obtained from the fit of the data can serve as a basis to model, at first order, the SPM concentration at any position in the river cross section. This approach, combined with acoustic Doppler current profiler (ADCP) water velocity transects, allows us to accurately estimate the depth‐integrated instantaneous sediment flux. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

19.
Field tracer experiments and model calibrations indicate that the magnitude of dispersivity increases as a function of the scale at which observations are made. Calculations presented in this study suggest that some part of this scaling may be explained as an artifact of the models used. Specifically, a scaling-up of dispersivity will occur whenever an (n − 1)-dimensional model is calibrated or otherwise employed to describe an n-dimensional system. The calibrated coefficients for such models will depend not only on size of the contaminant plume or tracer experiment at the time of calibration, but will exhibit a size dependency beyond the calibration period. The magnitude of scaling appears to be sufficient to encompass the range of differences between laboratory measurements of dispersivity and model calibrations.  相似文献   

20.
Large-scale advective transport through highly heterogeneous 3D formations is investigated using highly resolved numerical simulations and simple analytic models. Investigations are focused on impacts of two types of contaminant injection on transport through isotropic formations where flow conditions are uniform in the average. Transport is quantified by analyzing breakthrough curves for control planes at various distances from the injection zone. In flux-proportional injection mode local mass in injection zone is proportional to local groundwater flux; this setup models many practical cases such as contaminant injection through wells. In resident concentration mode local concentration in injection zone is constant. Results show that impacts of injection mode on breakthrough curves and their moments are strong and they persist for hundreds of correlation scales. The resident concentration mode leads to a fatter tails of the breakthrough curves, while the peaks are generally underpredicted. For a synthetic porous medium with logconductivity variance of 8, dispersivity computed using resident concentration mode at control plane 100 integral scales away from the injection zone was about 10 times larger than corresponding one for flux-proportional mode. Hence, injection mode impacts on transport through highly heterogeneous formations are strong and they persist for large distances from the injection zone.  相似文献   

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

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