共查询到20条相似文献，搜索用时 62 毫秒
1.
Vertically Integrated Hydraulic Conductivity: A New Parameter for Groundwater-Surface Water Analysis
Numerical modeling of groundwater-surface water interactions provides vital information necessary for determining the extent of nutrient transport, quantifying water budgets, and delineating zones of ecological support. The hydrologic data that drive these models are often collected at disparate scales and subsequently incorporated into numerical models through upscaling techniques such as piecewise constancy or geostatistical methods. However, these techniques either use basic interpolation methods, which often simplifies the system of interest, or utilize complex statistical methods that are computationally expensive, time consuming, and generate complex subsurface configurations. We propose a bulk parameter termed “vertically integrated hydraulic conductivity” (KV), and defined as the depth-integrated resistance to fluid flow sensed at the groundwater-surface water interface, as an alternative to hydraulic conductivity when investigating vertical fluxes across the groundwater-surface water interface. This bulk parameter replaces complex subsurface configurations in situations dominated by vertical fluxes and where heterogeneity is not of primary importance. To demonstrate the utility of KV, we extracted synthetic temperature time series data from a forward numerical model under a variety of scenarios and used those data to quantify vertical fluxes using the amplitude ratio method. These quantified vertical fluxes and the applied hydraulic head gradient were subsequently input into Darcy's Law and used to quantify KV. This KV was then directly compared to the equivalent hydraulic conductivity (KT) assuming an infinitely extending layer. Vertically integrated hydraulic conductivity allows for more accurate and robust flow modeling across the groundwater-surface water interface in instances where complex heterogeneities are not of primary concern. 相似文献
2.
Jon Albretsen 《Ocean Dynamics》2007,57(4-5):287-304
We perform eddy-permitting to eddy-resolving simulations of the Skagerrak/northern North Sea with a terrain-following numerical
ocean model. We demonstrate that realistic representations of freshwater input are not required when the focus is on modelling
mesoscale structures such as meanders and eddies. To arrive at this conclusion, we analyze the results using a recently developed
energy diagnostic scheme to study the sensitivity to realistic representations of the lateral freshwater flux provided to
the area from the Baltic Sea and by the major rivers. The scheme is suitable for analysis of growth of instabilities, and
it has four basic instability processes prominent. We recognize both horizontal and vertical shear instabilities. There are
two processes where average potential energy is converted to eddy kinetic energy, and they are related to the mean gradient
in surface elevation and the mean lateral density gradient, respectively. The latter process is known as frontal instability.
We demonstrate that the change in the eddy kinetic energy field is small, despite the large variations in the hydrographic
properties from experiment to experiment. Moreover, generation of eddy activity appears at the same locations and with approximately
the same strength regardless of actual representations of freshwater input. Furthermore, we find that vertical shear instability
dominates the energy conversion processes in the Norwegian Coastal Current. Finally, we find that the areas off the northwest
coast of Denmark recognized with enhanced eddy kinetic energy level is not caused by instability processes but eddy–eddy interaction
rooted in variations in the sea level. 相似文献
3.
The final part of the study is devoted to the process of rapid reconfiguration of the geomagnetotail after the upset of the equilibrium, i.e., to substorm activation. Such a version of spontaneous formation of nonlinear kinetic small-scale structures at the site occupied by an initial current sheet (CS), which results in the formation of a specific equilibrium current structure with strong ion anisotropy (forced kinetic current sheet, FKCS), has been revealed in the course of numerical simulation. This reconfiguration “channel” is realized when the relative value of the magnetic field component normal to CS is smaller than a certain critical value. In this case the disturbance intensity is SPONTANEOUSLY settled at a certain final level at the late stages. The obtained results are compared with the previous study of the system evolution under the action of the external trigger. The physical causes and specific features of different CS evolution regimes are considered. The studied plasma mechanisms are responsible for the effects of rapid energy transformation—magnetic field “annihilation”—that took place in the geomagnetotail. 相似文献
4.
Canonical balanced dynamic equations involving vertically sheared horizontal flow with heat or mass sources have emerged recently in systematic multi-scale modeling of the equatorial wave guide on a wide range of spatio-temporal scales. Here, a new self-contained derivation of these equations is developed briefly in a context for potential applications to the hurricane embryo. These canonical balanced equations are studied through a combination of exact solutions and simple numerics. The results below include elementary exact solutions given by velocity fields that are linear in the spatial coordinates combined with an exact nonlinear stability analysis for vertical vorticity amplification in such a preconditioned environment. Other elementary solutions studied here include the evolution of radial eddies, which represent “hot towers” in the hurricane embryo in a suitable radial preconditioned background environment. 相似文献
5.
The 1D version of the Model for Applications at Regional Scale is used to parameterize the effects of sea surface waves in
2D in a horizontally homogeneous offshore zone of the Iroise sea. Here we present the first simulation of the Iroise sea including
sea surface waves forcing, and more generally, the first study of a boundary layer including the Hasselmann force with a tidal
wave. We use a single equation turbulence closure based on a non-local diagnosis for energetic and dissipation length scales.
The turbulent energy flux at the surface due to whitecaps and the Hasselmann force induced by Stokes drift are assessed using
the whole sea surface waves spectrum given by the Wave Watch Third generation model. The ability of the parameterization to
reproduce surface currents over a period of 1 year (2007) is tested with high frequency radar using spectral and time-frequency
analysis. One problem with 1D modelling, corresponding to overestimation of current oscillating at inertial frequency is illustrated
by comparing 1D and 3D simulations. We found an overall improvement by including the Hasselmann force mainly within the bandwidth
of less than one cycle per day to one cycle per day for surface currents. Turbulence is induced by whitecaps decaying rapidly
below the ocean surface but the mixed layer below 40 m is deeper due to waves breaking on the sea surface. 相似文献
6.
Santu Das 《地球物理与天体物理流体动力学》2014,108(2):147-167
Damping of water waves by a vertical porous structure placed at some distance from a vertical wall is investigated within the framework of linear water wave theory. The rectangular porous structure is placed on a small rectangular elevation. An incident wave of small amplitude propagates through the structure – some portion gets reflected back while some portion gets transmitted to a third region bounded vertically by a rigid wall which is considered to be at a distance near the porous structure, and also away from the wall at a large distance as a separate case. Boundary value problems are set up in all three regions and, by using the matching conditions along the vertical boundaries, a system of linear equations is deduced. The roots of the relevant dispersion relation are used in setting up the system of equations. The overall scattering phenomenon is studied with respect to different relevant parameters. The dependence of the coefficients on the thickness (width) of the porous structure is investigated for different numbers of modes and porosity. It is observed that, except for the case when the porous structure is thin, the reflection and transmission coefficients give rise to values as expected. In the case, when the rigid wall is nearer to the structure, the reflection coefficient decreases rapidly for a thin structure and converges for all numbers of evanescent modes afterwards. The transmission coefficient also decreases as the width increases, ultimately converging and vanishing for a wide structure. When the wall is at a large distance away from the structure, the behavior of both the reflection and transmission coefficients remain the same. For both cases of the wall being nearer and away from the structure, higher porosity gives rise to lower reflection coefficients and higher transmission coefficients. However, the transmission coefficients converge and vanish when the porous structure is very wide. We also discuss the energy loss against the width of the porous structure for different values of number of modes and porosity. Irrespective of the positioning of the rigid wall, we observe that for higher values of porosity, energy loss is more pronounced when the structure is not thin, whereas energy loss is same for all numbers of modes. All our observations are supported by graphs. Good agreement of our result with earlier results justifies our model. 相似文献
7.
Waves, circulation and vertical dependence 总被引：1，自引：1，他引：0
George Mellor 《Ocean Dynamics》2013,63(4):447-457
Longuet-Higgins and Stewart (J Fluid Mech 13:481–504, 1962; Deep-Sea Res 11:529–562, 1964) and later Phillips (1977) introduced the problem of waves incident on a beach, from deep to shallow water. From the wave energy equation and the vertically integrated continuity equation, they inferred velocities to be Stokes drift plus a return current so that the vertical integral of the combined velocities was nil. As a consequence, it can be shown that velocities of the order of Stokes drift rendered the advective term in the momentum equation negligible resulting in a simple balance between the horizontal gradients of the vertically integrated elevation and wave radiation stress terms; the latter was first derived by Longuet-Higgins and Stewart. Mellor (J Phys Oceanogr 33:1978–1989, 2003a), noting that vertically integrated continuity and momentum equations were not able to deal with three-dimensional numerical or analytical ocean models, derived a vertically dependent theory of wave–circulation interaction. It has since been partially revised and the revisions are reviewed here. The theory is comprised of the conventional, three-dimensional, continuity and momentum equations plus a vertically distributed, wave radiation stress term. When applied to the problem of waves incident on a beach with essentially zero turbulence momentum mixing, velocities are very large and the simple balance between elevation and radiation stress gradients no longer prevails. However, when turbulence mixing is reinstated, the vertically dependent radiation stresses produce vertical velocity gradients which then produce turbulent mixing; as a consequence, velocities are reduced, but are still larger by an order of magnitude compared to Stokes drift. Nevertheless, the velocity reduction is sufficient so that elevation set-down obtained from a balance between elevation gradient and radiation stress gradients is nearly coincident with that obtained by the aforementioned papers. This paper includes four appendices. The first appendix demonstrates the numerical process by which Stokes drift is excluded from the turbulence stress parameterization in the momentum equation. A second appendix determines a bottom slope criterion for the application of linear wave relations to the derivation of the wave radiation stress. The third appendix explores the possibility of generalizing results by non-dimensionalization. The final appendix applies the basic theory to a problem introduced by Bennis and Ardhuin (J Phys Oceanogr 41:2008–2012, 2011). 相似文献
8.
A 3D, two-time-level, σS-z-σB hybrid-coordinate Marine Science and Numerical Modeling numerical ocean circulation model (HyMOM) is developed in this paper. In HyMOM, the σ coordinate is employed in the surface and bottom regions, and the z coordinate is used in the intermediate layers. This method can overcome problems with vanishing surface cells and minimize the unwanted deviation in representing bottom topography. The connection between the σ and z layers vertically includes an expanded “ghost” method and the linear interpolation. The governing equations in the σS-z-σB hybrid coordinate based on the complete Reynolds-averaged Navier-Stokes equations are derived in detail. The two-level time staggered and Eulerian forward and backward schemes, which are of second-order of accuracy, are adopted for the temporal difference in internal and external mode, respectively. The computation of the baroclinic gradient force is tested in an analytic test problem; the errors for two methods in HyMOM, which are relatively large only in the bottom layers, are obviously smaller than those in the pure σ and z models in almost all of the vertical layers. A quasi-global climatologic numerical experiment is constructed to test the simulation performance of HyMOM. With the monthly mean Levitus climatology data as reference, the HyMOM can improve the simulating accuracy compared with its pure z or σ coordinate implementation. 相似文献
9.
基于块石静、动室内三轴试验确定的广义塑性模型参数,对直立墙结构振动台试验进行有限元数值模拟,并与试验结果进行对比分析,进一步探讨直立墙结构在地震荷载作用下的破坏过程和破坏特征。计算表明:该模型可较合理地模拟直立墙结构的地震反应特性和破坏特征,计算结果与试验现象基本相符。位于抛石基床上的直立墙结构破坏模式为直立墙向外海侧的滑移、倾斜和竖向沉降,其破坏过程为:当输入加速度较小时,直立墙处于稳定状态;随着输入加速度逐渐增大,直立墙在自身惯性力和墙后回填块石的动土压力作用下缓慢向外海侧水平滑移、倾斜和竖向沉降,墙后回填块石出现沉陷,但变形较小;当加速度达到一定值时,直立墙向外海侧移动和回填块石沉陷速率急剧增加,变形较大。 相似文献
10.
A computational method of energy evaluation is derived to study the elastic responses and energy distribution of actively controlled single‐degree‐of‐freedom (SDOF) structures during earthquakes. Contrary to the common perception that applying active control force pumps energy into the structure, the applied control force can actually reduce the energy in the structure by reducing the input energy from earthquakes to the structure. In addition, applying control force can dissipate a large amount of energy in the structure when this control force is applied in the direction opposite to the displacement and velocity responses. To demonstrate this energy mechanism in active controlled structures, the two most popular control algorithms, optimal linear control (OLC) and instantaneous optimal control (IOC) algorithms, are used to calculate the control response and energy spectra. One‐step time delay is incorporated into the algorithms to take into consideration the practical aspect of active control. The effects of different earthquakes and damping ratios on control energy and response spectra are studied. These studies show that both OLC and IOC are very effective in reducing the structural displacement and velocity responses by reducing the input earthquake energy as well as dissipating a large amount of energy in the structure. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
11.
12.
强烈地震震源区和其周边地域的介质,在力源作用下开始破裂,即微破裂,而当其在应力不断积累、逐渐发展并形成具有一定规模的“破裂链”时,即应力集中达到临界状态时,震源介质所积累的巨大能量以波动的形式瞬间释放,并冲出地表发生地震,且在地表产生一系列的破坏和大型断裂.在强烈地震发生后依据设定的初始模型和远、近地震台网记录,通过震源机制解计算,可求得震源介质的破裂长度和破裂过程.事实表明,一系列大地震的“孕育”、发生和发展,由震源深处到地表均呈现出强烈的破裂效应,其破裂长度可由几公里、几十公里、乃至千余公里长.有的地震在发生前确亦存在着破裂响应和迹象.为此,若在地表和深井中能进行破裂效应的观测,以捕捉初始微破裂和其形成“破裂链”的动力过程,这对于短临地震预测可能是一个十分重要的途径. 相似文献
13.
《水文科学杂志》2013,58(3):502-506
Abstract Currently, environmental modelling is frequently conducted with the aid of artificial neural networks (ANNs) in an effort to achieve greater accuracy in simulation and forecasting beyond that typically obtained when using solely linear models. For the design of an ANN, modellers must contend with two key issues: (a) the selection of model input and (b) the determination of the number of hidden neurons. A novel approach is introduced to address the optimal design of ANNs based on a multi-objective strategy that enables the user to find a set of feasible ANNs, determined as optimal trade-off solutions between model simplicity and accuracy. This is achieved in a multi-objective fashion by simultaneously minimizing three different cost functions: the model input dimension, the hidden neuron number and the generalization error computed on a validation set of data. The multi-objective approach is based on the Pareto dominance criterion and an evolutionary strategy has been employed to solve the combinatorial optimization problem. From a theoretical perspective, the choice of a multi-objective approach marks an attempt to account for, and overcome, the “curse of dimensionality” and to circumvent the drawbacks of “overfitting” that are inherent in ANNs. Moreover, it is demonstrated that the strategy renders the choice of the ANN more robust, as is evident by “unseen data” in the testing stage, since structure determination is not merely based on the statistical evaluation of the generalization performance. The methodology is tested and the results are reported in a case study relating groundwater level predictions to total monthly rainfall. 相似文献
14.
Field test of sub‐basalt hydrocarbon exploration with marine controlled source electromagnetic and magnetotelluric data 
下载免费PDF全文
下载免费PDF全文 David Myer David Alumbaugh Oliver Hermann Randall Hobbet 《Geophysical Prospecting》2015,63(5):1284-1310
The recent use of marine electromagnetic technology for exploration geophysics has primarily focused on applying the controlled source electromagnetic method for hydrocarbon mapping. However, this technology also has potential for structural mapping applications, particularly when the relative higher frequency controlled source electromagnetic data are combined with the lower frequencies of naturally occurring magnetotelluric data. This paper reports on an extensive test using data from 84 marine controlled source electromagnetic and magnetotelluric stations for imaging volcanic sections and underlying sediments on a 128‐km‐long profile. The profile extends across the trough between the Faroe and Shetland Islands in the North Sea. Here, we focus on how 2.5D inversion can best recover the volcanic and sedimentary sections. A synthetic test carried out with 3D anisotropic model responses shows that vertically transverse isotropy 2.5D inversion using controlled source electromagnetic and magnetotelluric data provides the most accurate prediction of the resistivity in both volcanic and sedimentary sections. We find the 2.5D inversion works well despite moderate 3D structure in the synthetic model. Triaxial inversion using the combination of controlled source electromagnetic and magnetotelluric data provided a constant resistivity contour that most closely matched the true base of the volcanic flows. For the field survey data, triaxial inversion of controlled source electromagnetic and magnetotelluric data provides the best overall tie to well logs with vertically transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data a close second. Vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data provided the best interpreted base of the volcanic horizon when compared with our best seismic interpretation. The structural boundaries estimated by the 20‐Ω·m contour of the vertical resistivity obtained by vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data gives a maximum geometric location error of 11% with a mean error of 1.2% compared with the interpreted base of the volcanic horizon. Both the model study and field data interpretation indicate that marine electromagnetic technology has the potential to discriminate between low‐resistivity prospective siliciclastic sediments and higher resistivity non‐prospective volcaniclastic sediments beneath the volcanic section. 相似文献
15.
《地震工程与结构动力学》2018,47(5):1250-1269
Reinforced concrete waffle‐flat plate (WFP) structures present 2 important drawbacks for use as a main seismic resisting system: low lateral stiffness and limited ductility. Yet the former can serve a positive purpose when, in parallel, the flexible WFP structure is combined with a stiff system lending high‐energy dissipation capacity, to form a “flexible‐stiff mixed structure.” This paper experimentally investigates the seismic performance of WFP structures (flexible system) equipped with hysteretic dampers (stiff system) through shake‐table tests conducted on a 2/5‐scale test specimen. The WFP structure was designed only for gravitational loads. The lateral strength and stiffness provided by the dampers at each story were, respectively, about 3 and 7 times greater than those of the bare WFP structure. The mixed system was subjected to a sequence of seismic simulations representing frequent to very rare ground motions. Under the seismic simulations associated with earthquakes having return periods ranging from 93 to 1894 years, the WFP structure performed in the level of “immediate occupancy,” with maximum interstory drifts up to about 1%. The dampers dissipated most (75%) of the energy input by the earthquake. 相似文献
16.
Vertical Discretization Impact in Numerical Modeling of Matrix Diffusion in Contaminated Groundwater
Shahla K. Farhat David T. Adamson Arun R. Gavaskar Sophia A. Lee Ronald W. Falta Charles J. Newell 《Ground Water Monitoring & Remediation》2020,40(2):52-64
Understanding the effects of contaminants that can diffuse into low-permeability (“low-k”) zones is crucial for effective groundwater remedial decision-making. Because low-k zones can serve as low-level sources of contamination to more transmissive zones over time, an accurate evaluation of the impacts of matrix diffusion at contaminated sites is vital. This study compared numerical groundwater flow and transport simulations using MODFLOW/RT3D at a hypothetical site using three cases, each with increasing discretization of the vertical 10-m thick domain: (1) a coarse multilayer heterogeneous grid based on one layer for each of four different hydrogeological units, (2) a “low-resolution” discretization approach where the low-k units were divided into several sublayers giving the model 10 layers, and (3) a “high-resolution” numerical model with 199 layers that are a few centimeters thick. When comparing the results of each case, significant differences were observed between the discretizations used, even though all other model input data were identical. The conventional grid models (Cases 1 and 2) appeared to underestimate groundwater plume concentrations by a factor ranging from 1.1 to 36 when compared to the high-resolution grid model (Case 3), and underestimated predicted cleanup times by more than a factor of 10 for some of the hypothetical sampling points in the modeling domain. These results validate the implication of Chapman et al. (2012), that conventional vertical discretization of numerical groundwater flow and transport models at contaminated sites (with layers that are greater than 1 m thick) can lead to significant errors when compared to more accurate high-resolution vertical discretization schemes (layers that are centimeters thick). 相似文献
17.
The purpose of this paper is to present an accurate and efficient model for use in explicit soil-structure interaction analyses for seismic excitation. The main body of the paper is concerned with the implementation of non-reflecting boundaries at the base and vertical faces of a two-dimensional finite element or finite difference soil mesh. The paper proposes a scheme for implementing viscous dashpots as energy absorbers at the base of the model where the seismic excitation is applied. To achieve this, the form of the seismic input is modified from the conventional accelerogram to a force time history that can be obtained from a standard deconvolution program with very minor modifications. The proposed scheme is also applicable to frequency-domain solutions. For the lateral boundaries, a superposition non-reflecting boundary formulation is recommended. The paper shows how easily the standard formulation can be modified to accommodate a seismic excitation that is assumed to be vertically propagating from the base of the model. An example is presented to demonstrate the accuracy of the proposed model. 相似文献
18.
Several field studies in bays and estuaries have revealed pronounced subsurface maxima in the vertical profiles of the current amplitude of the principal tidal harmonic, or of its vertical shear, over the water column. To gain fundamental understanding about these phenomena, a semi-analytical model is designed and analysed, with focus on the sensitivity of the vertical structure of the tidal current amplitude to formulations of the vertical shape of the eddy viscosity. The new analytical solutions for the tidal current amplitude are used to explore their dependence on the degree of surface mixing, the vertical shape of eddy viscosity in the upper part of the water column and the density stratification. Sources of surface mixing are wind and whitecapping. Results show three types of current amplitude profiles of tidal harmonics, characterised by monotonically decreasing shear towards the surface, “surface jumps” (vertical shear of tidal current amplitude has a subsurface maximum) and “subsurface jets” (maximum tidal current amplitude below the surface), respectively. The “surface jumps” and “subsurface jets” both occur for low turbulence near the surface, whilst additionally the surface jumps only occur if the eddy viscosity in the upper part of the water column decreases faster than linearly to the surface. Furthermore, “surface jumps” take place for low density stratification, while and “subsurface jets” occur for high density stratification. The physics causing the presence of surface jumps and subsurface jets is also discussed. 相似文献
19.
Near-field horizontal and vertical earthquake ground motions 总被引：8，自引：0，他引：8
Strong-motion attenuation relationships are presented for peak ground acceleration, spectral acceleration, energy density, maximum absolute input energy for horizontal and vertical directions and for the ratio of vertical to horizontal of these ground motion parameters. These equations were derived using a worldwide dataset of 186 strong-motion records recorded with 15 km of the surface projection of earthquakes between Ms=5.8 and 7.8. The effect of local site conditions and focal mechanism is included in some of these equations. 相似文献