共查询到20条相似文献，搜索用时 46 毫秒
1.
2.
Vegetation constitutes one of the fundamental types of land use on Earth. The presence of trees in urban areas can decrease local winds and exchange sensible and latent heat with the surrounding environments, thus exerting notable influences on the urban microenvironment. A better understanding of the turbulent transfer of momentum and scalars around vegetation canopy could significantly contribute to improvement of the urban environment. This work develops a large-eddy simulation(LES) method that is applicable to model the flow and scalar transport over the forest canopy. We study the atmospheric flow over complex forested areas under typical weather conditions by coupling LES to the mesoscale model. Models of radiation and energy balance have been developed with explicit treatment of the vegetation canopy. By examining the flow over a forest canopy under a range of stability conditions, we found that buoyancy enhances or suppresses turbulent mixing in unstable or stable atmosphere respectively, with decreasing or increasing wind shear, respectively. From the multiscale modeling of the Beijing Olympic Forest Park, the present coupling scheme proves to better resolve the diurnal variations in wind speed,temperature, and relative humidity over complex urban terrains. The coupling scheme is superior to the traditional mesoscale model in terms of wind field simulation. This is mainly because the coupling scheme not only takes the influences of external mesoscale flow into consideration, but also resolves the heterogeneous urban surface at a fine scale by downscaling, thus better reproducing the complex flow and turbulent transport in the urban roughness sublayer. 相似文献
3.
《Advances in water resources》2005,28(3):257-271
We present a new approach to reservoir simulation that gives accurate resolution of both large-scale and fine-scale flow patterns. The method uses a mixed multiscale finite-element method (MMsFEM) to solve the pressure equation on a coarse grid and a streamline-based technique to solve the fluid transport on a fine-scale subgrid. The MMsFEM is based on the construction of special approximation velocity spaces that are adaptive to the local properties of the differential operator. As such, MMsFEM produces a detailed subgrid velocity field that reflects the impact of the fine-scale heterogeneous structures. By combining MMsFEM with rapid streamline simulation of the fluid transport, we aim towards a numerical scheme that facilitates routine reservoir simulation of large heterogeneous geomodels without upscaling. The new method is applied to two different test cases. The first test case consists of two (strongly) heterogeneous quarter five-spot problems in 2D. The second test case is a 3D upscaling benchmark taken from the 10th SPE Comparative Solution Project, a project whose purpose is to compare and validate upscaling techniques. The test cases demonstrate that the combination of multiscale methods and streamlines is a robust and viable alternative to traditional upscaling-based reservoir simulation. 相似文献
4.
Characterization of shale reservoirs, which are typically of low permeability, is very difficult because of the presence of multiscale structures. While three-dimensional (3D) imaging can be an ultimate solution for revealing important complexities of such reservoirs, acquiring such images is costly and time consuming. On the other hand, high-quality 2D images, which are widely available, also reveal useful information about shales’ pore connectivity and size. Most of the current modeling methods that are based on 2D images use limited and insufficient extracted information. One remedy to the shortcoming is direct use of qualitative images, a concept that we introduce in this paper. We demonstrate that higher-order statistics (as opposed to the traditional two-point statistics, such as variograms) are necessary for developing an accurate model of shales, and describe an efficient method for using 2D images that is capable of utilizing qualitative and physical information within an image and generating stochastic realizations of shales. We then further refine the model by describing and utilizing several techniques, including an iterative framework, for removing some possible artifacts and better pattern reproduction. Next, we introduce a new histogram-matching algorithm that accounts for concealed nanostructures in shale samples. We also present two new multiresolution and multiscale approaches for dealing with distinct pore structures that are common in shale reservoirs. In the multiresolution method, the original high-quality image is upscaled in a pyramid-like manner in order to achieve more accurate global and long-range structures. The multiscale approach integrates two images, each containing diverse pore networks – the nano- and microscale pores – using a high-resolution image representing small-scale pores and, at the same time, reconstructing large pores using a low-quality image. Eventually, the results are integrated to generate a 3D model. The methods are tested on two shale samples for which full 3D samples are available. The quantitative accuracy of the models is demonstrated by computing their morphological and flow properties and comparing them with those of the actual 3D images. The success of the method hinges upon the use of very different low- and high-resolution images. 相似文献
5.
《Advances in water resources》2005,28(3):303-314
In this paper we combine a multiscale data integration technique introduced in [Lee SH, Malallah A, Datta-Gupta A, Hidgon D. Multiscale data integration using Markov Random Fields. SPE Reservoir Evaluat Eng 2002;5(1):68–78] with upscaling techniques for spatial modeling of permeability. The main goal of this paper is to find fine-scale permeability fields based on coarse-scale permeability measurements. The approach introduced in the paper is hierarchical and the conditional information from different length scales is incorporated into the posterior distribution using a Bayesian framework. Because of a complicated structure of the posterior distribution Markov chain Monte Carlo (MCMC) based approaches are used to draw samples of the fine-scale permeability field. 相似文献
6.
Genetic algorithms have been shown to be powerful tools for solving a wide variety of water resources optimization problems. Applying these approaches to complex, large-scale water resources applications can be difficult due to computational limitations, especially when a numerical model is needed to evaluate different solutions. This problem is particularly acute for solving field-scale groundwater remediation design problems, where fine spatial grids are often needed for accuracy. Finer grids usually improve the accuracy of the solutions, but they are also computationally expensive. In this paper we present multiscale island injection genetic algorithms (IIGAs), in which the optimization algorithms have different multiscale populations working on different islands (groups of processors) and periodically exchanging information. This new approach is tested using a field-scale pump-and-treat design problem at the Umatilla Army Depot in Oregon, USA. The performance of several variations of this approach is compared with the results of a simple genetic algorithm. The new approach found the same solution as much as 81% faster than the simple genetic algorithm and 9–53% faster than other previously formulated multiscale strategies. These findings indicate substantial promise for multiscale IIGA approaches to improve solution of complex water resources applications at the field scale. 相似文献
7.
本文针对时间推移地震本身包含不同时期的两次或者两次以上的勘探反问题,构造了一种快速有效的反演方法--多尺度全变分法.通过引入全变分正则化来代替传统的Tikhonov正则化,针对待反演参数不连续的情况,提高了算法精度.为了提高计算效率,引入了多尺度方法,从而构造了多尺度全变分方法.在数值模拟中,针对一个时间推移地震反演问题对多尺度-Tikhonov正则化法、单一尺度全变分法、以及本文所构造的多尺度全变分法进行了比较.结果表明,本文所提出的多尺度全变分法是一种稳定、快速和精确的反演方法. 相似文献
8.
We invert three-dimensional seismic data by a multiscale phase inversion scheme, a modified version of full waveform inversion, which applies higher order integrations to the input signal to produce low-boost signals. These low-boost signals are used as the input data for the early iterations, and lower order integrations are computed at the later iterations. The advantages of multiscale phase inversion are that it (1) is less dependent on the initial model compared to full waveform inversion, (2) is less sensitive to incorrectly modelled magnitudes and (3) employs a simple and natural frequency shaping filtering. For a layered model with a three-dimensional velocity anomaly, results with synthetic data show that multiscale phase inversion can sometimes provide a noticeably more accurate velocity profile than full waveform inversion. Results with the Society of Exploration Geophysicists/European Association of Geoscientists and Engineers overthrust model shows that multiscale phase inversion more clearly resolves meandering channels in the depth slices. However, the data and model misfit functions achieve about the same values after 50 iterations. The results with three-dimensional ocean-bottom cable data show that, compared to the full waveform inversion tomogram, the three-dimensional multiscale phase inversion tomogram provides a better match to the well log, and better flattens angle-domain common image gathers. The problem is that the tomograms at the well log provide an incomplete low-wavenumber estimate of the log's velocity profile. Therefore, a good low-wavenumber estimate of the velocity model is still needed for an accurate multiscale phase inversion tomogram. 相似文献
9.
Kuo-Chin Hsu Kuan-Chih Chen 《Stochastic Environmental Research and Risk Assessment (SERRA)》2010,24(7):1053-1065
This paper proposes a multiscale flow and transport model which can be used in three-dimensional fractal random fields. The
fractal random field effectively describes a field with a high degree of variability to satisfy the one-point statistics of
Levy-stable distribution and the two-point statistics of fractional Levy motion (fLm). To overcome the difficulty of using
infinite variance of Levy-stable distribution and to provide the physical meaning of a finite domain in real space, truncated
power variograms are utilized for the fLm fields. The fLm model is general in the sense that both stationary and commonly
used fractional Brownian motion (fBm) models are its special cases. When the upper cutoff of the truncated power variogram
is close to the lower cutoff, the stationary model is well approximated. The commonly used fBm model is recovered when the
Levy index of fLm is 2. Flow and solute transport were analyzed using the first-order perturbation method. Mean velocity,
velocity covariance, and effective hydraulic conductivity in a three-dimensional fractal random field were derived. Analytical
results for particle displacement covariance and macrodispersion coefficients are also presented. The results show that the
plume in an fLm field moves slower at early time and has more significant long-tailing behavior at late time than in fBm or
stationary exponential fields. The proposed fractal transport model has broader applications than those of stationary and
fBm models. Flow and solute transport can be simulated for various scenarios by adjusting the Levy index and cutoffs of fLm
to yield more accurate modeling results. 相似文献
10.
We present the first attempt to explain slow earthquakes as cascading thermal-mechanical instabilities. To attain this goal we investigate brittle-ductile coupled thermal-mechanical simulation on vastly different time scales. The largest scale model consists of a cross section of a randomly perturbed elasto-visco-plastic continental lithosphere on the order of 100 × 100 km scale with no other initial structures. The smallest scale model investigates a km-scale subsection of the large model and has a local resolution of 40 × 40 m. The model is subject to a constant extension velocity applied on either side. We assume a free top surface and with a zero tangential stress along the other boundaries. Extension is driven by velocity boundary conditions of 1 cm/a applied on either side of the model. This is the simplest boundary condition, and makes it an ideal starting point for understanding the behavior of a natural system with multiscale brittle-ductile coupling. Localization feedback is observed as faulting in the brittle upper crust and ductile shearing in an elasto-viscoplastic lower crust. In this process brittle faulting may rupture at seismogenic rates, e.g., at 102–103 ms?1, whereas viscous shear zones propagate at much slower rates, up to 3 × 10?9 ms?1. This sharp contrast in the strain rates leads to complex short-time-scale interactions at the brittle-ductile transition. We exploit the multiscale capabilities from our new simulations for understanding the underlying thermo-mechanics, spanning vastly different, time- and length-scales. 相似文献
11.
Zachary K. Curtis Hua-sheng Liao Shu-Guang Li Prasanna Venkatesh Sampath David P. Lusch 《Ground water》2019,57(5):784-806
Managing nonpoint-source (NPS) pollution of groundwater systems is a significant challenge because of the heterogeneous nature of the subsurface, high costs of data collection, and the multitude of scales involved. In this study, we assessed a particularly complex NPS groundwater pollution problem in Michigan, namely, the salinization of shallow aquifer systems due to natural upwelling of deep brines. We applied a system-based approach to characterize, across multiple scales, the integrated groundwater quantity–quality dynamics associated with the brine upwelling process, assimilating a variety of modeling tools and data—including statewide water well datasets scarcely used for larger scientific analysis. Specifically, we combined (1) data-driven modeling of massive amounts of groundwater/geologic information across multiple spatial scales with (2) detailed analysis of groundwater salinity dynamics and process-based flow modeling at local scales. Statewide “hotspots” were delineated and county-level severity rankings were developed based on dissolved chloride (Cl−) concentration percentiles. Within local hotspots, the relative impact of upwelling was determined to be controlled by: (1) streams—which act as “natural pumps” that bring deeper (more mineralized) groundwater to the surface; (2) the occurrence of nearly impervious geologic material at the surface—which restricts fresh water dilution of deeper, saline groundwater; and (3) the space–time evolution of water well withdrawals—which induces slow migration of saline groundwater from its natural course. This multiscale, data-intensive approach significantly improved our understanding of the brine upwelling processes in Michigan, and has applicability elsewhere given the growing availability of statewide water well databases. 相似文献
12.
The identification of groundwater parameters in heterogeneous systems is a major challenge in groundwater modeling. Flexible parameterization methods are needed to assess the complexity of the spatial distributions of these parameters in real aquifers. In this article, we introduce an adaptative parameterization to identify the distribution of hydraulic conductivity within the large‐scale (4400 km2) Upper Rhine aquifer. The method is based on adaptative multiscale triangulation (AMT) coupled with an inverse problem procedure that identifies the parameters' distributions by reducing the error between measured and simulated heads. The AMT method has the advantage of combining both zonation and interpolation approaches. The AMT method uses area‐based interpolation rather than an interpolation based on stochastic features. The method is applied to a standard 2D groundwater model that takes into account the interactions between the aquifer and surface water bodies, groundwater recharge, and pumping wells. The simulation period covers 204 months, from January 1986 to December 2002. Recordings at 109 piezometers are used for model calibration. The simulated heads are globally quite accurate and reproduce the main dynamics of the system. The local hydraulic conductivities resulting from the AMT method agree qualitatively with existing local experimental observations across the Rhine aquifer. 相似文献
13.
Reactive contaminant transport models are used by hydrologists to simulate and study the migration and fate of industrial waste in subsurface aquifers. Accurate transport modeling of such waste requires clear understanding of the system’s parameters, such as sorption and biodegradation. In this study, we present an efficient sequential data assimilation scheme that computes accurate estimates of aquifer contamination and spatially variable sorption coefficients. This assimilation scheme is based on a hybrid formulation of the ensemble Kalman filter (EnKF) and optimal interpolation (OI) in which solute concentration measurements are assimilated via a recursive dual estimation of sorption coefficients and contaminant state variables. This hybrid EnKF-OI scheme is used to mitigate background covariance limitations due to ensemble under-sampling and neglected model errors. Numerical experiments are conducted with a two-dimensional synthetic aquifer in which cobalt-60, a radioactive contaminant, is leached in a saturated heterogeneous clayey sandstone zone. Assimilation experiments are investigated under different settings and sources of model and observational errors. Simulation results demonstrate that the proposed hybrid EnKF-OI scheme successfully recovers both the contaminant and the sorption rate and reduces their uncertainties. Sensitivity analyses also suggest that the adaptive hybrid scheme remains effective with small ensembles, allowing to reduce the ensemble size by up to 80% with respect to the standard EnKF scheme. 相似文献
14.
High-resolution geologic models that incorporate observed state data are expected to effectively enhance the reliability of reservoir performance prediction. One of the major challenges faced is how to solve the large-scale inverse modeling problem, i.e., to infer high-resolution models from the given observations of state variables that are related to the model parameters according to some known physical rules, e.g., the flow and transport partial differential equations. There are typically two difficulties, one is the high-dimensional problem and the other is the inverse problem. A multiscale inverse method is presented in this work to attack these problems with the aid of a gradient-based optimization algorithm. In this method, the model responses (i.e., the simulated state data) can be efficiently computed from the high-resolution model using the multiscale finite-volume method. The mismatch between the observations and the multiscale solutions is then used to define a proper objective function, and the fine-scale sensitivity coefficients (i.e., the derivatives of the objective function with respect to each node’s attribute) are computed by a multiscale adjoint method for subsequent optimization. The difficult high-dimensional optimization problem is reduced to a one-dimensional one using the gradient-based gradual deformation method. A synthetic single-phase transient flow example problem is employed to illustrate the proposed method. Results demonstrate that the multiscale framework presented is not only computationally efficient but also can generate geologically consistent models. By preserving spatial structure for inverse modeling, the method presented overcomes the artifacts introduced by the multiscale simulation and may enhance the prediction ability of the inverse-conditional realizations generated. 相似文献
15.
1 INTRODUCTION In recent years, due to the increase in population and industrial developments, mankind has faced manyproblems associated with rivers, coastal waters and reservoirs. Some of these problems are flood control,water supply, power generation, and irrigation. In addition, making new hydraulic structures changesnatural conditions. Prediction of these changes is necessary for designing such constructions. For solutionof these problems usually an assessment of flow pattern, sedim… 相似文献
16.
The main problems in three-dimensional gravity inversion are the non-uniqueness of the solutions and the high computational cost of large data sets. To minimize the high computational cost, we propose a new sorting method to reduce fluctuations and the high frequency of the sensitivity matrix prior to applying the wavelet transform. Consequently, the sparsity and compression ratio of the sensitivity matrix are improved as well as the accuracy of the forward modeling. Furthermore, memory storage requirements are reduced and the forward modeling is accelerated compared with uncompressed forward modeling. The forward modeling results suggest that the compression ratio of the sensitivity matrix can be more than 300. Furthermore, multiscale inversion based on the wavelet transform is applied to gravity inversion. By decomposing the gravity inversion into subproblems of different scales, the non-uniqueness and stability of the gravity inversion are improved as multiscale data are considered. Finally, we applied conventional focusing inversion and multiscale inversion on simulated and measured data to demonstrate the effectiveness of the proposed gravity inversion method. 相似文献
17.
Realistic models of lithologic structure are critical for predicting flow and transport through heterogeneous volcanic aquifers. Existing models of lava flows based on physical processes are able to realistically simulate flow geometry and lithology, but the computational intensity limits applicability in generating entire aquifers. Fast surface-based models have been developed for hazard mapping, but these do not incorporate 3D geometry or lithology critical for hydrogeologic applications. Here we develop a hybrid modeling method (HMM) based on a combination of a process-based model (PBM) and a surface-based model. The methodologies are presented and compared to a known single flow and to each other in a full aquifer simulation. Results indicate that both the PBM and HMM simulations reasonably reproduce the flow geometry (length, branching, thickness) of the 1984 eruption of Mauna Loa in Hawai’i. Simulations of a volcanic aquifer built from 100 flows with the PBM and HMM are similar in spatial distribution and overall proportions of lithology (aa, transitional, pahoehoe, ash), flow geometry, and aquifer geometry. Thus, the hybrid method is an efficient method to generate geologically realistic models of volcanic aquifer structure. Model realism and parameterization can be improved as more field data become available. 相似文献
18.
Billie F.Spencer Jr. Chia-Ming Chang Thomas M.Frankie Daniel A.Kuchma Pedro F.Silva Adel E.Abdelnaby 《地震工程与工程振动(英文版)》2014,(Z1)
Hybrid simulation has been shown to be a cost-effective approach for assessing the seismic performance of structures. In hybrid simulation,critical parts of a structure are physically tested,while the remaining portions of the system are concurrently simulated computationally,typically using a finite element model. This combination is realized through a numerical time-integration scheme,which allows for investigation of full system-level responses of a structure in a cost-effective manner. However,conducting hybrid simulation of complex structures within large-scale testing facilities presents significant challenges. For example,the chosen modeling scheme may create numerical inaccuracies or even result in unstable simulations; the displacement and force capacity of the experimental system can be exceeded; and a hybrid test may be terminated due to poor communication between modules(e.g.,loading controllers,data acquisition systems,simulation coordinator). These problems can cause the simulation to stop suddenly,and in some cases can even result in damage to the experimental specimens; the end result can be failure of the entire experiment. This study proposes a phased approach to hybrid simulation that can validate all of the hybrid simulation components and ensure the integrity largescale hybrid simulation. In this approach,a series of hybrid simulations employing numerical components and small-scale experimental components are examined to establish this preparedness for the large-scale experiment. This validation program is incorporated into an existing,mature hybrid simulation framework,which is currently utilized in the Multi-Axial Full-Scale Sub-Structuring Testing and Simulation(MUST-SIM) facility of the George E. Brown Network for Earthquake Engineering Simulation(NEES) equipment site at the University of Illinois at Urbana-Champaign. A hybrid simulation of a four-span curved bridge is presented as an example,in which three piers are experimentally controlled in a total of 18 degrees of freedom(DOFs). This simulation illustrates the effectiveness of the phased approach presented in this paper. 相似文献
19.
Billie F. SpencerJr Chia-Ming Chang Thomas M. Frankie Daniel A. Kuchma Pedro F. Silva Adel E. Abdelnaby 《地震工程与工程振动(英文版)》2014,13(1):63-77
Hybrid simulation has been shown to be a cost-effective approach for assessing the seismic performance of structures. In hybrid simulation, critical parts of a structure are physically tested, while the remaining portions of the system are concurrently simulated computationally, typically using a finite element model. This combination is realized through a numerical time-integration scheme, which allows for investigation of full system-level responses of a structure in a cost-effective manner. However, conducting hybrid simulation of complex structures within large-scale testing facilities presents significant challenges. For example, the chosen modeling scheme may create numerical inaccuracies or even result in unstable simulations; the displacement and force capacity of the experimental system can be exceeded; and a hybrid test may be terminated due to poor communication between modules (e.g., loading controllers, data acquisition systems, simulation coordinator). These problems can cause the simulation to stop suddenly, and in some cases can even result in damage to the experimental specimens; the end result can be failure of the entire experiment. This study proposes a phased approach to hybrid simulation that can validate all of the hybrid simulation components and ensure the integrity large-scale hybrid simulation. In this approach, a series of hybrid simulations employing numerical components and small-scale experimental components are examined to establish this preparedness for the large-scale experiment. This validation program is incorporated into an existing, mature hybrid simulation framework, which is currently utilized in the Multi-Axial Full-Scale Sub-Structuring Testing and Simulation (MUST-SIM) facility of the George E. Brown Network for Earthquake Engineering Simulation (NEES) equipment site at the University of Illinois at Urbana-Champaign. A hybrid simulation of a four-span curved bridge is presented as an example, in which three piers are experimentally controlled in a total of 18 degrees of freedom (DOFs). This simulation illustrates the effectiveness of the phased approach presented in this paper. 相似文献
20.
为了研究地球物理观测信号在多尺度空间中的复杂性,本文提出了一种新的广义熵谱（General Entropy Spectral）概念,其在零频率下的值正好给出了传统意义的信息熵。同时,还提出了一种新的局部多尺度熵（Local Multiscale Entropy）分析方法,该方法结合了信息熵和多尺度方法,可用于对数据特征、信息复杂度进行多尺度表达。上述这些新方法适用于分析各种复杂的地球物理信号,以便深入挖掘更多的信号特征和信息。 相似文献