首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 468 毫秒
1.
To assess the contribution of accumulated winter precipitation and glacial meltwater to the recharge of deep ground water flow systems in fracture crystalline rocks, measurements of environmental isotope ratios, hydrochemical composition, and in situ parameters of ground water were performed in a deep tunnel. The measurements demonstrate the significance of these ground water recharge components for deep ground water flow systems in fractured granites of a high alpine catchment in the Central Alps, Switzerland. Hydrochemical and in situ parameters, as well as delta(18)O in ground water samples collected in the tunnel, show only small temporal variations. The precipitation record of delta(18)O shows seasonal variations of approximately 14% and a decrease of 0.23% +/- 0.03% per 100 m elevation gain. delta(2)H and delta(18)O in precipitation are well correlated and plot close to the meteoric water line, as well as delta(2)H and delta(18)O in ground water samples, reflecting the meteoric origin of the latter. The depletion of 18O in ground water compared to 18O content in precipitation during the ground water recharge period indicates significant contributions from accumulated depleted winter precipitation to ground water recharge. The hydrochemical composition of the encountered ground water, Na-Ca-HCO3-SO4(-F), reflects an evolution of the ground water along the flowpath through the granite body. Observed tritium concentrations in ground water range from 2.6 to 16.6 TU, with the lowest values associated with a local negative temperature anomaly and anomalous depleted 18O in ground water. This demonstrates the effect of local ground water recharge from meltwater of submodern glacial ice. Such localized recharge from glaciated areas occurs along preferential flowpaths within the granite body that are mainly controlled by observed hydraulic active shear fractures and cataclastic faults.  相似文献   

2.
Adverse geologies are often encountered during tunnel construction, which could seriously endanger the construction. To ensure the safety, it is essential to detect adverse geologies and their water‐bearing situation ahead the tunnel face. Ground‐penetrating radar is a suitable instrument, but the accurate interpretation of its detection results is difficult. In this paper, at first, an improved back projection imaging algorithm is proposed, which can make reflection waves closer to the real geological boundaries with few artificial clutters. And then, forward modelling of ground‐penetrating radar is carried out for typical adverse geologies, such as karst caves, faults, fractured rock masses, fracture network, and water‐bearing body. Their corresponding response features are obtained, accumulating experience for geological interpretation. The above two methods provide the basis for target identification and geological interpretation. In the last part, the application of the above two methods in several engineering cases are given, and their effectiveness is verified.  相似文献   

3.
The problem of water flow around a tunnel cavity located in the saturated capillary fringe on top of a very permeable, freely draining substratum is considered for the critical non-leakage condition when there is uniform vertical downward flow through the upper surface of the saturated region. In this critical condition the soil–water pressure is equal to zero everywhere on the cavity wall that is also a streamline. The conditions at the upper fringe boundary are that the soil–water pressure is equal to the air-entry value of the soil and the flux through this surface is the uniform infiltration rate. The cavity surface and the fringe boundary which is elevated above the cavity position, are found through conformal mapping and the use of integral representations of non-standard mixed boundary-value problems. They are calculated for a range of infiltration rates and compared with those obtained by assuming the upper fringe boundary to be horizontal. The exact analysis given here gives larger tunnel cavities than those given by the approximate treatment of the problem. The results have application in the design of underground repositories against entry of seepage water, the construction of protective capillary barriers and in the design of interceptor drainage systems.  相似文献   

4.
Calculation of ground water ages--a comparative analysis   总被引:1,自引:0,他引:1  
Castro MC  Goblet P 《Ground water》2005,43(3):368-380
Ground water age is a fundamental, yet complex, concept in ground water hydrology. Discrepancies between results obtained through different modeling approaches for ground water age calculation have been reported, in particular, between ground water ages modeled by advection and direct simulation of ground water ages (e.g., age-mass approach), which includes effects of advection and dispersion. Here, through a series of two-dimensional (2D) simulations, the impact of water mixing through advection and dispersion on modeled 14C and directly simulated ground water ages is assessed. Impact of dispersion on modeled ages is systematically stronger in areas where water velocities are smaller and far more pronounced on 14C ages. This effect is also observed in one-dimensional models. 2D simulations show that longitudinal dispersion generally acts as a "source" of 14C, while vertical dispersion acts as a "sink," leading to apparent younger or older modeled 14C ages as compared to advective and directly simulated ground water ages. The presence of permeable and impermeable faults provides an equally important source for discrepancies, leading to major differences in modeled ages among the three methods considered. Overall, our results show that a 14C modeling approach using a solute transport model for calculating ground water age appears to be more reliable in ground water systems without faults and where water velocities are relatively high than in systems that are relatively more heterogeneous and those where faults are present. Among the three modeling approaches considered here, direct simulation of ground water age seems to yield the most consistent results in complex, heterogeneous ground water flow systems, giving a vertical age structure consistent with ages expected from consideration of the flow system.  相似文献   

5.
Loess tunnels are a common geo-hazard in the Loess Plateau and not only cause considerable soil and water loss, but also aggravate and even induce the occurrence and development of other disasters such as ground fissures, mudflows, collapses, and landslides. To date, research on the hydrological characteristics and erosion behaviour of loess tunnel systems has focused on field investigation data and limited river basin observation data, whereas field test information and data are very scarce. In this study, field surveys, observations, field scouring experiments, and laboratory-based sediment percentage tests were conducted to analyse the erosion characteristics, spatial distribution, and hydrological characteristics of a large-scale loess tunnel system in the hilly Loess Plateau southeast region of northern China. The results showed that the loess tunnel erosion exhibited periodicity. Tunnel erosion in each period shows a similar erosion process, that is, thin-layer water flow erosion and lateral expansion, tunnel wall collapse and deposit due to the loss of support from the lateral erosion, and erosion and transport of deposits by water flow. Waterfall erosion, lateral erosion, headward erosion, and the resulting collapses were the main forms of tunnel erosion. Besides this, the base level of erosion significantly affects the erosion characteristics of the loess tunnel. The hydrological characteristics during field scouring experiments exhibited three different stages: a lag effect, attributed to the temporary loss of water velocity in the first stage; small water flow fluctuations in the second stage; and an increase in total seepage loss with increased water injection flow in the third stage. The erosion rate was positively correlated with the flow quantity. The results of this study not only provide valuable reference data for research on the mechanism and velocity of erosion events in loess, but also provide a theoretical basis for the prevention of loess tunnel disasters in engineering construction. © 2020 John Wiley & Sons, Ltd.  相似文献   

6.
Electro-osmosis (EO), the movement of water through porous media in response to an electric field, offers a means for extracting contaminated ground water from fine-grained sediments, such as clays, that are not easily amenable to conventional pump-and-treat approaches. The EO-induced water flux is proportional to the voltage gradient in a manner analogous to the flux dependence on the hydraulic gradient under Darcy's law. The proportionality constant, the soil electro-osmotic conductivity or keo, is most easily measured in soil cores using bench-top tests, where flow is one-dimensional and interfering effects attributable to Darcy's law can be directly accounted for. In contrast, quantification of EO fluxes and keo in the field under deployment conditions can be difficult because electrodes are placed in ground water wells that may be screened across a heterogeneous mixture of lithologies. As a result, EO-induced water fluxes constitute an approximate radial flow system that is superimposed upon a Darcy flow regime through permeable pathways that may or may not be coupled with hydraulic head differences created by the EO-induced water fluxes. A single well comparative tracer test, which indirectly measures EO fluxes by comparing wellbore tracer dilution rates between background and EO-induced water fluxes, may provide a means for routinely quantifying the efficacy of EO systems in such settings. EO fluxes measured in field tests through this technique at a ground water contamination site were used to estimate a mean keo value through a semianalytic line source model of the electric field. The resulting estimate agrees well with values reported in the literature and with values obtained with bench-top tests conducted on a soil core collected in the test area.  相似文献   

7.
以隧道掘进机(Tunnel Boring Machine,TBM)为例模拟了隧道强干扰环境下,瞬变电磁超前探测的响应曲线,系统分析了异常体(以直立充水断层为例)与掌子面距离、围岩电阻率差异、TBM长度、异常体规模等条件下的曲线特征和影响规律,发现TBM干扰源表现为低电阻率目标特征,其影响主要集中在早期,对于电性差异较大或目标规模较大的低电阻率异常(充水断层)模型能够明显地通过衰减曲线区分.根据电磁场叠加原理,将隧道腔体中包含TBM模型的响应减去纯隧道腔体响应可以获得TBM的响应信号,以此作为干扰背景,从实际包含TBM和充水断层的隧道模型总响应中减除,获得去除TBM干扰的响应信号.通过8组算例进行对比,发现经过校正的衰减曲线与模型计算曲线吻合较好,视电阻率曲线差异相对较小,能够表现探测区域的电性分布情况,确认该方法在不同情况下的适用性.即使在TBM响应计算时给定背景电阻率与实际电阻率差异达到100%的情况下,依然能够通过校正获得合理的响应信号和视电阻率曲线.该方法不仅仅适用于隧道环境,对于其他诸如地面、航空、半航空、海洋瞬变电磁勘探同样适用.  相似文献   

8.
There are few studies on the hydrogeology of sedimentary rock aquitards although they are important controls in regional ground water flow systems. We formulate and test a three-dimensional (3D) conceptual model of ground water flow and hydrochemistry in a fractured sedimentary rock aquitard to show that flow dynamics within the aquitard are more complex than previously believed. Similar conceptual models, based on regional observations and recently emerging principles of mechanical stratigraphy in heterogeneous sedimentary rocks, have previously been applied only to aquifers, but we show that they are potentially applicable to aquitards. The major elements of this conceptual model, which is based on detailed information from two sites in the Maquoketa Formation in southeastern Wisconsin, include orders of magnitude contrast between hydraulic diffusivity (K/S(s)) of fractured zones and relatively intact aquitard rock matrix, laterally extensive bedding-plane fracture zones extending over distances of over 10 km, very low vertical hydraulic conductivity of thick shale-rich intervals of the aquitard, and a vertical hydraulic head profile controlled by a lateral boundary at the aquitard subcrop, where numerous surface water bodies dominate the shallow aquifer system. Results from a 3D numerical flow model based on this conceptual model are consistent with field observations, which did not fit the typical conceptual model of strictly vertical flow through an aquitard. The 3D flow through an aquitard has implications for predicting ground water flow and for planning and protecting water supplies.  相似文献   

9.
通过大型振动台模型试验并采用Midas-GTS有限元软件进行模拟计算,研究黄土隧道洞口段在地震作用下的动力响应特征、破坏过程和地震波在模型中的传递规律,分析影响黄土隧道洞口段地震动力响应的主要因素。结果表明:边坡沿弧形开裂面的垮塌受坡脚剪切和坡顶拉裂的共同作用;边坡会对其卓越频率内的地震波产生明显放大效应,且在1/2坡高以上放大效应出现饱和现象;隧道临空面是影响隧道洞口段地震动力响应的主要因素。考虑进洞高程效应时隧道洞口段抗震设防长度可取距洞口5倍洞径范围。振动台模型试验与数值计算在位移、加速度、应力三个响应特征上吻合较好,证明二者结果合理可靠。研究成果可为隧道工程设计和地下结构抗震理论研究提供有益参考。  相似文献   

10.
Heat as a ground water tracer   总被引:40,自引:0,他引:40  
Anderson MP 《Ground water》2005,43(6):951-968
Heat carried by ground water serves as a tracer to identify surface water infiltration, flow through fractures, and flow patterns in ground water basins. Temperature measurements can be analyzed for recharge and discharge rates, the effects of surface warming, interchange with surface water, hydraulic conductivity of streambed sediments, and basin-scale permeability. Temperature data are also used in formal solutions of the inverse problem to estimate ground water flow and hydraulic conductivity. The fundamentals of using heat as a ground water tracer were published in the 1960s, but recent work has significantly expanded the application to a variety of hydrogeological settings. In recent work, temperature is used to delineate flows in the hyporheic zone, estimate submarine ground water discharge and depth to the salt-water interface, and in parameter estimation with coupled ground water and heat-flow models. While short reviews of selected work on heat as a ground water tracer can be found in a number of research papers, there is no critical synthesis of the larger body of work found in the hydrogeological literature. The purpose of this review paper is to fill that void and to show that ground water temperature data and associated analytical tools are currently underused and have not yet realized their full potential.  相似文献   

11.
This study is aimed at investigating the vertical velocity profile of flow passing over a vegetal area by an analytical approach. The soil ground is considered as pervious and thus non-zero velocity at the ground surface can be estimated. The soil and vegetation layers are regarded as homogeneous and isotropic porous media. Therefore the solution of the flow can be obtained by applying the theory of turbulent flow and Biot’s theory of poroelasticity after dividing the flow field into three layers: homogenous water, vegetation and pervious soil. The velocity distribution is compared with the experimental data of [Rowiński PM, Kubrak J. A mixing-length model for predicting vertical velocity distribution on flows through emergent vegetation. J Hydrol Sci 2002;47(6):893–904] to show its validity. In addition, five dimensionless parameters denoting the variation of slope, permeability of soil, Reynolds stress, density of vegetation, and relative height of vegetation are proposed to reveal their effects on the surface water flow. The analytical solutions of flow velocity can also be simplified into simpler expressions to describe the flow passing over a non-vegetated area.  相似文献   

12.
《Continental Shelf Research》2007,27(10-11):1379-1398
Water masses in Monterey Bay are determined from the CTD casts of the Monterey Ocean Observing System (MOOS) Upper-water-column Science Experiment (MUSE) August 2000 dataset. It is shown through cluster analysis that the MUSE 2000 CTD dataset contains 5 water masses. These five water masses are: bay surface water (BSW), bay warm water (BWW), bay intermediate water (BIW), sub arctic upper water (SUW) and North Pacific deep water (NPDW). The BWW is a new water mass that exists in one area and is attributed to the effects of solar heating. The volumes occupied by each of the water masses are obtained. The BIW water is the most dominant water mass and occupies 68.8% of the volume. The statistical means and standard deviations for each water parameter, including spiciness and oxygen concentration, are calculated during separate upwelling and relaxed periods.The water mass content and structure are analyzed and studied during upwelling and a relaxed period. During upwelling, along a CTD track off Pt. Ano Nuevo, the water mass T, S distribution tended to be organized along three branches. Off Pt. Ano Nuevo the innovative coastal observation network (ICON) model showed the formation of a cyclonic eddy during the analyzed upwelling period. In time the eddy moved southwest and became absorbed into the southerly flow during the initial phases of the following wind-relaxed period.  相似文献   

13.
针对黄土边坡与隧道洞口段衬砌的相互作用问题,运用数值模拟的方法分析了以不同进洞高程进洞时黄土隧道洞口段衬砌的动力响应特征和洞口仰坡的动力稳定性。结果表明:进洞高程越大,洞口段隧道衬砌的位移响应与内力响应越大;随着进洞高程的增大,坡面位移放大系数在减小,不同进洞高程进洞时坡面位移放大系数均呈先增大后减小再增大的变化趋势。在0.2~0.6H时变化最为剧烈,0.4H左右时位移放大系数达到了最大值;不同进洞高程进洞时坡面中心和水平方向距离隧道结构1.5D处的坡面位移放大系数变化趋势基本一致,其大小关系为:纯边坡位移放大系数 < 有隧道结构中面位移放大系数 < 距隧道1.5D位移放大系数;随着进洞高程的增大,剪应变增量和坡面位移均在减小,坡面的稳定性在增强。该研究可为黄土地区隧道进洞高程的选择提供一定的参考。  相似文献   

14.
Rains MC  Mount JF 《Ground water》2002,40(5):552-563
In this study, we identify the origin of shallow ground water that supports regionally unique plant and wildlife habitats in a riparian and reservoir-fringe system using isotopic and chemical procedures. This study was conducted where Little Stony Creek flows into East Park Reservoir on the east front of the Coast Range, northern California. Little Stony Creek water, Hyphus Creek water, Franciscan Complex regional ground water, Great Valley Group regional ground water, and local shallow ground water were collected during wet and dry seasons and were analyzed for deuterium, oxygen-18, temperature, pH, redox potential, conductivity, and major cation and anion concentrations. Turnover in the local flow system is rapid indicating that local shallow ground water is dependent on recent recharge. Local shallow ground water is recharged primarily by Little Stony Creek water and Franciscan Complex ground water. In the wet season, Little Stony Creek is the more prominent source of local shallow ground water, and the ratio of Little Stony Creek water to Franciscan Complex ground water decreases with distance from the channel. In the dry season, Franciscan Complex ground water is the more prominent source of local shallow ground water, and the ratio of Little Stony Creek water to Franciscan Complex ground water decreases with distance down the valley. Franciscan Complex ground water discharges to local shallow ground water throughout the year, primarily because the local flow system is a regional low that lies perpendicular to the Franciscan Complex ground water flowpath. Little Stony Creek is a more prominent source of ground water in the wet season than in the dry season because Little Stony Creek flows continuously through the alluvial reach in the wet season and intermittently through the alluvial reach in the dry season. Extensive ground water withdrawals from the Franciscan Complex flow system could reduce the amount of water available to the local flow system, particularly during the dry season, and could substantially reduce the geographic extent of the regionally unique plant and wildlife habitats.  相似文献   

15.
Transient recharge to the water table is often not well understood or quantified. Two approaches for simulating transient recharge in a ground water flow model were investigated using the Trout Lake watershed in north-central Wisconsin: (1) a traditional approach of adding recharge directly to the water table and (2) routing the same volume of water through an unsaturated zone column to the water table. Areas with thin (less than 1 m) unsaturated zones showed little difference in timing of recharge between the two approaches; when water was routed through the unsaturated zone, however, less recharge was delivered to the water table and more discharge occurred to the surface because recharge direction and magnitude changed when the water table rose to the land surface. Areas with a thick (15 to 26 m) unsaturated zone were characterized by multimonth lags between infiltration and recharge, and, in some cases, wetting fronts from precipitation events during the fall overtook and mixed with infiltration from the previous spring snowmelt. Thus, in thicker unsaturated zones, the volume of water infiltrated was properly simulated using the traditional approach, but the timing was different from simulations that included unsaturated zone flow. Routing of rejected recharge and ground water discharge at land surface to surface water features also provided a better simulation of the observed flow regime in a stream at the basin outlet. These results demonstrate that consideration of flow through the unsaturated zone may be important when simulating transient ground water flow in humid climates with shallow water tables.  相似文献   

16.
Important observations and parameters for a salt water intrusion model   总被引:1,自引:0,他引:1  
Shoemaker WB 《Ground water》2004,42(6-7):829-840
Sensitivity analysis with a density-dependent ground water flow simulator can provide insight and understanding of salt water intrusion calibration problems far beyond what is possible through intuitive analysis alone. Five simple experimental simulations presented here demonstrate this point. Results show that dispersivity is a very important parameter for reproducing a steady-state distribution of hydraulic head, salinity, and flow in the transition zone between fresh water and salt water in a coastal aquifer system. When estimating dispersivity, the following conclusions can be drawn about the data types and locations considered. (1) The "toe" of the transition zone is the most effective location for hydraulic head and salinity observations. (2) Areas near the coastline where submarine ground water discharge occurs are the most effective locations for flow observations. (3) Salinity observations are more effective than hydraulic head observations. (4) The importance of flow observations aligned perpendicular to the shoreline varies dramatically depending on distance seaward from the shoreline. Extreme parameter correlation can prohibit unique estimation of permeability parameters such as hydraulic conductivity and flow parameters such as recharge in a density-dependent ground water flow model when using hydraulic head and salinity observations. Adding flow observations perpendicular to the shoreline in areas where ground water is exchanged with the ocean body can reduce the correlation, potentially resulting in unique estimates of these parameter values. Results are expected to be directly applicable to many complex situations, and have implications for model development whether or not formal optimization methods are used in model calibration.  相似文献   

17.
Surface-water and ground-water flow are coupled in the central Everglades, although the remoteness of this system has hindered many previous attempts to quantify interactions between surface water and ground water. We modeled flow through a 43,000 ha basin in the central Everglades called Water Conservation Area 2A. The purpose of the model was to quantify recharge and discharge in the basin's vast interior areas. The presence and distribution of tritium in ground water was the principal constraint on the modeling, based on measurements in 25 research wells ranging in depth from 2 to 37 m. In addition to average characteristics of surface-water flow, the model parameters included depth of the layer of ‘interactive’ ground water that is actively exchanged with surface water, average residence time of interactive ground water, and the associated recharge and discharge fluxes across the wetland ground surface. Results indicated that only a relatively thin (8 m) layer of the 60 m deep surfical aquifer actively exchanges surface water and ground water on a decadal timescale. The calculated storage depth of interactive ground water was 3.1 m after adjustment for the porosity of peat and sandy limestone. Modeling of the tritium data yielded an average residence time of 90 years in interactive ground water, with associated recharge and discharge fluxes equal to 0.01 cm d−1. 3H/3He isotopic ratio measurements (which correct for effects of vertical mixing in the aquifer with deeper, tritium-dead water) were available from several wells, and these indicated an average residence time of 25 years, suggesting that residence time was overestimated using tritium measurements alone. Indeed, both residence time and storage depth would be expected to be overestimated due to vertical mixing. The estimate of recharge and discharge (0.01 cm d−1) that resulted from tritium modeling therefore is still considered reliable, because the ratio of residence time and storage depth (used to calculated recharge and discharge) is much less sensitive to vertical mixing compared with residence time alone. We conclude that a small but potentially significant component of flow through the Everglades is recharged to the aquifer and stored there for years to decades before discharged back to surface water. Long-term storage of water and solutes in the ground-water system beneath the wetlands has implications for restoration of Everglades water quality.  相似文献   

18.
Ground water recharge and flow characterization using multiple isotopes   总被引:2,自引:0,他引:2  
Stable isotopes of delta(18)O, delta(2)H, and (13)C, radiogenic isotopes of (14)C and (3)H, and ground water chemical compositions were used to distinguish ground water, recharge areas, and possible recharge processes in an arid zone, fault-bounded alluvial aquifer. Recharge mainly occurs through exposed stream channel beds as opposed to subsurface inflow along mountain fronts. This recharge distribution pattern may also occur in other fault-bounded aquifers, with important implications for conceptualization of ground water flow systems, development of ground water models, and ground water resource management. Ground water along the mountain front near the basin margins contains low delta(18)O, (14)C (percent modern carbon [pmC]), and (3)H (tritium units [TU]), suggesting older recharge. In addition, water levels lie at greater depths, and basin-bounding faults that locally act as a flow barrier may further reduce subsurface inflow into the aquifer along the mountain front. Chemical differences in ground water composition, attributed to varying aquifer mineralogy and recharge processes, further discriminate the basin-margin and the basin-center water. Direct recharge through the indurated sandstones and mudstones in the basin center is minimal. Modern recharge in the aquifer is mainly through the broad, exposed stream channel beds containing coarse sand and gravel where ground water contains higher delta(18)O, (14)C (pmC), and (3)H (TU). Spatial differences in delta(18)O, (14)C (pmC), and (3)H (TU) and occurrences of extensive mudstones in the basin center suggest sluggish ground water movement, including local compartmentalization of the flow system.  相似文献   

19.
Gu A  Gray F  Eastoe CJ  Norman LM  Duarte O  Long A 《Ground water》2008,46(3):502-509
Sulfate (S and O) isotopes used in conjunction with sulfate concentration provide a tracer for ground water contributions to base flow. They are particularly useful in areas where rock sources of contrasting S isotope character are juxtaposed, where water chemistry or H and O isotopes fail to distinguish water sources, and in arid areas where rain water contributions to base flow are minimal. Sonoita Creek basin in southern Arizona, where evaporite and igneous sources of sulfur are commonly juxtaposed, serves as an example. Base flow in Sonoita Creek is a mixture of three ground water sources: A, basin ground water with sulfate resembling that from Permian evaporite; B, ground water from the Patagonia Mountains; and C, ground water associated with Temporal Gulch. B and C contain sulfate like that of acid rock drainage in the region but differ in sulfate content. Source A contributes 50% to 70%, with the remainder equally divided between B and C during the base flow seasons. The proportion of B generally increases downstream. The proportion of A is greatest under drought conditions.  相似文献   

20.
The combination of detailed multilevel ground water geochemistry samples, a natural-gradient tracer test, minislug tests, and a numerical flow and transport model was used to examine flow through a zero-valent iron permeable reactive barrier (PRB) installed to remove explosives from ground water. After 20 months of operation, the PRB continued to completely remove explosives from the ground water flowing through it. However, the data indicate that a portion of ground water flow was being diverted beneath the PRB. Ground water geochemistry was significantly altered by the PRB, and concentrations of some ions, including sulfate, carbonate, and calcium, were substantially reduced due to precipitation. Field data and numerical model results indicate that, after 20 months of operation, flow through the PRB was reduced to approximately one-third of its expected value.  相似文献   

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

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