首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The surface nuclear magnetic resonance (SNMR) method has been tested at a site in Haldensleben, northern Germany, to assess the suitability of this new method for groundwater exploration and environmental investigations. More information is obtained by SNMR, particularly with respect to aquifer parameters, than with other geophysical techniques. SNMR measurements were carried out at three borehole locations, together with 2D and 1D direct current geoelectrics, as well as ground-penetrating radar, and well logging (induction log, gamma-ray log and pulsed neutron-gamma log). Permeabilities were calculated from the grain-size distributions of core material determined in the laboratory. It is demonstrated that the SNMR method is able to detect groundwater and the results are in good agreement with other geophysical and hydrogeological data. Using the SNMR method, the water content of the unsaturated and saturated zones (i.e. porosity of an aquifer) can be reliably determined. This information and resistivity data permit in situ determination of other aquifer parameters. Comparison of the SNMR results with borehole data clearly shows that the water content determined by SNMR is the free or mobile water in the pores. The permeabilities estimated from the SNMR decay times are similar to those derived from sieve analysis of core material. Thus, the combination of SNMR with geoelectric methods promises to be a powerful tool for studying aquifer properties.  相似文献   

2.
Our study focuses on the potential usefulness of surface geophysical data to constrain the water content within an alluvial aquifer. On a study area where two wells have been drilled, we have performed several geophysical measurements, including ground penetrating radar, DC resistivity prospecting, seismic refraction survey and magnetic resonance soundings. From these data, we estimated several parameters, namely, the water height in the deposits, the effective porosity, the water content, the permeability, and the transmissivity of alluvial deposits. These physical parameters allow us to characterize the alluvial deposits in order to constrain the estimation of the potential water flow. The lithology and water flow rate known from the wells enabled us to compare geophysical results obtained in a high water flow rate zone to those in a low water flow rate zone. Correlation has been found between the water flow rate observed in both wells and the geophysical data obtained in the vicinity of these wells.  相似文献   

3.
Several field surveys of a waste rock pile were carried out during the summers of 2002 and 2003 using ground-penetrating radar, electromagnetic conductivity and DC resistivity imaging. The waste rock deposit is prone to generate acid mine drainage (AMD) due to the oxidation of sulphidic minerals. One of the most critical factors that lead to the production of AMD is unsaturated water flow and the ensuing moisture distribution in the waste rock. This geophysical characterization study, performed over a 30 m × 30 m test zone, was designed to image the internal structure controlling the water flux at shallow depth. The subsurface was found to consist of three zones for the first 6 m of the pile, mainly based on electrical resistivities: a thin superficial conductive material, an intermediate 2 to 3 m thick highly resistive zone, and a lower, more conductive medium. With the help of hydrogeological tests, chemical analyses and two 2.5 m-deep trenches, it is shown that the two conductive zones are correlated with fine-grained waste rock and the resistive zone correlates with a coarser material. In the two deeper zones, the contact between the two types of waste rock is typically highlighted by a sharp resistive/conductive boundary. An increase of conductance in the relatively thin upper layer towards the edge of the pile appears to be caused by an increase in thickness of the fine-grained material. Additional geophysical surveys carried out on a profile along the flank of the upper bench of the pile show that the main features of the internal structure are sub-parallel to the slope, at least for the first 3 m in depth. The data also show an increase in resistivity from the top to bottom of the slope, in accordance with expected particle segregation, from fine-grained material at the top to coarser material at the bottom. Wide-angle reflection GPR monitoring during large scale infiltration tests seems to indicate preferential flow paths towards the direction of coarser, more pervious material (which also appears to be less oxidized). Water preferentially flows through the coarse-grained material, but it is stored by capillary forces in the fine-grained material. Apart from the deposition methods, the results strongly suggest that factors such as machinery-induced mechanical alteration, construction history of the pile, and increased oxidization near the edges could explain the resistivity model. The model interpreted from geophysical imaging agrees well with the conceptual model of the rock pile. The resistivity and GPR methods appear to be efficient geophysical methods to characterize the internal structure and preferential flow patterns within unsaturated waste rock piles.  相似文献   

4.
The U.S. Environmental Protection Agency (EPA) recently proposed to amend federal regulations to require vadose zone monitoring at certain hazardous waste facilities. To support this proposal, EPA evaluated previous policy on vadose zone monitoring and examined advances in vadose zone monitoring technology. Changes in EPA vadose zone monitoring policy were driven by demonstrated advances in the available monitoring technology and improvements in understanding of vadose zone processes/When used under the appropriate conditions, currently available direct and indirect monitoring methods can effectively detect contamination that may leak from hazardous waste facilities into the vadose zone. Direct techniques examined include soil-core monitoring and soil-pore liquid monitoring. Indirect techniques examined include soil-gas monitoring, neutron moderation, complex resistivity, ground-penetrating radar, and electrical resistivity. Properly designed vadose zone monitoring networks can act as a complement to saturated zone monitoring networks at numerous hazardous waste facilities. At certain facilities, particularly those in arid climates where the saturated zone is relatively deep, effective vadose zone monitoring may allow a reduction in the scope of saturated zone monitoring programs.  相似文献   

5.
Electrical and electromagnetic methods are well suited for coastal aquifer studies because of the large contrast in resistivity between fresh water-bearing and salt water-bearing formations. Interpretation models for these aquifers typically contain four layers: a highly resistive unsaturated zone; a surficial fresh water aquifer of intermediate resistivity; an underlying conductive, salt water saturated aquifer; and resistive substratum. Additional layers may be added to allow for variations in lithology within the fresh water and salt water layers. Two methods are evaluated: direct current resistivity and time domain electromagnetic soundings. Use of each method alone produces nonunique solutions for resistivities and/or thicknesses of the different layers. We show that joint inversion of vertical electric and time domain electromagnetic soundings produces a more tightly constrained interpretation model at three test sites than is produced by inversion methods applied to each data set independently.  相似文献   

6.
A small‐diameter nuclear magnetic resonance (NMR) logging tool has been developed and field tested at various sites in the United States and Australia. A novel design approach has produced relatively inexpensive, small‐diameter probes that can be run in open or PVC‐cased boreholes as small as 2 inches in diameter. The complete system, including surface electronics and various downhole probes, has been successfully tested in small‐diameter monitoring wells in a range of hydrogeological settings. A variant of the probe that can be deployed by a direct‐push machine has also been developed and tested in the field. The new NMR logging tool provides reliable, direct, and high‐resolution information that is of importance for groundwater studies. Specifically, the technology provides direct measurement of total water content (total porosity in the saturated zone or moisture content in the unsaturated zone), and estimates of relative pore‐size distribution (bound vs. mobile water content) and hydraulic conductivity. The NMR measurements show good agreement with ancillary data from lithologic logs, geophysical logs, and hydrogeologic measurements, and provide valuable information for groundwater investigations.  相似文献   

7.
Relationships between porosity and hydraulic conductivity tend to be strongly scale- and site-dependent and are thus very difficult to establish. As a result, hydraulic conductivity distributions inferred from geophysically derived porosity models must be calibrated using some measurement of aquifer response. This type of calibration is potentially very valuable as it may allow for transport predictions within the considered hydrological unit at locations where only geophysical measurements are available, thus reducing the number of well tests required and thereby the costs of management and remediation. Here, we explore this concept through a series of numerical experiments. Considering the case of porosity characterization in saturated heterogeneous aquifers using crosshole ground-penetrating radar and borehole porosity log data, we use tracer test measurements to calibrate a relationship between porosity and hydraulic conductivity that allows the best prediction of the observed hydrological behavior. To examine the validity and effectiveness of the obtained relationship, we examine its performance at alternate locations not used in the calibration procedure. Our results indicate that this methodology allows us to obtain remarkably reliable hydrological predictions throughout the considered hydrological unit based on the geophysical data only. This was also found to be the case when significant uncertainty was considered in the underlying relationship between porosity and hydraulic conductivity.  相似文献   

8.
The travel time and amplitude of ground-penetrating radar (GPR) waves are closely related to medium parameters such as water content, porosity, and dielectric permittivity. However, conventional estimation methods, which are mostly based on wave velocity, are not suitable for real complex media because of limited resolution. Impedance inversion uses the reflection coefficient of radar waves to directly calculate GPR impedance and other parameters of subsurface media. We construct a 3D multiscale stochastic medium model and use the mixed Gaussian and exponential autocorrelation function to describe the distribution of parameters in real subsurface media. We introduce an elliptical Gaussian function to describe local random anomalies. The tapering function is also introduced to reduce calculation errors caused by the numerical simulation of discrete grids. We derive the impedance inversion workflow and test the calculation precision in complex media. Finally, we use impedance inversion to process GPR field data in a polluted site in Mongolia. The inversion results were constrained using borehole data and validated by resistivity data.  相似文献   

9.
An improved, iteratively re‐weighted factor analysis procedure is presented to interpret engineering geophysical sounding logs in shallow unsaturated sediments. We simultaneously process cone resistance, electric resistivity, and nuclear data acquired by direct‐push tools to give robust estimates of factor variables and water content in unconsolidated heterogeneous formations. The statistical procedure is based on the iterative re‐weighting of the deviations between the measured and calculated data using the most frequent value method famous for its robustness and high statistical efficiency. The iterative approach improves the result of factor analysis for not normally distributed data and extremely noisy measurements. By detecting a strong regression relation between one of the extracted factors and the fractional volume of water, we establish an independent method for water content estimation along the penetration hole. We verify the estimated values of water volume by using a highly over‐determined, quality‐checked interval inversion procedure. The multidimensional extension of the statistical method allows the estimation of water content distribution along both the vertical and the horizontal coordinates. Numerical tests using engineering geophysical sounding data measured in a Hungarian loessy–sandy formation demonstrate the feasibility of the most frequent value‐based factor analysis, which can be efficiently used for a more reliable hydrogeophysical characterisation of the unsaturated zone.  相似文献   

10.
The determination of clay content in near‐surface formations is crucial for geotechnical, hydrogeological and oil‐contamination studies. We have developed a technique for estimating clay content that consists of the minimization of the difference between the theoretically calculated and measured soil resistivities as a function of water salinity. To calculate the resistivity, we used a model that takes into account the electrochemical processes in the clay micropores. The experimental measurements of soil resistivity were performed on soil samples, completely saturated by brines at different concentrations of NaCl salt in the range 0.6–100 g/l, to obtain the resistivity versus salinity curve. The parameters obtained with this curve inversion are the clay content, the total porosity and the cation exchange capacity. To verify the new technique, we determined clay concentrations of artificial mixtures of calibrated sand and clay. The relative mean error in the clay content does not exceed 20% for a 5% fitting error of the resistivity versus salinity curves. Such evaluations allow the correct separation of the main lithological groups (sand, sandy loam, loam, and light, medium and heavy clay). We applied this technique to estimate the petrophysical parameters of soils (clay content, porosity and cation exchange capacity) at various sites in Mexico. The results improved the interpretation of the vertical electrical soundings, the lithological soil characterization and the delineation of oil‐contaminated areas.  相似文献   

11.
An important quantity in groundwater protection is the residence time of water in an aquifer. It relates to both the travel time of a pollutant to arrive at a well and the time span required for self-purification of a polluted aquifer after removal of pollutant inputs. Time scales for aquifers can be gained from artificial tracer experiments or from environmental tracer data, the latter offering the only realistic alternative if time scales of years or decades have to be taken into account.

Different tracers show different time scales due to their different transport mechanisms especially in the unsaturated zone. While solute tracers are moved advectively with the seepage water, gas tracers pass the unsaturated zone diffusively through the air phase. Depending on the properties of the unsaturated zone (hydraulic properties, thickness) this difference in behavior can be used to separate the subsurface transport process into the unsaturated and the saturated parts.

In a field study in Germany, SF6 and 3H were used as environmental tracers. Both have a relatively well-known input function. Interpretation of data from observation wells by a box model approach led to spatially and temporally varying residence times. This was an indication that the influence of the unsaturated zone could not be neglected. While the gas tracer SF6 shows only residence times in the saturated zone, the tracer 3H reflects the whole travel time of water including both the unsaturated and saturated zones. Using a one-dimensional plug-flow model for the unsaturated zone combined with a detailed two-dimensional flow and transport model for the saturated zone leads to a holistic and consistent interpretation of the measured tracer concentrations. The observed pattern of old water under thick loess cover and younger water under areas where the fractured basalt aquifer crops out is reproduced after adjusting only two parameters: the effective porosity of the saturated aquifer and the product of field capacity and thickness of the unsaturated zone. While the effective porosity of the saturated zone is adjusted by means of the SF6 data, the field capacity of the loess layer is adjusted by means of the 3H observations. The thickness of the unsaturated zone is deduced from geological and pedological maps. All flow data are obtained from a calibrated flow model, which is based on geological data, observed heads and pumping tests only.

The transport model for the saturated zone was calibrated by fitting the porosity by means of gaseous tracer concentrations (SF6). The combined saturated–unsaturated zone model was then calibrated by fitting the field capacity of the unsaturated zone by means of 3H concentrations. With this model it was possible to verify the observed NO3 concentrations at the drinking water wells and to develop predictions for their future development under various scenarios of fertilizer input reduction in specific areas.  相似文献   


12.
Methods for estimating the parameter distributions necessary for modeling fluid flow and contaminant transport in the shallow subsurface are in great demand. Soil properties such as permeability, porosity, and water retention are typically estimated through the inversion of hydrological data (e.g., measurements of capillary pressure and water saturation). However, ill-posedness and non-uniqueness commonly arise in such non-linear inverse problems making their solutions elusive. Incorporating additional types of data, such as from geophysical methods, may greatly improve the success of inverse modeling. In particular, ground-penetrating radar (GPR) methods have proven sensitive to subsurface fluid flow processes and appear promising for such applications. In the present work, an inverse technique is presented which allows for the estimation of flow parameter distributions and the prediction of flow phenomena using GPR and hydrological measurements collected during a transient flow experiment. Specifically, concepts from the pilot point method were implemented in a maximum a posteriori (MAP) framework to allow for the generation of permeability distributions that are conditional to permeability point measurements, that maintain specified patterns of spatial correlation, and that are consistent with geophysical and hydrological data. The current implementation of the approach allows for additional flow parameters to be estimated concurrently if they are assumed uniform and uncorrelated with the permeability distribution. (The method itself allows for heterogeneity in these parameters to be considered, and it allows for parameters of the petrophysical and semivariogram models to be estimated as well.) Through a synthetic example, performance of the method is evaluated under various conditions, and some conclusions are made regarding the joint use of transient GPR and hydrological measurements in estimating fluid flow parameters in the vadose zone.  相似文献   

13.
Archaeoseismology can provide a useful chronological tool for constraining earthquakes and documenting significant evidence that would otherwise be lost. In this paper, we report a case of surface faulting on ancient man-made structures belonging to the archaeological site of Santa Venera al Pozzo situated along the eastern flank of Mt. Etna volcano in eastern Sicily (southern Italy), which is affected by well-developed tectonic faults. Geological surveys highlight a set of fractures affecting the archaeological ruins, suggesting the occurrence of a capable fault zone across the area. An integrated geophysical survey was carried out in order to identify the main subsurface tectonic discontinuity ascribable to the fault zone. The information derived from different geophysical techniques, such as electrical resistivity tomography, seismic refraction tomography, ground-penetrating radar, and magnetic surveys allowed us to infer that the fractures observed at the surface could have been produced by coseismic rupture. They are conceivably linked to a strong earthquake that probably occurred in the Roman period, around mid-end of the third-century AD; time constraints are inferred through the dating of buildings of the archaeological site.  相似文献   

14.
The relationship between aquifer hydraulic conductivity and aquifer resistivity, either measured on the ground surface by vertical electrical sounding (VES) or from resistivity logs, or measured in core samples have been published for different types of aquifers in different locations. Generally, these relationships are empirical and semi-empirical, and confined in few locations. This relation has a positive correlation in some studies and negative in others. So far, there is no potentially physical law controlling this relation, which is not completely understood. Electric current follows the path of least resistance, as does water. Within and around pores, the model of conduction of electricity is ionic and thus the resistivity of the medium is controlled more by porosity and water conductivity than by the resistivity of the rock matrix. Thus, at the pore level, the electrical path is similar to the hydraulic path and the resistivity should reflect hydraulic conductivity. We tried in this paper to study the effect of degree of groundwater saturation in the relation between hydraulic conductivity and bulk resistivity via a simple numerical analysis of Archie’s second law and a simplified Kozeny-Carmen equation. The study reached three characteristic non-linear relations between hydraulic conductivity and resistivity depending on the degree of saturation. These relations are: (1) An inverse power relation in fully saturated aquifers and when porosity equals water saturation, (2) An inverse polynomial relation in unsaturated aquifers, when water saturation is higher than 50%, higher than porosity, and (3) A direct polynomial relation in poorly saturated aquifers, when water saturation is lower than 50%, lower than porosity. These results are supported by some field scale relationships.  相似文献   

15.
Electromagnetic geophysical methods, such as ground-penetrating radar (GPR), have proved to be optimal tools for detecting and mapping near-surface contaminants. GPR has the capability of mapping the location of hydrocarbon pools on the basis of contrasts in the effective permittivity and conductivity of the subsoil. At radar frequencies (50 MHz to 1 GHz), hydrocarbons have a relative permittivity ranging from 2 to 30, compared with a permittivity for water of 80. Moreover, their conductivity ranges from zero to 10 mS/m, against values of 200 mS/m and more for salt water. These differences indicate that water/hydrocarbon interfaces in a porous medium are electromagnetically 'visible'. In order to quantify the hydrocarbon saturation we developed a model for the electromagnetic properties of a subsoil composed of sand and clay/silt, and partially saturated with air, water and hydrocarbon. A self-similar theory is used for the sandy component and a transversely isotropic constitutive equation for the shaly component, which is assumed to possess a laminated structure. The model is first verified with experimental data and then used to obtain the properties of soils partially saturated with methanol and aviation gasoline. Finally, a GPR forward-modelling method computes the radargrams of a typical hydrocarbon spill, illustrating the sensitivity of the technique to the type of pore-fluid. The model and the simulation algorithm provide an interpretation methodology to distinguish different pore-fluids and to quantify their degree of saturation.  相似文献   

16.
Three controlled experiments were conducted at the Oregon Graduate Institute (OGI) with the purpose of evaluating electrical resistance tomography for imaging underground processes associated with in-situ site assessment and remediation. The OGI facilities are unique: a double-wall tank 10 m square and 5 m deep, filled with river bottom sediments and instrumented for geophysical and hydrological studies. At this facility, liquid contaminants could be released into the confined soil at a scale sufficiently large to represent real-world physical phenomena.In the first test, images of electrical resistivity were made before and during a controlled spill of gasoline into a sandy soil. The primary purpose was to determine if electrical resistivity images could detect the hydrocarbon in either the vadose or saturated zone. Definite changes in electrical resistivity were observed in both the vadose and saturated soils. The effects were an increase in resistivity of as much as 10% above pre-release values. A single resistive anomaly was imaged, directly below the release point, principally within the vadose zone but extending below the phreatic surface. The anomaly remained identifiable in tomograms taken two days after the release ended with clear indications of lateral spreading along the water table.The second test involved electrical resistance measurements before, during, and after air sparging in a saturated soil. The primary purpose was to determine if the electrical images could be used to detect and delineate the extent of the zone influenced by sparging. The images showed an increase of about 20% in resistivity over background values within the sparged zone and the extent of the imaged zone agreed with that inferred from other information.Electrical resistivity tomography measurements were made under a simulated oil storage tank in the third test. Comparison of images taken before and during separate releases of brine and water showed effects of changes induced by the water or brine. The simulated leak and its location were imaged as a conductive anomaly centered near the point of origin and were observed to spread with time during the release.  相似文献   

17.
Abstract

The resolution of the freshwater and saline water aquifers in a coastal terrain (Mahanadi Basin, India) is updated. We analysed electrical borehole log data at four sites and compared the water resistivity regime of the freshwater and saline water zones obtained from electrical borehole logging, with the resistivity regime obtained by interpreting vertical electrical sounding (VES) data. The multilayer VES data interpretation is modified to a simple model, containing only the freshwater zone and the saline water zone. The composite geophysical parameters of the freshwater and saline water zones, in particular the resistivity and longitudinal unit conductance regime, are identified. The resolution obtained from the composite geophysical data analyses is very clear and convincing. The composite longitudinal unit conductance regime of the saline water zones is very high compared to that of the freshwater zones. This makes the identification of the two aquifers easy and increases its reliability. A technique which enables analysis of composite geophysical data of freshwater and saline water zones at VES sites in the vicinity of the borehole log sites is proposed. The significance of longitudinal unit conductance in resolving the freshwater and saline water aquifers is illustrated graphically. The proposed technique is validated by correlating the longitudinal unit conductance and resistivity with the total dissolved solids. The efficiency of the technique is validated by carrying out discriminant function analysis.

Citation Hodlur, G. K., Dhakate, R., Sirisha, T. & Panaskar, D. B. (2010) Resolution of freshwater and saline water aquifers by composite geophysical data analysis methods. Hydrol. Sci. J. 55(3), 414–434.  相似文献   

18.
本文运用高密度电法对宁夏西吉县西南山区典型的黄土地震滑坡进行了探测,并结合钻探资料进行验证分析,目的是查明滑坡区域的地层结构、黄土厚度、基岩埋深、富水地段以及空间展布等特征.结果表明:电阻率参数能够很好的反应研究区的地层分布特征,表层疏松干燥的黄土为高阻反映,随着埋深增加,含水量较大的黄土呈现出低阻反映,在黄土与泥岩接触带饱水区呈现低阻反映;滥泥河流域典型的黄土梁斜坡具有阳坡黄土沉积薄、富水层薄、基岩埋深浅的特征,而在阴坡则表现出黄土沉积厚、富水层厚、基岩埋深大的特征.探测结果可有效的应用于黄土地震滑坡的勘察,为进一步开展黄土地震滑坡的成因机理和防治研究提供参考数据.  相似文献   

19.
Integrated geophysical and chemical study of saline water intrusion   总被引：3，自引：0，他引：3  
Choudhury K  Saha DK 《Ground water》2004,42(5):671-677
Surface geophysical surveys provide an effective way to image the subsurface and the ground water zone without a large number of observation wells. DC resistivity sounding generally identifies the subsurface formations-the aquifer zone as well as the formations saturated with saline/brackish water. However, the method has serious ambiguities in distinguishing the geological formations of similar resistivities such as saline sand and saline clay, or water quality such as fresh or saline, in a low resistivity formation. In order to minimize the ambiguity and ascertain the efficacy of data integration techniques in ground water and saline contamination studies, a combined geophysical survey and periodic chemical analysis of ground water were carried out employing DC resistivity profiling, resistivity sounding, and shallow seismic refraction methods. By constraining resistivity interpretation with inputs from seismic refraction and chemical analysis, the data integration study proved to be a powerful method for identification of the subsurface formations, ground water zones, the subsurface saline/brackish water zones, and the probable mode and cause of saline water intrusion in an inland aquifer. A case study presented here illustrates these principles. Resistivity sounding alone had earlier failed to identify the different formations in the saline environment. Data integration and resistivity interpretation constrained by water quality analysis led to a new concept of minimum resistivity for ground water-bearing zones, which is the optimum value of resistivity of a subsurface formation in an area below which ground water contained in it is saline/brackish and unsuitable for drinking.  相似文献   

20.
Hydrogeophysical methods are presented that support the siting and monitoring of aquifer storage and recovery (ASR) systems. These methods are presented as numerical simulations in the context of a proposed ASR experiment in Kuwait, although the techniques are applicable to numerous ASR projects. Bulk geophysical properties are calculated directly from ASR flow and solute transport simulations using standard petrophysical relationships and are used to simulate the dynamic geophysical response to ASR. This strategy provides a quantitative framework for determining site‐specific geophysical methods and data acquisition geometries that can provide the most useful information about the ASR implementation. An axisymmetric, coupled fluid flow and solute transport model simulates injection, storage, and withdrawal of fresh water (salinity ～500 ppm) into the Dammam aquifer, a tertiary carbonate formation with native salinity approximately 6000 ppm. Sensitivity of the flow simulations to the correlation length of aquifer heterogeneity, aquifer dispersivity, and hydraulic permeability of the confining layer are investigated. The geophysical response using electrical resistivity, time‐domain electromagnetic (TEM), and seismic methods is computed at regular intervals during the ASR simulation to investigate the sensitivity of these different techniques to changes in subsurface properties. For the electrical and electromagnetic methods, fluid electric conductivity is derived from the modeled salinity and is combined with an assumed porosity model to compute a bulk electrical resistivity structure. The seismic response is computed from the porosity model and changes in effective stress due to fluid pressure variations during injection/recovery, while changes in fluid properties are introduced through Gassmann fluid substitution.  相似文献   

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

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