首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
 Experiments on degassing of water-saturated granite melts with a pressure drop from 100 and 450 MPa to 40 and 120 MPa, respectively, at temperatures close to feldspar liquidus (750–700  °C), were carried out to determine the modality of water exsolution and vesicle formation at the liquidus temperature. Pressure-drop rates as small as approximately 100 bar/day were used. Uniform space distributions of bubbles of exsolved water were obtained with starting glass containing a small fraction (≈0.5 vol.%) of trapped air bubbles. Volume crystallization of feldspar was observed in degassed melts supplied with seeds. Bubble size distributions (BSD) measured in granite glasses after degassing are presented. Data on vesicle characteristics (number, radius, area, elongation) were acquired on images digitized with standard software, while the reconstruction of size distributions was performed with the Schwartz-Saltikov "unfolding" procedure. Bubble size distributions of size classes in the range 5–1000 μm were acquired with proper magnification and satisfactory statistical reliability of determined number densities. The BSDs of the experimental samples are compared with the results of measurements of rapidly degassed products of Mt. Etna and Vulcano Island. Many particular features of the bubble nucleation and growth can be distinguished in an individual BSD. However, the general BSD of the whole data set, including natural ones, can be relatively well described with linear regression in bilogarithmic coordinates. The slope of this regression is approximately 2.8±0.1. This dependence is in striking contrast with distributions theoretically predicted with classical nucleation models based on homogeneous nucleation of vesicles. The theoretical distribution requires the occurrence of strong maxima that are not observed in our experimental and natural samples, thus arguing for heterogeneous nucleation mechanisms. Received: 1 October 1998 / Accepted: 25 June 1999  相似文献   

2.
Saturated hydraulic conductivity (K) is one of the most important parameters determining groundwater flow and contaminant transport in both unsaturated and saturated porous media. Although several well‐established laboratory methods exist for determining K, in situ measurements of this parameter remain very complex and scale dependent. Often, the limited accessibility of subsurface sediments for sampling means an additional impediment to our ability to quantify subsurface K heterogeneity. One potential solution is the use of outcrops as analogues for subsurface sediments. This paper investigates the use of air permeameter measurements on outcrops of unconsolidated sediments to quantify K and its spatial heterogeneity on a broad range of sediment types. The Neogene aquifer in northern Belgium is used as a case study for this purpose. To characterize the variability in K, 511 small‐scale air permeability measurements were performed on outcrop sediments representative over five of the aquifer's lithostratigraphic units. From these measurements, outcrop‐scale equivalent K tensors were calculated using numerical upscaling techniques. Validation of the air permeameter‐based K values by comparison with laboratory constant head K measurements reveals a correlation of 0.93. Overall, the results indicate that hand‐held air permeameters are very efficient and accurate tools to characterize saturated K, as well as its small‐scale variability and anisotropy on a broad range of unconsolidated sediments. The studied outcrops further provided a qualitative understanding of aquifer hydrostratigraphy and quantitative estimates about K variability at the centimetre‐scale to metre‐scale. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

3.
In this paper we have developed a new method for measuring in situ soil permeability, which is based on the theory of radial gas advection through an isotropic porous medium. The method was tested in the laboratory and at several locations on the island of Vulcano (Aeolian Islands, Italy). It consists of a special device which generates a gas source at a depth of 50 cm and it permits measurement of the relative induced pressure in nearby soil at different depths. The characteristic error of the method was less than 10%. Furthermore, soil permeability measurements were carried out in the island of Vulcano during different periods of the year (between May 2000 and June 2001). A strong decrease in permeability in the upper layers of the soil during and after rainfall was noted, with very poor correlations between the spatial distributions of soil CO2 flux and shallow soil permeability.  相似文献   

4.
The gas permeability of volcanic rocks may influence various eruptive processes. The transition from a quiescent degassing dome to rock failure (fragmentation) may, for example, be controlled by the rocks permeability, in as much as it affects the speed by which a gas overpressure in vesicles is reduced in response to decompression. Using a modified shock-tube-based fragmentation bomb (Alidibirov and Dingwell 1996a,b; Spieler et al. 2003a), we have measured unsteady-state permeability at a high initial pressure differential. Following sudden decompression above the rock cylinder, pressurized gas flows through the sample. Two pressure transducers record the pressure signals above and below the sample. A transient 1D filtration code has been developed to calculate permeability using the experimental decay curve of the lower pressure transducer. Additionally an analytical steady-state method to achieve permeability is presented as an alternative to swiftly predict the sample permeability in a sufficiently precise manner. Over 100 permeability measurements have been performed on samples covering a wide range of porosity. The results show a general positive relationship between porosity and permeability with a high data scatter. Our preferred interpretation of the results is a combination of two different, but overlapping effects. We propose that at low porosities, gas escape occurs predominantly through microcracks or elongated micropores and therefore could be described by simplified forms of Kozeny–Carman relations (Carman 1956) and fracture flow models. At higher porosities, the influence of vesicles becomes progressively stronger as they form an increasingly connected network. Therefore, a model based on the percolation theory of fully penetrable spheres is used, as a first approximation, to describe the permeability-porosity trend. In the data acquired to date it is evident, that in addition to the porosity control, the samples bubble size, shape and distribution strongly influence the permeability. This leads to a range of permeability values up to 2.5 orders of magnitude at a given porosity.  相似文献   

5.
Cold room physical modelling of periglacial solifluction processes on an experimental slope of 12° is described, and data on soil temperatures, surface frost heave, thaw consolidation, downslope soil movement and porewater pressures over seven freeze–thaw cycles are presented. These data are analyzed in the context of laboratory determination of the rheometry of the experimental soils at high moisture contents. It is concluded that the observed thaw-induced solifluction represents pre-failure soil shear strain and results from loss of strength due to the combined effects of raised porewater pressures during thaw consolidation and upward seepage pressures as water flows towards the surface away from the thaw front. An investigation of the rheometry of thawing soils offers the prospect of an analytical model to predict rates and depths of periglacial solifluction. © 1997 by John Wiley & Sons, Ltd.  相似文献   

6.
Quantitative relationships between major fault parameters from geological observations and laboratory experiments are compared. Relationships are established between fault length, number, depth of fault penetration, amplitude of displacement and other characteristics. The width of destruction zones is estimated. Spacing between parallel faults of compatible length is evaluated. It is shown that there is a stable correlation between fault length and number, which is independent from the mode of loading the material under destruction. Destruction of the lithosphere is believed to occur according to the laws of deformation of Maxwell elasto-viscous body.  相似文献   

7.
Four runs of experimental landform development, with the same uplift rate, different rainfall intensity, and the same material of different permeability adjusted by the degree of compaction, showed complicated effects of rainfall and mound-forming material. In the run with more rainfall on less permeable material, low separated ridges developed in the uplifted area, because abundant overland flow promoted valley erosion and slope processes from early stages. In the run with less rainfall on less permeable material, valley incision proceeded mostly in major valleys where surface water converges. Canyons developed during early stages and later a high massive mountain emerged. The effect of rainfall difference, however, appeared completely opposite on more permeable material accompanied by lower shear strength. In the run with more rainfall on more permeable material, a massive mountain similar to that with less rainfall on less permeable material appeared, and low separated ridges appeared in the run with less rainfall on more permeable material as in the run with more rainfall on less permeable material. In the former case, similar amount of water available for Hortonian overland flow in early stages estimated from rainfall rate and permeability can explain the development of similar landforms. In the latter case, while abundant surface water with more rainfall on less permeable material made fluvial erosion active from early stages, the deficiency in surface water with less rainfall on more permeable material apparently attenuated fluvial erosion but possibly accentuated slope processes and slope failures by seepage water flow through more permeable material of low shear strength. The active erosion from early stages apparently resulted in the development of enduring similar low landforms later in the dynamic equilibrium stage. These experimental results indicate that similar landforms can emerge from different environmental and lithologic controls, and that process does not necessarily follow from form.  相似文献   

8.
This petrologic analysis of the Negra Muerta Volcanic Complex (NMVC) contributes to understanding the magmatic evolution of eruptive centres associated with prominent NW-striking fault zones in the southern Central Andes. Specifically, the geochemical characteristics and magmatic evolution of the two eruptive episodes of this Complex are analysed. The first one occurred as an explosive eruption at 9 Ma and is represented by a strongly welded, fiamme-rich, andesitic to dacitic ignimbrite deposit. The second commenced with an eruption of a rhyolitic ignimbrite at 7.6 Ma followed by effusive discharge of hybrid lavas at 7.3 Ma and by emplacement of andesitic to rhyodacitic dykes and domes. Both explosive and effusive eruptions of the second episode occurred within a short time span, but geochemical interpretations permit consideration of the existence of different magmas interacting in the same magma chamber. Our model involves an andesitic recharge into a partially cooled rhyolitic magma chamber, pressurising the magmatic system and triggering explosive eruption of rhyolitic magma. Chemical or mechanical evidence for interaction between the rhyolitic and andesitic magma in the initial stages are not obvious because of their difference in composition, which could have been strong enough to inhibit the interaction between the two magmas. After the initial explosive stages of the eruption at 7.6 Ma, the magma chamber become more depressurised and the most mafic magma settled in compositional layers by fractional crystallisation. Restricted hybridisation occurred and was effective between adjacent and thermally equivalent layers close to the top of the magma chamber. At 7.3 Ma, increments of caldera formation were accompanied by effusive discharge of hybrid lavas through radially disposed dykes whereby andesitic magma gained in importance toward the end of this effusive episode in the central portion of the caldera. Assimilation during turbulent ascent (ATA) is invoked to explain a conspicuous reversed isotopic signature (87Sr/86Sr and 143Nd/144Nd) in the entire volcanic series. Therefore, the 7.6 to 7.3 Ma volcanic rocks of the NMVC resulted from synchronous and mutually interacting petrological processes such as recharge, fractional crystallization, hybridisation, and Assimilation during Turbulent Ascent (ATA).Geochemical characteristics of both volcanic episodes show diverse type and/or depth in the sources and variable influence of upper crustal processes, and indicate a recurrence in the magma-forming conditions. Similarly, other minor volcanic centres in the transversal volcanic belts of the Central Andes repeated their geochemical signatures throughout the Miocene.  相似文献   

9.
The relationship between permeability and vesicularity in volcanic rocks has been used to infer the degassing behavior of hydrous magma. Recent data on natural samples from various eruptions show a wide variation, fitting a power–law relationship of the percolation models with low (< 30%) critical vesicularity (ФC). In this study, we present data on permeability and pore-connectivity of juvenile rhyolitic pumice clasts in a pyroclastic flow around Onikobe volcano, NE Japan, and investigate their relationship with vesicularity developed in a single eruption event. The permeability of the pumices having a relatively low abundance of microlites and microphenocrysts shows a trend increasing by 4 orders of magnitude (from 10− 13.8 to 10− 10.1 m2) in a high and narrow vesicularity range (from 72 to 80%). This trend intersects at a high angle with the fit to the permeability–vesicularity data in the previous studies that has a low ФC, and is located on the extension of the trend for the products of isotropic decompression experiments. The two-dimensional (2D) connectivities of pores for the pumices were also measured from thin sections. From the point of view of percolation theory, connectivity provides information about the probability of percolation. They showed a steep increase from ca. 0 to 0.7 in an almost similar vesicularity range, as compared to their permeabilities. We attribute the increase in 2D connectivity to the increasing amount of ruptured bubble walls, which might have provided less-tortuous paths through larger apertures for gas flow. This, in turn, would cause an effective increase in the permeability. Aggregates of bubble-wall-shaped glass shards were found in the pumices, and their amount and degree of welding are higher in the pumices that have a higher abundance of microlites and microphenocrysts. These pumices have relatively high permeability and 2D connectivity at low vesicularity, which is accounted for by the existence of large irregularly shaped pores. These textural characteristics suggest that a series of partial fragmentation processes, including local rupturing of bubble walls and subsequent foam-collapse with permeable gas flow, might have occurred before the ultimate bulk fragmentation, thus resulting in the increase in permeability. We suggest that the 2D connectivity of pores is a useful parameter to quantify the degree of fragmentation of bubble walls and has the potential for use to assess their permeability.  相似文献   

10.
Colors of plinian pumices were measured by spectrocolorimetry, and their quantitative color parameters in the L*a*b* color space were determined. A series of heating experiments of obsidian was conducted to simulate the color-change processes of rhyolitic glasses. In these experiments, following three stages of color-change processes were observed. Stage I showed a rapid b* (yellowishness) increase associated with fast dehydration controlled by water diffusivity (D water). In stage II, a* (reddishness) increase was accompanied by Fe2+ decrease. Both a* increase and Fe2+ decrease can be simulated by a diffusion model. Obtained diffusivity D oxidation were about two orders of magnitude smaller than D water . The a*-value increase after the oxidation in stage III appeared to be quasi-linear with time, indicating the zeroth order reaction corresponding to the formation of hematite-like structures in rhyolitic glasses. The diffusion-limited a* increase model in stage II was applied to a natural plinian pumice fall unit to evaluate time periods of color-change processes through oxidation by air of fragmented rhyolitic materials.  相似文献   

11.
The distribution of the metal ions Zn and Pb between particulate and dissolved phase in river Glatt was studied by field measurements and compared with calculated simulations, using parameters obtained by adsorption experiments with natural suspended particulate material. Differences in distribution coefficients obtained from field data are observed in function of the sampling locations and of the composition of the particulate matter.Experiments in which metal ion solutions are titrated with a suspension of natural particles and analyzed by anodic stripping voltammetry, are interpreted in terms of binding capacities and conditional stability constants of Zn and Pb with the surface sites. Binding constants of a particular metal ion varied very little for all samples. We obtained mean values for the conditional average complex formation constants at pH 8 of: logcond K Pb = 9.44 ± 0.18 and logcond K Zn = 8.17 ± 0.20. At this pH, binding capacities of 5 10–3 – 1.7 10–2 mol/kg of particles were obtained for samples collected at different locations and times; organic material, iron and manganese oxides are considered to be the main components that control the adsorption to the particles.Distribution coefficients are calculated from the experimentally obtained binding capacities and conditional stability constants. Calculated distribution coefficients for Zn agree with those obtained from the field data and are not very sensitive to changes in the composition of the solution. Good agreement was obtained for lead as well; for some samples it was important to take two types of sites with different affinity into consideration.  相似文献   

12.
Magnetotelluric (MT) measurements were conducted at Iwate volcano, across the entirety of the mountain, in 1997, 1999, 2003, 2006, and 2007. The survey line was 18 km in length and oriented E–W, comprising 38 measurements sites. Following 2D inversion, we obtained the resistivity structure to a depth of 4 km. The surface resistive layer (~ several hundreds of meters thick) is underlain by extensive highly conductive zones. Based on drilling data, the bottom of the highly conductive zone is interpreted to represent the 200 °C isotherm, below which (i.e., at higher temperatures) conductive clay minerals (smectite) are rare. The high conductivity is therefore mainly attributed to the presence of hydrothermally altered clay. The focus of this study is a resistive body beneath the Onigajo (West-Iwate) caldera at depths of 0.5–3 km. This body appears to have impeded magmatic fluid ascent during the 1998 volcanic unrest, as inferred from geodetic data. Both tectonic and low-frequency earthquakes are sparsely distributed throughout this resistive body. We interpret this resistive body as a zone of old, solidified intrusive magma with temperatures in excess of 200 °C. Given that a similar relationship between a resistive body and subsurface volcanic activity has been suggested for Asama volcano, structural controls on subsurface magmatic fluid movement may be a common phenomenon at shallow levels beneath volcanoes.  相似文献   

13.
14.
The effects of afforestation on water resources are still controversial. The aim of this paper is to (i) analyse the hydrological response of an afforested area in the Central Pyrenees and (ii) compare the hydrological response of an afforested area with the response observed in a natural undisturbed forest. The Araguás catchment was cultivated until the 1950s, and then afforested with pines in an effort to control the active degradation processes. The hydrological response was variable and complex, because the discharge was generated by a combination of distinct runoff processes. The hydrological response showed that (i) afforestation produced moderate peak discharges, stormflows and recession limbs, and long rising limbs; (ii) no one single variable was able to explain the hydrological response: rainfall volume and intensity did not explain the hydrological response and antecedent rainfall and initial discharge (indicating antecedent moisture conditions) did play an important role; (iii) seasonal differences were observed suggesting different runoff generating processes; and (iv) the effect of forest cover on peak discharges became less important as the size of the hydrological event increased. The stormflow coefficient showed a clear seasonal pattern with an alternation between a wet period, when the catchment was hydrologically responsive, and a dry summer period when the catchment rarely responded to any event, and two transitional periods (wetting and drying). Compared with a natural forest, the afforested area recorded greater flows and peak discharges, faster response times and shorter recession limbs. Afforestation reduces the water yield and the number of floods compared with non‐vegetated areas and abandoned lands. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

15.
In this study we performed three categories of steady- and unsteady-state core-flooding experiments to investigate capillary trapping, relative permeability, and capillary pressure, in a scCO2 + SO2/brine/limestone system at elevated temperature and pressure conditions, i.e., 60 °C and 19.16 MPa. We used a Madison limestone core sample acquired from the Rock Springs Uplift in southwest Wyoming. We carried out two sets of steady-state drainage-imbibition relative permeability experiments with different initial brine saturations to study hysteresis. We found that the final scCO2 + SO2 drainage relative permeability was very low, i.e., 0.04. We also observed a rapid reduction in the scCO2-rich phase imbibition relative permeability curve, which resulted in a high residual trapping. The results showed that between 62.8% and more than 76% of the initial scCO2 + SO2 at the end of drainage was trapped by capillary trapping mechanism (trapping efficiency). We found that at higher initial brine saturations, the trapping efficiency was higher. The maximum initial and residual scCO2-rich phase saturations at the end of primary drainage and imbibition were 0.525 and 0.329, respectively. Each drainage-imbibition cycle was followed by a dissolution process to re-establish Sw = 1. The dissolution brine relative permeabilities for both cycles were also obtained. We characterized the scCO2 + SO2/brine capillary pressure hysteresis behavior through unsteady-state primary drainage, imbibition, and secondary drainage experiments. We observed negative imbibition capillary pressure curve indicative of possible wettability alteration throughout the experiments due to contact with scCO2 + SO2/brine fluid system. The trapping results were compared to those reported in literature for other carbonate core samples. We noticed slightly more residual trapping in our sample, which might be attributed to heterogeneity, different viscosity ratio, and pore-space topologies. The impact of dynamic effects, i.e., high brine flow rate imbibition tests, on trapping of the scCO2-rich phase was also explored. We performed two imbibition experiments with relatively high brine flow rates. The residual scCO2 saturation dropped to 0.291 and 0.262 at the end of the first and second imbibition tests, i.e., 11.5% and 20.4%, respectively, compared to 0.329 under capillary-dominated regime.  相似文献   

16.
Mountainous river basins are one of the main sources of sediment. Over long time scales, sediment production is sustained by the persistent dissection of river basins, which is promoted by tectonic activity. The response or adjustment of rivers to forcing factors such as uplift is based on the concept of the graded river and a feedback mechanism between the incision and uplift. Although the development of graded rivers under natural circumstances has been discussed for a long time, knowledge about the transition of river basins under heterogeneous uplift is not enough. To understand the development of a river basin with a non‐uniform uplift rate, two simple cases are examined: landward and seaward tilting uplift, where the uplift rate varies linearly in space. For our study, laboratory experiments were conducted and the results were compared with those of natural river basins; two river basins in Yakushima Island were selected for this purpose. In both the laboratory and Yakushima, the longitudinal profile of the river basin under landward‐tilting uplift has a convex‐up zone and a specific knickpoint is formed at the upstream end of this zone. This knickpoint is inactive with respect to migration and incision owing to the insufficient cumulative uplift to the equilibrium state. It was also observed in both the experimental and natural cases that the profile of the river basin under seaward tilting is unlikely to have such a convex‐up zone in the long term, and will instead have a smooth concave profile. Therefore, the spatiotemporal pattern of dissection differs depending on the type of tilting uplift, which suggests that sediment production also varies in time and space according to the type of uplift.  相似文献   

17.
The first stage of alteration of a basalt glass by seawater has been studied experimentally between 200° and 320°C under vapor pressure, following both the chemical evolution of the reacting solutions during the experiments and the chemical and mineralogical transformations of the altered glass surface via microprobe, scanning transmission electron microscope and resonant nuclear reactions.The alteration is controlled by diffusion of dissolved species through an altered layer and not by surface chemical reaction. Selective removal of cations with respect to SiO2 leads to the formation of a thick, porous protonated surface layer enriched in SiO2. Aqueous diffusion of dissolved species from the reacting glass and from solution through the pores of this layer results in the formation of an amorphous silicate gel which is the precursor of normal crystalline clays (saponites) appearing in a later stage.Computer simulations of the alteration generated with the EQ3/6 computer software package account for the chemical composition of the rim and the sequence of secondary minerals found in the experiments.  相似文献   

18.
The present work investigates the influence of small geometrical defects on the behavior of slender rigid blocks. A comprehensive experimental campaign was carried out on one of the shake tables of CEA/Saclay in France. The tested model was a massive steel block with standard manufacturing quality. Release, free oscillations tests as well as shake table tests revealed a non‐negligible out‐of‐plane motion even in the case of apparently plane initial conditions or excitations. This motion exhibits a highly reproducible part for a short duration that was used to calibrate a numerical geometrically asymmetrical model. The stability of this model when subjected to 2000 artificial seismic horizontal bidirectional signals was compared with the stability of a symmetrical one. This study showed that the geometrical imperfections slightly increase the rocking and overturning probabilities for earthquake signals in a narrow range of peak ground acceleration. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

19.
This paper aims to identify the spatial distribution of exchangeable base cations in soils on an acid hillslope and to investigate possible cation release processes from slope soils to the stream. The basic assumption underlying this research is that the amount of exchangeable cations in soils reflects the nutrient stores and cation leaching processes across the slope where vegetation and parent materials are similar. The distribution of exchangeable Ca2+, Mg2+, K+ and Na+ has been investigated on a three-dimensional hillslope on the Quantock Hills, Somerset, UK. A two-way ANOVA shows that soil depth is predominant in explaining the total variance of exchangeable bases, despite the steep slope gradient and clear podzolic catena development. Major nutrient base cations, such as Ca2+, Mg2+ and K+, display homogeneous topsoil storage right across the slope. This spatial pattern may indicate that the spatial distribution of major nutrient cations is tightly controlled by the soil–vegetation system in nutrient-poor heathland environments. Na+ is an exception to this vegetation-controlled spatial distribution, because of its small involvement in the soil–vegetation and soil exchangeable systems. In subsurface soils, cations liberated from the soil–vegetation system are subject to redistribution over the slope according to the hydrological flowpaths operating on the slope, with some eventually released into the stream. The saturated wedge developed at the base of the slope plays a key role in the storage and release processes of base cations from slope soils to the stream. Ca2+, Mg2+ and Na+ carried by throughflow are stored in the saturated wedge and gradually released into the stream at times of high flow. K+, however, shows an apparently different spatial behaviour, being deficient in the saturated wedge. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

20.
Knowledge and understanding of shore platform erosion and tidal notch development in the tropics and subtropics relies mainly on short‐term studies conducted on recently deposited carbonate rocks, predominantly Holocene and Quaternary reef limestones and aeolianites. This paper presents erosion rates, measured over a 10 year period on notches and platforms developed on the Permian, Ratburi limestone at Phang Nga Bay, Thailand. In so doing it contributes to informing a particular knowledge gap in our understanding of the erosion dynamics of shore platform and tidal notch development in the tropics and subtropics – notch erosion rates on relatively hard, ancient limestones measured directly on the rock surface using a micro‐erosion meter (MEM) over time periods of a decade or more. The average intertidal erosion rate of 0.231 mm/yr is lower than erosion rates measured over 2–3 years on recent, weaker carbonate rocks. Average erosion rates at Phang Nga vary according to location and site and are, in rank order from highest to lowest: Mid‐platform (0.324 mm/yr) > Notch floor (0.289 mm/yr) > Rear notch wall (0.228 mm/yr) > Lower platform (0.140 mm/yr) > Notch roof (0.107 mm/yr) and Supratidal (0.095 mm/yr). The micro‐relief of the eroding rock surfaces in each of these positions exhibits marked differences that are seemingly associated with differences in dominant physical and bio‐erosion processes. The results begin to help inform knowledge of longer term shore platform erosion dynamics, models of marine notch development and have implications for the use of marine notches as indicators of changes in sea level and the duration of past sea levels. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

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

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