首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
  相似文献   

2.
Unlike the studies in small parcels by systematic measurements, the spatial variability of soil properties is expected to increase in those over relatively large areas or scales. Spatial variability of soil hydraulic conductivity (K h) is of significance for the environmental processes, such as soil erosion, plant growth, transport of the plant nutrients in a soil profile and ground water levels. However, its variability is not much and sufficiently known at basin scale. A study of testing the performance of cokriging of K h compared with that of kriging was conducted in the catchment area of Sarayköy II Irrigation Dam in Cank?r?, Turkey. A total of 300 soil surface samples (0–10 cm) were collected from the catchment with irregular intervals. Of the selected soil properties, because the water-stable aggregates (WSA) indicated the highest relationship with the hydraulic conductivity by the Pearson correlation analysis, it is used as an auxiliary variable to predict K h by the cokriging procedure. In addition, the sampling density was reduced randomly to n = 175, n = 150, n = 75 and n = 50 for K h to determine if the superiority of cokriging over kriging would exist. Statistically, the results showed that all reduced K h was as good as the complete K h when its auxiliary relations with WSA were used in cokriging. Particularly, the results of the “Relative Reduction in MSE” (RMSE) revealed that the reduced data set of n = 75 produced the most accurate map than the others. In this basin-scaled study, there was a clear superiority of the cokriging procedure by the reduction in data although a very undulating topography and topographically different aspects, two different land uses with non-uniform vegetation density, different parent materials and soil textures were present in the area. Hence, using the statistically significant auxiliary relationship between K h and WSA might bring about a very useful data set for watershed hydrological researches.  相似文献   

3.
The main objective of this paper was to investigate the dewatering behaviour of a clayey uranium ore slurry. The slurry (containing 28% clay size) exhibited moderate water adsorption (w l  = 83% and w p  = 30%). Primarily composed of muscovite (46%) and quartz (30%), the clay minerals included illite (8%), chlorite (5%) and kaolinite (2%) alongside a CEC of 41 (cmol(+)/kg) with Ca2+ and Mg2+ as the dominant cations. Likewise, the high EC (17,600 μS/cm) and ionic strength (1.15 mol/L) indicated a flocculated microstructure due to the presence of SO4 2? (22,600 mg/L) and Mg2+ (1340 mg/L) in the slurry water. Settling included sedimentation and consolidation at low initial solids condition (25–35%) whereas only consolidation was observed at high initial solids contents (40–50%). The average k reduced from 1.2 × 10?6 m/s (initial s = 25%) to 5.3 × 10?8 m/s (initial s = 50%) along with a void ratio reduction from 7.4 to 2.6. Due to thixotropic strength, volume compressibility during consolidation showed apparent pre-consolidation at low effective stress (0.3–2 kPa) with a reduction in void ratio from 2.6 to 2.5. The e s was found to be 2.46 at σ′ = 2 kPa and was followed by a steeper slope with the void ratio reducing to 2.1 at σ′ = 31 kPa. Likewise, the hydraulic conductivity during consolidation decreased from 2.6 × 10?9 m/s (at e = 2.6) to 2.0 × 10?10 m/s (at e = 2.1).  相似文献   

4.
Non-Darcian flow to a partially penetrating well in a confined aquifer with a finite-thickness skin was investigated. The Izbash equation is used to describe the non-Darcian flow in the horizontal direction, and the vertical flow is described as Darcian. The solution for the newly developed non-Darcian flow model can be obtained by applying the linearization procedure in conjunction with the Laplace transform and the finite Fourier cosine transform. The flow model combines the effects of the non-Darcian flow, partial penetration of the well, and the finite thickness of the well skin. The results show that the depression cone spread is larger for the Darcian flow than for the non-Darcian flow. The drawdowns within the skin zone for a fully penetrating well are smaller than those for the partially penetrating well. The skin type and skin thickness have great impact on the drawdown in the skin zone, while they have little influence on drawdown in the formation zone. The sensitivity analysis indicates that the drawdown in the formation zone is sensitive to the power index (n), the length of well screen (w), the apparent radial hydraulic conductivity of the formation zone (K r2), and the specific storage of the formation zone (S s2) at early times, and it is very sensitive to the parameters n, w and K r2 at late times, especially to n, while it is not sensitive to the skin thickness (r s).  相似文献   

5.
Soil saturated hydraulic conductivity (Ks) is considered as soil basic hydraulic property, and its precision estimation is a key element in modeling water flow and solute transport processes both in the saturated and vadose zones. Although some predictive methods (e.g., pedotransfer functions, PTFs) have been proposed to indirectly predict Ks, the accuracy of these methods still needs to be improved. In this study, some easily available soil properties (e.g., particle size distribution, organic carbon, calcium carbonate content, electrical conductivity, and soil bulk density) are employed as input variables to predict Ks using a fuzzy inference system (FIS) trained by two different optimization techniques: particle swarm optimization (PSO) and genetic algorithm (GA). To verify the derived FIS, 113 soil samples were taken, and their required physical properties were measured (113 sample points?×?7 factors?=?791 input data). The initial FIS is compared with two methods: FIS trained by PSO (PSO-FIS) and FIS trained by GA (GA-FIS). Based on experimental results, all three methods are compared according to some evaluation criteria including correlation coefficient (r), modeling efficiency (EF), coefficient of determination (CD), root mean square error (RMSE), and maximum error (ME) statistics. The results showed that the PSO-FIS model achieved a higher level of modeling efficiency and coefficient of determination (R2) in comparison with the initial FIS and the GA-FIS model. EF and R2 values obtained by the developed PSO-FIS model were 0.69 and 0.72, whereas they were 0.63 and 0.54 for the GA-FIS model. Moreover, the results of ME and RMSE indices showed that the PSO-FIS model can estimate soil saturated hydraulic conductivity more accurate than the GA-FIS model with ME?=?10.4 versus 11.5 and RMSE?=?5.2 versus 5.5 for PSO-FIS and GA-FIS, respectively.  相似文献   

6.
Effective soil thermal conductivity (λ eff) describes the ability of a multiphase soil to transmit heat by conduction under unit temperature gradient. It is a critical parameter for environmental science, earth and planetary science, and engineering applications. Numerous models are available in the literature, but their applicability is generally restricted to certain soil types or water contents (θ). The objective of this study was to develop a new model in the similar form of the Johansen 1975 model to simulate the λ eff(θ) relationship of soils of various soil textures and water contents. An exponential type model with two parameters is developed and a new function for calculating dry soil thermal conductivity is presented. Performance of the new model and six other normalized models were evaluated with published datasets. The results show that the new model is able to well mimic λ eff(θ) relationship of soils from sand to silt loam and from oven dry to full saturation. In addition, it has the best performance among the seven models under test (with root-mean-square error of 0.059 W m?1 °C?1, average deviations of 0.0009 W m?1 °C?1, and Nash–Sutcliffe efficiency of 0.994). The new model has potential to improve the reliability of soil thermal conductivity estimation and be incorporated into numerical modeling for environmental, earth and engineering studies.  相似文献   

7.
DEM simulations of sandstone under true triaxial compressive tests   总被引：1，自引：0，他引：1  
Numerically simulated true triaxial compression tests (σ 1 ≥ σ 2 ≥ σ 3) are conducted in this study to elucidate the failure mechanism of sandstone using 3D discrete element method (DEM), in particular the effect of the intermediate principal stress (σ 2). Eight series of tests (σ 3 = 0, 10, 20, 30, 40, 50, 70, and 100 MPa) are conducted. Within each series, σ 2 is varied from σ 2 = σ 3 to σ 2 = σ 1 from test to test. For each test, σ 1 is raised monotonically to failure while keeping σ 2 and σ 3 constant. The DEM simulations reveal the effect of σ 2 on the variations of peak stress, Young’s modulus, failure plane angles, the brittle–ductile transition, and the evolution of failure modes, the effect beyond the well-understood effect of σ 3. The simulation is in qualitative agreement with the results obtained experimentally. Detailed analyses performed on the particle-scale responses further the understanding of the microscopic mechanisms. The distribution of contact force becomes more homogeneous with the increase of σ 3, which leads to the resulting damage being more localized rather than diffused. The interaction between contact force distribution and coalescence of cracks determines the processes and patterns of fracturing in the sample scale. σ 2 is found to affect the microscopic stress distribution as well as structure evolution, and this effect weakens with the increase of σ 3.  相似文献   

8.
Understanding the changes in permeability of soil, when soil is subjected to high confining pressure and flow pressure, which may alter the textural and geomechanical characteristics of soil, is of great importance to many geo-engineering activities such as, construction of high-rise buildings near the coast or the water bodies, earthen dams, pavement subgrades, reservoir, and shallow repositories. It is now possible to evaluate the changes in permeability of soil samples under varying conditions of confining pressure and flow pressure using flexible wall permeameter (FWP). In the present study, investigation was carried out on a cylindrical sample of granular soil employing FWP under varied conditions of confining pressure (σ3)—50–300 kPa, which can simulate the stress conditions equivalent to depth of about 20 m under the earth’s crust, and a flow pressure (fp)—20–120 kPa, which is mainly present near the small earthen embankment dams, landfill liners, and slurry walls near the soft granular soil with high groundwater table. The obtained results indicate a linear relationship between hydraulic conductivity (k) with effective confining pressure (σeff.), k, decreasing linearly with an incremental change in σeff.. Further, k increases significantly with an increase in fp corresponding to each σeff., and q increases significantly with increase in the fp corresponding to each (σ3). It was also observed that corresponding to the low fp of 20 kPa, the reduction in k is nonlinear with σ3. The percentage reduction in k is observed to be 9, 13, and 27% corresponding to σ3 of 50–100, 100–200, and 200-300 kPa, respectively.  相似文献   

9.
  相似文献   

10.
Hydraulic Conductivity of Fly Ash-Amended Mine Tailings   总被引：1，自引：1，他引：0  
The objective of this study was to evaluate the effect of fly ash addition on hydraulic conductivity (k) of mine tailings. Mine tailings used in this study were categorized as synthetic tailings and natural tailings; two synthetic tailings were developed via blending commercially-available soils and natural tailings were collected from a garnet mine located in the U.S. Two fly ashes were used that had sufficient calcium oxide (CaO) content (17 and 18.9 %) to generate pozzolanic activity. Hydraulic conductivity was measured on pure tailings and fly ash-amended tailings in flexible-wall permeameters. Fly ash was added to mine tailings to constitute 10 % dry mass of the mixture, and specimens were cured for 7 and 28 days. The influence of fly ash-amendment on k of mine tailings was attributed to (1) molding water content and (2) plasticity of the mine tailings. Tailings that classified as low-plasticity silts with clay contents less than 15 % exhibited a decrease in k when amended with fly ash and prepared wet of optimum water content (w opt ). Tailings that classified as low-plasticity clay exhibited a one-order magnitude increase in k with addition of fly ash for materials prepared dry or near w opt . The decrease in k for silty tailings was attributed to formation of cementitious bonds that obstructed flow paths, whereas the increase in k for clayey tailings was attributed to agglomeration of clay particles and an overall increase in average pore size. The results also indicated that the effect of curing time on k is more pronounced during the early stages of curing (≤7 days), as there was negligible difference between k for 7 and 28-days cured specimens.  相似文献   

11.
Recognizing the heterogeneity of hydraulic conductivity and hyporheic flow is critical for understanding contaminant transfer and biogeochemical and hydrological processes involving streams and aquifers. In this study, the heterogeneity of hydraulic conductivity and Darcian flux in a submerged streambed and its adjacent exposed stream banks were investigated in the Beiluo River, northwest China. In the submerged streambed, Darcian flux was estimated by measurement of vertical hydraulic conductivity (K v) and vertical head gradient (VHG) using in-situ permeameter tests. On exposed stream banks, both horizontal hydraulic conductivity (K h) and K v were measured by on-site permeameter tests. In the submerged streambed, K v values gradually decreased with depth and the higher values were concentrated in the center and close to the erosional bank. Compared to the exposed stream banks, the K v values were higher in the streambed. From stream stage to the topmost layer of tested sediment, through increasing elevation, the K h values increased on the erosional bank, while they decreased on the depositional bank. The values of VHG along the thalweg illustrate that downwelling flux occurred in the deepest area while upwelling flux appeared in the other areas, which might result from the change of streambed elevation. The higher value of the Darcian flux in the submerged streambed existed near the erosional bank.  相似文献   

12.
The crystal chemistry and the elastic behavior under isothermal conditions up to 9 GPa of a natural, and extremely rare, 3T-phlogopite from Traversella (Valchiusella, Turin, Western Alps) [(K0.99Na0.05Ba0.01)(Mg2.60Al0.20Fe 0.21 2+ )[Si2.71Al1.29O10](OH)2, space group P3112, with a = 5.3167(4), c = 30.440(2) Å, and V = 745.16(9) ų] have been investigated by electron microprobe analysis in wavelength dispersion mode, single-crystal X-ray diffraction at 100 K, and in situ high-pressure synchrotron radiation powder diffraction (at room temperature) with a diamond anvil cell. The single-crystal refinement confirms the general structure features expected for trioctahedral micas, with the inter-layer site partially occupied by potassium and sodium, iron almost homogeneously distributed over the three independent octahedral sites, and the average bond distances of the two unique tetrahedra suggesting a disordered Si/Al-distribution (i.e., 〈T1-O〉 ~ 1.658 and 〈T2-O〉 ~ 1.656 Å). The location of the H-site confirms the orientation of the O–H vector nearly perpendicular to (0001). The refinement converged with R 1(F) = 0.0382, 846 unique reflections with F O > 4σ(F O) and 61 refined parameters, and not significant residuals in the final difference-Fourier map of the electron density (+0.77/?0.37 e ?/Å3). The high-pressure experiments showed no phase transition within the pressure range investigated. The P–V data were fitted with a Murnaghan (M-EoS) and a third-order Birch-Murnaghan equation of state (BM-EoS), yielding: (1) M-EoS, V 0 = 747.0(3) Å3, K T0 = 44.5(24) GPa, and K′ = 8.0(9); (2) BM-EoS, V 0 = 747.0(3) Å3, K T0 = 42.8(29) GPa, and K′ = 9.9(17). A comparison between the elastic behavior in response to pressure observed in 1M- and 3T-phlogopite is made.  相似文献   

13.
The high-pressure behavior of a vanadinite (Pb10(VO4)6Cl2, a = b = 10.3254(5), c = 7.3450(4) Å, space group P63/m), a natural microporous mineral, has been investigated using in-situ HP-synchrotron X-ray powder diffraction up to 7.67 GPa with a diamond anvil cell under hydrostatic conditions. No phase transition has been observed within the pressure range investigated. Axial and volume isothermal Equations of State (EoS) of vanadinite were determined. Fitting the P–V data with a third-order Birch-Murnaghan (BM) EoS, using the data weighted by the uncertainties in P and V, we obtained: V 0 = 681(1) Å3, K 0 = 41(5) GPa, and K′ = 12.5(2.5). The evolution of the lattice constants with P shows a strong anisotropic compression pattern. The axial bulk moduli were calculated with a third-order “linearized” BM-EoS. The EoS parameters are: a 0 = 10.3302(2) Å, K 0(a) = 35(2) GPa and K′(a) = 10(1) for the a-axis; c 0 = 7.3520(3) Å, K 0(c) = 98(4) GPa, and K′(c) = 9(2) for the c-axis (K 0(a):K 0(c) = 1:2.80). Axial and volume Eulerian-finite strain (fe) at different normalized stress (Fe) were calculated. The weighted linear regression through the data points yields the following intercept values: Fe a (0) = 35(2) GPa for the a-axis, Fe c (0) = 98(4) GPa for the c-axis and Fe V (0) = 45(2) GPa for the unit-cell volume. The slope of the regression lines gives rise to K′ values of 10(1) for the a-axis, 9(2) for the c-axis and 11(1) for the unit cell-volume. A comparison between the HP-elastic response of vanadinite and the iso-structural apatite is carried out. The possible reasons of the elastic anisotropy are discussed.  相似文献   

14.
High pressure in situ synchrotron X-ray diffraction experiment of strontium orthophosphate Sr3(PO4)2 has been carried out to 20.0 GPa at room temperature using multianvil apparatus. Fitting a third-order Birch–Murnaghan equation of state to the P–V data yields a volume of V 0 = 498.0 ± 0.1 Å3, an isothermal bulk modulus of K T  = 89.5 ± 1.7 GPa, and first pressure derivative of K T ′ = 6.57 ± 0.34. If K T ′ is fixed at 4, K T is obtained as 104.4 ± 1.2 GPa. Analysis of axial compressible modulus shows that the a-axis (K a  = 79.6 ± 3.2 GPa) is more compressible than the c-axis (K c  = 116.4 ± 4.3 GPa). Based on the high pressure Raman spectroscopic results, the mode Grüneisen parameters are determined and the average mode Grüneisen parameter of PO4 vibrations of Sr3(PO4)2 is calculated to be 0.30(2).  相似文献   

15.
Thermal inertia (I) is an important parameter in the Earth’s thermal study. There is no doubt that correct and up-to-date knowledge of thermal inertia particularly as it is affected by the land use/cover will provide good and useful information to agriculturists and environmental scientists. In this work, thermal inertia of Abeokuta City of Ogun State, southwestern Nigeria, was determined. Map of the study area was gridded using 2-min resolution which gave 5 × 5 sampling points from where core samples were collected. Bulk density (ρ s ), thermal conductivity (λ s ), and heat capacity (C s ) of the samples were determined. The pattern of variation of the city view of I (in the order of ×103 Jm?2 s?1/2 K?1) showed that the main urban built-up part of the study area, Abeokuta South, had the highest I mean, 0.76160 with standard deviation of 0.032547 and standard error of 0.01455541. The trend also varied along each sampling latitudinal line. It ranges between 0.509 and 0.756 on latitude 7° 14′, 0.557 and 0.768 on latitude 7° 12′, 0.642 and 0.782 on latitude 7° 10′, 0.7 and 0.794 on latitude 7° 08′, and between 0.642 and 0.728 on latitude 7° 06′. Hence, we have estimated thermal inertia of Abeokuta using the thermophysical properties of the study area. Thermal inertia tends to gradually decrease with distance from the city center (Abeokuta South). This result is promising in possible future consideration of urban ground heat energy conversion to other forms of energy.  相似文献   

16.
Subcritical crack growth plays an important role in evaluating the long-term stability of structures in rocks. By applying the constant-displacement-relax method, two groups of test specimens that one immersed in groundwater and the other in air were tested to get the parameters of subcritical crack growth in double torsion test. The relations of the stress intensity factor K I versus the subcritical crack growth velocity V were obtained under the two different environments, and the behavior of subcritical crack growth was also analyzed. The results showed: the relations of lg K I  ? lg V accorded with linear rules, which is in good agreement with Charles theory; Compared with specimens in nature state, the lg K I  ? lg V curves of saturated water specimens locate top left comer of those of air specimens. The slope of curve is smaller, and the intercept is bigger, which shows that the water–rock interaction speeds up the subcritical crack growth. And A increases 2.9 × 1018 folds but n decreases from 85.12 to 40.83 because of the water–rock interaction. Meanwhile, the fracture toughness K IC also decreases obviously from 2.55 in air to 2.26 in water due to water rock interaction. The testing results provide a basis for time-dependence of rock engineering stability.  相似文献   

17.
The electrical conductivity of aqueous fluids containing 0.01, 0.1, and 1 M NaCl was measured in an externally heated diamond cell to 600 °C and 1 GPa. These measurements therefore more than double the pressure range of previous data and extend it to higher NaCl concentrations relevant for crustal and mantle fluids. Electrical conductivity was generally found to increase with pressure and fluid salinity. The conductivity increase observed upon variation of NaCl concentration from 0.1 to 1 M was smaller than from 0.01 to 0.1 M, which reflects the reduced degree of dissociation at high NaCl concentration. Measured conductivities can be reproduced (R 2 = 0.96) by a numerical model with log \(\sigma\) = ?1.7060– 93.78/T + 0.8075 log c + 3.0781 log \(\rho\) + log \(\varLambda\) 0(T, \(\rho\)), where \(\sigma\) is the conductivity in S m?1, T is temperature in K, c is NaCl concentration in wt%, \(\rho\) is the density of pure water (in g/cm3) at given pressure and temperature, and \(\varLambda\) 0 (T, \(\rho\)) is the molar conductivity of NaCl in water at infinite dilution (in S cm2 mol?1), \(\varLambda\) 0 = 1573–1212 \(\rho\) + 537 062/T–208 122 721/T 2. This model allows accurate predictions of the conductivity of saline fluids throughout most of the crust and upper mantle; it should not be used at temperatures below 100 °C. In general, the data show that already a very small fraction of NaCl-bearing aqueous fluid in the deep crust is sufficient to enhance bulk conductivities to values that would be expected for a high degree of partial melting. Accordingly, aqueous fluids may be distinguished from hydrous melts by comparing magnetotelluric and seismic data. H2O–NaCl fluids may enhance electrical conductivities in the deep crust with little disturbance of v p or v p/v s ratios. However, at the high temperatures in the mantle wedge above subduction zones, the conductivity of hydrous basaltic melts and saline aqueous fluids is rather similar, so that distinguishing these two phases from conductivity data alone is difficult. Observed conductivities in forearc regions, where temperatures are too low to allow melting, may be accounted for by not more than 1 wt% of an aqueous fluid with 5 wt% NaCl, if this fluid forms a continuous film or fills interconnected tubes.  相似文献   

18.
In this study, the seismicity rate changes that can represent an earthquake precursor were investigated along the Sagaing Fault Zone (SFZ), Central Myanmar, using the Z value technique. After statistical improvement of the existing seismicity data (the instrumental earthquake records) by removal of the foreshocks and aftershocks and man-made seismicity changes and standardization of the reported magnitude scales, 3574 earthquake events with a M w ≥ 4.2 reported during 1977–2015 were found to directly represent the seismotectonic activities of the SFZ. To find the characteristic parameters specifically suitable for the SFZ, seven known events of M w ≥ 6.0 earthquakes were recognized and used for retrospective tests. As a result, utilizing the conditions of 25 fixed earthquake events considered (N) and a 2-year time window (T w), a significantly high Z value was found to precede most of the M w ≥ 6.0 earthquakes. Therefore, to evaluate the prospective areas of upcoming earthquakes, these conditions (N = 25 and T w = 2) were applied with the most up-to-date seismicity data of 2010–2015. The results illustrate that the vicinity of Myitkyina and Naypyidaw (Z = 4.2–5.1) cities might be subject to strong or major earthquakes in the future.  相似文献   

19.
Study of intact rock failure criteria is an important topic in rock mechanics. In this study, applicability of nine different intact rock failure criteria is investigated for intact coal strength data. PFC3D modeling was used to simulate the laboratory polyaxial tests for cubic intact coal blocks of side dimension 110 mm under different confining stress combinations. A modified grid search procedure is proposed and used to find the best-fitting parameter values and to calculate the coefficient of determination (R 2) values for each criterion. Detailed comparisons of the nine criteria are made using the following aspects: R 2 values, σ 1 ? σ 2 plots for different σ 3, shapes on the deviatoric plane, linearity or nonlinearity on the meridian planes. Through the comparisons of R 2 values, σ 1 ? σ 2 plots and meridian lines, the modified Wiebols–Cook and modified Lade criteria were found to fit the intact coal strength data best. The nine failure criteria are categorized into three types based on the appearances on the deviatoric plane.  相似文献   

20.
Synchrotron-based in situ angle-dispersive X-ray diffraction experiments were conducted on a natural uvite-dominated tourmaline sample by using an external-heating diamond anvil cell at simultaneously high pressures and temperatures up to 18 GPa and 723 K, respectively. The angle-dispersive X-ray diffraction data reveal no indication of a structural phase transition over the P–T range of the current experiment in this study. The pressure–volume–temperature data were fitted by the high-temperature Birch–Murnaghan equation of state. Isothermal bulk modulus of K 0 = 96.6 (9) GPa, pressure derivative of the bulk modulus of \(K_{0}^{\prime } = 12.5 \;(4)\), thermal expansion coefficient of α 0 = 4.39 (27) × 10?5 K?1 and temperature derivative of the bulk modulus (?K/?T) P  = ?0.009 (6) GPa K?1 were obtained. The axial thermoelastic properties were also obtained with K a0 = 139 (2) GPa, \(K_{a0}^{\prime }\) = 11.5 (7) and α a0 = 1.00 (11) × 10?5 K?1 for the a-axis, and K c0 = 59 (1) GPa, \(K_{c0}^{\prime }\) = 11.4 (5) and α c0 = 2.41 (24) × 10?5 K?1 for the c-axis. Both of axial compression and thermal expansion exhibit large anisotropic behavior. Thermoelastic parameters of tourmaline in this study were also compared with that of the other two ring silicates of beryl and cordierite.  相似文献   

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

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