首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
In the assessment of potentially contaminated land, the number of samples and the uncertainty of the measurements (including that from sampling) are both important factors in the planning and implementation of an investigation. Both parameters also effect the interpretation of the measurements produced, and the process of making decisions based upon those measurements. However, despite their importance, previously there has been no method for assessing if an investigation is fit‐for‐purpose with respect to both of these parameters. The Whole Site Optimised Contaminated Land Investigation (WSOCLI) method has been developed to address this issue, and to allow the optimisation of an investigation with respect to both the number of samples and the measurement uncertainty, using an economic loss function. This function was developed to calculate an ‘expectation of (financial) loss’, incorporating costs of the investigation itself, subsequent land remediation, and potential consequential costs. To allow the evaluation of the WSOCLI method a computer program ‘OCLISIM’ has been developed to produce sample data from simulated contaminated land investigations. One advantage of such an approach is that as the ‘true’ contaminant concentrations are created by the program, these values are known, which is not the case in a real contaminated land investigation. This enables direct comparisons between functions of the ‘true’ concentrations and functions of the simulated measurements. A second advantage of simulation for this purpose is that the WSOCLI method can be tested on many different patterns and intensities of contamination. The WSOCLI method performed particularly well at high sampling densities producing expectations of financial loss that approximated to the true costs, which were also calculated by the program. WSOCLI was shown to produce notable trends in the relationship between the overall cost (i.e., expectation of loss) and both the number of samples and the measurement uncertainty, which are: (a) low measurement uncertainty was optimal when the decision threshold was between the mean background and the mean hot spot concentrations. (b) When the hot spot mean concentration is equal to or near the decision threshold, then mid‐range measurement uncertainties were optimal. (c) When the decision threshold exceeds the mean of the hot spot, mid‐range measurement uncertainties were optimal. The trends indicate that the uncertainty may continue to rise if the difference between hot spot mean and the decision threshold increases further. (d) In any of the above scenarios, the optimal measurement uncertainty was lower if there is a large geochemical variance (i.e., heterogeneity) within the hot spot. (e) The optimal number of samples for each scenario was indicated by the WSOCLI method, and was between 50 and 100 for the scenarios considered generally; although there was significant noise in the predictions, which needs to be addressed in future work to allow such conclusions to be clearer.  相似文献   

2.
多数情况下,评估的研究区缺乏多年积累的地下水位及开采量资料。因此,通常采用土力学模型,根据太沙基固结理论,进行地面沉降量预测。在地面沉降量计算中,土层变形参数值的选取直接影响计算结果的精确性。深层含水层埋深通常大于100m。通过现场勘探,采取土样进行岩土物理力学性质试验获取参数值,成本过于昂贵。通过参考区域土层物理力学性质资料,结合工程经验获取参数值,预测结果的精度得不到保证。针对上述情况,作者提出用非线性优化算法-进化策略算法,根据实测地面沉降量反演土层变形参数值。通过这种方法确定参数值,既节约成本,又保证了计算结果的可靠性。进化策略算法通过模拟生物遗传及进化过程,利用转移概率来帮助指导搜索。搜索结果不依赖于初始点的选择,对于求解全局最优解有很强的鲁棒性。作者将进化策略算法用于某一工程实例土层变形参数的反演,结果表明了该算法的可行性及稳健性,值得在工程实践中推广应用。  相似文献   

3.
Soil loss evaluation by means of radioisotopic content measurements represents a promising technique, half-way between field surveys and theoretical models, which still suffers from its practical limits when applied at basin scale. A main limit is represented by the costs of field investigations and laboratory investigations on the high number of soil samples required by a traditional sampling scheme, such as regular grid or transects. A non-conventional soil sampling scheme was tested with the aim to consider a relative scanty number of soil samples and to verify the feasibility of this technique on large areas. This scheme was based on the hypothesis that land analysis and classification could point out areas characterised by homogeneous behaviour with respect to the 137Cs deposition and transmission model. A Geographic Information System (GIS)-aided procedure allowed to classify the selected basin area in pedo-morphological units, representative of the different pedologic, morphologic and land-use conditions, to locate few sampling points for each unit. Outcomes pointed out a low correlation between 137Cs contents and soil physical and compositional characteristics. Nevertheless, the isotopic methodology allowed to estimate a total soil loss value at basin scale almost consistent with both observed data, given by reservoir sedimentation measurements and estimates from the application of the RUSLE model. Thus, the results can be considered encouraging and they allow to deem that the isotopic methodology can be refined in order to account for erosion and deposition processes even at river basin scale and with a limited number of soil samples.  相似文献   

4.
Measurements taken to characterise environmental contamination contain uncertainty, which is generated by both field sampling and chemical analyses. Recently devised techniques have been applied for the first time to estimate this uncertainty in the commercial monitoring and assessment of contaminated land. The uncertainty reduces the reliability of the classification of the land that is made following a site investigation. The possible misclassification of areas of land, as a result of measurement uncertainty, can lead to substantial financial penalties, resulting from litigation or unnecessary remediation. Previous studies have developed methods for the estimation and financial optimisation of measurement uncertainty. These methods have now been applied to a series of six contrasting site investigations, which were conducted by various commercial organisations. The previous uses of these sites included a gas works, a tin mine and railway sidings. The measurement uncertainty was successfully estimated for each of the six investigations, showing its applicability to a wide range of different sampling methods, such as trial pits, window sampling and augering. The measurement uncertainty ranged widely between sites from 25% to 158%, indicating that investigations can differ widely in their reliability. The field sampling tended to generate the largest component of the measurement uncertainty when compared to the contribution from the chemical analysis. The Optimised Contaminated Land Investigation (OCLI) method was applied to each site, with the initial aim of estimating the financial losses that could be incurred as a result of misclassifying the land, due to the uncertainty. It showed that the expectation of loss value per sampling location ranged from only £58 at one site to over £ 11 000 at another. The optimal level of uncertainty that produced the minimal financial loss was then calculated for each site. It provided a reduction in the expectation of loss for the whole site of over £ 10 000 at two of the sites and over £90 000 at two others. These findings demonstrate that implementing concepts of uncertainty can have practical benefits in environmental monitoring, and can enable improvements to be made in the quality of sampling and hence of measurements in general.  相似文献   

5.
Soil texture is a key variable that reflect a number of soil properties such as soil permeability, water holding capacity, nutrient storage and availability, and soil erosion. The main objective of this study was to produce the kriged maps of soils of the Shahrekord region, central Iran. One hundred four soil samples were collected on a 375-m2 sampling grid from the depths of 0–30, 30–60, and 60–100 centimeter, and their particle sizes were determined using hydrometer method. The results showed a moderately spatial correlation in the soil particles among sampling soil layers and across the study area. Moreover, increasing clay and therewith observation of heavier soil textures is evident from surface to subsurface layers of the soils in the studied area due to rainfall and/or irrigation agriculture. These findings indicated that study of the soil texture variation with depth can be used as a clue for site-specific management and precision agriculture. Moreover, we suggest further analysis by using other data layers like topographical parameters, land use, parent material, soil erosion, and any other information which might influence the spatial distribution of soil texture.  相似文献   

6.
As part of an effort to estimate estuarine habitat values with respect to ecological indicators of benthic macrofaunal community condition, an optimal (effective and least costly) sampling protocol (sample unit size [area x depth], sieve mesh size, and sample number [n]) was determined. The goal was to use four ecological indicators (number of species, abundance, biomass, and fish and crab prey abundance) to detect differences among four intertidal habitats in Willapa Bay, Washington, United States. The four habitats were eelgrass (Zostera marina), Atlantic cordgrass (Spartina alterniflora), mud shrimp (Upogebia pugettensis), and ghost shrimp (Neotrypaea californiensis). Four sample unit areas (0.005, 0.010, 0.015, and 0.020 m2), two sample unit depths (0–5 and 0–10 cm), and two sieve mesh sizes (1.0 and 0.5 mm) were evaluated. The optimal sampling protocol was defined as the least costly protocol capable of reliably (statistical power, 1?β≥0.80) detecting significant (α=0.05) differences among ≥4 of the 6 pairwise habitat contrasts by ANOVA on all four ecological indicators. The relative cost of each sampling protocol was estimated as a direct function of the sample unit size and number and the cost-in-processing-time ratios of 1 (5 cm deep):1.7 (10 cm deep) and 1 (≥1.0 mm macrofauna size fraction); 2.5 (≥0.5 mm macrofauna size fraction), which were taken from previous studies. The optimal sampling protocol was 15–20, 0.01-m2×5-cm deep, 0.5-mm mesh samples per habitat.  相似文献   

7.
Drastically disturbed soils caused by opencast mining can result in the severe loss of soil structure and increase in soil compactness. To assess the effects of mining activities on reconstructed soils and to track the changes in reclaimed soil properties, the variability of soil properties (soil particle distribution, penetration resistance (PR), pH, and total dissolved salt (TDS)) in the Shanxi Pingshuo Antaibao opencast coal-mine inner dump after dumping and before reclamation was analyzed using a geostatistics method, and the number of soil monitoring points after mined land reclamation was determined. Soil samples were equally collected at 78 sampling sites in the study area with an area of 0.44 km2. Soil particle distribution had moderate variability, except for silt content at the depth of 0–20 cm with a low variability and sand content at the depth of 20–40 cm with a high variability. The pH showed a low variability, and TDS had moderate variability at all depths. The variability of PR was high at the depth of 0–20 cm and moderate at the depth of 20–40 cm. There was no clear trend in the variance with increasing depth for the soil properties. Interpolation using kriging displayed a high heterogeneity of the reconstructed soil properties, and the spatial structure of the original landform was partially or completely destroyed. The root-mean-square error (RMSE) can be used to determine the number of sampling points for soil properties, and 40 is the ideal sampling number for the study site based on cross-validation.  相似文献   

8.
Fluoride in drinking water has both beneficial and detrimental effects on public health, and a narrow range between .6 and 1.5 mg/L is optimal for consumption. However, natural groundwater sources exceed these guidelines affecting the entire population. This study aims to assess the distribution and controlling factors of fluoride concentration in the Tamiraparani River basin, South India. A total of 124 groundwater samples were analyzed for their fluoride content and other hydrogeochemical parameters. The fluoride concentration in the study area varied from .01 to 1.67 mg/L, and the highest concentrations were measured in the northern and central parts of the study area, which is underlain by charnockites and hornblende biotite gneiss. The sampling indicated (as per the Bureau of Indian Standards) that 53.9% of the area has fluoride concentrations below levels that are protective of teeth from dental caries (<.6 mg/L). .1% of the area is considered to be at risk of dental fluorosis, and the remaining 46% of the area is considered to have fluoride levels at desirable to permissible limit in groundwater. The groundwater in the study area belongs to Ca–Mg–Cl–SO4 and Ca–Mg–HCO3 types. A positive correlation between fluoride and TDS, Na+, K+ and HCO3 ? indicates its geogenic origin, and positive loading between pH and fluoride shows that alkaline environment enhances the dissolution of fluoride-bearing minerals into the groundwater. An empirical Bayesian kriging model was applied to interpolate the fluoride concentration in the study area. This geostatistical model is found to be better than other kriging methods, and it yielded an average standard error of .332 and root-mean-square standardized value of .986.  相似文献   

9.
The glacial isostatic adjustment (GIA) of the British Isles is of interest due to the constraints that can be provided on key model parameters such as the global meltwater signal, local ice sheet history and viscoelastic earth structure. A number of recent studies have modelled relative sea‐level (RSL) data from this region to constrain model parameters. As indicated in these studies, the sensitivity of these data to both local and global parameters results in a highly non‐unique problem. This study aims to address this inherent non‐uniqueness by combining a previously published British–Irish ice model that is based on the most recent geomorphological data with a new global ice sheet model that provides an accurate prediction of eustatic sea‐level change. In addition, constraints from Global Positioning System (GPS) measurements of present‐day vertical land motion are considered alongside the entirety of RSL data from both Great Britain and Ireland. A model solution is found that provides a high‐quality fit to both the RSL data and the GPS data. Within the range of earth viscosity values considered, the optimal data model fits were achieved with a relatively thin lithosphere (71 km), upper mantle viscosities in the range 4–6 × 1020 Pa s and lower mantle viscosities ≥ 3 × 1022 Pa s. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

10.
In the Daqing region of China there are 34 groundwater well fields with a groundwater withdrawal of 81.9×104 m3/d. Due to over-abstraction of the groundwater resources from the 1960s to present, a cone of depression up to 4,000 km2 has formed in the area. To monitor the change in the groundwater environment, it is necessary to design an effective groundwater-monitoring network. The sites for monitoring groundwater level were selected by applying the finite-element method coupled with Kalman filtering to the area in which the groundwater resources have been extensively exploited. The criterion is a threshold value of the standard deviation of estimation error. This threshold value is determined by the tradeoff between maximum information and minimum cost, in which the maximum information is characterized by the standard deviation and the minimum cost is equivalent to the number of observation wells. The groundwater flow model was calibrated by an optimal algorithm coupled the finite-element method with Kalman filtering by using the data from 16 observation wells from 1986 to 1993. A simulation algorithm coupled with the finite-element method with Kalman filtering analyzed the location data obtained from the existing 38 observation wells in the same region. The spatial distribution of standard deviation of estimation error is computed and the locations that have the maximum standard deviation are selected as additional sites for augmenting the existing observational well network at a given threshold value of the standard deviation surface. Based on the proposed method for selecting a groundwater level monitoring network, an optimal monitoring network with 88 observation wells with the measurement frequency of 12 times per year is selected in the Daqing region of China.  相似文献   

11.
《China Geology》2019,2(4):501-511
With the adjustment of industrial structure, many high pollution enterprises will gradually shut down. This process will produce a large number of high-risk pollution plots and cause a series of environmental problems. In this study, geophysical detection and direct push technology were applied to investigate soil and groundwater pollution in a typical brownfield in Zhenjiang, Jiangsu province, China. The results showed that Chlorotoluene and its ramification were the main pollutants for the brownfield. Pollution decreased with depth and was quite uneven in the brownfield. Production and wastewater treatment areas were heavily polluted, where volatile organic compounds contaminated and semi-volatile organic compounds were found 12 m and 15 m below the land surface, respectively. About 18880 m3 groundwater was contaminated. This study is significant to develop a geological survey work mode for the postindustrial era.  相似文献   

12.
This study proposes an interactive sampling strategy for locating the hot spot or maximum regions of a concerned attribute in a given area of survey. In the proposed strategy, information analysis is performed based on the ordinary kriging from the existing sample data to suggest a new batch of samples under the criterion of the highest information free energy. The information free energy (F) is a function of information energy (U) and information entropy (S) through F = U - TS, where T is information temperature and is used to coordinate the contribution of U and S to F. Information energy is the value of the concerned attribute, and information entropy is the transformed error variance of kriging and therefore measures the evenness and density of coverage of samples over the area under survey. At early sampling batches, information temperature is high and information entropy dominates the information free energy, and samples are suggested to give an even and dense enough coverage of the whole area under investigation. As samples accumulate, information temperature decreases to enlarge the contribution of information energy, and future samples are taken toward the locations with high attribute values. Two examples demonstrate the efficiency and effectiveness of the proposed sampling strategy in locating the hot spot regions of various fields: (1) a heavy metal contaminated site reproduced by modeling on 55 real field data; (2) a simulated two-dimensional field by the random phase volume (RPV) model. The results show that the proposed strategy, a robust interactive sampling procedure, is able to locate hot spot regions without compromising with the overall profile of an under-survey area.  相似文献   

13.
Sampling design optimization for spatial functions   总被引：4，自引：0，他引：4  
A new procedure is presented for minimizing the sampling requirements necessary to estimate a mappable spatial function at a specified level of accuracy. The technique is based on universal kriging, an estimation method within the theory of regionalized variables. Neither actual implementation of the sampling nor universal kriging estimations are necessary to make an optimal design. The average standard errorand maximum standard error of estimationover the sampling domain are used as global indices of sampling efficiency. The procedure optimally selects those parameters controlling the magnitude of the indices, including the density and spatial pattern of the sample elements and the number of nearest sample elements used in the estimation. As an illustration, the network of observation wells used to monitor the water table in the Equus Beds of Kansas is analyzed and an improved sampling pattern suggested. This example demonstrates the practical utility of the procedure, which can be applied equally well to other spatial sampling problems, as the procedure is not limited by the nature of the spatial function.  相似文献   

14.
Watershed degradation due to soil erosion and sedimentation is considered to be one of the major environmental problems in Iran. In order to address the critical conditions of watershed degradation in arid and semiarid regions, a study based on the Modified Pacific Southwest Inter-Agency Committee (MPSIAC) model was carried out at Golestan watershed, northeast of Iran. The model information layers comprising nine effective factors in erosion and sedimentation at the watershed site were obtained by digitalization and spatial interpolation of the basic information data in a GIS program. These factors are geology, soil, climate, runoff, topography, land cover, land use, channel, and upland erosion. The source data for the model were obtained from available records on rainfall and river discharge and sediment, topography, land use, geology, and soil maps as well as field surveys and laboratory analysis. The results of the MPSIAC model indicated that 60.75 % (194.4 km2) and 54.97 % (175.9 km2) of the total watershed area were classified in the heavy sedimentation and erosion classes, and the total basin sediment yield and erosion were calculated as 4,171.1 and 17,813.4 m3 km?2 year?1, respectively. In the sensitivity analysis, it was found that the most sensitive parameters of the model in order of importance were topography (slope), land cover and use, runoff, and channel erosion (R 2?=?0.92–0.94), while geology, climate (rainfall), soil, and upland erosion factors were found to have moderate effect to the model output (R 2?=?0.74–0.59).  相似文献   

15.
选择吉林省公主岭市大岭地区作为东北平原区典型代表区域开展土壤元素空间变异性、经典统计学合理取样数及不同采用密度数据空间插值对比研究。结果表明:(1)受地形平坦及成土母质相对单一等因素影响,研究区土壤元素空间变异性总体较小,大部分以轻中度变异为主(变异系数<15%),受人为因素影响较大的Cd、Hg变异系数分别为35.3%、136.6%,属于高度变异。(2)经典统计学确定的研究区合理采样数为80,该样本量可在95%的置信区间及允许误差为30%的条件下反应区内土壤元素含量的均值与方差,但因未考虑样本的空间属性,不足以反应区内土壤元素空间变异特征,具有一定的局限性。(3)通过对均匀抽稀后4种不同采样密度数据与实测数据空间插值对比研究,在定量评估空间插值相对误差、地块预测值相对误差及预测等级与实测等级一致性的基础上,结合土地质量地球化学调查工作精度要求,提出研究区地块尺度地球化学评价工作合理采样密度为8个点/km2,该密度可在确保评价精度的前提下,大幅减少采样数量和工作成本。上述结论为东北平原及类似地区大面积开展地块尺度土地质量地球化学评价工作提供了关键的技术支撑,对进一步完善土地质量地球化学评价方法技术具有重要意义。  相似文献   

16.
Improved knowledge of processes determining groundwater quality is an important precondition for the solution of various ecological and water management problems. In areas with highly fluctuating groundwater levels, time-limited access, local pollution sources or temporary interactions between surface water and groundwater, a temporary groundwater sampling technique could be of advantage. Furthermore, depth-specific sampling is of high value for investigating groundwater pollution related to seepage or surface water infiltration. A stainless steel core probe has been developed to obtain groundwater samples and to measure the hydraulic head distribution at various defined depths. The sampling technique is applicable only for non-volatile water constituents. An advantage of the core probe is that it can be driven into soil or sediments using ordinary low cost percussion equipment. The probe enables hydraulic head measurements and water sampling over vertical intervals of 0.3 m. Results from field experiments using the stainless steel core probe were in good correspondence with results from groundwater sampling at nearby observation wells. In the upper layer of the aquifer, the intrinsic spatial change in concentrations of sulphate, chloride and other water constituents is a function of distance between observation points and groundwater surface. Results indicate strong effects of a fluctuating groundwater level on groundwater quality at certain depths.  相似文献   

17.
A pot experiment was conducted to monitor the dynamic response of photosynthesis of Amorpha fruticosa seedlings to different concentrations of petroleum-contaminated soils from April to September. The results showed that the photosynthetic rates, stomatal conductance and transpiration rate of seedlings significantly decreased in 5–20 g kg?1 petroleum-contaminated soil during the three given sampling period of July 31 (early), August 30 (mid-term) and September 29 (late). However, the intercellular CO2 concentration significantly increased in 10 g kg?1 contaminated soil, while declined in 20 g kg?1 contaminated soil during the early sampling period as well as in 20 g kg?1 contaminated soil during the late sampling period. The leaf relative water content of seedlings significantly increased in 20 g kg?1 contaminated soil during the early sampling period, while it dropped dramatically in 15–20 g kg?1 contaminated soil during the late sampling period. The contents of chlorophyll a, chlorophyll b and the total chlorophyll of seedlings showed a sharp decline during the three sampling periods in contaminated soil. Comprehensively, considering the negative effects of petroleum on the photosynthesis, growth performance and remediation effect on petroleum of A. fruticosa seedlings, this plant was tolerant of petroleum-contaminated soil and was potentially useful for the phytoremediation of petroleum-contaminated sites in northern Shaanxi, China.  相似文献   

18.
The groundwater in the upper Kodaganar basin is contaminated due to the discharge of effluents from tannery industries. The water in the wells, whose physico-chemical characteristics are altered due to the influence of the effluents, is statistically analyzed. The physico-chemical variables such as EC, Na+, K+, Ca2+, Mg2+, F?, Cl?, HCO3 ?,CO3 2?, NO3 ?, SO4 2?, pH, and Crtotal were used for this study. An attempt was made to identify the contaminated wells based on suitability for drinking, suitability for industrial requirements, and through principal component analysis (PCA). Classification based on suitability helped in identifying the contaminated wells. However, this resulted in failure when identifying the wells that are contaminated by tanneries. PCA has proved to be effective in the segregation of contaminated wells influenced by tannery industries. The physico-chemical variables that are 13 in number are transformed into two orthogonal components and Eigen values based on the variance. The Eigen values are used to select the first two principal components PC1 (7.26) and PC2 (2.24) that accounted for 73.04% variance in the data. The components of the variables and the wells are plotted in a biplot to isolate the contaminated samples. The contaminated samples are analyzed in the spatial domain in geographic information system and found to be clustered around the tannery belt. The study reveals that 35% of the samples are contaminated due to discharge from tannery industries.  相似文献   

19.
土壤污染防治工作已成为提升耕地质量、保护国土生态安全的重要任务之一。为了科学预测我国大宗农作物(如水稻、小麦)与土壤重金属含量的关系,减少安全利用类农用地的大量农产品与土壤的协同监测,实现重金属污染农用地的安全利用,本研究以重金属Cd为例,选取对水稻、小麦Cd含量影响较大的土壤Cd含量、土壤pH值、土壤阳离子交换量(CEC)和土壤有机碳(OC)含量作为输入因子,水稻、小麦Cd含量作为输出因子,分别建立了多元回归模型与神经网络模型。结果表明:水稻、小麦Cd含量与土壤Cd含量呈现正相关关系;模拟出的水稻、小麦与土壤Cd的多元线性回归模型的预测能力分别为67.8%和83.8%;利用神经网络构建了水稻、小麦Cd含量预测模型,在训练集、验证集和测试集中都表现出很好的预测能力,R值均大于多元线性回归模型,且MSE(均方误差)值较小,神经网络对水稻、小麦Cd含量预测具有很好的适用性,模拟精度总体优于多元回归预测模型。研究结果可为污染农用地的安全利用评价及优化配置提供一定的理论依据和参考。  相似文献   

20.
The environmental issues associated with mining have damaged the industry’s substantial global economic value. In particular, the mining industry has a negative legacy of contaminated land. The effective reclamation of contaminated soil is therefore required before former mining land can be further developed for residential and commercial purposes. The objective of this study was to technically evaluate the feasibility of reclamation techniques for agricultural soils contaminated with toxic elements (As, Cd, Cu, Pb, and Zn) associated with metal mining. The reclamation methods investigated were covering without stabilization, covering with stabilization, and exchange with stabilization. The thickness of the soil layer used in covering and exchange was in the range of 30–50 cm. Limestone, furnace slag, and a mixture of limestone and furnace slag were applied as soil amendments. After reclamation, the contamination level in surface tillage soils and crops was monitored regularly. Four years of monitoring data revealed that surface soil contamination levels could be maintained at acceptable levels, although at some sites, the metal levels in crops exceeded legislative limits. Soil reclamation at former mining sites in Korea has not yet been perfected, but the results of this study show that there is potential for safe agricultural operations on large sites in a cost-effective manner, as long as the appropriate control of surface soil contamination and adequate agronomic management is undertaken.  相似文献   

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

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