首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Studies investigating the effects of inland recharge on coastal groundwater dynamics were carried out typically in unconfined aquifers, with few in confined aquifers. This study focused on the groundwater dynamics in confined aquifers with seasonally sinusoidally fluctuated inland groundwater head and constant sea level by numerical simulations. It is known that the mixing zone (MZ) of saltwater wedge in response to the seasonal oscillations of inland groundwater head swings around the steady-state MZ. However, our simulation results indicate that even the most landward freshwater-saltwater interface over a year is seaward from the steady-state location when the hydraulic conductivity K is ≤10−4 m/s under certain boundary conditions with given parameter values. That is, seasonal oscillations of inland groundwater head may reduce seawater intrusion in confined coastal aquifers when K ≤ 10−4 m/s. Sensitivity analysis indicates that for aquifers of K ≤ 10−4 m/s, the larger the inland head fluctuation amplitude is, the less the seawater intrudes. This is probably due to the reason that the seawater intrusion time decreases with the increase of fluctuation amplitude when K ≤ 10−4 m/s. Numerical simulations demonstrate that seasonal inland groundwater head oscillations promote the annual averaged recirculated seawater discharge across the seaward boundary.  相似文献   

2.
High transmissivity aquifers typically have low hydraulic gradients (i.e., a flat water table). Measuring low gradients using water levels can be problematic because measurement error may be greater than the true difference in water levels (i.e., a low signal-to-noise ratio). In this study, the feasibility of measuring a hydraulic gradient in the range of 10−6 to 10−5 m/m was demonstrated. The study was performed at a site where the depth to water from land surface ranged from 40.1 to 94.2 m and the aquifer transmissivity was estimated at 41,300 m2/d (hydraulic conductivity of 18,800 m/d). The goals of the study were to reduce measurement error as much as practicable and assess the importance of factors affecting water level measurement accuracy. Well verticality was the largest source of error (0.000 to 0.168 m; median of 0.014 m), and geodetic survey of casing elevations was the next most important source of error (0.002 to 0.013 m; median of 0.005 m). Variability due to barometric pressure fluctuations was not an important factor at the site. Hydraulic heads were measured to an accuracy of ±0.0065 m, and the average hydraulic gradient was estimated to be 8.0 × 10−6 (±0.9 × 10−6) m/m. The improvement in accuracy allowed for two reversals in the groundwater flow direction to be identified, after which the gradient averaged 2.5 × 10−5 (±0.4 × 10−5) m/m. This study showed it is possible to sufficiently control sources of error to measure hydraulic gradients in the 10−6 to 10−5 m/m range.  相似文献   

3.
Glacial aquifers are an important source of groundwater in the United States and require accurate characterization to make informed management decisions. One parameter that is crucial for understanding the movement of groundwater is hydraulic conductivity, K. Nuclear magnetic resonance (NMR) logging measures the NMR response associated with the water in geological materials. By utilizing an external magnetic field to manipulate the nuclear spins associated with 1H, the time-varying decay of the nuclear magnetization is measured. This logging method could provide an effective way to estimate K at submeter vertical resolution, but the models that relate NMR measurements to K require calibration. At two field sites in a glacial aquifer in central Wisconsin, we collected a total of four NMR logs and obtained measurements of K in their immediate vicinity with a direct-push permeameter (DPP). Using a bootstrap algorithm to calibrate the Schlumberger-Doll Research (SDR) NMR-K model, we estimated K to within a factor of 5 of the DPP measurements. The lowest levels of accuracy occurred in the lower-K (K < 10−4 m/s) intervals. We also evaluated the applicability of prior SDR model calibrations. We found the NMR calibration parameters varied with K, suggesting the SDR model does not incorporate all the properties of the pore space that control K. Thus, the expected range of K in an aquifer may need to be considered during calibration of NMR-K models. This study is the first step toward establishing NMR logging as an effective method for estimating K in glacial aquifers.  相似文献   

4.
《Continental Shelf Research》2006,26(17-18):2125-2140
Sediment delivered to coastal systems by rivers (15×109 tons) plays a key role in the global carbon and nutrient cycles, as deltas and continental shelves are considered to be the main repositories of organic matter in marine sediments. The Mississippi River, delivering more than 60% of the total dissolved and suspended materials from the conterminous US, dominates coastal and margin processes in the northern Gulf of Mexico. Draining approximately 41% of the conterminous US, the Mississippi and Atchafalaya river system deliver approximately 2×108 tons of suspended matter to the northern Gulf shelf each year. Unlike previous work, this study provides a comprehensive evaluation of sediment accumulation covering majority of the shelf (<150 m water depth) west of the Mississippi Delta from 92 cores collected throughout the last 15 years. This provides a unique and invaluable data set of the spatial and modern temporal variations of the sediment accumulation in this dynamic coastal environment.Three types of 210Pb profiles were observed from short cores (15–45 cm) collected on the shelf. Proximal to Southwest Pass in 30–100 m water depths, non-steady-state profiles were observed indicating rapid accumulation. Sediment accumulation rates in this area are typically >2.5 cm yr−1 (>1.8 g cm−2 yr−1). Kasten cores (∼200 cm in length) collected near Southwest Pass also indicate rapid deposition (>4 cm yr−1; >3 g cm−2 yr−1) on a longer timescale than that captured in the box cores. Near shore (<20 m), profiles are dominated by sediments reworked by waves and currents with no accumulation (the exception is an area just south of Barataria Bay where accumulation occurs). The remainder of the shelf (distal of Southwest Pass) is dominated by steady-state accumulation beneath a ∼10-cm thick mixed layer. Sediment accumulation rates for the distal shelf are typically <0.7 cm yr−1 (<0.5 g cm−2 yr−1). A preliminary sediment budget based on the distribution of 210Pb accumulation rates indicates that 40–50% of the sediment delivered by the river is transported out of the study region. Sediment is moved to distal regions of the shelf/slope through two different mechanisms. Along-isobath sediment movement occurs by normal resuspension processes west of the delta, whereas delivery of sediments south and southwest of the delta may be also be influenced by mass movement events on varying timescales.  相似文献   

5.
Fluid flow exerts a critical impact on the convection of thermal energy in geological media, whereas heat transport in turn affects fluid properties, including fluid dynamic viscosity and density. The interplay of flow and heat transport also affects solute transport. To unravel these complex coupled flow, heat, and solute transport processes, here, we present a theory for the idealized scale‐dependent Poiseuille flow model considering a constant temperature gradient (?T) along a single fracture, where fluid dynamic viscosity connects with temperature via an exponential function. The idealized scale‐dependent model is validated based on the solutions from direct numerical simulations. We find that the hydraulic conductivity (K) of the Poiseuille flow either increases or decreases with scales depending on ?T > 0°C/m or ?T < 0°C/m, respectively. Indeed, the degree of changes in K depends on the magnitude of ?T and fracture length. The scale‐dependent model provides an alternative explanation for the well‐known scale‐dependent transport problem, for example, the dispersion coefficient increases with travel distance when ?T > 0°C/m according to the Taylor dispersion theory, because K (or equivalently flux through fractures) scales with fracture length. The proposed theory unravels intertwined interactions between flow and transport processes, which might shed light on understanding many practical geophysical problems, for example, geothermal energy exploration.  相似文献   

6.
An important constraint on the reliability of cosmogenic nuclide exposure dating is the rigorous determination of production rates. We present a new dataset for 10Be production rate calibration from Mount Billingen, southern Sweden, the site of the final drainage of the Baltic Ice Lake, an event dated to 11,620 ± 100 cal yr BP. Five samples of flood-scoured bedrock surfaces (58.5°N, 13.7°E, 105–120 m a.s.l.) unambiguously connected to the drainage event yield a reference 10Be production rate of 4.19 ± 0.20 atoms g−1 yr−1 for the CRONUS-Earth online calculator Lm scaling and 4.02 ± 0.18 atoms g−1 yr−1 for the nuclide specific LSDn scaling. We also recalibrate the reference 10Be production rates for four sites in Norway and combine three of these with the Billingen results to derive a tightly clustered Scandinavian reference 10Be production rate of 4.13 ± 0.11 atoms g−1 yr−1 for the CRONUS Lm scaling and 3.95 ± 0.10 atoms g−1 yr−1 for the LSDn scaling scheme.  相似文献   

7.
Volume measurements for magnesiowüstite (Mg0.6Fe0.4)O, were carried out up to pressures of 10.1 GPa in the temperature range 300–1273 K, using energy-dispersive synchrotron X-ray diffraction. These data allow reliable determination of the temperature dependence of the bulk modulus and good constraint on the thermal expansitivity at ambient pressure which was previously not known for magnesiowüstite. From these data, thermal and elastic parameters were derived from various approaches based on the Birch–Murnaghan equation of state (EOS) and on the relevant thermodynamic relations. The results from three different equations of state are remarkably consistent. With (∂KT/∂P)T fixed at 4, we obtained K0=158(2) GPa, (∂KT/∂T)P=−0.029(3) GPa K−1, (∂KT/∂T)V=−3.9(±2.3)×10−3 GPa K−1, and αT=3.45(18)×10−5+1.14(28)×10−8T. The K0, (∂KT/∂T)P, and (∂KT/∂T)V values are in agreement with those of Fei et al. (1992) and are similar to previously determined values for MgO. The zero pressure thermal expansitivity of (Mg0.6Fe0.4)O is found to be similar to that for MgO (Suzuki, 1975). These results indicate that, for the compositional range x=0–0.4 in (Mg1−xFex)O, the thermal and elastic properties of magnesiowüstite exhibit a dependence on the iron content that is negligibly small, within uncertainties of the experiments. They are consequently insensitive to the Fe–Mg partitioning between (Mg, Fe)SiO3 perovskite and magnesiowüstite when applied to compositional models of the lower mantle. With the assumption that (Mg0.6Fe0.4)O is a Debye-like solid, a modified equation of heat capacity at constant pressure is proposed and thermodynamic properties of geophysically importance are calculated and tabulated at high temperatures.  相似文献   

8.
The distribution of extraterrestrial chromite grains (> 63 μm) has been studied in the resurge deposit of the early Late Ordovician (458 Ma) Lockne impact structure in central Sweden. A 1-kg-sample of resurge sediment from the rim of the crater is extremely rich in extraterrestrial chromite. In total the sample contains > 125 chromium spinel grains kg 1, most grains being chromite. Due to post-depositional alterations, the origin of several of the grains is dubious, although it is clear that the major part of the chromites is of extraterrestrial origin (> 75 grains kg 1), most likely derived from the impactor. The alterations are primarily due to the hydrothermal system induced by the impactor, with zinc enrichment in the chromites as the most common characteristic. The element chemistry of the least altered chromite grains indicates that the impactor was an ordinary L chondrite. This concurs with a suggested increase in the flux of L chondritic asteroids to Earth, following the disruption of the L chondrite parent body at ca. 470 Ma.More than 170 impact craters are known from around the Earth, but only few of the impactors have been identified, mainly due to the low survival rate of projectile material. Further studies of impact craters will reveal if survival of extraterrestrial chromite is unique for the Lockne structure.  相似文献   

9.
《Continental Shelf Research》2005,25(9):1081-1095
The mesoscale distribution and seasonal variation of the size structure of phytoplankton biomass, as measured by chlorophyll a (chl a), was studied in the Ebro shelf area (NW Mediterranean) during three different seasons: autumn, winter and summer. In autumn and summer, when the water column was, respectively, slightly or strongly stratified and nutrient concentrations were low at surface, average total chl a values were 0.31 and 0.29 mg m−3, respectively. In winter, the intrusion of nutrients into the photic zone by intense vertical mixing and strong riverine inputs, produced an increase of the total autotrophic biomass (0.76 mg m−3). In the three seasons, the main contributor to total chl a was the picoplanktonic (<2 μm) size fraction (42% in winter and around 60% in autumn and summer). The nanophytoplankton (2–20 μm) contribution to total chl a showed the lowest variability amongst seasons (between 29% and 39%). The microplanktonic (>20 μm) chl a size fraction was higher in winter (27%) than in the other seasons (less than 13%). The maximum total chl a concentrations were found at surface in winter, at depths of 40 m in autumn and between 50 and 80 m in summer. The relative contribution of the <2 μm size fraction at these levels of the water column tended to be higher than at other depths in autumn and winter and lower in summer. In autumn and winter, nutrient inputs from Ebro river discharge and mixing processes resulted in an increase on the >2 μm contribution to total chl a in the coastal zone near the Ebro Delta area. In summer, the contribution of the <2 and >2 μm chl a size fractions was homogeneously distributed through the sampling area. In autumn and summer, when deep chl a maxima were observed, the total amount of the autotrophic biomass in the superficial waters (down to 10 m) of most offshore stations was less than 10% of the whole integrated chl a (down to 100 m or to the bottom). In winter, this percentage increased until 20% or 40%. The >2 μm chl a increased linearly with total chl a values. However, the <2 μm chl a showed a similar linear relationship only at total chl a values lower than 1 mg m−3 (in autumn and summer) or 2 mg m−3 (winter). At higher values of total chl a, the contribution of the <2 μm size fraction remained below an upper limit of roughly 0.5 mg m−3. Our results indicate that the picoplankton fraction of phytoplankton may show higher seasonal and mesoscale variability than is usually acknowledged.  相似文献   

10.
《Geofísica Internacional》2014,53(4):411-423
An investigation of the effect of some human activities on the magnetic susceptibility and frequency dependent susceptibility was conducted on top soil samples from, a commercial area, a motor park and a school environment in Jalingo, Taraba State, N-E Nigeria. The purpose was to assess the variation of magnetic susceptibility with different land use, detect pollution hotspots using magnetic proxy parameters and evaluate the contribution of superparamagnetic (SP) grain size contribution to the magnetic susceptibility from calculation of the frequency dependence of magnetic susceptibility (MS). The results of the mass specific low frequency magnetic susceptibility measurements showed significant enhancement with values ranging from 67.8  495.3 × 10−8 m3kg−1 with a mean value of 191.61 × 10−8 × m3kg−1 for the Jalingo College of Education (JCOE) data; 520.1  1612.8 × 10−8 m3kg−1 with a mean value of 901.34 × 10−8 m3kg−1 for the Jalingo main Market (JMM) and 188.5  1203.6 × 10−8m3kg−1 with an average value of 574 92 × 10−6 m3kg−1 for the Jalingo Motor Park (JMP). The significant magnetic enhancement indicates high concentration of ferrimagnetic minerals in the soil and hence increased pollution. The magnetic susceptibility of the different land use studied decreased in the order commercial area (market) > motor park > school premises. The results of the percentage frequency dependence susceptibility showed that most of the samples had a mixture of SP and coarse multi domain grains or SP grains < 0.05 µm. The value of χfd% range from 2.68 to 13.80% with an average value of 8.67% in the JCOE samples, 0.49 to 10.04% with an a-verage of 5.05% in the JMM samples and 0.56 to 13.04% with an average value of 5.86% in the JMP samples.  相似文献   

11.
12.
《Continental Shelf Research》2007,27(10-11):1477-1500
Constraining the exchange of water from the shoreline to the mid-shelf is necessary for the development of accurate and predictive models of nearshore circulation. Ra isotopes, which emanate from sediments and have a variety of half-lives, may be useful in measuring cross-shelf mixing rates. The distributions of Ra isotopes were measured in transects extending perpendicular from the shoreline at Sunset Beach and Huntington Beach, CA. The average inventory at Sunset Beach was four times greater than at Huntington Beach. Building on previous research on Ra inputs and circulation in San Pedro Bay, a two-dimensional model for surface water Ra was developed to identify the importance of onshore flow and cross-shelf mixing near Huntington Beach. For the mean summertime conditions, the eddy diffusivity (Kh) was 1.4±0.4 m2 s−1, with 8% of the water from Sunset Beach moving down the coast. The remaining water must be low-Ra water that has moved onshore. At time scales greater than a week, the short-lived Ra inventory at Huntington Beach varied by 50%, which reflects changes in the fractions of water moving down-coast and/or in the longshore advection rate. The shoreline Ra concentration varied on time scales of hours, which may be generated by tidal changes in the Ra input at the shoreline and short-period fluctuations in the mixing rate. The low Kh observed in this study in comparison to higher values measured further offshore is evidence that Kh increases with distance offshore. When scale-dependent mixing beyond 455 m offshore is incorporated into the model, the results are consistent with the observed data for 223Ra, 224Ra, and 228Ra. Using the model, the 228Ra input flux to the summertime mixed layer was between 3.4×106 and 4.0×106 atoms s−1 (m shoreline)−1.  相似文献   

13.
At longer periods, scattered ScS waves sometimes dominate over coda waves at large lapse times. Examining recordings of seismic envelopes at 9 IRIS seismic stations of regional earthquakes with focal depths deeper than 150 km in periods from 1 to 20 s for a wide lapse time range up to 2000 s, we found significant frequency dependence. The coda decay gradient at short periods is steeper than that at longer periods; however, the change of coda gradient associated with the ScS arrival becomes distinct as the period becomes longer. In particular, a clear offset of coda amplitude appears in central Asia for 10 and 15 s period bands. The multiple isotropic scattering process of S-waves in the heterogeneous mantle can be simply simulated by using the Monte Carlo simulation method based on the radiative transfer theory in scattering media. Assuming a two-plane-layer attenuation structure and smoothed velocity model of the PREM, we estimated the average total scattering coefficients of S-waves such as 7.52 × 10 4∼1.32 × 10 3 km 1 and 2.08 × 10 4∼6.23 × 10 4 km 1 at 4 s, and 4.51 × 10 4∼7.37 × 10 4 km 1 and 2.80 × 10 5∼2.71 × 10 4 km 1 at 10 s, for the lithosphere and the upper mantle and for the lower mantle, respectively. Our results indicate that scattering occurs mostly in the lithosphere and the upper mantle and support that medium heterogeneity spreads over the whole mantle though its scattering power is small. Strong scattering occurs beneath central Asia and Papua New Guinea, whereas the scattering beneath Italy and regions of east Russia is much weaker. The numerical calculation enables us to confirm that much stronger scattering than intrinsic attenuation causes the offset behavior with coda decay gradient change after the ScS arrival for 4 and 10 s period bands in some regions.  相似文献   

14.
Pressure–volume–temperature relations have been measured to 32 GPa and 2073 K for natural magnesite (Mg0.975Fe0.015Mn0.006Ca0.004CO3) using synchrotron X-ray diffraction with a multianvil apparatus at the SPring-8 facility. A least-squares fit of the room-temperature compression data to a third-order Birch–Murnaghan equation of state (EOS) yielded K0 = 97.1 ± 0.5 GPa and K′ = 5.44 ± 0.07, with fixed V0 = 279.55 ± 0.02 Å3. Further analysis of the high-temperature compression data yielded the temperature derivative of the bulk modulus (∂KT/∂T)P = −0.013 ± 0.001 GPa/K and zero-pressure thermal expansion α = a0 + a1T with a0 = 4.03 (7) × 10−5 K−1 and a1 = 0.49 (10) × 10−8 K−2. The Anderson–Grüneisen parameter is estimated to be δT = 3.3. The analysis of axial compressibility and thermal expansivity indicates that the c-axis is over three times more compressible (KTc = 47 ± 1 GPa) than the a-axis (KTc = 157 ± 1 GPa), whereas the thermal expansion of the c-axis (a0 = 6.8 (2) × 10−5 K−1 and a1 = 2.2 (4) × 10−8 K−2) is greater than that of the a-axis (a0 = 2.7 (4) × 10−5 K−1 and a1 = −0.2 (2) × 10−8 K−2). The present thermal EOS enables us to accurately calculate the density of magnesite to the deep mantle conditions. Decarbonation of a subducting oceanic crust containing 2 wt.% magnesite would result in a 0.6% density reduction at 30 GPa and 1273 K. Using the new EOS parameters we performed thermodynamic calculations for magnesite decarbonation reactions at pressures to 20 GPa. We also estimated stability of magnesite-bearing assemblages in the lower mantle.  相似文献   

15.
Analyses (n = 525) of chloride (Cl), bromide (Br), nitrate as nitrogen (NO3-N), sodium (Na+), calcium (Ca2+) and potassium (K+) in stream water, tile-drain water and groundwater were conducted in an urban-agricultural watershed (10% urban/impervious, 87% agriculture) to explore potential differences in the signature of Cl originating from an urban source as compared with an agricultural source. Only during winter recharge events did measured Cl concentrations exceed the 230 mg/L chronic threshold. At base flow, nearly all surface water and tile water samples had Cl concentrations above the calculated background threshold of 18 mg/L. Mann–Whitney U tests revealed ratios of Cl to Br (p = .045), to NO3-N (p < .0001), to Ca2+ (p < .0001), and to Na+ (p < .0001) to be significantly different between urban and agricultural waters. While Cl ratios indicate that road salt was the dominant source of Cl in the watershed, potassium chloride fertilizer contributed as an important secondary source. Deicing in watersheds where urban land use is minimal had a profound impact on Cl dynamics; however, agricultural practices contributed Cl year-round, elevating stream base flow Cl concentrations above the background level.  相似文献   

16.
We present the first demonstration of hydraulic tomography (HT) to estimate the three-dimensional (3D) hydraulic conductivity (K) distribution of a fractured aquifer at high-resolution field scale (HRFS), including the fracture network and connectivity through it. We invert drawdown data collected from packer-isolated borehole intervals during 42 pumping tests in a wellfield at the former Naval Air Warfare Center, West Trenton, New Jersey, in the Newark Basin. Five additional tests were reserved for a quality check of HT results. We used an equivalent porous medium forward model and geostatistical inversion to estimate 3D K at high resolution (K blocks <1 m3), using no strict assumptions about K variability or fracture statistics. The resulting 3D K estimate ranges from approximately 0.1 (highest-K fractures) to approximately 10−13 m/s (unfractured mudstone). Important estimated features include: (1) a highly fractured zone (HFZ) consisting of a sequence of high-K bedding-plane fractures; (2) a low-K zone that disrupts the HFZ; (3) several secondary fractures of limited extent; and (4) regions of very low-K rock matrix. The 3D K estimate explains complex drawdown behavior observed in the field. Drawdown tracing and particle tracking simulations reveal a 3D fracture network within the estimated K distribution, and connectivity routes through the network. Model fit is best in the shallower part of the wellfield, with high density of observations and tests. The capabilities of HT demonstrated for 3D fractured aquifer characterization at HRFS may support improved in situ remediation for contaminant source zones, and applications in mining, repository assessment, or geotechnical engineering.  相似文献   

17.
The Huancané II moraines deposited by the Quelccaya Ice Cap in southern Peru were selected by the CRONUS-Earth Project as a primary site for evaluating cosmogenic-nuclide scaling methods and for calibrating production rates. The CRONUS-Earth Project is an effort to improve the state of the art for applications of cosmogenic nuclides to earth-surface chronology and processes. The Huancané II moraines are situated in the southern Peruvian Andes at about 4850 m and ∼13.9°S, 70.9°W. They are favorable for cosmogenic-nuclide calibration because of their low-latitude and high-elevation setting, because their age is very well constrained to 12.3 ± 0.1 ka by 34 radiocarbon ages on peat bracketing the moraines, and because boulder coverage by snow or soil is thought to be very unlikely. However, boulder-surface erosion by granular disintegration is observed and a ∼4% correction was applied to measured concentrations to compensate. Samples from 10 boulders were analyzed for 10Be, 26Al, and 36Cl. Interlaboratory bias at the ∼5% level was the largest contributor to variability of the 10Be samples, which were prepared by three laboratories (the other two nuclides were only prepared by one laboratory). Other than this issue, variability for all three nuclides was very low, with standard deviations of the analyses only slightly larger than the analytical uncertainties. The site production rates (corrected for topographic shielding, erosion, and radionuclide decay) at the mean site elevation of 4857 m were 45.5 ± 1.6 atoms 10Be (g quartz)−1 yr−1, 303 ± 15 atoms 26Al (g quartz)−1 yr−1, and 1690 ± 100 atoms 36Cl (g K)−1 yr−1. The nuclide data from this site, along with data from other primary sites, were used to calibrate the production rates of these three nuclides using seven global scaling methods. The traditional Lal formulation and the new Lifton-Sato-Dunai calibrations yield average ages for the Huancané samples that are in excellent-to-good agreement with the radiocarbon age control (within 0.7% for 10Be and 36Cl and 6% for 26Al). However, all of the neutron-monitor-based methods yielded ages that were too young by about 20%. The nuclide production ratios at this site are 6.74 ± 0.34 for 26Al/10Be in quartz and 37.8 ± 2.3 (atoms 36Cl (g K)−1) (atom 10Be (g SiO2)−1)−1 for 36Cl/10Be, in sanidine and quartz, respectively.  相似文献   

18.
Although the bulk moduli (KT0) of silicate melts have a relatively narrow range of values, the pressure derivatives of the isothermal bulk modulus (KT0) can assume a broad range of values and have an important influence on the compositional dependence of the melt compressibility at high pressure. Based on the melt density data from sink/float experiments at high pressures in the literature, we calculate KT0 using an isothermal equation of state (EOS) (e.g., Birch–Murnaghan EOS and Vinet EOS) with the previously determined values of room-pressure density (ρ0) and room-pressure bulk modulus (KT0). The results show that best estimates of KT0 vary considerably from ~ 3 to ~ 7 for different compositions. KT0 is nearly independent of Mg # (molar Mg/(Mg + Fe)), but decreases with SiO2 content. Hydrous melts have anomalously small KT0 leading to a high degree of compression at high pressures. For anhydrous melts, KT0 is ~ 7 for peridotitic melts, ~ 6 for picritic melts, ~ 5 for komatiitic melts, and ~ 4 for basaltic melts.  相似文献   

19.
《Journal of Hydrology》2006,316(1-4):13-27
A linearized approach to quantifying predictive uncertainty in a 2-D model of shallow water flow in response to uncertainty in friction parameterization is presented. The resulting uncertain finite volume (UFV) method is tested against Monte Carlo simulations for uncertain models over channel only, floodplain only and channel and floodplain meshes. The results show that the UFV model performs well in predicting mean and standard deviations of water depths, for problems with two independent Manning's n values, with standard deviations of up to 0.02 m1/3 s−1 with a mean value of 0.03 m1/3 s−1. For depth averaged velocities, mean values are well represented, but standard deviations are poorly predicted by UFV. UFV also performs well when modelling flow over an uneven fractal topography and for a distributed (11 degrees of freedom) parameterization. A computation time advantage of >50 when compared to the Monte Carlo method is observed.  相似文献   

20.
Storage in deep geological formations is a potential solution for the management of high-level radioactive wastes. In this context, different types of rocks such as argillite are extensively studied. In the Mont Terri underground laboratory (Switzerland), several experiments have been performed in order to characterize the properties of the opalinus clay. One of these experiments, called HE-C, has consisted in measuring in situ the time evolution of the rock temperature submitted to a heating source. Experimental measurements have shown that the thermal behaviour of the clay was not homogeneous around the borehole where the heater was installed. Furthermore, 3D direct numerical simulations of this experiment performed with the code Cast3M have proved that it was necessary to introduce a new parameter α to model the amount of electric power lost in cables and by air convection inside the metallic tube containing the heater. A numerical simulation–optimization technique has been used to estimate the thermal longitudinal and transverse conductivities (λ// and λ) of the host rock. It consists in minimizing an objective function that is the sum of the squared differences between measured and calculated temperatures. But this method induced a lot of Cast3M simulations. In order to drastically reduce the CPU time, we used a neural network approximation built from a sample training of 1100 Cast3M simulations. It allowed us to calculate the objective function for 500 000 different values of the triplet (λ//,λ,α).Finally, we obtained the following values for the thermal conductivities
  • –on one side of the borehole, λ// = 1.84 ± 0.04 W m− 1 K− 1 and λ = 0.55 ± 0.03 W m− 1 K− 1;
  • –on the other side, λ// = 1.90 ± 0.07 W m− 1 K− 1 and λ = 1.07 ± 0.09 W m− 1 K− 1.
The estimated thermal conductivities λ perpendicular to the bedding plane are quite different. It is perhaps caused by the presence of an intensive fractured zone on one side of the borehole, due to bentonite swelling. It can also be due to the presence of a bed of carbonated rock.  相似文献   

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

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