首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 250 毫秒
1.
Flooding of the Flint River in July 1994 triggered the collapse of at least 312 sinkholes in the karstic Dougherty Plain at Albany, Georgia. We examined the distribution and morphology of these new sinkholes to evaluate the mode of formation, to characterize early stages of the evolution of sinkhole form, and to estimate the lowering of the surface associated with the development of new sinkholes.Eighty-eight percent of sinkholes occur inside the limits of flooding, especially in areas of sandy overburden, and they often follow joint-controlled linear trends. Sinkhole dimensions are log-normally distributed with median values of circumference = 5.7 m, length = 1.8 m, width = 1.6 m, and depth = 0.7 m; asymmetry (L:W) = 1.2. Cross-sectional forms range from narrow cylinders to large bowls, with many sinkholes having undercut sides.Flooding triggered the formation of sinkholes by saturating and liquefying overburden, which caused soil arches to collapse and flow into cavities in bedrock. The prevalence of sinkholes near the periphery of flooding suggests that drainage and loss of buoyant support as flood waters subsided may also have contributed to failure. A volume ratio index is used to quantify the three dimensional geometric form of sinkholes. Initially, small cylindrical shafts open over a bedrock joint, followed by progressive slumping that leads to widening and increases in volume to a final bowl form. Estimates of the aggregate volume of overburden transported underground in flooded areas range from 7,990 to 11,130 m3. Averaged over flooded areas, this accounts for 0.26 to 0.37 mm/km2 lowering of the surface. Based on a 500 year recurrence interval for the flood event, values for lowering of the surface range from 0.52 to 0.74 mm per 1,000 years. These values are an order of magnitude less than estimates of carbonate dissolution and suggest that transport of overburden underground is limited by triggering events.  相似文献   

2.
As closed topographic depressions, sinkholes effectively divide karst regions into a series of subwatersheds, where swallets within each sinkhole define the downstream end of each subwatershed. Water and sediment are temporarily stored in sinkholes before continuing to underground streams. While pathways and travel times of water in karst terrain are commonly found using tracing techniques, such techniques do not provide information about the movement and travel time of sediment through these systems. This paper reviews the literature regarding sediment storage within sinkholes and presents a case study, the objective of which is to estimate the magnitude and timing of sediment deposition within three sinkholes in an urbanized area in Tennessee, United States of America. Sedimentation rates in sinkholes were estimated based on age determination of cultural artifacts (buried trash and coal fragments) and by 14C dating. For the three study sinkholes, minimum deposition rates (maximum residence times) were 15?cm?yr ?1(6 yrs), 1.5?cm?yr ?1(81?yrs), and 0.2?cm?yr?1(640 yrs). The analysis revealed that sinkholes undergo fill and flush cycles at different timescales caused by changes in upland sediment supply and morphologic changes at swallet openings. Sinkholes that lack swallet openings, or with swallets that have become blocked by debris, are effectively disconnected from master streams with regard to sediment. Sediment stored in sinkholes may be released downstream when swallet conditions change.  相似文献   

3.
Recent karst forms in gypsiferous beds are reported from Mathiesondalen in Central Spitsbergen, where several sinkholes (dolines) and swallow holes have been formed after the deposition of Holocene raised beach sediments. The bedrock of the area is mostly interbedded gypsum/anhydrite and dolomites of the Carboniferous Ebbadalen Formation.  相似文献   

4.
易顺民  卢薇  周心经 《热带地理》2021,41(4):801-811
通过对广州夏茅村岩溶塌陷进行地质勘察,系统探讨了当地的地质环境条件、岩溶塌陷的灾害特征、岩溶发育条件和岩溶塌陷的形成原因。结果表明,区内第四系覆盖土层厚度大且分布不均,土洞大量发育,隐伏可溶岩的富水性好,溶洞发育;同时,覆盖层内的松散岩类孔隙水与下伏灰岩裂隙溶洞水之间的水力联系密切,含水砂层内地下水位的改变会直接引起下伏灰岩裂隙溶洞水的水位波动,构成岩溶塌陷的水动力因素;钻孔桩施工时钻孔揭穿覆盖层内土洞顶部的黏性土层盖板,造成盖板上部的细砂及淤泥流入土洞,是岩溶塌陷的直接触发因素。在此基础上,提出建立岩溶分布区的地下水动态长期监测网,对岩溶塌陷地段进行灌浆处理,是根治夏茅村岩溶塌陷的主要有效途径。  相似文献   

5.
This study presented a quantitative comparison of cockpit and doline karst by examining the numbers and characteristics of typical types of landform entities that are developed in Guilin(Guangxi, China), La Alianza(PR, USA), Avalton(KY, USA), and Oolitic(IN, USA). Five types of landform entities were defined: isolated hill(IH), clustered hills(CHs), isolated sinkhole(IS), clustered sinkholes(CSs), and clustered hills with sinkholes(CHSs). An algorithm was developed to automatically identify these types of landform entities by examining the contour lines on topographic maps of two cockpit karst areas(Guilin and La Alianza) and two doline karst areas(Oolitic and Avalton). Within each specific study area, the CHSs is the least developed type yet with a larger size and higher relief. The IH and IS entities are smaller in size, lower in relief, and outnumber their clustered counterparts. The total numbers of these types of entities are quite different in cockpit and doline karst areas. Doline karst is characterized by more negative(IS and CSs) than positive(IH and IHs) landforms and vice versa for cockpit karst. For example, the Guilin study area has 1192 positive landform entities in total, which occupy 9.81% of the total study area. It has only 622 negative landform entities occupying only 3.91% of the total study area. By contrast, the doline karst in Oolitic has 130 negative while only 10 positive landform entities. The positive and negative landforms in Oolitic occupy 12.68% and 2.61% of the total study area, respectively. Furthermore, average relief and slope of the landform entities are much higher and steeper in the cockpit karst than the doline karst areas. For instance, the average slope of CHs in Alvaton is 3.90 degrees while it is 19.78 degrees in La Alianza. The average relief of CSs is 4.07 m and 34.29 m in Oolitic and Guilin respectively. Such a difference within a specific area or between the cockpit and doline karst may reveal different controls on the development of karst landscape.  相似文献   

6.
This study presented a quantitative comparison of cockpit and doline karst by examining the numbers and characteristics of typical types of landform entities that are developed in Guilin（Guangxi, China）, La Alianza（PR, USA）, Avalton（KY, USA）, and Oolitic（IN, USA）. Five types of landform entities were defined： isolated hill（IH）, clustered hills（CHs）, isolated sinkhole（IS）, clustered sinkholes（CSs）, and clustered hills with sinkholes（CHSs）. An algorithm was developed to automatically identify these types of landform entities by examining the contour lines on topographic maps of two cockpit karst areas（Guilin and La Alianza） and two doline karst areas（Oolitic and Avalton）. Within each specific study area, the CHSs is the least developed type yet with a larger size and higher relief. The IH and IS entities are smaller in size, lower in relief, and outnumber their clustered counterparts. The total numbers of these types of entities are quite different in cockpit and doline karst areas. Doline karst is characterized by more negative（IS and CSs） than positive（IH and IHs） landforms and vice versa for cockpit karst. For example, the Guilin study area has 1192 positive landform entities in total, which occupy 9.81% of the total study area. It has only 622 negative landform entities occupying only 3.91% of the total study area. By contrast, the doline karst in Oolitic has 130 negative while only 10 positive landform entities. The positive and negative landforms in Oolitic occupy 12.68% and 2.61% of the total study area, respectively. Furthermore, average relief and slope of the landform entities are much higher and steeper in the cockpit karst than the doline karst areas. For instance, the average slope of CHs in Alvaton is 3.90 degrees while it is 19.78 degrees in La Alianza. The average relief of CSs is 4.07 m and 34.29 m in Oolitic and Guilin respectively. Such a difference within a specific area or between the cockpit and doline karst may reveal different controls on  相似文献   

7.
A landslide susceptibility map is proposed for the Pays de Herve (E Belgium), where large landslides affect Cretaceous clay outcrop areas. Based on a Bayesian approach, this GIS-supported probabilistic map identifies the areas most susceptible to deep landslides. The database is comprised of the source areas of ten pre-existing landslides (i.e. a sample of 154 grid cells) and of six environmental data layers, namely lithology, proximity to active faults, slope angle and aspect, elevation and distance to the nearest valley-floor. A 30-m-resolution DEM from the Belgian National Geographical Institute is used for the analysis. Owing to the small size of the sample, a special cross-validation procedure of the susceptibility map is performed, which uses in an iterative way each of the landslides to test the predictive power of the map derived from the other landslides. Four different sets of variables are used to produce four susceptibility maps, whose prediction curves are compared. While the prediction rates associated with the models not involving the “proximity to active fault” criterion are comparable to those of the models considering this variable, strong weaknesses inherent in the fault data on which the latter rely suggest that the final susceptibility map should be based on a model that excludes any reference to fault. This highlights the difference between a triggering factor and determining factors, and in the same time broadens the scope of the produced map. A single reactivated slide is also used to test the possibility of predicting future reactivation of existing landslides in the area. Finally, the need for geomorphological control over the mathematical treatment is underlined in order to obtain realistic prediction maps.  相似文献   

8.
闽东北沿海罗源县土地利用空间分布格局的 多尺度分析   总被引：4，自引：0，他引：4  
在不同的空间尺度上,制约土地利用空间分布的影响因子及其影响程度并不相同,因此区域空间分布格局 分析应优先考虑制约空间分布的影响因子识别及其影响程度的尺度依赖性的研究。本文以闽东沿海的罗源县为研 究区域,采样统计方法与GIS 技术,选取20 个候选影响因子,研究了该县主要土地利用空间分布格局的影响因子 及其空间尺度相关性。研究表明模型的解释能力、影响因子及其影响系数均会随研究尺度发生不同程度的变化,回 归模型的解释能力以及主要影响因子的制约程度总体上均随研究尺度增大呈增强趋势。除受坡度、海拔高程等地 形条件的严格制约外,罗源县主要地类均在一定程度上受到人口因素以及若干可达性因素的影响。  相似文献   

9.
In a spatio-temporal data set, identifying spatio-temporal clusters is difficult because of the coupling of time and space and the interference of noise. Previous methods employ either the window scanning technique or the spatio-temporal distance technique to identify spatio-temporal clusters. Although easily implemented, they suffer from the subjectivity in the choice of parameters for classification. In this article, we use the windowed kth nearest (WKN) distance (the geographic distance between an event and its kth geographical nearest neighbour among those events from which to the event the temporal distances are no larger than the half of a specified time window width [TWW]) to differentiate clusters from noise in spatio-temporal data. The windowed nearest neighbour (WNN) method is composed of four steps. The first is to construct a sequence of TWW factors, with which the WKN distances of events can be computed at different temporal scales. Second, the appropriate values of TWW (i.e. the appropriate temporal scales, at which the number of false positives may reach the lowest value when classifying the events) are indicated by the local maximum values of densities of identified clustered events, which are calculated over varying TWW by using the expectation-maximization algorithm. Third, the thresholds of the WKN distance for classification are then derived with the determined TWW. In the fourth step, clustered events identified at the determined TWW are connected into clusters according to their density connectivity in geographic–temporal space. Results of simulated data and a seismic case study showed that the WNN method is efficient in identifying spatio-temporal clusters. The novelty of WNN is that it can not only identify spatio-temporal clusters with arbitrary shapes and different spatio-temporal densities but also significantly reduce the subjectivity in the classification process.  相似文献   

10.
Double Fourier series models of two doline and two cockpit landscapes in northern Jamaica, each 2 × 2 km in area, explained 92% and 90%, and 73% and 58% of the variance in the topographic data. The ten most significant waves accounted for 74% and 76%, and 61% and 58% of the variance in each model. The greater importance of frequency pairs of longer wavelength in the doline karst models suggests that there are fewer horizontally and vertically persistent bedrock fractures controlling topographic development in the doline areas than in the cockpit terrains—a fact confirmed by fracture-trace mapping. Frequency pair orientations and powers indicate that north-west- and northeast-trending fractures exert a major influence on topographic development in the cockpit terrains and that east-trending fractures are relatively more important in the doline areas. Further studies are needed to determine if Fourier models of the doline, cockpit, and tower karst styles differ in a consistent fashion and to find out to what extent these differences are related to the bedding, fracture, and relief characteristics of the karst bedrock. [Key words: Jamaica, karst, morphometry, double Fourier series, dolines, cockpits.]  相似文献   

11.
Landslide hazard assessment, effected by means of geostatistical methods, is based on the analysis of the relationships between landslides and the spatial distributions of some instability factors. Frequently such analyses are based on landslide inventories in which each record represents the entire unstable area and is managed as a single instability landform. In this research, landslide susceptibility is evaluated through the study of a variety of instability landforms: landslides, scarps and areas uphill from crown. The instability factors selected were: bedrock lithology, steepness, topographic wetness index and stream power index. The instability landform densities computed for all the factors, which were arranged in Unique Condition Unit, allowed us to derive a total of three prediction images for each landslide typology. The role of the instability factors and the effects generated by the use of different landforms were analyzed by means of: a) bivariate analysis of the relationships between factors and landslide density; b) predictive power validations of the prediction images, based on a random partition strategy.The test area was the Iato River Basin (North-Western Sicily), whose slopes are moderately involved in flow and rotational slide landslides (219 and 28, respectively). The area is mainly made up of the following complexes: Numidian Flysch clays (19%, 1%), Terravecchia sandy clays (5%, 1%), Terravecchia clayey sands (3%, 0.3%) and San Cipirello marly clays (9%, 0%). The steepness parameter shows the highest landslide density in the [11–19°] class for both the typologies (8%, 1%), even if the density distributions for rotational slides are right-asymmetric and right-shifted. We obtained significant differences in shape when we used different instability landforms. Unlike scarps and areas uphill from crowns, landslide areas produce left-asymmetric and left-shifted density distributions for both the typologies. As far as the topographic wetness index is concerned, much more pronounced differences were detected among the instability landforms of rotational slides. In contrast, the flow landslides produce normal-like density distributions. The latter and the rotational slide landslide areas produce the highest density values in the class [5.5–6.7], despite an abrupt decreasing trend starting from the first class [3.2–4.4], which is generated by the density values of the rotational slide scarps and areas uphill from crowns. The stream power index at the foot of the slopes, which was automatically derived using a GIS-procedure, shows a positive correlation with the landslide densities marked by the maximum classes: [4.8–6.0] for flows, and [6.0–7.2] for rotational slides. The validation procedure results confirmed that the choice of instability landform influences the results of the susceptibility analysis. Furthermore, the validation procedure indicates that: a) the predictive models are generally satisfactory; b) scarps and zones uphill from crown areas are the most diagnostically unstable landforms, for flow and rotational slide landslides respectively.  相似文献   

12.
The surfaces of salt diapirs in the Zagros Mountains are mostly covered by surficial deposits, which significantly affect erosion rates, salt karst evolution, land use and the density of the vegetation cover. Eleven salt diapirs were selected for the study of surficial deposits in order to cover variability in the geology, morphology and climate in a majority of the diapirs in the Zagros Mountains and Persian Gulf Platform. The chemical and mineralogical compositions of 80 selected samples were studied mainly by X-ray powder diffraction and X-ray fluorescence. Changes in salinity along selected vertical profiles were studied together with the halite and gypsum distribution. The subaerial residuum formed from minerals and rock detritus released from the dissolved rock salt is by far the most abundant material on the diapirs. Fluvial sediments derived from this type of residuum are the second most common deposits found, while submarine residuum and marine sediments have only local importance. The mineralogical/chemical composition of surficial deposits varies amongst the three end members: evaporite minerals (gypsum/anhydrite and minor halite), carbonates (dolomite and calcite) and silicates-oxides (mainly quartz, phyllosilicates, and hematite). Based on infiltration tests on different types of surficial deposits, most of the rainwater will infiltrate, while overland flow predominates on rock salt exposures. Recharge concentration and thick accumulations of fine sediment support relatively rich vegetation cover in some places and even enable local agricultural activity. The source material, diapir relief, climatic conditions and vegetation cover were found to be the main factors affecting the development and erosion of surficial deposits. A difference was found in residuum type and landscape morphology between the relatively humid NW part of the studied area and the arid Persian Gulf coast: In the NW, the medium and thick residuum seems to be stable under current climatic conditions. Large sinkholes and blind valleys with sinking streams are common. On other diapirs, the original thick residuum is undergoing erosion and the new morphology is currently represented by salt exposures and badland-like landscapes or by fields of small sinkholes developed in the thin residuum. Models for evolution of the subaerial residuum and the diapir landscape/morphology are described in this paper. While the thick residuum with vegetation has very low erosion rates, the salt exposures and thin residuum are eroded rapidly. During wet periods (e.g. early Holocene), the diapirs rose and salt glaciers expanded as the influx of salt mass was much faster compared to erosion. After the onset of an arid climate, c. 6 ka BP, the rising of the some diapir surfaces decreased or even reversed due to acceleration of erosion thanks to vegetation degradation and changes in the residuum type and thickness.  相似文献   

13.
The influence of land use and land cover on ecological environment is a focus of global change research. The paper chooses an industrial city-Shuicheng in Guizhou Province-as a study area because the karst water quality around the city is deteriorating with land use and land cover change. The natural susceptibility of karst water system is an important factor leading to karst water pollution. But land use and land cover change is also a main factor according to the chemical analysis of karst water quality and land use change. So it is a good way to protect karst water through rational planning and managing of land use and land cover.  相似文献   

14.
土地利用变化对水城盆地岩溶水水质的影响   总被引：1，自引：0，他引：1  
The influence of land use and land cover on ecological environment is a focus of global change research. The paper chooses an industrial city-Shuicheng in Guizhou Province-as a study area because the karst water quality around the city is deteriorating with land use and land cover change.The natural susceptibility of karst water system is an important factor leading to karst water pollution.But land use and land cover change is also a main factor according to the chemical analysis of karst water quality and land use change. So it is a good way to protect karst water through rational planning and managing of land use and land cover.  相似文献   

15.
土地利用变化对水城盆地岩溶水水质的影响   总被引：38，自引：3，他引：35  
贾亚男  袁道先 《地理学报》2003,58(6):831-838
贵州水城盆地是我国裸露型岩溶分布区的一个新兴的工业城市,几十年来社会经济发展迅速,土地利用状况发生了很大的变化,同时对岩溶地下水产生了很大的负面影响。尽管岩溶水文系统具有天然的脆弱性,很容易遭受地表各种污染物质的进入,但水城盆地20年来水质变化,以及水质污染区分布与土地利用分布的一致性,无疑可以认为由于经济发展所造成的各种土地利用矛盾,引起的土地利用的变化,是本地区水质变化的主要原因。从土地利用/土地覆被的角度研究人类活动对岩溶水的影响对合理规划土地利用,保护岩溶地下水,无疑具有重要的指导意义。  相似文献   

16.
This paper presents a new derivative-free search method for finding models of acceptable data fit in a multidimensional parameter space. It falls into the same class of method as simulated annealing and genetic algorithms, which are commonly used for global optimization problems. The objective here is to find an ensemble of models that preferentially sample the good data-fitting regions of parameter space, rather than seeking a single optimal model. (A related paper deals with the quantitative appraisal of the ensemble.)
  The new search algorithm makes use of the geometrical constructs known as Voronoi cells to derive the search in parameter space. These are nearest neighbour regions defined under a suitable distance norm. The algorithm is conceptually simple, requires just two 'tuning parameters', and makes use of only the rank of a data fit criterion rather than the numerical value. In this way all difficulties associated with the scaling of a data misfit function are avoided, and any combination of data fit criteria can be used. It is also shown how Voronoi cells can be used to enhance any existing direct search algorithm, by intermittently replacing the forward modelling calculations with nearest neighbour calculations.
  The new direct search algorithm is illustrated with an application to a synthetic problem involving the inversion of receiver functions for crustal seismic structure. This is known to be a non-linear problem, where linearized inversion techniques suffer from a strong dependence on the starting solution. It is shown that the new algorithm produces a sophisticated type of 'self-adaptive' search behaviour, which to our knowledge has not been demonstrated in any previous technique of this kind.  相似文献   

17.
Sinkholes up to 50 m diameter and 15 m deep are developed in the lateritised Mullaman (Lower Cretaceous) quartzitic and argillaceous beds of the western Sturt Plateau, in the monsoonal north of the Northern Territory, Australia. Some sinkholes are degraded and of some antiquity: though their age is not precisely known, they formed during the late Cenozoic following the uplift and desiccation of the latertte. Others are fresh and a few have formed in living memory. Some sinkholes may be due to collapse of the Mullaman beds into voids in the underlying Cambrian limestone, but most have evolved wholly within the highly siliceous Mullaman beds, which have been dissolved or altered, creating either voids or compartments filled with uncosolidated material. Subsequent water-table lowering and collapse of the overlying beds into the low-density zones caused the formation of sinkholes.Fractures have significantly influenced the location, pattern and detailed plan shape of sinkholes. Sinkholes also occur within and at the margins of old drainage lines. The water table stood higher in the late Pleistocene, bringing alkaline groundwaters into contact with the country rock. A late Cenozoic fall in regional water table caused slacking or disaggregation. Litter derived from a Eucalyptus-dominated woodland facilitated iron mobilisation. Biogenic agents may have accelerated the weathering of the silica, silicates and iron oxides of which the country rock and laterite are largely composed.  相似文献   

18.
Geomorphological information can be combined with decision-support tools to assess landslide hazard and risk. A heuristic model was applied to a rural municipality in eastern Cuba. The study is based on a terrain mapping units (TMU) map, generated at 1:50,000 scale by interpretation of aerial photos, satellite images and field data. Information describing 603 terrain units was collected in a database. Landslide areas were mapped in detail to classify the different failure types and parts. Three major landslide regions are recognized in the study area: coastal hills with rockfalls, shallow debris flows and old rotational rockslides denudational slopes in limestone, with very large deep-seated rockslides related to tectonic activity and the Sierra de Caujerí scarp, with large rockslides. The Caujerí scarp presents the highest hazard, with recent landslides and various signs of active processes. The different landforms and the causative factors for landslides were analyzed and used to develop the heuristic model. The model is based on weights assigned by expert judgment and organized in a number of components such as slope angle, internal relief, slope shape, geological formation, active faults, distance to drainage, distance to springs, geomorphological subunits and existing landslide zones. From these variables a hierarchical heuristic model was applied in which three levels of weights were designed for classes, variables, and criteria. The model combines all weights into a single hazard value for each pixel of the landslide hazard map. The hazard map was then divided by two scales, one with three classes for disaster managers and one with 10 detailed hazard classes for technical staff. The range of weight values and the number of existing landslides is registered for each class. The resulting increasing landslide density with higher hazard classes indicates that the output map is reliable. The landslide hazard map was used in combination with existing information on buildings and infrastructure to prepare a qualitative risk map. The complete lack of historical landslide information and geotechnical data precludes the development of quantitative deterministic or probabilistic models.  相似文献   

19.
In August 1983, 18 samples were collected in Norwegian karst environments: seepage waters, surface flow upstream sinkholes or caves, allogenic and authigenie karst springs. Hydrochemical data (ion contents and computed variables for CO-CaCO system) are described by means of a factor analysis, the principal components analysis. Four factors, responsible for about 90 % of the total variance, are extracted and interpreted. Five clusters of samples are then identified:

– Authigenic karst waters (seepage and only carbonate aquifer);

– Waters of very short residence time (surface or karst flow);

– CO-rich waters (seepage and possibly organic-rich surface waters), CO2 part is discussed in more detail;

– CO-poor waters (surface waters above tree-line);

– a central group, combining 2 or more of these characters.  相似文献   

20.
以湖北省鄂州程潮铁矿和黄石大冶铁矿为例,利用GIS空间分析功能对研究区数据进行提取分级、赋值统计及归一化等处理,构建了包括高程、坡度、地层、地下开采点的分布密度、相距最近地下开采点的距离、开采厚度与深度比值、蚀变接触带缓冲区、地下水深度以及地表地物类型的矿区采空塌陷易发性评价指标数据集;借助IDL语言调用Matlab神经网络工具箱,将研究区2011和2012年的指标数据集作为输入数据,塌陷易发性作为期望输出,建立基于BP神经网络的矿区采空塌陷易发性预测模型;通过选取并优化训练样本,实现对2013年矿山塌陷易发性的预测。结果表明,高易发区及以上的区域包含89.91%的采空塌陷,随着易发等级的提高,采空塌陷面积占易发等级面积比也随之增大;采空塌陷的分布具有明显的地带性,高易发区基本沿着岩体与围岩的接触带分布。模型解决了塌陷预测中的非线性映射问题,预测结果与实际调查情况基本吻合。BP神经网络模型与GIS技术相结合预测矿区采空塌陷的易发性具有可行性。  相似文献   

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

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