The stability of density-driven flows in saturated heterogeneous porous media     

Jude L. Musuuza  Florin A. Radu  Sabine Attinger 《Advances in water resources》2011,34(11):1464-1482

This work concludes the investigations into the stability of haline flows in saturated porous media. In the first part [33] a stability criterion for density-driven flow in a saturated homogeneous medium was derived excluding dispersion. In the second part [34], the effects of dispersion were included. The latter criterion made reasonable predictions of the stability regimes (indicated by the number of fingers present) as a function of density and dispersivity variations. We found out that destabilising variables caused an increase in the number of fingers and vice versa. The investigation is extended here for the effects of the medium heterogeneity. The cell problem derived via homogenization theory [20] is solved and its solution used to evaluate the elements of the macrodispersion tensor as functions of time for flow aligned parallel to gravity. The longitudinal coefficient exhibits asymptotic behaviour for favourable and moderately unfavourable density contrasts while it grows indefinitely for higher density contrasts. The range of densities stabilised by medium heterogeneities can thus be estimated from the behaviour of the coefficient. The d³f software program is used for the numerical simulations. The code uses the cell-centred finite volume and the implicit Euler techniques for the spatial and temporal discretisations respectively.  相似文献   

Probabilistic collocation and lagrangian sampling for advective tracer transport in randomly heterogeneous porous media     

Florian Müller  Patrick JennyDaniel W. Meyer 《Advances in water resources》2011,34(12):1527-1538

The Karhunen-Loeve (KL) decomposition and the polynomial chaos (PC) expansion are elegant and efficient tools for uncertainty propagation in porous media. Over recent years, KL/PC-based frameworks have successfully been applied in several contributions for the flow problem in the subsurface context. It was also shown, however, that the accurate solution of the transport problem with KL/PC techniques is more challenging. We propose a framework that utilizes KL/PC in combination with sparse Smolyak quadrature for the flow problem only. In a subsequent step, a Lagrangian sampling technique is used for transport. The flow field samples are calculated based on a PC expansion derived from the solutions at relatively few quadrature points. To increase the computational efficiency of the PC-based flow field sampling, a new reduction method is applied. For advection dominated transport scenarios, where a Lagrangian approach is applicable, the proposed PC/Monte Carlo method (PCMCM) is very efficient and avoids accuracy problems that arise when applying KL/PC techniques to both flow and transport. The applicability of PCMCM is demonstrated for transport simulations in multivariate Gaussian log-conductivity fields that are unconditional and conditional on conductivity measurements.  相似文献   

各向异性双相介质弹性波场褶积算法数值模拟     

宫猛  李信富 《地震工程学报》2011,33(4):313-318

从各向同性介质中波场数值模拟的褶积微分算子法出发,推导出了各向异性双相介质中波场传播数值计算的褶积新算法.将常见的二阶微分Biot波动方程用等效的一阶速度—应力双曲方程表示,其中未知的波场向量包括固相和流体的速度分量和应力分量,由此对方程的时间项使用交错网格差分方法计算,而对空间项则采用褶积微分算法进行求解.对各向异性双相介质在单层介质模型和双层介质模型中的波场特征进行了研究.研究的结果显示,在两层介质分界面上当地震波产生反射时能观测到两类纵波和横波,并且在衰减系数大的介质里慢纵波很难见到.  相似文献   

Luminescence dating of Holocene sediment cores from a wave-dominated and mountainous river delta in central Vietnam     

《Quaternary Geochronology》2022

Sediments of river deltas provide valuable records of past coastal environments. Optically-stimulated luminescence (OSL) dating has become an alternative to radiocarbon dating for constraining the sediment chronology in large deltas that allow for sufficient sunlight bleaching of sediments during the fluvial transport. However, its applicability to smaller deltas with mountainous riverine systems has not been confirmed yet. To check this, we examine multiple signals from two Holocene sediment cores in the wave-dominated Thu Bon River delta in central Vietnam. Two cores were collected, respectively, 3.9 km and 1.2 km inland from the present shoreline and both show a >-25-m thick succession of coarsening-upward mud to sand deposits. Coarse grains (180–250 μm in diameter) of quartz and K-feldspar, and fine grains (4–11 μm in diameter) of quartz and polymineral were extracted from the upper and lower parts of the cores for multi-grain measurements of quartz OSL, and of feldspar infrared-stimulated luminescence (IRSL) at 50 °C (IR₅₀) and post-IR IRSL at 175 °C (pIRIR₁₇₅) to determine burial ages. In addition, facies analysis and radiocarbon dating were conducted. The landward core consists of transgressive to early regressive estuarine and prodelta facies, which is overlain by a sandy beach-shoreface facies. The seaward core consists of a relatively simple shallowing-upward succession from muddy prodelta facies to sandy beach-shoreface facies. All luminescence ages except for pIRIR₁₇₅ of fine grains are mostly consistent with the radiocarbon ages. Instead, pIRIR₁₇₅ ages of fine grains are significantly overestimated with variable offsets. OSL and IR₅₀ of fine grains provide reasonable age estimates, as these grains were likely well bleached during the transport even along a short and steep mountainous river. Consistent age estimates of all signals from sand indicate that sand was well-bleached in the beach and shoreface owing to the frequent sediment reworking by waves and currents. These results support the hypothesis that luminescence dating is applicable to Holocene wave-dominated deltas and reiterate that comparing different luminescence signals is an effective way to check reliability of the age estimates in environments where the sunlight bleaching is not ensured.  相似文献   

Appraisal of fracture sampling methods and a new workflow to characterise heterogeneous fracture networks at outcrop     

《Journal of Structural Geology》2015

Characterising fractures at outcrop for use as analogues to fractured reservoirs can use several methods. Four important fracture data collection methods are linear scanline sampling, areal sampling, window sampling and circular scanline sampling. In regions of homogeneous fracture networks these methods are adequate to characterise fracture patterns for use as outcrop analogues, however where fractures are heterogeneous, it is more difficult to characterise fracture networks and a different approach is needed.We develop a workflow for fracture data collection in a region of heterogeneous fractures in a fold and thrust belt, which we believe has applicability to a wide variety of fracture networks in different tectonic settings. We use an augmented circular scanline method, along with areal sampling to collect a range of fracture attribute data, including orientation, length, aperture, spatial distribution and intensity. This augmented circular scanline method more than halves the time taken for data collection, provides accurate, unbiased data that is representative of local fracture network attributes and involves data collection of a wider range of fracture attributes than other sampling techniques alone.  相似文献   

非均质含水层中渗流与溶质运移研究进展简   总被引：1，自引：0，他引：1  

Qin Ronggao  Cao Guangzhu  Wu Yanqing 《地球科学进展》2014,29(1):30-41

Natural aquifer heterogeneity controls the groundwater flow and solute transport, and how to accurately quantify the flow and solute transport in heterogeneous aquifers has received wide attention by many scholars, and has become a hot research topic in earth science. Theoretically, a systematic review is given by the following aspect: flow and solute transport model, moment analysis, multi scale analysis. The resolved and remained issues for scale conversion in hydrogeology research are pointed out. Secondly, recent advances of heterogeneous field test, uncertainty and velocity connectivity are analyzed. Finally, the geophysical inversion of aquifer heterogeneity, stochastic theory and development of stochastic simulation software, scale conversion and uncertainty of velocity connectivity, and the relationship between heterogeneity and hydrogeological condition on the major four aspects of the future research direction is pointed out.  相似文献   

Degradation of Textile Dyeing Wastewater by a Modified Solar Photo‐Fenton Process Using Steel Scrap/H2O2     

Ganesan Ramachandran  Thanasekaran Kumarasamy 《洁净——土壤、空气、水》2013,41(3):267-274

This experimental research deals with using steel scrap as a heterogeneous catalyst. This catalyzes the oxidation reaction of real textile dye wastewater based on a modified solar photo‐Fenton oxidation process. Morphologic analysis and mapping of the elementary composition of the steel scrap have been carried out by scanning electron microscopy. The effects of concentration of H₂O₂, the pH of the solution and the catalyst loading on the degradation of textile dye wastewater are elucidated. Kinetic studies have been performed for the decolorization of wastewater under optimum conditions. It could be concluded that the steel scrap is a potential substitute for ferrous salts as a catalyst for the solar photo‐Fenton reaction.  相似文献   

Water quality and coral bleaching thresholds: Formalising the linkage for the inshore reefs of the Great Barrier Reef, Australia   总被引：1，自引：0，他引：1  

Scott A. Wooldridge 《Marine pollution bulletin》2009,58(5):745-751

The threats of wide-scale coral bleaching and reef demise associated with anthropogenic climate change are widely known. Here, the additional role of poor water quality in lowering the thermal tolerance (i.e. bleaching ‘resistance’) of symbiotic reef corals is considered. In particular, a quantitative linkage is established between terrestrially-sourced dissolved inorganic nitrogen (DIN) loading and the upper thermal bleaching thresholds of inshore reefs on the Great Barrier Reef, Australia. Significantly, this biophysical linkage provides concrete evidence for the oft-expressed belief that improved coral reef management will increase the regional-scale survival prospects of corals reefs to global climate change. Indeed, for inshore reef areas with a high runoff exposure risk, it is shown that the potential benefit of this ‘local’ management imperative is equivalent to ∼2.0-2.5 °C in relation to the upper thermal bleaching limit; though in this case, a potentially cost-prohibitive reduction in end-of-river DIN of >50-80% would be required. An integrated socio-economic modelling framework is outlined that will assist future efforts to understand (optimise) the alternate tradeoffs that the water quality/coral bleaching linkage presents.  相似文献   

Impact of log-transmissivity variogram structure on groundwater flow and transport predictions     

M. Riva  M. Willmann 《Advances in water resources》2009

We analyze the impact of the choice of the variogram model adopted to characterize the spatial variability of natural log-transmissivity on the evaluation of leading (statistical) moments of hydraulic heads and contaminant travel times and trajectories within mildly (randomly) heterogeneous two-dimensional porous systems. The study is motivated by the fact that in several practical situations the differences between various variogram types and a typical noisy sample variogram are small enough to suggest that one would often have a hard time deciding which of the tested models provides the best fit. Likewise, choosing amongst a set of seemingly likely variogram models estimated by means of geostatistical inverse models of flow equations can be difficult due to lack of sensitivity of available model discrimination criteria. We tackle the problem within the framework of numerical Monte Carlo simulations for mean uniform and radial flow scenarios. The effect of three commonly used isotropic variogram models, i.e., Gaussian, Exponential and Spherical, is analyzed. Our analysis clearly shows that (ensemble) mean values of the quantities of interest are not considerably influenced by the variogram shape for the range of parameters examined. Contrariwise, prediction variances of the quantities examined are significantly affected by the choice of the variogram model of the log-transmissivity field. The spatial distribution of the largest/lowest values of the relative differences observed amongst the tested models depends on a combination of variogram shape and parameters and relative distance from internal sources and the outer domain boundary. Our findings suggest the need of developing robust techniques to discriminate amongst a set of seemingly equally likely alternative variogram models in order to provide reliable uncertainty estimates of state variables.  相似文献   

从湍流经典理论到大气湍流非平衡态热力学理论   总被引：2，自引：0，他引：2  

胡隐樵  陈晋北  吕世华 《高原气象》2012,31(1):1-27

湍流是日常生活中一种普遍的自然现象,也是经典物理学仍未完全解决的难题。湍流更是大气运动的最基本特征。本文系统地回顾了大气湍流经典理论发展简史,进一步详细介绍了大气湍流非平衡态热力学理论。大气湍流非平衡态热力学理论在熵平衡方程中引入动力过程,进而统一推导出大气湍流输送的Fourier定律、Flick定律和Newton定律,证明了Dufour效应、Soret效应、可逆动力过程与热力不可逆湍流输送过程之间的交叉耦合效应,以及湍流强度定理。这些定律和定理中得到了观测的事实验证,同时它们的唯象系数也由观测资料所确定。湍流强度定理揭示,湍流发展的宏观原因是速度和温度的剪切效应,Reynolds湍流和Rayleigh-Bénard湍流共存于大气湍流中。热力过程和动力过程间耦合效应现象的发现突破了传统湍流输送理论,即Fourier定律、Flick定律和Newton定律的观点——一个宏观量的输送通量等价于这个宏观量的梯度湍流输送通量。热力和动力过程间的耦合原理认为,一个宏观量的输送通量包括这个量的梯度湍流输送通量和速度耦合输送通量两部分。因此,能量和物质的垂直输送通量除了相应物理量梯度造成的湍流输送外,还应包括垂直速度耦合效应,即辐散或辐合运动造成的耦合效应。在一个很宽的尺度范围内,地表面的空间特征是非均匀的。下垫面非均匀性造成的对流运动将引起大气的辐散或辐合运动。这可能是导致地表能量收支不平衡的重要原因之一。垂直速度对垂直湍流输送的交叉耦合效应为非均匀下垫面大气边界层理论的发展,并为克服地表能量收支不平衡问题及非均匀下垫面大气边界层参数化遇到的困难提供了可能的线索。  相似文献