首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 156 毫秒
1.
Numerical Modelling in Geo-Electromagnetics: Advances and Challenges   总被引:11,自引:0,他引:11  
During the last decade, tremendous advances have been observed in the broad field of numerical modelling for geo-electromagnetic applications. This trend received support due to increasing industrial needs, mainly caused by hydrocarbon and ore exploration industry. On the other hand, the increasing reliability and accuracy of data acquisition techniques further spurs this development. In this review, we will focus on advances and challenges in numerical modelling in geo-electromagnetics. We review recent developments in the discrete solution of the 3-D induction problem in the time and frequency domains. Particularly, advantages and disadvantages of the common numerical techniques for solving partial differential equations such as the Finite Difference and Finite Element methods will be considered.  相似文献   

2.
Multiple intertidal bars are common features of wave‐dominated sandy beaches, yet their short‐term (<1 month) and small‐scale (<1 km) morphology and dynamics remain poorly understood. This study describes the morphodynamics of multiple intertidal bars in North Lincolnshire, England, during single and lunar tidal cycles under two contrasting conditions – first when significant wave height was <0·5 m and second when significant wave height frequently exceeded 1 m. The relative importance of swash, surf and shoaling processes in determining morphological change was examined using detailed field observations and a numerical model. The beach featured four intertidal bars and both cross‐shore and longshore bar morphology evolved during the field investigation, particularly under medium to high wave‐energy conditions. Numerical modelling suggests shoaling processes are most common on the seaward two bars under calm wave conditions (Hs < 0·5 m) and that surf zone processes become more common during neap tides and under more energetic (Hs < 0·5 m) conditions. Surf processes dominate the inner two bars, though swash influence increases in a landward direction. The numerical modelling results combined with low tide survey data and high‐resolution morphological measurements strongly suggest changes in the intertidal bar morphology are accomplished by surf zone processes rather than by shoaling wave or swash processes. This is because shoaling waves do not induce significant sediment transport to have any morphological effect, whereas swash action generally does not have enough scope to act as the swash zone is much narrower than the surf zone. It was found, however, that the absolute rate of morphological change under swash action and surfzone processes are of similar magnitudes and that swash action may induce a significant amount of local morphological change when the high tide mark is located on the upper bar, making this process important for bar morphodynamics. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

3.
Salinization of land and water is a significant challenge in most continents and particularly in arid and semi-arid regions. The need to accurately forecast surface and groundwater interactions has promoted the use of physically-based numerical modelling approaches in many studies. In this regard, two issues can be considered as the main research challenges. First, in contrast with surface water, there is generally less observed level and salinity data available for groundwater systems. These data are critical in the validation and verification of numerical models. The second challenge is to develop an integrated surface-groundwater numerical model that is capable of salt mobilization modelling but which can be validated and verified against accurate observed data. This paper reviews the current state of understanding of groundwater and surface water interactions with particular respect to the numerical modelling of salt mobilization. 3D physically-based fully coupled surface-subsurface numerical model with the capability of modelling density-dependent, saturated-unsaturated solute transport is an ideal tool for groundwater-surface water interaction studies. It is concluded that there is a clear need to develop modelling capabilities for the movement of salt to, from, and within wetlands to provide temporal predictions of wetland salinity which can be used to assess ecosystem outcomes.  相似文献   

4.
With the development of computing technology, numerical models are often employed to simulate flow and water quality processes in coastal environments. However, the emphasis has conventionally been placed on algorithmic procedures to solve specific problems. These numerical models, being insufficiently user-friendly, lack knowledge transfers in model interpretation. This results in significant constraints on model uses and large gaps between model developers and practitioners. It is a difficult task for novice application users to select an appropriate numerical model. It is desirable to incorporate the existing heuristic knowledge about model manipulation and to furnish intelligent manipulation of calibration parameters. The advancement in artificial intelligence (AI) during the past decade rendered it possible to integrate the technologies into numerical modelling systems in order to bridge the gaps. The objective of this paper is to review the current state-of-the-art of the integration of AI into water quality modelling. Algorithms and methods studied include knowledge-based system, genetic algorithm, artificial neural network, and fuzzy inference system. These techniques can contribute to the integrated model in different aspects and may not be mutually exclusive to one another. Some future directions for further development and their potentials are explored and presented.  相似文献   

5.
A combination of photogeologic mapping, analysis of Viking Orbiter thermal inertia data, and numerical modelling of eruption conditions has permitted us to construct a new model for the evolution of the martian volcano Alba Patera. Numerous digitate channel networks on the flanks of the volcano are interpreted to be carved by sapping due to the release of non-juvenile water from unconsolidated flank deposits. Using the thermal inertia measurements, we estimate the particle size of these deposits to be 3–10 µm, which, together with theoretical modelling of the disperison of explosively derived volcanic materials, leads us to conclude that the flank deposits on Alba Patera are low-relief pyroclastic flows. The recognition of numerous late-stage summit and sub-terminal lava flows thus makes Alba Patera a unique martian volcano that is transitional between the older pyroclastic-dominated highland paterae and the more recent effusive central-vent volcanoes such as the Tharsis Montes.  相似文献   

6.
Summary A three-dimensional numerical model of the divided-bar thermal conductivity measuring device that can include both solid and fluid sample portions has been developed. The model has been employed to investigate the effects of porosity in divided-bar measurements and the effects of the distribution of solid sample and contained fluids in porous material on the thermal conductivity calculated from the bar temperature measurements. The results show that the positions of the temperature sensors and thus the nature of the vertical column between the temperature measuring points is of prime importance, and that the sizes of the solid and fluid portions of the sample affect the sample conductivity estimates. The divided bar measures the thermal conductivity parallel to the axis of the bar. Thermal conductivity variations in the radial direction have little effect on the calculated conductivity. Whether or not the sides of the bar are insulated strongly affects the calculated conductivity.Presented at the International Meeting on Terrestrial Heat Flow and the Structure of Lithosphere, Bechyn Castle, Czech Republic, September 2 – 7, 1991.  相似文献   

7.
We compute the transfer of oceanic lithosphere material from the surface of the model to the inner convective mantle at successive stages of the supercontinental cycle, in the time interval from the beginning of convergence of the continents to their complete dispersal. The sequence of stages of a supercontinental cycle (Wilson cycle) is calculated with a two-dimensional numerical model of assembling and dispersing continents driven by mantle flows; in turn, the flows themselves are forming under thermal and mechanical influence of continents. We obtain that during the time of the order of 300 Myr the complete stirring of oceanic lithosphere through whole mantle does not occur. This agrees with current ideas on the circulation of oceanic crust material. Former oceanic crust material appears again at the Earth’s surface in the areas of mantle upstreams. The numerical simulation demonstrates that the supercontinental cycle is a factor which intensifies stirring of the material, especially in the region beneath the supercontinent. The reasons are a recurring formation of plumes in that region as well as a global restructuring of mantle flow pattern due to the process of joining and separation of continents. The computations of viscous shear stresses are also carried out in the mantle as a function of spatial coordinates and time. With a simplified model of uniform mantle viscosity, the numerical experiment shows that the typical maximal shear stresses in the major portion of the mantle measure about 5 MPa (50 bar). The typical maximal shear stresses located in the uppermost part of mantle downgoing streams (in a zone that measures roughly 200 × 200 km) are approximately 8 times greater and equal to 40 MPa (400 bar).  相似文献   

8.
In this paper we consider the concept of modelling dynamical systems using numerical–experimental substructuring. This type of modelling is applicable to large or complex systems, where some part of the system is difficult to model numerically. The substructured model is formed via the adaptive minimal control synthesis (MCS) algorithm. The aim of this paper is to demonstrate that substructuring can be carried out in real time, using the MCS algorithm. Thus, we reformulate the MCS algorithm into a substructuring form. We introduce the concepts of a transfer system, and carry out numerical simulations of the substructuring process using a coupled three mass example. These simulations are compared with direct simulations of a three mass system. In addition we consider the stability of the substructuring algorithm, which we discuss in detail for a class of second‐order transfer systems. A numerical–experimental system is considered, using a small‐scale experimental system, for which the substructuring algorithm is implemented in real time. Finally we discuss these results, with particular reference to the future application of this method to modelling large‐scale structures subject to earthquake excitation. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

9.
Bars are key morphological units in river systems, fashioning the sediment regime and bedload transport processes within a reach. Reworking of these features underpins channel adjustment at larger scales, thereby acting as a key determinant of channel stability. Despite their importance to channel evolution, few investigations have acquired spatially continuous data on bar morphology and sediment-size to investigate bar reworking. To this end, four bars along a 10 km reach of a wandering gravel-bed river were surveyed with terrestrial laser scanning (TLS), comparing downstream changes in slope, bed material size and channel planform. Detrended standard deviations (σz) were extracted from TLS point clouds and correlated to underlying physically measured median grain-size (D50), across a greater range of σz values than have hitherto been reported. The resulting linear regression model was used to create a 1 m resolution median grain-size map. A fusion of airborne LiDAR and optical-empirical bathymetric mapping was used to develop reach-scale digital elevation models (DEMs) for rapid two-dimensional hydraulic modelling using JFlow® software. The ratio of dimensionless shear stress over critical shear stress was calculated for each raster cell to calculate the effectiveness of a range of flood events (2.33–100 year recurrence intervals) to entrain sediment and rework bar units. Results show that multiple bar forming discharges exist, whereby frequent flood flows rework tail and back channel areas, while much larger, less frequent floods are required to mobilise the coarser sediment fraction on bar heads. Valley confinement is shown to exert a primary influence on patterns of bar reworking. Historical aerial photography, hyperscale DEMs and hydraulic modelling are used to explain channel adjustment at the reach scale. The proportion of the bar comprised of more frequently entrained units (tail, back channel, supra-platform) relative to more static units (bar head) exerts a direct influence upon geomorphic sensitivity. © 2018 John Wiley & Sons, Ltd.  相似文献   

10.
In the Hanford Reach of the Columbia River, a thin layer of recent alluvium overlies the sedimentary formations that comprise the unconfined groundwater aquifer. Experimental and modelling studies have demonstrated that this alluvial layer exerts significant control on the exchange of groundwater and surface water (hydrologic exchange flux), and is associated with elevated levels of biogeochemical activity. This layer is also observed to be strongly heterogeneous, and quantifying the spatial distribution of properties over the range of scales of interest is challenging. Facies are elements of a sediment classification scheme that groups complex geologic materials into a set of discrete classes according to distinguishing features. Facies classifications have been used as a framework for assigning heterogeneous material properties to grid cells of numerical models of flow and reactive transport in subsurface media. The usefulness of such an approach hinges on being able to relate facies to quantitative properties needed for flow and reactive transport modelling, and on being able to map facies over the domain of interest using readily available information. Although aquifer facies have been used in various modelling contexts, application of this concept to riverbed sediments is relatively new. Here, we describe an approach for categorizing and mapping recent alluvial (riverbed) sediments based on the integration of diverse observations with numerical simulations of river hydrodynamics. The facies have distinct distributions of sediment texture that correspond to variations in hydraulic properties, and therefore provide a useful framework for assigning heterogeneous properties in numerical simulations of hydrologic exchange flows and biogeochemical processes.  相似文献   

11.
A non-standard new code to solve multiphase viscous thermo–mechanical problems applied to geophysics is presented. Two numerical methodologies employed in the code are described: A level set technique to track the position of the materials and an enrichment of the solution to allow the strain rate to be discontinuous across the interface. These techniques have low computational cost and can be used in standard desktop PCs. Examples of phase tracking with level set are presented in two and three dimensions to study slab detachment in subduction processes and Rayleigh–Taylor instabilities, respectively. The modelling of slab detachment processes includes realistic rheology with viscosity depending on temperature, pressure and strain rate; shear and adiabatic heating mechanisms; density including mineral phase changes and varying thermal conductivity. Detachment models show a first prolonged period of thermal diffusion until a fast necking of the subducting slab results in the break–off. The influence of several numerical and physical parameters on the detachment process is analyzed: The shear heating exerts a major influence accelerating the detachment process, reducing the onset time to one half and lubricating the sinking of the detached slab. The adiabatic heating term acts as a thermal stabilizer. If the mantle temperature follows an adiabatic gradient, neglecting this heating term must be included, otherwise all temperature contrasts are overestimated. As expected, the phase change at 410 km depth (olivine–spinel transition) facilitates the detachment process due to the increase in negative buoyancy. Finally, simple plume simulations are used to show how the presented numerical methodologies can be extended to three dimensions.  相似文献   

12.
Since the late 1970s, numerical modelling has become established as an important technique for the understanding of ice sheet and glacier dynamics, and several models have been developed over the years. Ice sheet models are particularly relevant for predicting the possible response of ice sheets to climate change. Recent observations suggest that ice dynamics could play a crucial role for the contribution of ice sheets to future sea level rise under global warming conditions, and the need for further research into the matter was explicitly stated in the Fourth Assessment Report (AR4) of the United Nations Intergovernmental Panel on Climate Change (IPCC). In this paper, we review the state of the art and current problems of ice sheet and glacier modelling. An outline of the underlying theory is given, and crucial processes (basal sliding, calving, interaction with the solid Earth) are discussed. We summarise recent progress in the development of ice sheet and glacier system models and their coupling to climate models, and point out directions for future work.  相似文献   

13.
The whole subject of three-dimensional (3-D) electromagnetic (EM) modelling and inversion has experienced a tremendous progress in the last decade. Accordingly there is an increased need for reviewing the recent, and not so recent, achievements in the field. In the first part of this review paper I consider the finite-difference, finite-element and integral equation approaches that are presently applied for the rigorous numerical solution of fully 3-D EM forward problems. I mention the merits and drawbacks of these approaches, and focus on the most essential aspects of numerical implementations, such as preconditioning and solving the resulting systems of linear equations. I refer to some of the most advanced, state-of-the-art, solvers that are today available for such important geophysical applications as induction logging, airborne and controlled-source EM, magnetotellurics, and global induction studies. Then, in the second part of the paper, I review some of the methods that are commonly used to solve 3-D EM inverse problems and analyse current implementations of the methods available. In particular, I also address the important aspects of nonlinear Newton-type optimisation techniques and computation of gradients and sensitivities associated with these problems.  相似文献   

14.
— Today, wavelets are recognized to have a wide range of useful properties that allow them to treat effectively multifacet problems, such as data compression, scale-localization analysis, feature extraction, statistics, numerical simulation, visualization, and communication. Second-generation wavelets represent a recent improvement of the wavelet algorithm, allowing for greater flexibility in the spatial domain and other computational advantages. We will show how these wavelets can be employed to extract large-scale coherent structures from (1) three-dimensional turbulent flows and (2) high Rayleigh number thermal convection. We will discuss the concept of modelling via decomposition into coherent and incoherent fields, taking into account the effect of the incoherent field via statistical modelling. We will discuss wavelet properties and how they can be utilized and integrated in handling large-scale problems in earthquake physics and other nonlinear phenomena in the geosciences.  相似文献   

15.
富营养深水水库底栖动物群落与浮游生物相关性分析   总被引:1,自引:0,他引:1  
底栖动物和浮游生物通过食物网发生复杂的交互作用,是水库生态系统物质循环的重要组分.南湾水库和汤浦水库均为存在富营养化现象的深水水库,具有较长时间的水温分层期.本文选取这两座供水水库为对象,比较其水体营养状态、底栖动物的群落结构、现存量和多样性.通过群落状态指数(IICS)定量评估底栖动物群落状态,并采用PLS模型、Spearman秩相关以及回归图探讨浮游生物各类群密度与IICS的关系.研究结果显示,南湾水库富营养程度高于汤浦水库,寡毛类在南湾水库占据优势地位,摇蚊类在汤浦水库数量上占优.南湾水库以收集者为主,其数量占比远高于汤浦水库.深水水库通常存在较强的热分层,富营养化会加重水库底部的缺氧或厌氧状态,导致底栖动物遭受周期性的缺氧扰动,影响种类多样性和现存量,改变浮游生物类群与底栖动物群落的相关性,因此底栖动物群落与浮游生物的相关性能够提供判断水库富营养化程度的重要信号,进而为水库管理方提供更多有价值的参考信息.  相似文献   

16.
The subject of environmental engineering is currently of great interest. Field experiments as well as numerical models have proven their worth in this research field. An introduction to hydrodynamic modelling, coupled to the modelling of vegetation biomass is described. The developed Strive (STream RIVer Ecosystem) model is set up in the Femme (‘Flexible Environment for Mathematically Modelling the Environment’) environment and has already proven its worth in a large number of calculations (De Doncker et al., 2006 , 2008b ). Discharges and water levels are modelled together with modelling of electrical conductivity (EC). Extensive measurement campaigns are carried out to collect a large number of observations and calibration of the model is based on this data set. Furthermore, calibration methods and the discussion of this process are displayed. As a result, it is seen that the developed Strive model can model both, hydrodynamic and ecological processes, in an accurate way. The work highlights the importance of detailed determination of Manning's coefficient, dependent on discharge and amount of biomass, as an important calibration parameter for accurate modelling. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

17.
The thermal regime of the middle atmosphere is determined to a great extent by the balance between the incoming solar and outgoing infrared radiation. To account for these processes in numerical models of the middle atmosphere, parameterizations that are capable of quickly and accurately calculating infrared cooling and solar heating rates are required. These parameterizations should include the breakdown of local thermodynamic equilibrium (LTE) conditions and allow for feedbacks by ensuring that dependencies on all input parameters are accounted for. This paper discusses the major mechanisms responsible for maintaining the radiative energy budget of the middle atmosphere and presents a brief review of approaches and numerical schemes currently available for use in general circulation models. The main focus of the paper is on the approaches and schemes designed for non-LTE treatment.  相似文献   

18.
We derive a governing second-order acoustic wave equation in the time domain with a perfectly matched layer absorbing boundary condition for general inhomogeneous media. Besides, a new scheme to solve the perfectly matched layer equation for absorbing reflections from the model boundaries based on the rapid expansion method is proposed. The suggested scheme can be easily applied to a wide class of wave equations and numerical methods for seismic modelling. The absorbing boundary condition method is formulated based on the split perfectly matched layer method and we employ the rapid expansion method to solve the derived new perfectly matched layer equation. The use of the rapid expansion method allows us to extrapolate wavefields with a time step larger than the ones commonly used by traditional finite-difference schemes in a stable way and free of dispersion noise. Furthermore, in order to demonstrate the efficiency and applicability of the proposed perfectly matched layer scheme, numerical modelling examples are also presented. The numerical results obtained with the put forward perfectly matched layer scheme are compared with results from traditional attenuation absorbing boundary conditions and enlarged models as well. The analysis of the numerical results indicates that the proposed perfectly matched layer scheme is significantly effective and more efficient in absorbing spurious reflections from the model boundaries.  相似文献   

19.
With respect to radioactive waste confinement indurated clays are potential host rocks as well as barrier rock formations. Regarding the tight coupling between fluid flow processes and time-dependent deformation in argillaceous rock mass this paper presents recent experimental findings obtained in the Mont Terri URL which are interpreted in terms of potential driving mechanisms for the time-dependent deformation behaviour of the rock. Within the interpretation process, numerical modelling work, reflecting the actual properties of the rock in situ, is performed and the simulation results are compared to the experimental results. The discussion also touches on the work which is in progress on validating the Finite Element Code for physical modelling and numerical simulation (CODE_BRIGHT).  相似文献   

20.
Assessment of the seismic performance of existing structures requires due consideration of both aleatory and epistemic sources of uncertainty; the former being typically associated with the randomness in ground motion records and the latter with the uncertainty in numerical modelling. Using a numerical modelling approach calibrated to available experimental test data collected from the literature, the uncertainty associated with different modelling parameters for existing reinforced concrete frames in Italy was quantified via an extensive numerical study. This was done to quantify the propagation of modelling parameter type uncertainty to the overall dispersion of the demand parameters typically used in seismic assessment, namely peak storey drift and peak floor accelerations. In addition, the impact of such modelling uncertainty on the median intensity and dispersion of the collapse fragility function was also examined. From the results of this study, empirical values of modelling parameter uncertainty were quantified with a view to being used in the assessment of existing reinforced concrete frames with masonry infill designed prior to the introduction of seismic design provisions in Italy during the 1970s. Comparing these empirical values to those available in the literature, it is seen how the fundamental behaviour of the frames differs from more modern frames with ductile detailing to the extent that values available in guidelines such as FEMA P58 cannot be reasonably adopted for these structural typologies.  相似文献   

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

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