首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 968 毫秒
1.
The atmospheric chloride mass balance (CMB) method allows spatial evaluations of the average diffuse aquifer recharge by rainfall () in large and varied territories when long‐term steady conditions can be assumed. Often, the distributed average CMB variables necessary to calculate have to be estimated from the available variable‐length data series, which may be of suboptimal quality and spatial coverage. This paper explains the use of these data and the reliability of the results in continental Spain, chosen as a large and varied territory. The CMB variables have been regionalized by ordinary kriging at the same 4976 nodes of a 10 km × 10 km grid. Nodal values vary from 14 to 810 mm year–1, 90% ranging from 30 to 300 mm year–1. The recharge‐to‐precipitation ratios vary from 0.03 in low‐permeability formations and semiarid areas to 0.65 in some carbonate massifs. Integrated average results for the whole of continental Spain yield a potential aquifer recharge of 64 km3 year?1, the net recharge over permeable formations (40% of the territory) being 32 km3 year?1. Two main sources of uncertainty affecting (given by the coefficient of variation, CV), induced by the inherent natural variability of the variables (CVR) and from mapping (), have been segregated. The average CVR is 0.13 and could be improved with longer data series. The average is 0.07 and may be decreased with better data coverage. The estimates were compared with other regional and local recharge estimates, being 4% and 1% higher, respectively. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

2.
Release of nitrogen compounds into groundwater, particularly those compounds from excessive agricultural fertilization, is a major concern in an aquifer recharge. Among the nitrogen compounds, ammonium ( ) is a common one. In order to assess the risk of agricultural fertilizer contamination to an aquifer through infiltration, adsorption onto a loamy agricultural soil profile (0–0.60 m depth) was studied using a soil column experiment and modelling simulation. The soil used in the experiment was drawn from an agricultural field in Xinzhen, Fangshan district, Beijing, China, and reconstituted in laboratory soil columns. Column experiments were conducted using bromide (conservative tracer) and ‐bearing aqueous solutions. The ammonium concentrations in the soil water samples were measured, and their values were plotted as the breakthrough curves. The chemical's soil–water distribution coefficients (Kd) were calculated using breakthrough curves. Then the retardation factor (R) in saturated soil was calculated. For the ‐bearing aqueous solutions, the strongest adsorption occurred at the soil depth of 0.30–0.45 m. The convection–dispersion equation model and chemical non‐equilibrium model in Hydrus‐1D were used to simulate transport in the loamy soil. The two‐site chemical non‐equilibrium model in Hydrus‐1D was best to simulate transport through the soil column. Parameter sensitivity study was conducted to investigate the influences of solute transport by Kd, the fraction of exchange sites assuming to be in equilibrium with the solution phase (f), the longitudinal dispersivity (λ), and the first‐order rate coefficients (ω). The sensitivity analysis results indicate Kd is the most critical parameter.  相似文献   

3.
D. Markovic  M. Koch 《水文研究》2015,29(7):1806-1816
Hydrological processes commonly exhibit long‐term persistence, also known as the ‘Hurst phenomenon’. Here, we examine long‐term precipitation and streamflow time series from the Elbe River Basin to quantify differences in the spectral properties and in the Hurst parameter estimates () of the individual hydrological cycle components. Precipitation‐runoff modelling is performed for the Elbe River sub‐catchment Striegis using the Soil and Water Assessment Tool (SWAT). For 38 daily 50 years long streamflow time series from the Elbe River Basin, baseflow separation and spectral analysis is performed. The results show a spectral shift towards low‐frequency scales (>2 years) from precipitation to baseflow, with a parallel increase of from 0.52 (precipitation) to 0.65 (baseflow). The SWAT model is able to reproduce both, the main low‐frequency mode (≈7 yr.) and the (0.62) of the observed Striegis River flow time series. The baseflow appears to be the main component which shapes the low‐frequency response and of streamflow in the Elbe River Basin to the input precipitation. This conclusion is further confirmed through PMWIN‐MODFLOW groundwater modelling of a hypothetic phreatic stream‐connected aquifer system that consists of various soils (sand, loamy sand and silt). A power shift towards lower frequencies and an increase of for the hydraulic heads is obtained, as the aquifer's lateral dimensions increase and its hydraulic conductivity decreases. The average of the groundwater heads is 0.80, 0.90 and 1.0 for sand, loamy sand and silt aquifers, respectively. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

4.
Evapotranspiration (ET) can cause diel fluctuations in the elevation of the water table and the stage in adjacent streams. The diel fluctuations of water levels change head gradients throughout the day, causing specific discharge through near‐stream sediment to fluctuate at the same time scale. In a previous study, we showed that specific discharge controls the residence time of groundwater in streambed sediment that, in turn, exerted the primary control on removal from groundwater passing through the streambed. In this study, we examine the magnitude of diel specific discharge patterns through the streambed driven by ET in the riparian zone with a transient numerical saturated–unsaturated groundwater flow model. On the basis of a first‐order kinetic model for removal, we predicted diel fluctuations in stream concentrations. Model results indicated that ET drove a diel pattern in specific discharge through the streambed and riparian zone (the removal zones). Because specific discharge is inversely proportional to groundwater travel time through the removal zones and travel time determines the extent of removal, diel changes in ET can result in a diel pattern in concentration in the stream. The model predictions generally matched observations made during summertime base‐flow conditions in a small coastal plain stream in Virginia. A more complicated pattern was observed following a seasonal drawdown period, where source components to the stream changed during the receding limb of the hydrograph and resulted in diel fluctuations being superimposed over a multi‐day trend in concentrations. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

5.
Although changes in rainfall characteristics have been noted across the world, few studies have reported those in mountainous areas. This study was undertaken to clarify spatial and temporal variations in rainfall characteristics such as annual rainfall amount (Pr), mean daily rainfall intensity (η), and ratio of rain days (λ) in mountainous and lowland areas in Taiwan. To this aim, we examined spatial and year‐to‐year variations and marginal long‐term trends in Pr, η, and λ, based on rainfall data from 120 stations during the period 1978–2008. The period mean rainfall () at the lowland stations had strong relationships with the period mean daily rainfall intensity () and the period mean ratio of rain days () during those 31 years. Meanwhile, was only strongly related to at mountainous stations, indicating that influences on spatial variations in were different between lowland and mountainous stations. Year‐to‐year variations in Pr at each station were primarily determined from the variation in η at most stations for both lowland and mountainous stations. Long‐term trend analysis showed that Pr and η increased significantly at 10% and 31% of the total 120 stations, respectively, and λ decreased significantly at 6% of the total. The increases in Pr were mostly accompanied by increases in η. Although stations with significant η increases were slightly biased toward the western lowland area, increases or decreases in Pr and λ were not common. These results contribute to understanding the impacts of possible climate changes on terrestrial hydrological cycles. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

6.
A study of the changes in the ionic loads of NO, NH, SO and H+ in a boreal forest snowpack at Lake Laflamme, Québec was carried out using hydrological and chemical data from field lysimeters. The results showed that depletion of the N-containing species occurs periodically in the snowpack during meltwater discharge. Rain-on-snow events led to in-pack losses of NO and NH at a rate of 130 μeq m?2 day?1 and 101·3 μeq m?2day?1 respectively. On dry days, however, dry deposition and deposition of organic debris from the canopy resulted in increases of 183·3 μeq m?2day?1 for NO and 4·5 μeq m?2day?1 for NH in the pack. In contrast, SO42? showed continual in-pack increases due to deposition of 5·0 μeq m?2day?1 for wet days and 92·6 μeq m?2day?1 for dry days. The depletion of NO and NH is due to microbiological uptake of these nutrients during periods when the free water content of the pack is high. Controlled melts in a laboratory snowmelt simulator containing snow and organic matter from the forest canopy at Lake Laflamme showed losses of NO and NH similar to those observed in the field. As the microbiological uptake proceeds at a rate comparable to that of ionic load increases in the pack by dry deposition, models of the chemical dynamics of snowmelt should take the former into account in any system where organic content of the snowpack is appreciable.  相似文献   

7.
  相似文献   

8.
Errors in the kinematic wave and diffusion wave approximations for time-independent (or steady-state) cases of channel flow were derived for three types of boundary conditions: zero flow at the upstream end, and critical flow depth and zero depth gradient at the downstream end. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors in the range 1–2% for values of KF (? 7.5), where K is the kinematic wave number and F0 is the Froude number. Even for small values of KF (e.g. KF20 = 0.75), the errors were typically less than 15%. The accuracy of the diffusion wave approximation was greatly influenced by the downstream boundary condition. The error of the kinematic wave approximation was found to be less than 13% in the region 0.1 ? x ? 0.95 for KF = 7.5 and was greater than 30% for smaller values of KF (? 0.75). This error increased with strong downstream boundary control.  相似文献   

9.
Soils in post‐wildfire environments are often characterized by a low infiltration capacity with a high degree of spatial heterogeneity relative to unburned areas. Debris flows are frequently initiated by run‐off in recently burned steeplands, making it critical to develop and test methods for incorporating spatial variability in infiltration capacity into hydrologic models. We use Monte Carlo simulations of run‐off generation over a soil with a spatially heterogenous saturated hydraulic conductivity (Ks) to derive an expression for an aerially averaged saturated hydraulic conductivity ( ) that depends on the rainfall rate, the statistical properties of Ks, and the spatial correlation length scale associated with Ks. The proposed method for determining is tested by simulating run‐off on synthetic topography over a wide range of spatial scales. Results provide a simplified expression for an effective saturated hydraulic conductivity that can be used to relate a distribution of small‐scale Ks measurements to infiltration and run‐off generation over larger spatial scales. Finally, we use a hydrologic model based on to simulate run‐off and debris flow initiation at a recently burned catchment in the Santa Ana Mountains, CA, USA, and compare results to those obtained using an infiltration model based on the Soil Conservation Service Curve Number.  相似文献   

10.
We investigate fracture‐induced attenuation anisotropy in a cluster of events from a microseismic dataset acquired during hydraulic fracture stimulation. The dataset contains 888 events of magnitude ?3.0 to 0.0. We use a log‐spectral‐amplitude‐ratio method to estimate change in over a half‐hour time period where fluid is being injected and an increase in fracturing from S‐wave splitting analysis has been previously inferred. A Pearson's correlation analysis is used to assess whether or not changes in attenuation with time are statistically significant. P‐waves show no systematic change in during this time. In contrast, S‐waves polarised perpendicular to the fractures show a clear and statistically significant increase with time, whereas S‐waves polarised parallel to the fractures show a weak negative trend. We also compare between the two S‐waves, finding an increase in with time. A poroelastic rock physics model of fracture‐induced attenuation anisotropy is used to interpret the results. This model suggests that the observed changes in t* are related to an increase in fracture density of up to 0.04. This is much higher than previous estimates of 0.025 ± 0.002 based on S‐wave velocity anisotropy, but there is considerably more scatter in the attenuation measurements. This could be due to the added sensitivity of attenuation measurement to non‐aligned fractures, fracture shape, and fluid properties. Nevertheless, this pilot study shows that attenuation measurements are sensitive to fracture properties such as fracture density and aspect ratio.  相似文献   

11.
Small bounded plots were used to investigate dissolved nutrient losses in overland flow in eucalyptus and pine forests in northern Portugal following understorey fires. Losses of NON, POCa2+, Mg2+ and K+ were measured over 19 months during the first 2–3 years after fire. Solute losses in overland flow increased after fire in eucalyptus and pine terrain due to the mineralization of litter and vegetation and increased overland flow. The data suggest that (i) elevated losses persist for at least 2 and 3 years at the pine and eucalyptus sites respectively and (ii) soluble losses of POK+ in a post‐fire cycle may exceed those adsorbed to eroded material. Losses of POalthough relatively small, are potentially the most detrimental to soil fertility and forest productivity because of the limited opportunity for P replenishment and correspondingly high P losses adsorbed to eroded sediment. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

12.
Processes driving carbonate diagenesis in islands of the northern Bahamas are investigated using major ion, dissolved oxygen and dissolved organic carbon analyses of water samples from surface and ground waters, and measurements of soil gas P. Meteoric waters equilibrate with aragonite, but reactions are water controlled rather than mineral‐controlled and drive dissolution rather than concurrent precipitation of calcite. Surface runoff waters equilibrate with atmospheric P and rapidly recharge the vadose zone, limiting subaerial bedrock dissolution to only 6·6–15 mg l?1 Ca. P of soil gas measured in the summer wet season ((7·4 ± 3·7) × 10?3 atm) is elevated compared with that of the atmosphere, despite the thin skeletal organic nature of the soil and the discontinuous soil cover. Soil waters retained in surface pockets are equilibrated with respect to aragonite and have dissolved 51 ± 19 mg l?1 Ca. This is substantially less than the 93 ± 18 mg l?1 Ca in samples from pumping boreholes that sample meteoric waters from the freshwater lens. The high P of the freshwater lens ((16 ± 8·3) × 10?3 atm for pumping boreholes) suggests that significant additional CO2 may be derived by oxidation of soil‐ and surface‐derived organic carbon within the lens. The suboxic nature of the majority of the freshwater lens and the observed depletion in sulphate support this suggestion, and indicate that both aerobic and anaerobic oxidation may take place. Shallow lens samples from observation boreholes are calcite supersaturated and have a lower P than deeper lens waters, indicating that CO2 degasses from the water table, driving precipitation of calcite cements. We suggest that the geochemical evolution of waters in the vadose zone and upper part of the freshwater lens may be determined by the presence of a body of ground air with P controlled by production in the freshwater lens and soil and by degassing to the atmosphere. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

13.
Our objective was to discover the effect of variations in fluid properties and fracture geometry on the velocity of seismic wave propagation in fluid‐saturated media with parallel planar fractures. We used numerical models calculated by analytical solutions to examine the behaviour of P‐wave phase velocity dispersion in the normal direction to layering, in non‐porous and porous media with planar fractures. We also examined the anisotropy of low frequency phase and group velocities of fast and slow P‐waves and angular‐dependent reflection coefficients in media with planar fractures, under conditions of saturation by fluids with varying bulk moduli, densities, and fracture apertures. We defined several parametre , , and characterising dispersion, characterising anisotropy, characterising the difference between fast and slow modes, and R0 and characterising reflection. Our results show that the behaviour of dispersion shows wider stopbands in the case of gas saturation. Concavity indicator of dispersion for gas saturation was greater than that for liquid saturation and is usually greater than one. Anisotropy is more sensitive to bulk modulus contrast than to density contrast between the solid and the fluid, and is more sensitive to density contrast than to bulk modulus contrast. The case of gas saturation usually had a greater negative R0 and a greater value of compared with those of brine and heavy and light oil saturations. Our results are helpful in distinguishing fluid types saturating geophysical fractures and estimating the aperture and spacing of planar fractures. In seismic exploration, bulk modulus and fluid density can provide useful information in distinguishing among brine, oil, and gas; fracture geometry is important to estimate the permeability of reservoirs.  相似文献   

14.
Es wird eine Fermentationsvorrichtung für die Messung des Sauerstoffverbrauchs bei wählbaren Sauerstoffpartialdrucken vorgestellt, die in einem geschlossenen Gaskreislauf einen Laboratoriumsfermentor mit Meß- und Regeleinrichtungen sowie Möglichkeiten der Probenahme aufweist. Die Druckregression durch die Kohlendioxidadsorption steuert die synchron registrierte Nachlieferung von Sauerstoffgas. Durch diskontinuierliche Kultivation eines adaptierten Belebtschlammes wurden von einem Braunkohlenpyrolyseabwasser die Sauerstoffverbrauchskurven aufgenommen. Es wurde ein polyauxischer Kurvenverlauf erhalten, wobei auf den Abbau der niederen Fettsauren, im wesentlichen Azetat, 70 … 80% des Gesamtsauerstoffverbrauchs entfielen. Hierfür werden biokinetische Daten bestimmt; ; ; μmax = 0,56 h–1 (ausgewählte Werte). Die bei der Erfassung und Errechnung der Daten auftretenden Probleme wurden diskutiert.  相似文献   

15.
T. H. Brikowski 《水文研究》2015,29(7):1746-1756
Adaptation and mitigation efforts related to global trends in climate and water scarcity must often be implemented at the local, single‐catchment scale. A key requirement is understanding the impact of local climate and watershed characteristics coupled with these regional trends. For surface water, determination of multi‐parameter runoff elasticities is a promising tool for achieving such understanding, as explored here for two surface‐water dependent basins in Texas. The first basin is the water supply for Dallas‐Ft. Worth (DFW), and exhibits relatively high precipitation elasticity (proportional change in runoff to change in precipitation) εP = 2.64, and temperature elasticity εT = ? 0.41. Standard precipitation–temperature elasticity diagrams exhibit unusual concave contours of runoff change, indicating influence of additional parameters, which can be isolated using multi‐parameter approaches. The most influential local parameter in DFW is unexpected reduced runoff fraction in cooler wetter years. Those years exhibit increased summer (JJA) precipitation fraction, but predominant cracking soils in DFW minimize JJA runoff, yielding negative . A comparative basin near Houston shows positive , reflecting the local impact of tropical cyclones and lesser abundance of cracking soils. Both basins exhibit positive elasticity to 1‐year previous precipitation (e.g. DFW εP ? 1 = 1.24), reflecting the influence of soil moisture storage. Only DFW exhibits negative elasticity to 2‐year previous precipitation (εP ? 2 = ? 0.65), reflecting multi‐year influence of vegetation growth and increased evapotranspiration. Using these elasticities, analysis of historical multi‐decadal climate departures for DFW indicates the 80% decrease in runoff during the 1950–1957 drought of record was primarily caused by reduced precipitation. Runoff 56% above‐normal during an unprecedented 1986–1998 wet period was primarily caused by increased precipitation. Since 2000, despite precipitation slightly above normal, runoff has decreased 20%, primarily in response to ～ 1°C warming. Future precipitation droughts superimposed on this new drier normal are likely to be much more severe than historical experience would indicate. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

16.
The nonlinearity of the seismic amplitude‐variation‐with‐offset response is investigated with physical modelling data. Nonlinearity in amplitude‐variation‐with‐offset becomes important in the presence of large relative changes in acoustic and elastic medium properties. A procedure for pre‐processing physical modelling reflection data is enacted on the reflection from a water‐plexiglas boundary. The resulting picked and processed amplitudes are compared with the exact solutions of the plane‐wave Zoeppritz equations, as well as approximations that are first, second, and third order in , , and . In the low angle range of 0°–20°, the third‐order plane‐wave approximation is sufficient to capture the nonlinearity of the amplitude‐variation‐with‐offset response of a liquid‐solid boundary with , , and ρ contrasts of 1485–2745 m/s, 0–1380 m/s, and 1.00–1.19 gm/cc respectively, to an accuracy value of roughly 1%. This is in contrast to the linear Aki–Richards approximation, which is in error by as much as 25% in the same angle range. Even‐order nonlinear corrective terms are observed to be primarily involved in correcting the angle dependence of , whereas the odd‐order nonlinear terms are involved in determining the absolute amplitude‐variation‐with‐offset magnitudes.  相似文献   

17.
A simple relationship is proposed in this paper to construct damage‐based inelastic response spectra including the effect of ground motion duration that it can be used for damage control in seismic design of structures. This relation is established for three groups of ground motions with short‐duration, moderate‐duration, and long‐duration ranges. To develop the model, the duration effect is included in the cyclic ductility of structures by an energy‐based method, and then strength reduction factors are computed based on this modified ductility (named ). The strength reduction factors were calculated for 44 stiffness‐degrading oscillators having vibration periods between 0.05 and 4.0 s, four ultimate ductility capacities, and five damage levels subjected to 296 earthquake records. The results showed that ductility capacity, damage level, and ground motion duration are effective parameters in the energy dissipation of structures, which affect the spectra. The values of short‐period oscillators (e.g., low‐rise structures) under short‐duration records are generally greater than those under moderate‐duration and long‐duration records. Residual analysis has been made in terms of magnitude and distance to examine the validity of the proposed simple expression. Finally, the introduced spectra were compared with three previously published proposals. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

18.
N. Subba Rao 《水文研究》2012,26(22):3344-3350
A pollution index of groundwater (PIG) is proposed for quantification of water contamination. PIG quantifies the status of concentrations of water quality measures with respect to their drinking water quality standards. The validity of the proposed index is verified by choosing the data of groundwater quality of the Varaha River Basin (Visakhapatnam District, Andhra Pradesh, India) as a case study. The computed index from the study area varies from 0.83 to 2.55. The index disseminates the area into zones of insignificant (PIG <1.0), low (PIG: 1.0 to 1.5), moderate (PIG: 1.5 to 2.0), high (PIG 2.0 to 2.5) and very high (PIG >2.5) pollution. Insignificant pollution zone is observed from the upstream area, where the groundwater is dominated by , and very high pollution zone from the downstream area, where the groundwater is associated with Cl?. This indicates that the quality of groundwater in the study area is mainly influenced by the source of geogenic origin, but it is subsequently modified by the effects of anthropogenic and marine sources. Geochemical ratios (Na+ : Cl?, : Cl?, Na+ : Ca2+ and Mg2+ : Ca2+) also form the quantitative basis of the index. The present study paves the way to implement appropriate management strategies at a specific site to circumvent the pollution. As the classification of the pollution zones with PIG depends upon the drinking water quality standards, it becomes a universal assessment tool for groundwater contamination at any test area. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

19.
Linear prediction filters are an effective tool for reducing random noise from seismic records. Unfortunately, the ability of prediction filters to enhance seismic records deteriorates when the data are contaminated by erratic noise. Erratic noise in this article designates non‐Gaussian noise that consists of large isolated events with known or unknown distribution. We propose a robust f‐x projection filtering scheme for simultaneous erratic noise and Gaussian random noise attenuation. Instead of adopting the ?2‐norm, as commonly used in the conventional design of f‐x filters, we utilize the hybrid ‐norm to penalize the energy of the additive noise. The estimation of the prediction error filter and the additive noise sequence are performed in an alternating fashion. First, the additive noise sequence is fixed, and the prediction error filter is estimated via the least‐squares solution of a system of linear equations. Then, the prediction error filter is fixed, and the additive noise sequence is estimated through a cost function containing a hybrid ‐norm that prevents erratic noise to influence the final solution. In other words, we proposed and designed a robust M‐estimate of a special autoregressive moving‐average model in the f‐x domain. Synthetic and field data examples are used to evaluate the performance of the proposed algorithm.  相似文献   

20.
Under certain circumstances, seismic propagation within porous media may be associated to the conversion of mechanical energy to electromagnetic energy, which is known as a seismo‐electromagnetic phenomenon. The propagation of fast compressional P‐waves is more specifically associated to the manifestations of a seismoelectric field linked to the fluid flows within the pores. The analysis of seismoelectric phenomena, which requires the combination of the theory of electrokinetics and Biot's theory of poroelasticity, provides us with transfer function that links the coseismic seismoelectric field E to the seismic acceleration . To measure the transfer function, we have developed an experimental setup enabling seismoelectric laboratory observation in unconsolidated quartz sand within the kilohertz range. The investigation focused on the impact of fluid conductivity and water saturation over the coseismic seismoelectric field. During the experiment, special attention was given to the accuracy of electric field measurements. We concluded that, to obtain a reliable estimate of the electric field amplitude, the dipole from which the potential differences are measured should be of much smaller length than the wavelength of the propagating seismic field. Time‐lapse monitoring of the seismic velocities and seismoelectric transfer functions were performed during imbibition and drainage experiments. In all cases, the quantitative analysis of the seismoelectric transfer function was in good agreement with theoretical predictions. While investigating saturation variations from full to residual water saturation, we showed that the ratio undergoes a switch in polarity at a particular saturation , which also implies a sign change of the filtration, traducing a reversal of the relative fluid displacement with respect to the frame. This sign change at critical saturation stresses a particular behaviour of the poroelastic medium: the dropping of the coseismic electric field to zero traduces the absence of relative pore/fluid displacements representative of a Biot dynamically compatible medium. We concluded from our experimental study in loose sand that the measurements of the coseismic seismoelectric coupling may provide information on fluid distribution within the pores and that the reversal of the seismoelectric field may be used as an indicator of the dynamically compatible state of the medium.  相似文献   

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

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