首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Power spectrum analysis of the Bouger gravity values in the Eastern Alps suggests that the gravity field may be separated into long and short wavelength components. The long wavelength component is assumed to be caused by Alpine crustal thickening. This long wavelength component was subjected to a gravimetric single density-interface inversion procedure, giving a gravimetric Mohorovi?i? model which is generally in good agreement with Moho-depths derived by refraction and reflection seismology. The residual high-frequency gravity component correlates well with the main surface geological units in the Eastern Alps. Apparent density mapping by applying an inverse density deconvolution filter to the short wavelength gravity component gives density values for the upper crust which correspond well with averaged density values from rock samples.  相似文献   

2.
This study suggested a numerical model using the Tabu search algorithm along with the Adjoint State method to identify the hydrogeological characteristics of an anisotropic groundwater aquifer. The Tabu search algorithm was applied to identify the anisotropic transmissivity components to avoid a local optimum. Then, the Adjoint State method was used to calculate the sensitivity of the parameters in order to increase the efficiency of the optimization. For an anisotropic and homogeneous aquifer, results showed that the optimal procedure presented combining the Tabu search algorithm and the Adjoint State method might successfully identify the values of the transmissivity components. If the duration of the pumping test was long enough (12‐h pumping test), the value of the transmissivity components could be optimized with type‐curve, straight‐line, and Tabu search methods, along with the Adjoint State methods. If the duration of the pumping test was short (0·5‐h pumping test), the Tabu search method, along with the Adjoint State method proposed herein, might successfully optimize the transmissivity components. For an anisotropic but heterogeneous aquifer, results showed that the suggested optimal procedure still successfully identified the values of the transmissivity components. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

3.
A two‐parameter transfer function with an infinite characteristic time is proposed for conceptual rainfall–runoff models. The large time behaviour of the unit response is an inverse power function of time. The infinite characteristic time allows long‐term memory effects to be accounted for. Such effects are observed in mountainous and karst catchments. The governing equation of the model is a fractional differential equation in the limit of long times. Although linear, the proposed transfer function yields discharge signals that can usually be obtained only using non‐linear models. The model is applied successfully to two catchments, the Dud Koshi mountainous catchment in the Himalayas and the Durzon karst catchment in France. It compares favourably to the linear, non‐linear single reservoir models and to the GR4J model. With a single reservoir and a single transfer function, the model is capable of reproducing hysteretic behaviours identified as typical of long‐term memory effects. Computational efficiency is enhanced by approximating the infinite characteristic time transfer function with a sum of simpler, exponential transfer functions. This amounts to partitioning the reservoir into several linear sub‐reservoirs, the output discharges of which are easy to compute. An efficient partitioning strategy is presented to facilitate the practical implementation of the model. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

4.
Migration velocity analysis, a method for determining long wavelength velocity structure, is a critical step in prestack imaging. Solution of this inverse problem is made difficult by a multimodal objective function; a parameter space often vast in extent; and an evaluation procedure for candidate solutions, involving the calculation of depth-migrated image gathers, that can be prohibitively expensive. Recognizing the global nature of the problem, we employ a genetic algorithm (GA) in the search for the optimum velocity model. In order to describe a model efficiently, regions of smooth variation are identified and sparsely parametrized. Region boundaries are obtained via map migration of events picked on the zero-offset time section. Within a region, which may contain several reflectors, separate components describe long and short wavelength variations, eliminating from the parameter space, models with large velocity fluctuations. Vital to the success of the method is rapid model evaluation, achieved by generating image gathers only in the neighbourhood of specific reflectors. Probability of a model, which we seek to maximize, is derived from the flatness of imaged events. Except for an initial interpretation of the zero-offset time section, our method is automatic in that it requires no picking of residual moveout on migrated gathers. Using an example data set from the North Sea, we show that it is feasible to solve for all velocity parameters in the model simultaneously: the method is global in this respect also.  相似文献   

5.
Hydrological trend analysis with innovative and over-whitening procedures   总被引:1,自引:0,他引:1  
Different statistical methodologies can be employed to identify possible trend components in any hydro-meteorological time series. A pre-whitening (P-W) procedure has been suggested to reduce the serial correlation effect on Mann-Kendall (M-K) trend analysis. In this paper, instead of P-W, an over-whitening (O-W) procedure is suggested, which generates serially independent series with the same trend slope value. Analytically necessary formulations for O-W are presented with a non-parametric but simple innovative trend assessment procedure, which are supported by extensive simulation studies. The applications of the methodology are presented for eight factual time series records from tropical, temperate and arid regions including temperature, rainfall, streamflow, relative humidity and CO2 concentrations for different short and long durations. Relative humidity and CO2 records are monthly time series and, hence, there are trend and periodicity components. It is noticed in all cases that the natural trends remain as they were after the O-W procedure, thus providing an opportunity to determine reliably the trends embedded even in the serially dependent series. The O-W procedure is applicable even in the cases of periodicity in the original records.  相似文献   

6.
Introduction The history of research work on the response of buried pipeline under the fault movement hasbeen about 30 years. Several simplified design methods have been proposed to obtain the maxi-mum stress or strain in pipe. These methods include the theoretical method and the finite elementmethod (FEM). In the theoretical method, the pipe is usually modeled as a cable (Newmark, Hall1975; Kennedy, et al, 1977) or a beam (Wang, Wang, 1995; LIU, ZHANG, 2002). These theoreti-cal me…  相似文献   

7.
This study aims to develop an improved time series model to overcome difficulties in modeling monthly short term stream flows. The periodic, serial dependent and independent components of the classical time series models are improved separately by information transfer from a surrounding long term gauging station to the considered flow section having short term records. Eventually, an improved model preserving the mathematical model structure of the classical time series model, while improving general and monthly statistics of the monthly stream flows, is derived by using the improved components instead of the short term model components in the time series modeling. The correlative relationships between the current short term and surrounding long term stations are used to improve periodic and serial dependent behaviors of monthly flows. Independent components (residuals) are improved via the parameters defining their theoretical probability distribution. The improved model approach is tested by using 50 year records of Göksu-Himmetli (1801) and Göksu-Gökdere (1805) flow monitoring stations located on the Ceyhan river basin, in south of Turkey. After 50 year records of the station 1801 are separated into five 10 year sub series, their improved and classical time series models are computed and compared with the real long-term (50 year) time series model of this station to reveal efficiencies of the improved models for each subseries (sub terms with 10 year observation). The comparisons are realized based on the model components, model estimates and general/monthly statistics of model estimates. Finally, some evaluations are made on the results compared to the regression method classically applied in the literature.  相似文献   

8.
The surface water quality monitoring is an important concern of public organizations due to its relevance to the public health. Statistical methods are taken as consistent and essential tools in the monitoring procedures in order to prevent and identify environmental problems. This work presents the study case of the hydrological basin of the river Vouga, in Portugal. The main goal is discriminate the water monitoring sites using the monthly dissolved oxygen concentration dataset between January 2002 and May 2013. This is achieved through the extraction of trend and seasonal components in a linear mixed-effect state space model. The parameters estimation is performed with both maximum likelihood method and distribution-free estimators in a two-step procedure. The application of the Kalman smoother algorithm allows to obtain predictions of the structural components as trend and seasonality. The water monitoring sites are discriminated through the structural components by a hierarchical agglomerative clustering procedure. This procedure identified different homogenous groups relatively to the trend and seasonality components and some characteristics of the hydrological basin are presented in order to support the results.  相似文献   

9.
A new shell finite element method (FEM) model with an equivalent boundary is presented for estimating the response of a buried pipeline under large fault movement. The length of affected pipeline under fault movement is usually too long for a shell-mode calculation because of the limitation of memory and time of computers. In this study, only the pipeline segment near fault is modeled with plastic shell elements to study the local buckling and the large section deformation in pipe. The material property of pipe segment far away from the fault is considered as elastic, and nonlinear spring elements at equivalent boundaries are obtained and applied to two ends of shell model. Compared with the fixed-boundary shell model, the shell model with an equivalent boundary proposed by the study can remarkably reduce the needed memory and calculating time. Foundation item: National Natural Sciences Foundation of China (50078049). Contribution No. 04FE1017, Institute of Geophysics, China Earthquake Administration.  相似文献   

10.
The use of a neuro‐fuzzy approach is proposed to model the dynamics of entrainment of a coarse particle by rolling. It is hypothesized that near‐bed turbulent flow structures of different magnitude and duration or frequency and energy content are responsible for the particle displacement. A number of Adaptive Neuro‐Fuzzy Inference System (ANFIS) architectures are proposed and developed to link the hydrodynamic forcing exerted on a solid particle to its response, and model the underlying nonlinear dynamics of the system. ANFIS combines the advantages of fuzzy inference (If‐Then) rules with the power of learning and adaptation of the neural networks. The model components and forecasting procedure are discussed in detail. To demonstrate the model's applicability for near‐threshold flow conditions an example is provided, where flow velocity and particle displacement data from flume experiments are used as input and output for the training and testing of the ANFIS models. In particular, a Laser Doppler velocimeter (LDV) is employed to obtain long records of local streamwise velocity components upstream of a mobile exposed particle. These measurements are acquired synchronously with the time history of the particle's position detected by a setup including a He‐Ne laser and a photodetector. The representation of the input signal in the time and frequency domain is implemented and the best performing models are found capable of reproducing the complex dynamics of particle response. Following a trial and error approach the different models are compared in terms of their efficiency and forecast accuracy using a number of performance indices. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

11.
Complexity in simulating the hydrological response in large watersheds over long times has prompted a significant need for procedures for automatic calibration. Such a procedure is implemented in the basin‐scale hydrological model (BSHM), a physically based distributed parameter watershed model. BSHM simulates the most important basin‐scale hydrological processes, such as overland flow, groundwater flow and stream–aquifer interaction in watersheds. Here, the emphasis is on estimating the groundwater parameters with water levels in wells and groundwater baseflows selected as the calibration targets. The best set of parameters is selected from within plausible ranges of parameters by adjusting the values of hydraulic conductivity, storativity, groundwater recharge and stream bed permeability. The baseflow is determined from stream flow hydrographs by using an empirical scheme validated using a chemical approach to hydrograph separation. Field studies determined that the specific conductance for components of the composite hydrograph were sufficiently unique to make the chemical approach feasible. The method was applied to the Big Darby Creek Watershed, Ohio. The parameter set selected for the groundwater system provides a good fit with the estimated baseflow and observed water well data. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

12.
Abstract

We present a procedure for estimating Q95 low flows in both gauged and ungauged catchments where Q95 is the flow that is exceeded 95% of the time. For each step of the estimation procedure, a number of alternative methods was tested on the Austrian data set by leave-one-out cross-validation, and the method that performed best was used in the final procedure. To maximise the accuracy of the estimates, we combined relevant sources of information including long streamflow records, short streamflow records, and catchment characteristics, according to data availability. Rather than deriving a single low flow estimate for each catchment, we estimated lower and upper confidence limits to allow local information to be incorporated in a practical application of the procedure. The components of the procedure consist of temporal (climate) adjustments for short records; grouping catchments into eight seasonality-based regions; regional regressions of low flows with catchment characteristics; spatial adjustments for exploiting local streamflow data; and uncertainty assessment. The results are maps of lower and upper confidence limits of low flow discharges for 21 000 sub-catchments in Austria.  相似文献   

13.
Sufficient low‐frequency information is essential for full‐waveform inversion to get the global optimal solution. Multi‐scale envelope inversion was proposed using a new Fréchet derivative to invert the long‐wavelength component of the model by directly using the low‐frequency components contained in an envelope of seismic data. Although the new method can recover the main structure of the model, the inversion quality of the model bottom still needs to be improved. Reflection waveform inversion reduces the dependence of inversion on low‐frequency and long‐offset data by using travel‐time information in reflected waves. However, when the underground medium contains strong contrast or the initial model is far away from the true model, it is hard to get reliable reference reflectors for the generation of reflected waves. Here, we propose a combination inversion algorithm, i.e., reflection multi‐scale envelope inversion, to overcome the limitations of multi‐scale envelope inversion and reflection waveform inversion. First, wavefield decomposition was introduced into the multi‐scale envelope inversion to improve the inversion quality of the long‐wavelength components of the model. Then, after the initial model had been established to be accurate enough, migration and de‐migration were introduced to achieve multi‐scale reflection waveform inversion. The numerical results of the salt‐layer model and the SEG/EAGE salt model verified the validity of the proposed approach and its potential.  相似文献   

14.
Stratification is incorporated into an unsteady model of shelf currents by splitting the dynamic response of the flow into two parts, each with its own time scale. The barotropic part of the response is independent of depth and varies rapidly on a short time scale, whereas the baroclinic part depends on depth and changes slowly with time on a long time scale.The three-dimensional model has a continental shelf sloping down from an eastern boundary to the deep ocean. The equations for the barotropic component of the pressure field contain forcing by the wind stress and feedback from the baroclinic field. An integral of the heat equation over the long time scale determines the slow changes in the temperature field and hence in the baroclinic component of the velocity distribution.The temperature field is specified at the start of the numerical calculation. Its subsequent development is controlled by the numerical procedure. It is found that significant changes in the temperature field require a long period of upwelling favourable winds, whereas the longshore currents react more quickly to changes in the wind stress.  相似文献   

15.
Wind in Ireland: long memory or seasonal effect?   总被引:1,自引:1,他引:0  
Since Haslett and Raftery’s paper Space-Time Modelling with Long-Memory Dependence: Assessing Ireland’s Wind Power Resource (1989), modelling meteorological time series with long memory processes, in particular the ARFIMA model has become very common. Haslett and Raftery fitted an ARFIMA model on Irish daily wind speeds. In this paper, we try to reproduce Haslett and Raftery’s results (focusing on the dynamic of the wind process, and not on cross-correlation and space dependencies), and show that an ARFIMA model does not properly capture the behaviour of the series (in Modelling daily windspeed in Ireland section). Indeed, the series show a periodic behaviour, that is not taken into account by the ARFIMA model. Removing this periodic behaviour yields no results either, we therefore try to fit a GARMA model that takes into account both seasonality and long memory (in Seasonality and long memory using GARMA models section). If a GARMA process can be fitted to the data to model Irish daily data, we will show that these models could also be used to model Dutch hourly data.
Arthur CharpentierEmail:
  相似文献   

16.
Integrating migration velocity analysis and full waveform inversion can help reduce the high non‐linearity of the classic full waveform inversion objective function. The combination of inverting for the long and short wavelength components of the velocity model using a dual objective function that is sensitive to both components is still very expensive and have produced mixed results. We develop an approach that includes both components integrated to complement each other. We specifically utilize the image to generate reflections in our synthetic data only when the velocity model is not capable of producing such reflections. As a result, we get the migration velocity analysis working when we need it, and we mitigate its influence when the velocity model produces accurate reflections (possibly first for the low frequencies). This is achieved using a novel objective function that includes both objectives. Applications to a layered model and the Marmousi model demonstrate the main features of the approach.  相似文献   

17.
山西断陷盆地带形成机制的初步数值模拟   总被引:2,自引:0,他引:2       下载免费PDF全文
本文在分析山西断陷盆地带的各种地质、地球物理资料的基础上,用三维粘-弹塑性有限元方法对该断陷带的形成机制进行了模拟计算。表明山西断陷盆地带的形成是地壳深部的上隆作用和北东—南西向挤压应力场在深部隐伏断裂上派生的右旋滑动共同作用的结果。地壳深部的上隆是一个主要的不可忽视的因素  相似文献   

18.
Global optimization methods such as simulated annealing, genetic algorithms and tabu search are being increasingly used to solve groundwater remediation design and parameter identification problems. While these methods enjoy some unique advantages over traditional gradient based methods, they typically require thousands to tens of thousands of forward simulation runs before reaching optimal or near-optimal solutions. Thus, one severe limitation associated with these global optimization methods is very long computation time. To mitigate this limitation, this paper presents a new approach for obtaining, repeatedly and efficiently, the solutions of a linear forward simulation model subject to successive perturbations. The proposed approach takes advantage of the fact that successive forward simulation runs, as required by a global optimization procedure, usually involve only slight changes in the coefficient matrices of the resultant linear equations. As a result, the new solution to a system of linear equations perturbed by the changes in aquifer properties and/or sinks/sources can be obtained as the sum of a non-perturbed base solution and the solution to the perturbed portion of the linear equations. The computational efficiency of the proposed approach arises from the fact that the perturbed solution can be derived directly without solving the linear equations again. A two-dimensional test problem with 20 by 30 nodes demonstrates that the proposed approach is much more efficient than repeatedly running the simulation model, by more than 15 times after a fixed number of model evaluations. The ratio of speedup increases with the number of model evaluations and also the size of simulation model. The main limitation of the proposed approach is the large amount of computer memory required to store the inverse matrix. Effective ways for limiting the storage requirement are briefly discussed.  相似文献   

19.
本文将随机振动的虚拟激励法与拱坝-地基动力相互作用FE-BE-IBE时域模型结合,发展了一个可以考虑多维随机地震动作用下的拱坝动力响应计算模型,并用Monte Garlo方法对模型进行了验证,计算结果表明,地震动分量的相关性对结构的动力响应存在一定影响,合理考虑地震动各方向分量的相关性可以更好地计算实际地震作用下的拱坝动力响应。  相似文献   

20.
Abstract

The automatic calibration is done not by a hill-top climbing method but by a trial and error method carried out automatically by a computer program. The feedback procedure is made by comparing some criteria obtained from the observed hydrograph and the calculated hydrograph output from the working tank model. The two criteria are discharge volume and the shape of the hydrograph. The feedbacks of these two criteria correspond to dispacement feedback and velocity feedback in automatic control. The output of the working tank model is composed of components, the outputs from each of the tanks. Correspondingly, the whole period is divided into subperiods, in each of which each of the components plays the main part. The volume and shape are calculated in each subperiod and are used for the adjustment of the respective tanks. The feedback procedure starts from some initial model and converges very quickly after several (usually less than 15) iterations, and the result obtained is very good.  相似文献   

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

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