首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
In this research, a light transmission visualization (LTV) method was used to quantify dense non-aqueous phase liquids (DNAPL) saturation in two-dimensional (2-D), two fluid phase systems. The method is an expansion of earlier LTV methods and takes into account both absorption and refraction light theories. Based on this method, DNAPL and water saturations can rapidly be obtained point wise across sand-packed 2-D flow chambers without the need to develop a calibration curve. A single point calibration step is, however, needed when dyed DNAPL is used to account for the change in the transmission factor at the dyed DNAPL–water interface. The method was applied to measure, for the first time, undyed DNAPL saturation in small 2-D chambers. Known amounts of DNAPL, modeled by tetrachloroethylene (PCE), were added to the chamber and these amounts were compared to results obtained by this LTV method. Strong correlation existed between results obtained based on this method and the known PCE amounts with an R2 value of 0.993. Similar experiments conducted using dyed PCE showed a stronger correlation between results obtained by this LTV method and the known amounts of dyed PCE added to the chamber with an R2 value of 0.999. The method was also used to measure dyed PCE saturation in a large 2-D model following sparging experiments. Results obtained from image analyses following each sparging event were compared to results obtained by two independent techniques, namely gas chromatography–mass spectrometry (GC/MS) analyses and carbon column extraction. There was a good agreement between the results obtained by this LTV method and those obtained by the two independent techniques when experiments were carried out under stable light source conditions and errors in mass balance were minor. The method presented here can be expanded to measure fluid contents in three fluid phase systems and provide a non-destructive, non-intrusive tool to investigate changes in DNAPL architecture and flow characteristics in laboratory experiments.  相似文献   

2.
Field investigations were carried out to determine the occurrence of tetrachloroethene (PCE) dense nonaqueous phase liquid (DNAPL), the source zone architecture and the aquitard integrity at a 30‐ to 50‐year old DNAPL release site. The DNAPL source zone is located in the clay till unit overlying a limestone aquifer. The DNAPL source zone architecture was investigated through a multiple‐lines‐of‐evidence approach using various characterization tools; the most favorable combination of tools for the DNAPL characterization was geophysical investigations, membrane interface probe, core subsampling with quantification of chlorinated solvents, hydrophobic dye test with Sudan IV, and Flexible Liner Underground Technologies (FLUTe) NAPL liners with activated carbon felt (FACT). While the occurrence of DNAPL was best determined by quantification of chlorinated solvents in soil samples supported by the hydrophobic dye tests (Sudan IV and NAPL FLUTe), the conceptual understanding of source zone architecture was greatly assisted by the indirect continuous characterization tools. Although mobile or high residual DNAPL (S t > 1%) only occurred in 11% of the source zone samples (intact cores), they comprised 86% of the total PCE mass. The dataset, and associated data analysis, supported vertical migration of DNAPL through fractures in the upper part of the clay till, horizontal migration along high permeability features around the redox boundary in the clay till, and to some extent vertical migration through the fractures in the reduced part of the clay till aquitard to the underlying limestone aquifer. The aquitard integrity to DNAPL migration was found to be compromised at a thickness of reduced clay till of less than 2 m.  相似文献   

3.
Porous aquifer materials are often characterized by layered heterogeneities that influence groundwater flow and present complexities in contaminant transport modeling. Such flow variations also have the potential to impact the dissolution flux from dense nonaqueous phase liquid (DNAPL) pools. This study examined how these heterogeneous flow conditions affected the dissolution of a tetrachloroethene (PCE) pool in a two-dimensional intermediate-scale flow cell containing coarse sand. A steady-state mass-balance approach was used to calculate the PCE dissolution rate at three different flow rates. As expected, aqueous PCE concentrations increased along the length of the PCE pool and higher flow rates decreased the aqueous PCE concentration in the effluent. Nonreactive tracer studies at two flow rates confirmed the presence of a vertical flow gradient, with the most rapid velocity located at the bottom of the tank. These results suggest that flow focusing occurred near the DNAPL pool. Effluent PCE concentrations and pool dissolution flux rates were compared to model predictions assuming local equilibrium (LE) conditions at the DNAPL pool/aqueous phase interface and a uniform distribution of flow. The LE model did not describe the data well, even over a wide range of PCE solubility and macroscopic transverse dispersivity values. Model predictions assuming nonequilibrium mass-transfer-limited conditions and accounting for vertical flow gradients, however, resulted in a better fit to the data. These results have important implications for evaluating DNAPL pool dissolution in the field where subsurface heterogeneities are likely to be present.  相似文献   

4.
The influence of aquifer property correlation on multiphase fluid migration, entrapment and recovery was explored by incorporating correlated and uncorrelated porosity, permeability, and capillary pressure-saturation (Pc-Sat) parameter fields in a cross-sectional numerical multiphase flow model. Comparison of two-dimensional entrapped organic saturation distributions for a simulated tetrachloroethylene (PCE) spill in ensembles of aquifer realizations suggests that the degree of spatial correlation in Pc-Sat parameters exerts a controlling influence on dense nonaqueous phase liquid (DNAPL) spreading and redistribution in saturated aquifers. The predicted evolution of DNAPL source zones and resultant remediation efficiency under surfactant enhanced aquifer remediation (SEAR) also appear to be strongly influenced by the spatial correlation of aquifer parameters and multiphase flow constitutive relationships. Results for a limited number of realizations selected from each ensemble showed that removal of 60% to 99% of entrapped PCE could reduce dissolved contaminant concentration and mass flux by approximately two orders of magnitude under natural gradient conditions. Aqueous phase contaminant mass flux did not vary uniformly as a function of % DNAPL removed, however, and notable differences in behavior were observed for models incorporating correlated versus uncorrelated Pc-Sat and permeability fields. Although these results must be confirmed through analysis of additional realizations, it is likely that similar or larger differences between correlated and uncorrelated system behavior will be observed in aquifers with greater spatially variability than that of the nonuniform, homogeneous sand aquifer studied here. Funding for this research was provided by the United States Environmental Protection Agency, Great Lakes and Mid-Atlantic Center for Hazardous Substance Research under Grant No. R-825540, the Michigan Department of Environmental Quality under Contract No. Y80011, and the Strategic Environmental Research and Development Program under Project No. CU-1293. The content of this publication does not necessarily represent the views of these agencies and has not been subject to agency review.  相似文献   

5.
Air sparging was evaluated for remediation of tetrachloroethylene (PCE) present as dense nonaqueous phase liquid (DNAPL) in aquifers. A two-dimensional laboratory tank with a transparent front wall allowed for visual observation of DNAPL mobilization. A DNAPL zone 50 cm high was created, with a PCE pool accumulating on an aquitard. Detailed process control and analysis yielded accurate mass balances and insight into the mass-transfer limitations during air sparging. Initial PCE recovery rates were high, corresponding to fast removal of residual DNAPL within the zone influenced directly by air channels. The vadose zone DNAPL was removed within a few days, and the recovery in the extracted soil vapors decreased to low values. Increasing the sparge rate and pulsing the air injection led to improved mass recovery, as the pulsing induced water circulation and increased the DNAPL dissolution rate. Dissolved PCE concentrations both within and outside the zone of air channels were affected by the pulsing. Inside the sparge zone, aqueous concentrations decreased rapidly, matching the declining effluent PCE flux. Outside the sparge zone, PCE concentrations increased because highly contaminated water was pushed away from the air injection point. This overall circulation of water may lead to limited spreading of the contaminant, but accelerated the time-weighted average mass removal by 40% to 600%, depending on the aggressiveness of the pulsing. For field applications, pulsing with a daily or diurnal cycling time may increase the average mass removal rate, thus reducing the treatment time and saving in the order of 40% to 80% of the energy cost used to run the blowers. However, air sparging will always fail to remove DNAPL pools located below the sparge point because the air will rise upward from the top of a screen, unless very localized geological layers force the air to migrate horizontally. Unrecognized presence of DNAPL at chlorinated solvent sites residual and pools could potentially hamper success of air sparging cleanups, since the presence of small DNAPL pools, ganglia or droplets can greatly extend the treatment time.  相似文献   

6.
A three-dimensional two-phase flow model is coupled to a non-linear reactive transport model to study the efficacy of potassium permanganate treatment on dense, non-aqueous phase liquid (DNAPL) source removal in porous media. A linear relationship between the soil permeability (k) and concentration of manganese dioxide precipitate ([MnO2(s)]), k = ko + Srind [MnO2(s)], is utilized to simulate nodal permeability reductions due to precipitate formation. Using published experimental column studies, an Srind = −5.5 × 10−16 m2 L/mg was determined for trichloroethylene (TCE) DNAPL. This Srind was then applied to treatment simulations on three-dimensional TCE DNAPL source zones comprising either DNAPL at residual saturations, or DNAPL at pooled saturations.  相似文献   

7.
In situ chemical oxidation (ISCO) followed by a bioremediation step is increasingly being considered as an effective biphasic technology. Information on the impact of chemical oxidants on organohalide respiring bacteria (OHRB), however, is largely lacking. Therefore, we used quantitative PCR (qPCR) to monitor the abundance of OHRB (Dehalococcoides mccartyi, Dehalobacter, Geobacter, and Desulfitobacterium) and reductive dehalogenase genes (rdh; tceA, vcrA, and bvcA) at a field location contaminated with chlorinated solvents prior to and following treatment with sodium persulfate. Natural attenuation of the contaminants tetrachloroethene (PCE) and trichloroethene (TCE) observed prior to ISCO was confirmed by the distribution of OHRB and rdh genes. In wells impacted by persulfate treatment, a 1 to 3 order of magnitude reduction in the abundances of OHRB and complete absence of rdh genes was observed 21 days after ISCO. Groundwater acidification (pH<3) and increase in the oxidation reduction potential (>500 mV) due to persulfate treatment were significant and contributed to disruption of the microbial community. In wells only mildly impacted by persulfate, a slight stimulation of the microbial community was observed, with more than 1 order of magnitude increase in the abundance of Geobacter and Desulfitobacterium 36 days after ISCO. After six months, regeneration of the OHRB community occurred, however, neither D. mccartyi nor any rdh genes were observed, indicating extended disruption of biological natural attenuation (NA) capacity following persulfate treatment. For full restoration of biological NA activity, additional time may prove sufficient; otherwise addition electron donor amendment or bioaugmentation may be required.  相似文献   

8.
Based on the three Dimensional Particle Tracking Velocimetry (3D PTV) system, the characteristics of motion of particles with four different diameters were investigated under the steady flow conditions The longitudinal average velocity profiles of these particles were in accordance with Log-law, while the vertical and transverse velocities remained very low with minimal fluctuation. The time-average velocity of particles in the bed load layer was 8.50u., close to Bagnold's assumptionUn -60. The vertical concentration distribution of particles in the suspension region agreed with the Rouse equation. When the diameter of particles was relatively large, there existed an evident concentration gradient in the bed load layer.  相似文献   

9.
We present a 3D network model with particle tracking to upscale 3D Brownian motion of non-reactive tracer particles subjected to a velocity field in the network bonds, representing both local diffusion and convection. At the intersections of the bonds (nodes) various jump conditions are implemented. Within the bonds, two different velocity profiles are used. At the network scale the longitudinal dispersion of the particles is quantified through the coefficient DL, for which we evaluate a number of methods already known in the literature. Additionally, we introduce a new method for derivation of DL based on the first-arrival times distribution (FTD). To validate our particle tracking method, we simulate Taylor’s classical experiments in a single tube. Subsequently, we carry out network simulations for a wide range of the characteristic Péclet number Pe? to assess the various methods for obtaining DL. Using the new method, additional simulations have been carried out to evaluate the choice of nodal jump conditions and velocity profile, in combination with varying network heterogeneity. In general, we conclude that the presented network model with particle tracking is a robust tool to obtain the macroscopic longitudinal dispersion coefficient. The new method to determine DL from the FTD statistics works for the full range of Pe?, provided that for large Pe? a sufficiently large number of particles is used. Nodal jump conditions should include molecular diffusion and allow jumps in the upstream direction, and a parabolic velocity profile in the tubes must be implemented. Then, good agreement with experimental evidence is found for the full range of Pe?, including increased DL for increased porous medium heterogeneity.  相似文献   

10.
Invasion percolation (IP) models of dense non‐aqueous phase liquid (DNAPL) invasion into saturated horizontal fractures typically neglect viscous and gravity forces, as it is assumed that capillarity dominates in many situations. An IP model simulating DNAPL invasion into saturated horizontal fractures was modified to include gravity as a local effect. The model was optimized using a genetic algorithm, and demonstrated that the inclusion of gravity is important for replicating the architecture of the DNAPL invasion pattern. The optimized gravity‐included simulation showed the DNAPL invasion pattern to be significantly more representative of the experimentally observed pattern (80% accuracy) than did the optimized gravity‐neglected simulation (70% accuracy). Additional simulations of DNAPL invasion in 360 randomly generated fractures were compared with and without gravity forces. These simulations showed that with increasing fracture roughness, the minimum difference between simulations with and without gravity increases to 35% for a standard deviation of the mid‐aperture elevation field (SDz) of 10 mm. Even for low roughness (SDz = 0.1 mm), the difference was as high as 30%. Furthermore, a scaled Bond Number is defined which includes data regarding DNAPL type, media type and statistical characteristics of the fracture. The value of this scaled Bond Number can be used to determine the conditions under which gravity should be considered when simulating DNAPL invasion in a macroscopically horizontal fracture. Finally, a set of equations defining the minimum and maximum absolute percentage difference between gravity‐included and gravity‐neglected simulations is presented based on the fracture and DNAPL characteristics.  相似文献   

11.
Large laboratory columns (15.2 cm diameter, 183 cm long) were fed with groundwater containing trichloroethylene (TCE), were biostimulated and bioaugmented, and were monitored for over 7.5 years. The objective of the study was to observe how the selection of the carbon and energy source, i.e., whey, Newman Zone® standard surfactant emulsified oil and Newman Zone nonionic surfactant emulsified oil, affected the rate and extent of dechlorination. Column effluent was monitored for TCE and its degradation products, redox indicators (nitrate‐N, Fe(II), sulfate), and changes in iron mineralogy. Total bacteria and Dehalococcoides mccartyi strains were quantified using q‐PCR. Complete dechlorination was only observed in the whey treated columns, occurring 1 year after bioaugmentation with addition of a culture known to dechlorinate TCE to ethene, and 3 years later in the non‐bioaugmented column. The addition of the emulsified oils with or without bioaugmentation resulted in dechlorination only through cis‐DCE and vinyl chloride. While Dehalococcoides mccartyi strains are the only known bacteria that can fully dechlorinate TCE, their presence, either natural or augmented, was not the sole determiner of complete dechlorination. The establishment of a supporting microbial community and biogeochemistry that developed with continuous feeding of whey, in addition to the presence of D. mccartyi, were necessary to support complete reductive dechlorination. Results confirm that careful selection of a biostimulant is critical to the success of TCE dechlorination in complex soil environments.  相似文献   

12.
Two-dimensional numerical simulations of two-phase (DNAPL-water) flow in spatially correlated random fields demonstrate the influence of nonwetting phase (NWP) relative permeability–saturation (kr,NSW) relationships correlated to porous media intrinsic permeability (k). Both the volume of porous media invaded by the NWP and the length of time during which the NWP is migrating are under predicted if kr,Nk correlation is not accounted for in the model formulation. Not accounting for the kr,Nk correlation resulted in under predicting the volume of porous media invaded by up to approximately 10%, which is likely not significant for many practical applications. However, not accounting for the kr,Nk correlation resulted in under predicting field scale migration times by up to a factor of 4, which is likely significant in that the migration times are on the order of years to several decades for the DNAPL (1,2-DCE) considered in this study. The under prediction of migration times was greater for lower permeability aquifers.  相似文献   

13.
Cone penctrometer test (CPT) based Raman spectroscopy was used to identify separate phase tetrachloroethylene (PCE) and trichlorocthylene (TCE) contamination in the subsurface at two locations during field tests conducted at the U.S. Department of Energy's (DOE) Savannah River site. Clear characteristic Raman spectral peaks for PCE and TCE were observed at two sites and several depths during CPT deployment. Because of the uniqueness of a Raman spectrum for a given compound, these data are compelling evidence of the presence of the two compounds. The Raman spectral results correlated with high PCE and TCE concentrations in soil samples collected from the same subsurface zones, confirming that the method is a viable dense nonaqueous phase liquid (DNAPL) characterization technique. The Raman spectroscopic identification of PCE and TCE in these tests represents the first time that DNAPLs have been unequivocally located in the subsurface by an in situ technique.
The detection limit of the Raman spectroscopy is related to the probability of contaminant droplets appearing on the optical window in the path of the probe light. Based on data from this fieldwork the Raman technique may require a threshold quantity of DNAPL to provide an adequate optical cross section for spectroscopic response. The low aqueous solubility of PCE and TCE and relatively weak optical intensity of the Raman signal precludes the detection of aqueous phase contaminants by this method, making it selective for DNAPL contaminants only.  相似文献   

14.
The F-region peak electron densities NmF2 measured during daytime quiet geomagnetic conditions at low solar activity on January 22, 2008, April 8, 1997, July 12, 1986, and October 26, 1995, are compared. Ionospheric parameters are measured by the ionosonde and incoherent scatter radar at Millstone Hill and calculated with the use of a 1D nonstationary ionosphere–plasmasphere model of number densities and temperatures of electrons and ions at middle geomagnetic latitudes. The formation of the semiannual anomaly of the midlatitudinal NmF2 under daytime quiet geomagnetic conditions at low solar activity is studied. The study shows that the semiannual NmF2 anomaly occurs due to the total impact of three main causes: seasonal variations in the velocity of plasma drift along the geomagnetic field due to the corresponding variations in the components of the neutral wind velocity; seasonal variations in the composition and temperature of the neutral atmosphere; and the dependence of the solar zenith angle on a number of the day in the year at the same solar local time.  相似文献   

15.
The structure and dynamics of the ionosphere and plasmasphere at low solar activity under quiet geomagnetic conditions on January 15–17, 1985, and July 10–13, 1986, over Millstone Hill station and Argentine Islands ionosonde, the locations of which are approximately magnetically conjugate, have been theoretically calculated. The detected correction of the model input parameters makes it possible to coordinate the measured and calculated anomalous variations in the electron density NmF2 at the height hmF2 of the ionospheric F2 layer over Argentine Islands ionosonde as well as the calculated and measured values of NmF2 and electron temperature at the hmF2 height over Millstone Hill station. It has been shown that vibrationally excited N2 and O2 molecules almost do not influence the formation of the winter anomaly under the conditions of low solar activity. A difference between the influence of electronically excited O+ on N e ions under winter and summer conditions forms not more than 11% of the N e winter anomaly event in the F 2 layer and topside ionosphere. The model without electronically excited O+ ions reduces the duration of the N e winter anomaly event. It has been shown that the seasonal variations in the composition of the neutral atmosphere form mainly the NmF2 winter anomaly event over the Millstone Hill radar at low solar activity.  相似文献   

16.
Fine particulate organic matter (FPOM) represents a major component of stream organic matter budgets, and its dynamics greatly affect the productivity and metabolism of a stream community. FPOM transport dynamics has been well documented in high-gradient streams with rocky substrates, but information from low-gradient, sandy-bottom streams has been lacking. We estimated FPOM retention patterns in Payne Creek, a 2nd order Coastal Plain stream (USA), under naturally varying hydraulic conditions (discharge and velocity). Corn pollen, as an FPOM analogue, was released along with a conservative solute tracer and the particle retention coefficient (k p) was calculated by fitting the ratio of total pollen remaining in the water column against the longitudinal transport distance to an exponential decay model. Pollen k p (n = 4) ranged from 0.034 to 0.214 /m, and particle transport distance (S p) ranged from 4.7 to 29.7 m. The S p measured in Payne Creek was in the lowest range of previously reported values, and such rapid particle retention was attributed to the low channel slope and slow current velocity. S p was significantly correlated to water velocity and the channel friction factor, but not to discharge (Q). Two summer experiments conducted in contiguous stream segments resulted in the shortest (4.7 m) and longest (29.7 m) S p, despite the similar Q. This was attributed to the segment-scale channel alterations that occurred during the previous winter, which led to very different hydraulic conditions in the two stream segments. In Payne Creek, seasonal changes in hydrology and segment-scale variation in channel morphology were the main factors controlling FPOM transport and retention.  相似文献   

17.
The equipment and methodical characteristics of determining the vertical component of the ionospheric plasma motion velocity Vz based on an incoherent scatter radar of Institute of Ionosphere, National Academy of Sciences and Ministry of Education and Science of Ukraine (Kharkiv), which is the only radar of such type in Central Europe, are described. Based on the radar data, the patterns of altitude and diurnal variations in Vz near the maximum of solar cycle 24 for the typical geophysical conditions (around the summer and winter solstices, the spring and fall equinoxes) at low geomagnetic activity and the specifics of these changes during ionospheric storms are presented. The results of modeling of the dynamic processes in ionospheric plasma under the conditions of the undisturbed ionosphere, including the determination of altitudetime variations in the thermospheric wind velocity, are presented. It has been established that this velocity can significantly differ from the thermospheric wind velocity calculated by the known empirical global models. This difference is likely related to the regional features of thermospheric wind that are not shown in the global models.  相似文献   

18.
Formation mechanism of the spring–autumn asymmetry of the F2-layer peak electron number density of the midlatitudinal ionosphere, NmF2, under daytime quiet geomagnetic conditions at low solar activity are studied. We used the ionospheric parameters measured by the ionosonde and incoherent scatter radar at Millstone Hill on March 3, 2007, March 29, 2007, September 12, 2007, and September 18, 1984. The altitudinal profiles of the electron density and temperature were calculated for the studied conditions using a one-dimensional, nonstationary, ionosphere–plasmasphere theoretical model for middle geomagnetic latitudes. The study has shown that there are two main factors contributing to the formation of the observed spring–autumn asymmetry of NmF2: first, the spring–autumn variations of the plasma drift along the geomagnetic field due to the corresponding variations in the components of the neutral wind velocity, and, second, the difference between the composition of the neutral atmosphere under the spring and autumn conditions at the same values of the universal time and the ionospheric F2-layer peak altitude. The seasonal variations of the rate of O+(4S) ion production, which are associated with chemical reactions with the participation of the electronically excited ions of atomic oxygen, does not significantly affect the studied NmF2 asymmetry. The difference in the degree of influence of O+(4S) ion reactions with vibrationally excited N2 and O2 on NmF2 under spring and autumn conditions does not significantly change the spring–autumn asymmetry of NmF2.  相似文献   

19.
Microbiological degradation of perchloroethylene (PCE) under anaerobic conditions follows a series of chain reactions, in which, sequentially, trichloroethylene (TCE), cis‐dichloroethylene (c‐DCE), vinylchloride (VC) and ethene are generated. First‐order degradation rate constants, partitioning coefficients and mass exchange rates for PCE, TCE, c‐DCE and VC were compiled from the literature. The parameters were used in a case study of pump‐and‐treat remediation of a PCE‐contaminated site near Tilburg, The Netherlands. Transport, non‐equilibrium sorption and biodegradation chain processes at the site were simulated using the CHAIN_2D code without further calibration. The modelled PCE compared reasonably well with observed PCE concentrations in the pumped water. We also performed a scenario analysis by applying several increased reductive dechlorination rates, reflecting different degradation conditions (e.g. addition of yeast extract and citrate). The scenario analysis predicted considerably higher concentrations of the degradation products as a result of enhanced reductive dechlorination of PCE. The predicted levels of the very toxic compound VC were now an order of magnitude above the maximum permissible concentration levels. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

20.
Growth, stem morphology and some biochemical parameters were studied of one completely submerged (Myriophyllum spicatum) and two floating leaved macrophytes (Nymphoides peltata and Trapa japonica) under different turbulence velocities. The root mean square velocities of the high, medium and low amount of turbulence that was generated for the experiment were 2.18 ± 0.66, 1.48 ± 0.26 and 0.70 ± 0.07 cm s−1, respectively, in the microcosm. All three experimental plants survived exposed to all turbulence conditions provided, although a decrease in shoot elongation rate was associated with an increase in turbulence. Acceleration of tissue H2O2 generation and MDA content increased during the study period in all plant species. Oxidative enzymatic activities (POD, IAA and CKX) increased with time in plants under medium and high turbulence velocities. The shoot elongation rate, stem and leaf diameter, chlorophyll content and carbohydrate fractionations were found to be affected by this abiotic stress. It is evident from this study that high turbulence velocity inhibits normal metabolic activities of all three plants, while low to medium turbulence does not harm the floating leaved plants. Moreover, floating leaved plants were found to possess highly capable strategies to cope with this mechanical stress than completely submerged species.  相似文献   

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

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