首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
Water and gas samples were collected from (1) nine shallow groundwater aquifers overlying Marcellus Shale in north‐central West Virginia before active shale gas drilling, (2) wells producing gas from Upper Devonian sands and Middle Devonian Marcellus Shale in southwestern Pennsylvania, (3) coal‐mine water discharges in southwestern Pennsylvania, and (4) streams in southwestern Pennsylvania and north‐central West Virginia. Our preliminary results demonstrate that the oxygen and hydrogen isotope composition of water, carbon isotope composition of dissolved inorganic carbon, and carbon and hydrogen isotope compositions of methane in Upper Devonian sands and Marcellus Shale are very different compared with shallow groundwater aquifers, coal‐mine waters, and stream waters of the region. Therefore, spatiotemporal stable isotope monitoring of the different sources of water before, during, and after hydraulic fracturing can be used to identify migrations of fluids and gas from deep formations that are coincident with shale gas drilling.  相似文献   

2.
Horizontal drilling and hydraulic fracturing have enhanced unconventional hydrocarbon recovery but raised environmental concerns related to water quality. Because most basins targeted for shale‐gas development in the USA have histories of both active and legacy petroleum extraction, confusion about the hydrogeological context of naturally occurring methane in shallow aquifers overlying shales remains. The Karoo Basin, located in South Africa, provides a near‐pristine setting to evaluate these processes, without a history of conventional or unconventional energy extraction. We conducted a comprehensive pre‐industrial evaluation of water quality and gas geochemistry in 22 groundwater samples across the Karoo Basin, including dissolved ions, water isotopes, hydrocarbon molecular and isotopic composition, and noble gases. Methane‐rich samples were associated with high‐salinity, NaCl‐type groundwater and elevated levels of ethane, 4He, and other noble gases produced by radioactive decay. This endmember displayed less negative δ13C‐CH4 and evidence of mixing between thermogenic natural gases and hydrogenotrophic methane. Atmospheric noble gases in the methane‐rich samples record a history of fractionation during gas‐phase migration from source rocks to shallow aquifers. Conversely, methane‐poor samples have a paucity of ethane and 4He, near saturation levels of atmospheric noble gases, and more negative δ13C‐CH4; methane in these samples is biogenic and produced by a mixture of hydrogenotrophic and acetoclastic sources. These geochemical observations are consistent with other basins targeted for unconventional energy extraction in the USA and contribute to a growing data base of naturally occurring methane in shallow aquifers globally, which provide a framework for evaluating environmental concerns related to unconventional energy development (e.g., stray gas).  相似文献   

3.
Environmental concerns regarding the potential for drinking water contamination in shallow aquifers have accompanied unconventional energy development in the northern Appalachian Basin. These activities have also raised several critical questions about the hydrogeological parameters that control the naturally occurring presence and migration of hydrocarbon gases in shallow aquifers within petroliferous basins. To interrogate these factors, we analyzed the noble gas, dissolved ion, and hydrocarbon gas (molecular and isotopic composition) geochemistry of 98 groundwater samples from south‐central New York. All samples were collected ?1km from unconventional drilling activities and sample locations were intentionally targeted based on their proximity to various types of documented fault systems. In agreement with studies from other petroliferous basins, our results show significant correlations between elevated levels of radiogenic [4He], thermogenic [CH4], and dissolved ions (e.g., Cl, Br, Sr, Ba). In combination, our data suggest that faults have facilitated the transport of exogenous hydrocarbon‐rich brines from Devonian source rocks into overlying Upper Devonian aquifer lithologies over geologic time. These data conflict with previous reports, which conclude that hydrodynamic focusing regulates the occurrence of methane and salt in shallow aquifers and leads to elevated levels of these species in restricted flow zones within valley bottoms. Instead, our data suggest that faults in Paleozoic rocks play a fundamental role in gas and brine transport from depth, regulate the distribution of their occurrence in shallow aquifers, and influence the geochemistry of shallow groundwater in this petroliferous basin.  相似文献   

4.
There is concern about adverse impacts of natural gas (primarily methane) production on groundwater quality; however, data on trace element concentrations are limited. The objective of this study was to compare the distribution of trace elements in groundwater samples with and without dissolved methane in aquifers overlying the Barnett Shale (Hood and Parker counties, 207 samples) and the Haynesville Shale (Panola County, 42 samples). Both shales have been subjected to intensive hydraulic fracturing for gas production. Well clusters with high dissolved methane were previously found in these counties and are thought to be of natural origin. Overall, groundwater in these counties is of excellent quality with typically low elemental concentrations. Several statistical analyses strongly suggest that most trace element concentrations, generally at low background levels, are no higher and even reduced when dissolved methane is present. In addition, trace element concentrations are not correlated with distance to gas wells. The reduction in trace element concentrations is attributed to anaerobic microbial degradation of methane, is associated with a higher pH (>8.5), and, likely, with precipitation of carbonates and pyrite and formation of clays. Trace and other elements are likely incorporated within the precipitating mineral crystalline network or sorbed. High pH values are found throughout these high‐methane clusters (e.g., Parker‐Hood cluster), even in subregions where methane is not present, which is consistent with a pervasive natural origin of dissolved methane rather than a limited gas well source.  相似文献   

5.
Assessing natural vs. anthropogenic sources of methane in drinking water aquifers is a critical issue in areas of shale oil and gas production. The objective of this study was to determine controls on methane occurrences in aquifers in the Eagle Ford Shale play footprint. A total of 110 water wells were tested for dissolved light alkanes, isotopes of methane, and major ions, mostly in the eastern section of the play. Multiple aquifers were sampled with approximately 47 samples from the Carrizo‐Wilcox Aquifer (250‐1200 m depth range) and Queen City‐Sparta Aquifer (150‐900 m depth range) and 63 samples from other shallow aquifers but mostly from the Catahoula Formation (depth <150 m). Besides three shallow wells with unambiguously microbial methane, only deeper wells show significant dissolved methane (22 samples >1 mg/L, 10 samples >10 mg/L). No dissolved methane samples exhibit thermogenic characteristics that would link them unequivocally to oil and gas sourced from the Eagle Ford Shale. In particular, the well water samples contain very little or no ethane and propane (C1/C2+C3 molar ratio >453), unlike what would be expected in an oil province, but they also display relatively heavier δ13Cmethane (>?55‰) and δDmethane (>?180‰). Samples from the deeper Carrizo and Queen City aquifers are consistent with microbial methane sourced from syndepositional organic matter mixed with thermogenic methane input, most likely originating from deeper oil reservoirs and migrating through fault zones. Active oxidation of methane pushes δ13Cmethane and δDmethane toward heavier values, whereas the thermogenic gas component is enriched with methane owing to a long migration path resulting in a higher C1/C2+C3 ratio than in the local reservoirs.  相似文献   

6.
Tom Myers 《Ground water》2012,50(6):872-882
Hydraulic fracturing of deep shale beds to develop natural gas has caused concern regarding the potential for various forms of water pollution. Two potential pathways—advective transport through bulk media and preferential flow through fractures—could allow the transport of contaminants from the fractured shale to aquifers. There is substantial geologic evidence that natural vertical flow drives contaminants, mostly brine, to near the surface from deep evaporite sources. Interpretative modeling shows that advective transport could require up to tens of thousands of years to move contaminants to the surface, but also that fracking the shale could reduce that transport time to tens or hundreds of years. Conductive faults or fracture zones, as found throughout the Marcellus shale region, could reduce the travel time further. Injection of up to 15,000,000 L of fluid into the shale generates high pressure at the well, which decreases with distance from the well and with time after injection as the fluid advects through the shale. The advection displaces native fluids, mostly brine, and fractures the bulk media widening existing fractures. Simulated pressure returns to pre‐injection levels in about 300 d. The overall system requires from 3 to 6 years to reach a new equilibrium reflecting the significant changes caused by fracking the shale, which could allow advective transport to aquifers in less than 10 years. The rapid expansion of hydraulic fracturing requires that monitoring systems be employed to track the movement of contaminants and that gas wells have a reasonable offset from faults.  相似文献   

7.
In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.  相似文献   

8.
Horizontal drilling and hydraulic fracturing have enabled hydrocarbon recovery from unconventional reservoirs, but led to natural gas contamination of shallow groundwaters. We describe and apply numerical models of gas‐phase migration associated with leaking natural gas wells. Three leakage scenarios are simulated: (1) high‐pressure natural gas pulse released into a fractured aquifer; (2) continuous slow leakage into a tilted fractured formation; and (3) continuous slow leakage into an unfractured aquifer with fluvial channels, to facilitate a generalized evaluation of natural gas transport from faulty natural gas wells. High‐pressure pulses of gas leakage into sparsely fractured media are needed to produce the extensive and rapid lateral spreading of free gas previously observed in field studies. Transport in fractures explains how methane can travel vastly different distances and directions laterally away from a leaking well, which leads to variable levels of methane contamination in nearby groundwater wells. Lower rates of methane leakage (≤1 Mcf/day) produce shorter length scales of gas transport than determined by the high‐pressure scenario or field studies, unless aquifers have low vertical permeabilities (≤1 millidarcy) and fractures and bedding planes have sufficient tilt (～10°) to allow a lateral buoyancy component. Similarly, in fractured rock aquifers or where permeability is controlled by channelized fluvial deposits, lateral flow is not sufficiently developed to explain fast‐developing gas contamination (0‐3 months) or large length scales (～1 km) documented in field studies. Thus, current efforts to evaluate the frequency, mechanism, and impacts of natural gas leakage from faulty natural gas wells likely underestimate contributions from small‐volume, low‐pressure leakage events.  相似文献   

9.
The expanding use of horizontal drilling and hydraulic fracturing technology to produce oil and gas from tight rock formations has increased public concern about potential impacts on the environment, especially on shallow drinking water aquifers. In eastern Kentucky, horizontal drilling and hydraulic fracturing have been used to develop the Berea Sandstone and the Rogersville Shale. To assess baseline groundwater chemistry and evaluate methane detected in groundwater overlying the Berea and Rogersville plays, we sampled 51 water wells and analyzed the samples for concentrations of major cations and anions, metals, dissolved methane, and other light hydrocarbon gases. In addition, the stable carbon and hydrogen isotopic composition of methane (δ13C‐CH4 and δ2H‐CH4) was analyzed for samples with methane concentration exceeding 1 mg/L. Our study indicates that methane is a relatively common constituent in shallow groundwater in eastern Kentucky, where methane was detected in 78% of the sampled wells (40 of 51 wells) with 51% of wells (26 of 51 wells) exhibiting methane concentrations above 1 mg/L. The δ13C‐CH4 and δ2H‐CH4 ranged from ?84.0‰ to ?58.3‰ and from ?246.5‰ to ?146.0‰, respectively. Isotopic analysis indicated that dissolved methane was primarily microbial in origin formed through CO2 reduction pathway. Results from this study provide a first assessment of methane in the shallow aquifers in the Berea and Rogersville play areas and can be used as a reference to evaluate potential impacts of future horizontal drilling and hydraulic fracturing activities on groundwater quality in the region.  相似文献   

10.
Gaining streams can provide an integrated signal of relatively large groundwater capture areas. In contrast to the point‐specific nature of monitoring wells, gaining streams coalesce multiple flow paths. Impacts on groundwater quality from unconventional gas development may be evaluated at the watershed scale by the sampling of dissolved methane (CH4) along such streams. This paper describes a method for using stream CH4 concentrations, along with measurements of groundwater inflow and gas transfer velocity interpreted by 1‐D stream transport modeling, to determine groundwater methane fluxes. While dissolved ionic tracers remain in the stream for long distances, the persistence of methane is not well documented. To test this method and evaluate CH4 persistence in a stream, a combined bromide (Br) and CH4 tracer injection was conducted on Nine‐Mile Creek, a gaining stream in a gas development area in central Utah. A 35% gain in streamflow was determined from dilution of the Br tracer. The injected CH4 resulted in a fivefold increase in stream CH4 immediately below the injection site. CH4 and δ13CCH4 sampling showed it was not immediately lost to the atmosphere, but remained in the stream for more than 2000 m. A 1‐D stream transport model simulating the decline in CH4 yielded an apparent gas transfer velocity of 4.5 m/d, describing the rate of loss to the atmosphere (possibly including some microbial consumption). The transport model was then calibrated to background stream CH4 in Nine‐Mile Creek (prior to CH4 injection) in order to evaluate groundwater CH4 contributions. The total estimated CH4 load discharging to the stream along the study reach was 190 g/d, although using geochemical fingerprinting to determine its source was beyond the scope of the current study. This demonstrates the utility of stream‐gas sampling as a reconnaissance tool for evaluating both natural and anthropogenic CH4 leakage from gas reservoirs into groundwater and surface water.  相似文献   

11.
The rapid expansion in unconventional gas development over the past two decades has led to concerns over the potential impacts on groundwater resources. Although numerical models are invaluable for assessing likelihood of impacts at particular sites, simpler analytical models are also useful because they help develop hydrological understanding. Analytical approaches are also valuable for preliminary assessments and to determine where more complex models are warranted. In this article, we present simple analytical solutions that can be used to predict: (1) the spatial extent of drawdown from horizontal wells drilled into the gas‐bearing formation, and rate of recovery after gas production ceases; (2) the potential for upward transport of contaminants from the gas‐bearing formation to shallow aquifers during hydraulic fracturing operations when pressures in the gas‐bearing formation are greatly increased; and (3) the potential downward leakage of water from shallow aquifers during depressurization of gas‐bearing formations. In particular, we show that the recovery of pressure after production ceases from gas‐bearing shale formations may take several hundred years, and we present critical hydraulic conductivity values for intervening aquitards, below which the impact on shallow aquifers will be negligible. The simplifying assumptions inherent in these solutions will limit their predictive accuracy for site‐specific assessments, compared to numerical models that incorporate knowledge of spatial variations in formation properties and which may include processes not considered in the simpler solutions.  相似文献   

12.
Enhanced production of unconventional hydrocarbons in the United States has driven interest in natural gas development globally, but simultaneously raised concerns regarding water quantity and quality impacts associated with hydrocarbon extraction. We conducted a pre‐development assessment of groundwater geochemistry in the critically water‐restricted Karoo Basin, South Africa. Twenty‐two springs and groundwater samples were analyzed for major dissolved ions, trace elements, water stable isotopes, strontium and boron isotopes, hydrocarbons and helium composition. The data revealed three end‐members: a deep, saline groundwater with a sodium‐chloride composition, an old, deep freshwater with a sodium‐bicarbonate‐chloride composition and a shallow, calcium‐bicarbonate freshwater. In a few cases, we identified direct mixing of the deep saline water and shallow groundwater. Stable water isotopes indicate that the shallow groundwater was controlled by evaporation in arid conditions, while the saline waters were diluted by apparently fossil meteoric water originated under wetter climatic conditions. These geochemical and isotopic data, in combination with elevated helium levels, suggest that exogenous fluids are the source of the saline groundwater and originated from remnant seawater prior to dilution by old meteoric water combined with further modification by water‐rock interactions. Samples with elevated methane concentrations (>14 ccSTP/kg) were strongly associated with the sodium‐chloride water located near dolerite intrusions, which likely provide a preferential pathway for vertical migration of deeply sourced hydrocarbon‐rich saline waters to the surface. This pre‐drill evaluation indicates that the natural migration of methane‐ and salt‐rich waters provides a source of geogenic contamination to shallow aquifers prior to shale gas development in the Karoo Basin.  相似文献   

13.
Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.  相似文献   

14.
A large‐scale groundwater flow and transport model is developed for a deep‐seated (100 to 300 m below ground surface) sedimentary aquifer system. The model is based on a three‐dimensional (3D) hydrostratigraphic model, building on a sequence stratigraphic approach. The flow model is calibrated against observations of hydraulic head and stream discharge while the credibility of the transport model is evaluated against measurements of 39Ar from deep wells using alternative parameterizations of dispersivity and effective porosity. The directly simulated 3D mean age distributions and vertical fluxes are used to visualize the two‐dimensional (2D)/3D age and flux distribution along transects and at the top plane of individual aquifers. The simulation results are used to assess the vulnerability of the aquifer system that generally has been assumed to be protected by thick overlaying clayey units and therefore proposed as future reservoirs for drinking water supply. The results indicate that on a regional scale these deep‐seated aquifers are not as protected from modern surface water contamination as expected because significant leakage to the deeper aquifers occurs. The complex distribution of local and intermediate groundwater flow systems controlled by the distribution of the river network as well as the topographical variation (Tóth 1963) provides the possibility for modern water to be found in even the deepest aquifers.  相似文献   

15.
Understanding the source of dissolved methane in drinking‐water aquifers is critical for assessing potential contributions from hydraulic fracturing in shale plays. Shallow groundwater in the Texas portion of the Haynesville Shale area (13,000 km2) was sampled (70 samples) for methane and other dissolved light alkanes. Most samples were derived from the fresh water bearing Wilcox formations and show little methane except in a localized cluster of 12 water wells (17% of total) in a approximately 30 × 30 km2 area in Southern Panola County with dissolved methane concentrations less than 10 mg/L. This zone of elevated methane is spatially associated with the termination of an active fault system affecting the entire sedimentary section, including the Haynesville Shale at a depth more than 3.5 km, and with shallow lignite seams of Lower Wilcox age at a depth of 100 to 230 m. The lignite spatial extension overlaps with the cluster. Gas wetness and methane isotope compositions suggest a mixed microbial and thermogenic origin with contribution from lignite beds and from deep thermogenic reservoirs that produce condensate in most of the cluster area. The pathway for methane from the lignite and deeper reservoirs is then provided by the fault system.  相似文献   

16.
Gas‐saturated groundwater forms bubbles when brought to atmospheric pressure, preventing precise determination of its in situ dissolved gas concentrations. To overcome this problem, a modeling approach called the atmospheric sampling method is suggested here to recover the in situ dissolved gas concentrations of groundwater collected ex situ under atmospheric conditions at the Horonobe Underground Research Laboratory, Japan. The results from this method were compared with results measured at the same locations using two special techniques, the sealed sampler and pre‐evacuated vial methods, that have been developed to collect groundwater under its in situ conditions. In gas‐saturated groundwater cases, dissolved methane and inorganic carbon concentrations derived using the atmospheric sampling method were mostly within ±4 and ±10%, respectively, of values from the sealed sampler and pre‐evacuated vial methods. In gas‐unsaturated groundwater, however, the atmospheric sampling method overestimated the in situ dissolved methane concentrations, because the groundwater pressure at which bubbles appear (Pcritical) was overestimated. The atmospheric sampling method is recommended for use where gas‐saturated groundwater can be collected only ex situ under atmospheric conditions.  相似文献   

17.
Evaluation of Methane Sources in Groundwater in Northeastern Pennsylvania   总被引：1，自引：0，他引：1  
Testing of 1701 water wells in northeastern Pennsylvania shows that methane is ubiquitous in groundwater, with higher concentrations observed in valleys vs. upland areas and in association with calcium‐sodium‐bicarbonate, sodium‐bicarbonate, and sodium‐chloride rich waters—indicating that, on a regional scale, methane concentrations are best correlated to topographic and hydrogeologic features, rather than shale‐gas extraction. In addition, our assessment of isotopic and molecular analyses of hydrocarbon gases in the Dimock Township suggest that gases present in local water wells are most consistent with Middle and Upper Devonian gases sampled in the annular spaces of local gas wells, as opposed to Marcellus Production gas. Combined, these findings suggest that the methane concentrations in Susquehanna County water wells can be explained without the migration of Marcellus shale gas through fractures, an observation that has important implications for understanding the nature of risks associated with shale‐gas extraction.  相似文献   

18.
Although temperature is an important determinant of many biogeochemical processes in groundwater, very few studies have attempted to forecast the response of groundwater temperature to future climate warming. Using a composite linear regression model based on the lagged relationship between historical groundwater and regional air temperature data, empirical forecasts were made of groundwater temperature in several aquifers in Switzerland up to the end of the current century. The model was fed with regional air temperature projections calculated for greenhouse‐gas emissions scenarios A2, A1B, and RCP3PD. Model evaluation revealed that the approach taken is adequate only when the data used to calibrate the models are sufficiently long and contain sufficient variability. These conditions were satisfied for three aquifers, all fed by riverbank infiltration. The forecasts suggest that with respect to the reference period 1980 to 2009, groundwater temperature in these aquifers will most likely increase by 1.1 to 3.8 K by the end of the current century, depending on the greenhouse‐gas emissions scenario employed.  相似文献   

19.
BIOSCREEN is a well‐known simple tool for evaluating the transport of dissolved contaminants in groundwater, ideal for rapid screening and teaching. This work extends the BIOSCREEN model for the calculation of stable isotope ratios in contaminants. A three‐dimensional exact solution of the reactive transport from a patch source, accounting for fractionation by first‐order decay and/or sorption, is used. The results match those from a previously published isotope model but are much simpler to obtain. Two different isotopes may be computed, and dual isotope plots can be viewed. The dual isotope assessment is a rapidly emerging new approach for identifying process mechanisms in aquifers. Furthermore, deviations of isotope ratios at specific reactive positions with respect to “bulk” ratios in the whole compound can be simulated. This model is named BIOSCREEN‐AT‐ISO and will be downloadable from the journal homepage.  相似文献   

20.
The aim of this research was to refine the actual conceptual model related to the activation of high‐altitude temporary springs within the carbonate Apennines in southern Italy. The research was carried out through geophysical, hydrogeological, hydrochemical and isotopic investigations at the Acqua dei Faggi experimental site during five hydrologic years. The research demonstrated that, in carbonate aquifers where low‐permeability faults cause the aquifer system to be compartmentalized, high‐altitude temporary springs may be recharged by groundwater. In such settings, neither surface water infiltration in karst systems nor perched temporary aquifers play a role of utmost importance. The rare (once or a few time a year) activation of such springs is due to the fact that groundwater unusually reach the threshold head that allows the spring to flow. The activation of the studied high‐altitude temporary spring also depended on relationships between a low‐permeability fault core and a karst system that locally interrupts the low‐permeability barrier. In fact, when the hydraulic head did not reach the karst system, the concentrated head loss within the fault core did not allow the spring to flow, because the groundwater entirely flowed through the fault towards the downgradient compartment. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

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

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