首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
This study proposes a solution to the problem of maintaining the performance and sustainability of district-scale, cooling-dominated ground coupled heat pump (GCHP) systems. These systems tend to overheat because heat dissipates slowly in relation to the size of the borefields. To demonstrate this problem, a 2000-borehole field is considered at a district-scale GCHP system in the Upper Midwest, US. The borefield’s ground and fluid temperature responses to its design heating and cooling loads are simulated using computational fluid dynamics implemented by applying the finite volume method. The ground temperature is predicted by applying the thermal loads uniformly over the borefield and simulating heat dissipation to the surrounding geology through conduction coupled with advection due to groundwater flow. The results show that a significant energy imbalance will develop in the ground after the first few years of GCHP operation, even with high rates of groundwater flow. The model presented in this study predicts that the temperature at the center of the borefield will reach 18 °C after 5 years and approximately 50 °C after 20 years of operation in the absence of any mitigation strategies. The fluid temperature in the boreholes is then simulated using a single borehole model to estimate the heat pump coefficient of performance, which decreases as the modeled system heats up. To balance the energy inputs/outputs to the ground—thus maintaining the system’s performance—an operating scheme utilizing cold-water circulation during the winter is evaluated. Under the simulated conditions, this mitigation strategy carries the excess energy out of the borefield. Therefore, the proposed mitigation strategy may be a viable measure to sustaining the operating efficiency of cooling-dominated, district-scale borefields in climates with cold winters.  相似文献   

2.
The numerical simulation of collapsible settlement in loess soil subjected to rising ground water table is presented in this paper. A coupled hydro-mechanical model is proposed. Comparisons between the results of numerical simulations and those of oedometer and in situ water immersion field test in Lanzhou, northwest China, reveal good agreement, which validates the proposed model formulation. Factors that influence the ground settlement of loess including initial elevation of ground water table, rising water height and velocity are then evaluated. The results of the analyses reveal that the most critical situation of largest possible ground settlement due to ground water rising in loess involves initial water table elevation of 10 m and rising water velocity of 0.5 m/year. Two upper bound lines of predicted maximum possible ground settlement are proposed to facilitate a preliminary quick evaluation of ground settlement due to rising water under different water table scenarios in loess.  相似文献   

3.
This article presents a sensitivity analysis investigating the impact of using high-resolution site conditions databases in portfolio earthquake loss estimation. This article also estimates the effects of variability in the site condition databases on probabilistic earthquake loss ratios and their geographical pattern with respect to structural characteristics of different building types. To perform the earthquake loss estimation here, the OpenQuake software developed by Global Earthquake Model is implemented in Clemson University’s supercomputer. The probabilistic event-based risk analysis is employed considering several notional portfolios of different building types in the San Francisco area as the inventory exposure. This analysis produces the stochastic event sets worth for 10,000 years including almost 8000 synthetically simulated earthquakes. Then, the ground motion prediction equations are used to calculate the ground motion per event and incorporate the effect of five site conditions, on amplifying or de-amplifying the ground motions on notional building exposure locations. Notional buildings are used to account for various building characteristics in conformance with the building taxonomy represented in HAZUS software. The HAZUS damage functions are applied to model the vulnerability of various structural types of buildings. Finally, the 50-year average mean loss and probabilistic loss for multiple values for probability of exceedance (2, 10, 20, and 40%) in 50 years are calculated, and the impact of different site condition databases on portfolio loss ratios is investigated for different structural types and heights of buildings. The results show the aggregated and geographical variation of loss and loss ratio throughout the region for various site conditions. Comparing the aggregated loss and loss ratio, while considering different databases, represents normalized differences that are limited to 6% for all building taxonomy with various heights and for all PoEs. However, site-specific loss ratio errors are significantly greater and in some cases are more than 20%.  相似文献   

4.
Thermokarst lake is a significant indicator of permafrost degradation. However, the evaluation of thermokarst lake development is very difficult through physical mechanism analysis because the development is influenced by many factors. In the present study, the factors influencing thermokarst lake development were investigated from the perspective of frozen soil and the geographical environment. The influences of six factors on thermokarst lake development in Beiluhe basin, located in the hinterland of the Qinghai–Tibet Plateau, were analyzed: permafrost type, ground temperature, vegetation type, soil type, hydrogeological type, and slope. Sensitivity coefficients were calculated for these factors using statistical methods. The results show that thermokarst lake development was influenced by the analyzed factors as follows: permafrost > soil type > ground temperature > hydrogeological type. Furthermore, 80.1% of the study area was highly sensitive to thermokarst lake development. Overall, thermokarst lake development in the Beiluhe basin was thoroughly evaluated based on sensitivity factors and an established statistical evaluation method. The method detailed in this paper provides a reference for evaluating the likelihood and severity of thermokarst lake development.  相似文献   

5.
A. Joshi  Sandeep  Kamal 《Natural Hazards》2014,71(1):587-609
Modification in the semi-empirical technique for the simulation of strong ground motion has been introduced to incorporate the strong motion generation areas (SMGA) in the modeled rupture plane. Strong motion generation areas identified within the rupture plane of the Tohoku earthquake of March 11, 2011 (M w = 9.0), have been modeled using this modified technique. Two different source models having four and five SMGAs, respectively, are considered for modeling purpose. Strong motion records using modified semi-empirical technique have been simulated at two near-field stations located at epicentral distance of 137 and 140 km, respectively, using two different source models. Comparison of the observed and simulated acceleration waveforms is made in terms of root mean square error (RMSE) at both stations. Minimum root mean square error of the waveform comparison has been obtained at both the stations for source model having five SMGAs. Simulations from same rupture model have been made at other four stations lying at epicentral distance between 154 and 249 km. Comparison of observed and simulated records has been made in terms of RMSE in acceleration records, velocity records and response spectra at each six station. Simulations have been made at six other stations to obtain distribution of peak ground acceleration and peak ground velocity with hypocentral distance. Peak ground acceleration and velocity from simulated and observed records are compared at twelve stations surrounding the source of Tohoku earthquake. Comparison of waveforms and parameters extracted from observed and simulated strong motion records confirms the efficacy of the developed modified technique to model earthquake characterized by SMGAs.  相似文献   

6.
Two polynomial constitutive equations of nonlinear stress–strain relationship were used to construct two nonlinear 1-D wave equations with external pressure (source term) applied. Nonlinear model 1 was a concave downward curve, and nonlinear model 2 was a concave upward curve. The time-dependent stress and strain of a 300-m length were calculated. The computation for nonlinear model 1 terminated at time t?=?5.0 s. The evaluated stress versus position at different times was mainly at position x?=?0 m. The stress versus position had a concave point at position x?=?0.3 m. Between x?=?0 m and x?=?0.3 m, the graph of stress versus position showed a small convex upward curve. Moreover, at the two sides of x?=?0.3 m, the strain position had distinctively different slopes. The distinctive difference in the slope of strain at position x?=?0.3 m can therefore be used to estimate the rupture position of a rock. The terminal evaluation time for nonlinear model 2 was at t?=?1.55 s. The stress versus position and the strain versus position to time change were within the 0–3 m interval from the pressure end. Time increase produced the phenomenon of stress and strain solitons. These stress and strain solitons moved forward and increased in peak value with time. During the compression process, the stress and strain soliton resulted in instability which rendered the rock situation more easily broken. The position of the rock rupture might have occurred away from the pressure side.  相似文献   

7.
The bacterium Pseudomonas aeruginosa BCH decolorized and degraded the sulphonated azo dye Remazol Orange in plain distilled water. The effects of different parameters, i.e. pH, temperature and cell mass concentration on the biodegradation of dye in aqueous phase was evaluated using response surface methodology. Optimization was carried out using three-level Box–Behnken design. Predicted values were found to be in good agreement with experimental values (R 2 0.9997 and pred R 2 0.9984), which indicated suitability of the employed model and the success of response surface methodology. Optimum dye decolorization was maximized and the favourable conditions were pH 7.43, temperature 29.39 °C and cell mass concentration 2.88 g l?1, which resulted in 96.01 % decolorization within 5 h. It was validated from the predicted response (97.37 %). According to the analysis of variance results, the proposed model can be used to navigate the design space. 3D plot analysis disclosed the significant interaction between all three tested factors on decolorization process. The combinations of the three variables predicted during response surface methodology were confirmed through confirmatory experiments. Observations indicated that higher cell mass accelerated the decolorization process. Degradation of the dye was verified through high performance liquid chromatography analysis. Phytotoxicity studies carried out with dye and dye metabolites using Phaseolus mungo, Triticum aestivum and Sorghum vulgare indicated the detoxification of dye.  相似文献   

8.
The specific energy (SE) is the most important parameter to estimate the energy consumption in tunnel boring machines (TBMs). It is defined as the amount of required energy to excavate a unit volume of rock mass which used to predict the performance of TBMs. Several models are used to estimate the SE based on different parameters such as the rock mass properties, disc cutter dimensions and cutting geometry. The aim of this work is to propose new relations between the SE and the strain energy of rock mass (W) using the geological mappings of rock mass and TBM operational parameters from Amir-Kabir Water Transferring Tunnel of Iran. W is an appropriate criterion to estimate SE because it is a function of different parameters such as rock mass behavior, pre and post failure properties and peak and residual strains. In this study, to increase the correlation coefficient of relation between the mentioned parameters, the rock mass is classified in two methods, in the first method according to the geological strength index (GSI) all data is classified in three classes such as weak, fair and good and in the second method using the drop to deformation modulus ratio (η) the classification of data is performed in three classes such as η < 0.05, 0.05 ≤ η < 10 and η ≥ 10. The results show that there are direct relations between both parameters. It is suggested to estimate SE in all rock mass classes using the proposed relations based on GSI classification.  相似文献   

9.
Domain configuration and several physical parameterization settings such as planetary boundary layer, cumulus convection, and ocean–atmosphere surface flux parameterizations can play significant roles in numerical prediction of tropical cyclones. The present study focuses to improve the prediction of the TC Gonu by investigating the sensitivity of simulations to mentioned configurations with the Advanced Hurricane WRF model. The experiments for domain design sensitivity with 27 km resolution has been shown moving the domains towards the east improve the results, due to better account for the large-scale process. The fixed and movable nests on a 9-km grid were considered separately within the coarse domain and their results showed that despite salient improvement in simulated intensity, an accuracy reduction in simulated track was observed. Increasing horizontal resolution to 3 km incredibly reduced the simulated intensity accuracy when compared to 27 km resolution. Thereafter, different initial conditions were experimented and the results have shown that the cyclone of 1000 hPa sea level pressure is the best simulation initial condition in predicting the track and intensity for cyclone Gonu. The sensitivity of simulations to ocean–atmosphere surface-flux parameterizations on a 9-km grid showed the combination of ‘Donelan scheme’ for momentum exchanges along with ‘Large and Pond scheme’ for heat and moisture exchanges provide the best prediction for cyclone Gonu intensity. The combination of YSU and MYJ PBL scheme with KF convection for prediction of track and the combination of YSU PBL scheme with KF convection for prediction of intensity are found to have better performance than the other combinations. These 22 sensitivity experiments also implicitly lead us to the conclusion that each particular forecast aspect of TC (e.g., track, intensity, etc.) will require its own special design.  相似文献   

10.
At 6:10 p.m. on September 23, 1991, a catastrophic rock avalanche occurred in Zhaotong, Yunnan, southwestern China. Over 216 people were killed when the Touzhai village was overwhelmed directly in the path of the landslide. The landslide involved the failure of about 12 Mm3 of jointed basaltic rock mass from the source area. The displaced materials ran out a horizontal distance of 3650 m over a vertical distance of 960 m, equivalent to a Fahrböschung of 14.7°, and covered an area of 1.38 km2. To provide information for hazard zonation of similar type of potential landslides in the same area, we used a dynamic model (DAN-W) with three alternative rheological models to simulate the runout behaviour of the displaced landslide materials and found that a combination of the frictional model and Voellmy model could provide the best performance in simulating this landslide. The simulated results indicated that the duration of the movement is estimated at about 175 s for a mean velocity 21 m/s.  相似文献   

11.
Aquifer thermal energy storage in urban and industrial areas can lead to an increase in subsurface temperature to 70 °C and more. Besides its impacts on mineral and sorption equilibria and chemical reaction kinetics in an aquifer, temperature sensitively influences microbial activity and thus redox processes, such as sulphate reduction. Microorganism species can only operate within limited temperature ranges and their adaptability to temperature is a crucial point for the assessment of the environmental consequences of subsurface heat storage. Column experiments with aquifer sediment and tap water at 10, 25, 40, and 70 °C showed that under the constant addition of acetate sulphate reduction could be initiated after 26–63 pore volumes exchanged at all temperatures. Fastest initiation of sulphate reduction with the highest reduction rates was found at 40 °C. Maximum rate constants during experimental run-time were 0.56 h?1 at 40 °C and 0.33, 0.36, and 0.25 h?1 at 10 and, 25, and 70 °C, respectively. Hence, microbial activity was enhanced by a temperature increase to 40 °C but was significantly lowered at 70 °C. At 25 °C methane was found in solution, indicating the presence of fermenting organisms; at 10, 40, and 70 °C no methane production was observed. It could be shown that redox processes in an aquifer generally can adapt to temperatures significantly higher than in situ temperature and that the efficiency of the reduction process can be enhanced by temperature increase to a certain limit. Enhancement of sulphate reduction in an aquifer due to temperature increase could also allow enhanced degradation of organic ground water contaminants such as BTEX, where sulphate is an important electron acceptor.  相似文献   

12.
To study the displacement caused by underground mining, the displacement field in the strata overlying the mine was simulated in nine similar-material models with coal seam dip angles of 0°–80°. Digital close-range photogrammetry was employed to observe the displacement in these models and produce displacement vector diagrams of observation points set on the model surfaces. The movement boundary of the rock mass was extracted and determined. According to the displacement vectors’ direction, the displacement field within the movement boundary was divided into five zones: two zones with displacement vectors pointing towards the goaf center, two zones with vectors pointing towards the coal pillar, and one zone with displacement vectors pointing vertically downward. A symmetry index was defined to analyze the symmetry of the displacement field, and the results show that as the coal seam dip angle increased from 0° to 80°, the displacement field in the strata changed according to the following pattern: symmetric → asymmetric → nearly symmetric, with a transition angle of about 45°. The percentage of the area of the above-defined five zones relative to the displacement field’s total area also changed, showing a regular pattern. This study also revealed that the movement boundary inside the overlaying rock layers is not a straight line, but an S-shaped curve. These findings can enhance our understanding of the internal mechanism of ground subsidence and contribute to more effective prediction of the deformation occurring inside a given rock mass.  相似文献   

13.
The Sanxicun landslide occurred on July 10, 2013, in Sanxicun Village, which is located in Dujiangyan City, Sichuan Province, China. It travelled up to 1200 m, destroyed 11 houses, and killed 166 people in the village. To explain how this landslide could travel such a large distance and cause such serious damage, this study used a thermo-poro-elastic approach coupled with the Savage–Hutter model to simulate the dynamic process of the Sanxicun landslide. The simulated results were compared with the actual results, as well as those of other researchers. It showed that the simulated results for the landslide profile and mass accumulation scope were basically consistent with the actual results. The simulated landslide runout was 1242 m, which was quite close to the actual value. The simulated maximum mass accumulation thickness was 16.4 m. The maximum velocity was 32.6 m/s, which was between those calculated by Yin et al. (J Eng Geol 22(2):309–318, 2014), Yin et al. (Landslide 13:9–23, 2016), and the various trends were found to be consistent. The temperature change and the pore water pressure evolution in the shear zone during sliding are also obtained by simulation. This study had recreated the Sanxicun landslide motion process from the view of thermo-poro-elastic coupling within the shear zone.  相似文献   

14.
Since the 1960s, an early explosive activity in the Mont-Dore Massif is associated with a major pyroclastic rhyolitic eruption (5–7 km3) known as the “Grande Nappe” (GN). This event, linked to the formation of a 6-km-diameter cryptic caldera named “Haute Dordogne”, was before our investigation dated by 40Ar/39Ar at 3.07 ± 0.04 Ma. Our new single-crystal laser fusion 40Ar/39Ar dates obtained on two outcrops of the GN (Rochefort-Montagne and Ludières) questioned several hypotheses made concerning this “landmark” event of the Mont-Dore Massif history. We demonstrate that: (1) the GN rhyolitic eruption has occurred much later than previously estimated (i.e. 2.77 ± 0.02–0.07 Ma full external uncertainties); (2) the correlation made between the Vendeix rhyolitic complexes (intra-caldera position) dated back to 2.74 ± 0.04 Ma and the GN is proposed; (3) xenocryst contamination could be very high (i.e. 70% for the Rochefort-Montagne GN outcrop) and explains the noticeable older age obtained previously; (4) a link between the GN eruption and the formation of a caldera is questionable; the hypothesis of a northward-oriented blast channeled eastward toward the paleo-Allier River is thus proposed.  相似文献   

15.
Evaporation capacity is an important factor that cannot be ignored when judging whether extreme precipitation events will produce groundwater recharge. The evaporation layer’s role in groundwater recharge was evaluated using a lysimeter simulation experiment in the desert area of Dunhuang, in the western part of the Hexi Corridor in northwestern China’s Gansu Province. The annual precipitation in the study area is extremely low, averaging 38.87 mm during the 60-year study period, and daily pan evaporation amounts to 2,486 mm. Three simulated precipitation regimes (normal, 10 mm; ordinary annual maximum, 21 mm; and extreme, 31 mm) were used in the lysimeter simulation to allow monitoring of water movement and weighing to detect evaporative losses. The differences in soil-water content to a depth of 50 cm in the soil profile significantly affected rainfall infiltration during the initial stages of rainfall events. It was found that the presence of a dry 50-cm-deep sand layer was the key factor for “potential recharge” after the three rainfall events. Daily precipitation events less than 20 mm did not produce groundwater recharge because of the barrier effect created by the dry sand. Infiltration totaled 0.68 mm and penetrated to a depth below 50 cm with 31 mm of rainfall, representing potential recharge equivalent to 1.7 % of the rainfall. This suggests that only extreme precipitation events offer the possibility of recharge of groundwater in this extremely arid area.  相似文献   

16.
The occurrence of rockfall incidents on the transportation network may cause injuries, and even casualties, as well as severe damage to infrastructure such as dwellings, railways, road corridors, etc. Passive protective measures (i.e., rockfall barriers, wire nets, etc.) are mainly deployed by operators of ground transport networks to minimize the impact of detrimental effects on these networks. In conjunction with these passive measures, active rockfall monitoring should ideally include the magnitude of each rockfall, its initial and final position, and the triggering mechanism that might have caused its detachment from the slope. In this work, the operational principle of a low-cost rockfall monitoring and alerting system is being presented. The system integrates measurements from a multi-channel seismograph and commercial cameras as the primary equipment for event detection. A series of algorithms analyze these measurements independently in order to reduce alarms originated by surrounding noise and sources other than rockfall events. The detection methodology employs two different sets of algorithms: Time–frequency analyses of the rockfall event’s seismic signature are performed using moving window pattern recognition algorithms, whereas image processing techniques are utilized to deliver object detection and localization. Training and validation of the proposed approach was performed through field tests that involved manually induced rockfall events and recording of sources (i.e., passing car, walking people) that may cause a false alarm. These validation tests revealed that the seismic monitoring algorithms produce a 4.17 % false alarm rate with an accuracy of 93 %. Finally, the results of a 34-day operational monitoring period are presented and the ability of the imaging system to identify and exclude false alarms is discussed. The entire processing cycle is 10–15 s. Thus, it can be considered as a near real-time system for early warning of rockfall events.  相似文献   

17.
Structural remains of the site refer to clear evidence of large historical earthquakes that strongly affected the site during the last 2,000 years. In this work, we present a study of dynamic behavior of Jabal Al-Qalaa site as well as the dynamic behavior of some archeological remaining structures on the site using the free oscillation of surface geology and structures. In this study, seven ground records were obtained along approximately east–west profile, crossing the top of the hill zone. Structural records were obtained at three selected structures following administrative requirements. All structural records were conducted at the top level of each structure using three components seismometer of 2 Hz free oscillation. Horizontal sensors were oriented with respect to the longitudinal and transverse direction of structural horizontal projection. Used orientations are essential for obtaining both structural fundamental modes in the longitudinal and transverse directions. Moreover, damping factors for each horizontal fundamental mode were calculated based on the obtained FFT spectrums of each orientation using half-width band method. Analysis shows that all inspected structures in this study are of high frequency. Interesting results of records obtained at the Umayyad Alcazar, relatively low damping factor as well as sharp fundamental modes in its both longitudinal and transverse directions are recorded at this well-preserved old historical structure. Free-field records were conducted along approximately E-W profile crossing the top of the hill. The purpose of this methodology is to identify the variation of H/V amplification along with, and to find out the topographical effect of this convex geological feature. Obtained H/V ground resonance dominant modes as well as the site amplification reflect well the topographical effect of the citadel hill zone. Amplification gradually increases to reach 4.5 at the top of the hill.  相似文献   

18.
In conventional seismic hazard analysis, uniform distribution over area and magnitude range is assumed for the evaluation of source seismicity which is not able to capture peculiar characteristic of near-fault ground motion well. For near-field hazard analysis, two important factors need to be considered: (1) rupture directivity effects and (2) occurrence of scenario characteristic ruptures in the nearby sources. This study proposed a simple framework to consider these two effects by modifying the predictions from the conventional ground motion model based on pulse occurrence probability and adjustment of the magnitude frequency distribution to account for the rupture characteristic of the fault. The results of proposed approach are compared with those of deterministic and probabilistic seismic hazard analyses. The results indicate that characteristic earthquake and directivity consideration both have significant effects on seismic hazard analysis estimates. The implemented approach leads to results close to deterministic seismic hazard analysis in the short period ranges (T < 1.0 s) and follows probabilistic seismic hazard analysis results in the long period ranges (T > 1.0 s). Finally, seismic hazard maps based on the proposed method could be developed and compared with other methods.  相似文献   

19.
Noblesse multi-collector noble gas mass spectrometer is specially designed for multi-collection of Ar isotopes with different beam sizes, especially for small ion beams, precisely, and hence is perfectly suitable for 40Ar/39Ar geochronology. We have analyzed widely used sanidine, muscovite, and biotite standards with recommended ages of ~ 1.2–133 Ma, with the aim to assess the reliability of Noblesse for 40Ar/39Ar dating. An ESI MIR10 30W CO2 laser was used for total fusion or incremental heating samples. Extracted gases were routinely purified by four SAES NP10 getters (one at ~ 400 °C and others at room temperature). A GP50 getter and a metal cold finger cooled by liquid N (? 196 °C) were also attached for additional purification if necessary. The Ar isotopes were then measured by Noblesse using Faraday or multiplier according to the signal intensities. Over a period of 1.5 months 337 air calibrations produced a weighted mean 40Ar/36Ar of 296.50 ± 0.08 (2σ, MSWD = 4.77). Fish Canyon sanidine is used to calculate J-values, which show good linear relationship with position in irradiation. The age of four mineral standards (Alder Creek sanidine, Brione muscovite, Yabachi sanidine, and Fangshan biotite) are within error of the accepted ages. Five Alder Creek sanidine aliquots yielded an age range of 1.174–1.181 ± 0.013 Ma (2σ) which broadly overlaps the established age of the standard and the uncertainty approaches those of the foremost Ar/Ar laboratories in the world. The weighted mean ages of four Brione muscovite aliquots (18.75 ± 0.16 Ma, 2σ), five Yabachi sanidine aliquots (29.50 ± 0.19 Ma, 2σ), and three Fangshan biotite aliquots (133.0 ± 0.76 Ma, 2σ) are consistent with the recommended values of these standards, and the uncertainties are typical of modern Ar/Ar laboratories world-wide.  相似文献   

20.
In order to solve the bulging deformation and fracture at the top of widely used gravel piles in treating ground consolidation, a new, optimized composite foundation form was proposed. The composite foundation was constructed using discrete materials and concrete piles. Additionally, various parameters of this new composite foundation were analyzed, including foundation forms, construction technologies, bearing mechanism and failure mode. By applying cavity expansion theory, the Vesic cavity spreading pressure of the discrete material-concrete pile is solved as a polar axis symmetric problem on the basis of Mohr–Coulomb yield criterion. Then the computing formula for the ultimate bearing capacity of the discrete materials-concrete pile is elicited when the internal friction angle of soil in the piles is φ = 0 and φ ≠ 0. Finally, the ultimate bearing capacity value of the composite foundation is acquired through analytic calculation and numerical simulation. Finally, it is found that the calculation result is 14.4% lower than that of the simulated result, which is within the acceptable accuracy range and therefore proves the accuracy of the analytic calculation method for bearing capacity of the new composite foundation.  相似文献   

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

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