首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Tulsishyam thermal springs are located in the Saurashtra region of Gujarat, India with discharge temperatures varying from 39 to 42 °C. The pH of these thermal springs varies from 7.1 to 7.4, indicating neutral character. Though these thermal springs propagate through the near surface layer of Deccan basalt, detailed geochemical analysis of the thermal waters using Piper diagram suggests that the water is interacting with the granitic basement rock. Silica and cation geothermometry estimates have reservoir temperature in the range of 138 to 207 °C categorizing it into a low to moderate enthalpy geothermal system. Furthermore, the area has high heat flow values of 53–90 mW/m2 because of shallow Moho depth. The prevailing conditions suggest that the geothermal energy can potentially be exploited through an enhanced geothermal system (EGS). The study also indicates different mineral phases that may precipitate out of water during exploitation of geothermal energy and it should be taken into account while designing an EGS for the area.  相似文献   

2.
This study focuses on the hydrochemical characteristics of 47 water samples collected from thermal and cold springs that emerge from the Hammam Righa geothermal field, located in north-central Algeria. The aquifer that feeds these springs is mainly situated in the deeply fractured Jurassic limestone and dolomite of the Zaccar Mount. Measured discharge temperatures of the cold waters range from 16.0 to 26.5 °C and the hot waters from 32.1 to 68.2 °C. All waters exhibited a near-neutral pH of 6.0–7.6. The thermal waters had a high total dissolved solids (TDS) content of up to 2527 mg/l, while the TDS for cold waters was 659.0–852.0 mg/l. Chemical analyses suggest that two main types of water exist: hot waters in the upflow area of the Ca–Na–SO4 type (Hammam Righa) and cold waters in the recharge zone of the Ca–Na–HCO3 type (Zaccar Mount). Reservoir temperatures were estimated using silica geothermometers and fluid/mineral equilibria at 78, 92, and 95 °C for HR4, HR2, and HR1, respectively. Stable isotopic analyses of the δ18O and δD composition of the waters suggest that the thermal waters of Hammam Righa are of meteoric origin. We conclude that meteoric recharge infiltrates through the fractured dolomitic limestones of the Zaccar Mount and is conductively heated at a depth of 2.1–2.2 km. The hot waters then interact at depth with Triassic evaporites located in the hydrothermal conduit (fault), giving rise to the Ca–Na–SO4 water type. As they ascend to the surface, the thermal waters mix with shallower Mg-rich groundwater, resulting in waters that plot in the immature water field in the Na–K–Mg diagram. The mixing trend between cold groundwaters from the recharge zone area (Zaccar Mount) and hot waters in the upflow area (Hammam Righa) is apparent via a chloride-enthalpy diagram that shows a mixing ratio of 22.6 < R < 29.2 %. We summarize these results with a geothermal conceptual model of the Hammam Righa geothermal field.  相似文献   

3.
Analysis of geothermal and hydrogeologic characteristics of Paleozoic interior basins of Brazil has identified an association between the geographic distribution of thermal springs and areas of occurrences of groundwater flow. Specifically, thermal springs are found to be absent in regions inferred to have lateral flows of groundwater. This trend is evident in the basins of the Amazon region, in the central parts of the Parnaíba basin and in the west-central parts of the Paraná basin. Model studies help to elucidate mutual exclusion of regions of thermal springs and sub-horizontal flows of groundwater. Numerical simulations indicate that groundwater flows with velocities ≥ 1 cm/year are capable of masking the occurrence of thermal anomalies. Also, down flow through distributed recharge zones can lead to development of large zones of relatively low temperature. The observational data sets of temperature gradients and Peclet numbers have been employed outlining advection-convection domains of subsurface strata in the sedimentary basins of the Amazon region, Parnaíba and Paraná. Results obtained indicate that thermal buoyancy forces are incapable of overcoming advective flows in basins of the Amazon region. Similar conditions are also found to prevail in the central parts of the Parnaíba and Paraná basins.  相似文献   

4.
The Gulf of Suez is characterized by the presence of many hot springs and deep thermal wells scattered around its coastal areas. So it is considered one of the promised geothermal areas in Egypt. In this study, the main emphasis is to investigate the geothermal potential around the Gulf of Suez using the available logging and geothermometer datasets. The temperature profiles and well logging data of some hot springs and deep wells around or within the coastal area of the Gulf of Suez are used in this study. The temperature profiles are analyzed and some important thermophysical properties are estimated (geothermal gradient, thermal conductivity, heat flow, and specific heat capacity). Such analysis revealed that a medium to high geothermal gradient (22.0–30°C/Km) is given for the Gulf of Suez as a whole, with some spots of much higher gradient in the order of 35.0–44°C/Km (Ras Fanar and Hammam Faraun areas). The compiled thermal plots show that the thick evaporites and rock salt lithology, which is a major constituent in this area, attain the highest thermal conductivity (>3.10 W/m/K) and heat flow (>90 mW/m2) and the lowest specific heat capacity (<0.30 J/kg/K). The available gamma ray and the natural gamma ray spectroscopy logs are used to conduct a radioactive-based heat generation study using the characteristic radioactive nature of some elements like; 238U, 235U, 232Th, and of the isotope of 40K. A good linearity is observed between the heat production (A in microwatt per cubic meter) and the gamma ray (API) along a wide range of datasets (0–150 API) in all wells. The heat production factor increases in the carbonate lithology (up to 3.20?μW/m3) and is proportional to the shale volume. A geothermometer-based study is used to estimate the subsurface formation temperature and heat flow from the geochemical analysis of some water samples collected from the studied hot springs. The estimated thermal parameters are in harmony with the regional thermal regime concluded form logging data. A thermal basin growth study, in relation to the clay diagenesis is conducted concerning the thermal effects that take place with depth giving rise to another clay mineral (illite). Furthermore, a number of 2D thermal–burial history diagrams are constructed for the complied sections of some of the studied areas to show the vertical distribution of the estimated petrothermal properties. A reserve evaluation study is carried out to estimate the economic geothermal capacity of these hot springs to be used as alternative clean source for possible energy production (electricity) and other low-temperature purposes.  相似文献   

5.
Jowshan geothermal system comprises 6 thermal springs with outlet temperatures ranging from 39.3 to 46.6°C. The thermal water of these springs is presently used for swimming and as a treatment for rheumatism, sinusitis and skin diseases. The pH value of these springs is slightly acidic to neutral and the electrical conductivities about 1500 μS/Cm. The presence of many faults in the area, the alignment of all springs along the Sirch Fault and the similar chemical and isotopic composition of all springs in combination with the hydrogeological setting and geochemistry of water samples indicate that these springs are associated with deep circulation of meteoric water. According to this heating mechanism, meteoric waters infiltrate through fault openings to depth and after heating by geothermal gradient rise to the ground surface due to the hydraulic and buoyancy forces, a mechanism which is common in the southern parts of Iran. The use of various chemical geothermometers and mineral equilibrium states suggests a range of temperature about 50–90°C for the reservoir of Jowshan geothermal system.  相似文献   

6.
Geothermal water is plentiful in Changbai Mountain region, northeastern China, due to the volcanic activities and widespread faults. For the exploration of geothermal resources, this study uses quartz and cation geothermometer to estimate the temperatures of the geothermal reservoir and uses the tubular models to evaluate the thermal gradient. The hydrogeochemical characteristics of the geothermal resources were also evaluated by hydrogeochemical analysis. The results showed that the geothermal reservoir temperatures of the four major thermal springs in Changbai Mountain region range from 72 to 169 °C. The average geothermal reservoir temperatures of Jinjiang hot springs, Changbai hot springs I, Xianrenqiao hot springs, and Changbai hot springs II are 129.25, 169, 89, and 73.67 °C, respectively. The geothermal gradient values of the four major thermal springs have different characteristics. The geothermal gradient values of Jinjiang hot springs and Changbai hot springs I are 4.6 and 3.1 °C/100 m, respectively. The geothermal gradient values of Xianrenqiao thermal springs and Changbai thermal springs II are both lower than 1.5 °C/100 m, with the values of 1.1 and 1.4 °C/100 m. And the geothermal gradients are influenced by Changbai Mountain Tianchi volcano. In addition, the water chemical analyses showed that the geothermal water types are HCO3-Na with higher concentrations of Na+, Cl?, SO4 2?, TDS, and HCO3 ? than the non-thermal waters, which suggested a deep and long water cycle of the thermal water, and therefore a sufficient water-rock interaction.  相似文献   

7.
In this paper, the hydrochemical isotopic characteristics of samples collected from geothermal springs in the Ilica geothermal field, Eastern Anatolia of Turkey, are examined and described. Low-temperature geothermal system of Ilica (Erzurum, Turkey) located along the Eastern Anatolian fault zone was investigated for hydrogeochemical and isotopic characteristics. The study of ionic and isotopic contents shows that the thermal water of Ilica is mainly, locally fed by groundwater, which changes chemically and isotopically during its circulation within the major fault zone reaching depths. The thermal spring has a temperature of 29–39 °C, with electrical conductivity ranging from 4,000 to 7,510 µS/cm and the thermal water is of Na–HCO3–Cl water type. The chemical geothermometers applied in the Ilica geothermal waters yielded a maximum reservoir temperature of 142 °C according to the silica geothermometers. The thermal waters are undersaturated with respect to gypsum, anhydrite and halite, and oversaturated with respect to dolomite. The dolomite mineral possibly caused scaling when obtaining the thermal waters in the study area. According to the enthalpy chloride-mixing model, cold water to the thermal water-mixing ratio is changing between 69.8 and 75 %. The δ18O–δ2H compositions obviously indicate meteoric origin of the waters. Thermal water springs derived from continental precipitation falling on to higher elevations in the study area. The δ13C ratio for dissolved inorganic carbonate in the waters lies between 4.63 and 6.48 ‰. In low-temperature waters carbon is considered as originating from volcanic (mantle) CO2.  相似文献   

8.
An integrated survey program involving geological, hydrogeological and geophysical techniques has been employed to characterize the aquifer geometry, recharge and circulation dynamics of thermal springs within a shallow aquifer system in Ethiopia. The selected springs for the case study are Sodere and Gergedi, which are situated within the tectonically active Main Ethiopian Rift (MER). Geologically, the studied springs are located on Plio-Quaternary volcanic rocks. The geophysical results indicate the presence of subsurface weak zones represented by extensional tectonics and weathering zones which are responsible for thermal water circulation and facilitate recharge from the adjacent surface-water bodies. The structures inferred by the resistivity survey, both sounding and electrical tomography, present contrasts in rock resistivity response. The anomalous zones in the magnetic data are in good agreement with the zones that are revealed by geological mapping and surface manifestation of the thermal water discharge zones. The shallow aquifer of the central MER is under the influence of thermal water, which increases the groundwater temperature and mineral content.  相似文献   

9.
Kangding geothermal area is located in the western Sichuan, belonging to southeastern margin of Tibetan Plateau. Similar to world-renowned south Tibetan and western Yunnan geothermal belt, western Sichuan has intensive surface thermal manifestations including boiling and hot springs. The emerging temperature of thermal waters ranges from 47 to 79 °C with total dissolved solids lying between 899 and 2550 mg/L. δ2H–δ18O isotopes indicate a meteoric source for the thermal waters and a significant positive oxygen-18 shift in the southern region. It is suggested that southern thermal waters experienced stronger water–rock interaction and are closer to thermodynamic equilibrium, which is also proved by the water type classification. The reservoir temperature calculated by empirical and theoretical chemical thermometry is 180–225 °C for the north and 225–310 °C for the south. Evidences of hydrogeochemistry, stable isotopes, geothermometry and radiocarbon dating indicate that southern region of Kangding area shows greater geothermal potential than the northern region. In addition, based on the hydrogeochemical modeling of mineral saturation, underlying problem of scaling is likely to occur in the study area. According to the results of reservoir temperature, south Kangding sub-district has greater potential in geothermal power generation and development than northern Kangding. Therefore, further exploration and drilling work should give priority to the south Kangding area.  相似文献   

10.
Subsurface thermal structure in Tohoku district are characterized by existing data such as geothermal resources maps, drill hole thermal gradients, Curie point depths and hypocenters distribution maps. The collected data are registered in a database system, then, compared in plan view, cross-section and bird's-eye pictures. The comparison indicates that subsurface temperatures extrapolated from drill hole thermal gradients are generally concordant to the Curie point depth, assumed to be 650 °C. Tohoku district is generally divided into 5 type areas; fore arc lowland, fore arc mountain country, Quaternary volcanic terrain, back arc lowland and back arc mountain country. The surface thermal manifestations in Quaternary volcanic terrain are mainly controlled by the magma chambers as heat sources, while, surface thermal features such as hot springs in non-volcanic areas are controlled by degrees of heat flows, and hydrothermal flows in permeable Cenozoic formations and along permeable fault zones.  相似文献   

11.
Three groups of thermal springs with temperatures close to 70 °C discharge both in the core (at Bagni di Vinadio and Terme di Valdieri) and on the external margin (at Berthemont-Les-Bains) of the Argentera Massif. Detailed structural field analysis carried out on the hydrothermal sites allows us to delineate both a model of Alpine tectonic evolution of the Argentera Massif and the patterns of hydrothermal circulation that were active during its final exhumation. The observed fault rock assemblages provide information relative to deformation that occurred in viscous, frictional-to-viscous and frictional crustal regimes. During the Early Miocene, the Bersezio Fault Zone and the Fremamorta Shear Zone, two main mylonitic shear zones, mainly accommodated regional transpression and provided pathways for fluid flow promoting mineral reactions in greenschist facies. During the Late Miocene–Early Pliocene, frictional-to-viscous deformation affected the massif, which underwent predominant transpression in the internal sectors and extension on the external margin. During the Plio-Pleistocene, deformation in frictional condition accompanied the final exhumation of the massif in a transpressive regime and resulted in the development of the NW–SE striking cataclastic zones. The hydraulic properties of these structures mainly influence the patterns of the active thermal circulations and the localization of the recharge and discharge zones. At Berthemont these faults represent conduits, whereas at Vinadio and Valdieri they form complex systems of conduits and barriers. In these two latter sites, the cataclastic faults compose flower structures that constrain laterally the thermal fluid flows while intensely fractured granites sited at depth constitute a highly-transmissive geothermal reservoir. Less permeable migmatitic gneisses overlaying the granites prevent a massive infiltration of the cold fluids at depth. This context favours within the high-permeability fractures granites the development of buoyancy-driven flows which combined with topographically-driven flows, provided the conditions for the upflow of the high-temperature waters.  相似文献   

12.
Uttarakhand geothermal area, located in the central belt of the Himalayan geothermal province, is one of the important high temperature geothermal fields in India. In this study, the chemical characteristics of the thermal waters are investigated to identify the main geochemical processes affecting the composition of thermal waters during its ascent toward the surface as well as to determine the subsurface temperature of the feeding reservoir. The thermal waters are mainly Ca–Mg–HCO3 type with moderate silica and TDS concentrations. Mineral saturation states calculated from PHREEQC geochemical code indicate that thermal waters are supersaturated with respect to calcite, dolomite, aragonite, chalcedony, quartz (SI > 0), and undersaturated with respect to gypsum, anhydrite, and amorphous silica (SI < 0). XRD study of the spring deposit samples fairly corroborates the predicted mineral saturation state of the thermal waters. Stable isotopes (δ18O, δ2H) data confirm the meteoric origin of the thermal waters with no oxygen-18 shift. The mixing phenomenon between thermal water with shallow ground water is substantiated using tritium (3H) and chemical data. The extent of dilution is quantified using tritium content of thermal springs and non-thermal waters. Classical geothermometers, mixing model, and multicomponent fluid geothermometry modeling (GeoT) have been applied to estimate the subsurface reservoir temperature. Among different classical geothermometers, only quartz geothermometer provide somewhat reliable estimation (96–140 °C) of the reservoir temperature. GeoT modeling results suggest that thermal waters have attained simultaneous equilibrium with respect to minerals like calcite, quartz, chalcedony, brucite, tridymite, cristobalite, talc, at the temperature 130 ± 5 °C which is in good agreement with the result obtained from the mixing model.  相似文献   

13.
为查明色达—松潘断块地热资源赋存状态及热源来源,以四川黑水县内3处温泉(热水塘、上达古、卡龙沟)为研究对象,采集温泉水样进行水化学分析和同位素测试,研究地热水的补给来源和热储温度。研究结果显示,热水塘温泉的地下水化学类型为HCO3-Na型,上达古温泉和卡龙沟温泉的地下水化学类型为HCO3-Ca型,补给水源主要为大气降水,补给高程分别为5 121 m、3 890 m、3 921 m。结合矿物饱和指数,采用SiO2地热温标计算3处温泉的热储温度,分别为119.036 ℃、49.034 ℃、30.215 ℃。综合分析认为研究区地下热水的成因主要为大气降水经高山补给区入渗至储集层,吸取地下深部向上传导的热量和放射性元素衰变释放的热量,并与围岩发生水-岩作用形成地下热水,在断裂发育部位热水沿断裂带向上运移,最后在地表出露形成温泉。  相似文献   

14.
Pamancalan is located in Lebak Regency, Province of Banten. This area is located in the western part of Java Island, Indonesia. The geothermal manifestation in this area is in the form of hot spring in Cipamancalan River. But, how the structure of the geothermal system, the reservoir depth, and the thermal source in the Pancamalan area has not been studied much. Therefore, there is a need to conduct a geophysics study by surveying the gravity and magnetic field. A study, which used magnetic and gravity data to discover reservoir, has been conducted in the Pamancalan geothermal area in Lebak, Banten. Topographic map for total magnetic and earth gravity anomaly shows that the anomaly is located in the center and southern part of the presumed reservoir. 2.5-D section model of magnetic anomaly shows that there is a rock formation which shaped the geothermal system in Pamancalan. The thermal source is diorite which is a by-product of Hanjawar Mountain; the reservoir rocks consist of sandstone, limestone, and breccia; and the cap rock is in the form of clay and tuff. It is predicted that there is an intrusive body which functions as thermal source in the depth between 1650 and 4000 m, the reservoir depth is around 700 m, and the depth of clay cap is around 0 to 700 m. The geothermal manifestation in Pamancalan area is controlled by Cigeledug fault from the southwest and Cipamancalan fault from the north and south.  相似文献   

15.
为研究四川省康定市二道桥地区地下热水稳定同位素特征和热储温度,对二道桥地区5个温泉(井)即二道桥温泉(SC107、SC107-2)、康巴人家温泉(SC107-3)、自流热水井(SC107-4)、自喷热水井(SC107-5)进行调查和分析。研究区温泉的分布及出露主要受雅拉沟断裂和雅拉河谷控制。温泉水温33.2~46℃,为中低温温泉,pH为6~6.5。水样的氢氧稳定同位素特征表明研究区地下热水的补给来源为大气降水。利用氢氧稳定同位素高程效应及温度效应估算区内地下热水补给区高程为3 000~4 500 m,补给区温度为-3.5~-0.3℃,表明地下热水有一部分补给源自附近山区的冰雪融水。Na-K-Mg三角图显示研究区热水均为未成熟水,不宜用阳离子地热温标计算热储温度。应用SiO2地热温标、多矿物饱和指数法以及用固定铝方法对部分温泉多矿物平衡图进行修正,得出研究区地下热水的热储温度为65~75℃。研究区温泉在东部跑马山以及西部农戈山附近接受大气降水补给,降水沿着大雪山—农戈山断裂和跑马山断裂下渗,地下水经历深循环,在此过程中获得大地热流加热,最终在雅拉河谷雅拉沟断裂附近出...  相似文献   

16.
The temperature of groundwater will be increasing gradually while it is recharged and infiltrated down to the depth because of the geothermal gradient. So the density of water is smaller and smaller. The author puts forward the heat-driving mechanism that the movement of the groundwater is caused by the density differences from different temperatures. Based on this idea, the author proposes the short range recharge model of geothermal water in flat terrain hilly area. The model explains the coexistence mechanism of recharge and discharge (hot springs) zones can be at the same altitude. Considering that the temperature is an important driving force of groundwater movement, the paper puts forward and demonstrates that the hot springs can expose at higher lands while the recharge zone being situated lower lying areas.  相似文献   

17.
Deep 3D thermal modelling for the city of Berlin (Germany)   总被引:1,自引:1,他引:0  
This study predicts the subsurface temperature distribution of Germany’s capital Berlin. For this purpose, a data-based lithosphere-scale 3D structural model is developed incorporating 21 individual geological units. This model shows a horizontal grid resolution of (500 × 500) m and provides the geometric base for two different approaches of 3D thermal simulations: (1) calculations of the steady-state purely conductive thermal field and (2) simulations of coupled fluid flow and heat transport. The results point out fundamentally different structural and thermal configurations for potential geothermal target units. The top of the Triassic Middle Buntsandstein strongly varies in depth (159–2,470 m below sea level) and predicted temperatures (15–95 °C), mostly because of the complex geometry of the underlying Permian Zechstein salt. The top of the sub-salt Sedimentary Rotliegend is rather flat (2,890–3,785 m below sea level) and reveals temperatures of 85–139 °C. The predicted 70 °C-isotherm is located at depths of about 1,500–2,200 m, cutting the Middle Buntsandstein over large parts of Berlin. The 110 °C-isotherm at 2,900–3,700 m depth widely crosscuts the Sedimentary Rotliegend. Groundwater flow results in subsurface cooling the extent of which is strongly controlled by the geometry and the distribution of the Tertiary Rupelian Clay. The cooling effect is strongest where this clay-rich aquitard is thinnest or missing, thus facilitating deep-reaching forced convective flow. The differences between the purely conductive and coupled models highlight the need for investigations of the complex interrelation of flow- and thermal fields to properly predict temperatures in sedimentary systems.  相似文献   

18.
Shallow geothermal prospecting ( < 700 m) has been performed in four zones in Morocco for which few deep data are available: northwestern basin, northeastern basin, Tadla Basin and Agadir Basin. These areas are different geologically and hydrogeologically. The temperature data from 250 wells at depths between 15 and 500 m have been analysed in order to estimate the natural geothermal gradient in these areas, to determine the principal thermal anomalies, to identify the main thermal indices and to characterise the recharge, discharge and potential mixing limits of the aquifers.The hydrostratigraphical study of each basin revealed several potential reservoir layers in which the Turonian carbonate aquifer (Tadal and Agadir Basins) and Liassic acquifer (Moroccan northwestern and northeastern basins) are the most important hot water reservoirs in Morocco. The recharge zones of each aquifer are characterised by high topography, high water potential, shallow cold water, low geothermal gradient and negative anomalies. The discharge zones are characterized by low topography, low piezometric level, high geothermal gradient, high temperature with hot springs and positive anomalies. The main thermal indices and the principal thermal anomalies that coincide with the artesian zones of the Turonian and Liassic aquifers have been identified.  相似文献   

19.
The groundwater flow regime at great depth within the Molasse Basin (SW Germany) was studied. Data relevant for a flow model at 600–1,600 m depth are sparse in the western part of the basin. However, temperature measurements are available covering much of the area at a wide range of depths. Therefore, a thermal 3D steady-state model was set up with the aim of comparing modeled with observed subsurface temperatures. Stratigraphic information from many boreholes was also available, but only a few values of rock thermal conductivity and heat-production rate could be obtained. Some strong thermal residual anomalies were identified with respect to the purely conductive model, especially along fault zones, and within stratigraphic layers with high hydraulic conductivity. These anomalies can be explained by various advective heat-transport mechanisms, yet most explanations can be eliminated. The most plausible constellation explaining the major positive thermal anomalies of 10 Kelvin and more is a fault zone of E–W strike, intersected by an aquifer with flow parallel to the fault zone. This concept was investigated by using a simplified type model. In spite of some shortcomings, the method presented here can be used to identify temperature anomalies, and to identify possible explanations.  相似文献   

20.
Hydrothermal simulations are used to provide insight into the subsurface thermal regime of the Perth metropolitan area (PMA) in Western Australia. High average permeabilities and estimated fluid flow rates in shallow aquifers of the PMA suggest that advection and convection may occur in these aquifers. These processes are simulated, using a new geological model of the PMA to constrain the geometry of aquifers, aquitards and faults. The results show that advection has a strong influence on subsurface temperature, especially in the north of the PMA, where aquifer recharge creates an area of anomalously low temperature. Convection may be important, depending on the permeability of the Yarragadee Aquifer. If convection occurs, it creates thermal highs and lows with a spacing of approximately 5 km. Some of these thermal anomalies migrate over geological time due to coupling between advection and convection, but they are stationary on human timescales. Fault permeability influences the pattern of convection. Advection and convection cause variations in the geothermal gradient which cannot be predicted by conductive models; therefore, these processes should be considered in any model that is used for assessment of geothermal resources in the PMA.  相似文献   

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

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