首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 296 毫秒
1.
In a special opportunity, detailed measurements of the flow in an overbank flow in the Flood Channel Facility at HR Wallingford were used in conjunction with tracer test data to assess the effectiveness of dispersion models based around random particle tracking (RPT). Ten different RPT models based on different assumptions and levels of information about the nature of the Lagrangian velocity field were investigated. Multiple simulations were used to calibrate variable parameters controlling the average magnitude of the perturbations for each model by comparison with observed concentrations at one cross-section. The calibrated models were then used to predict concentration distributions further downstream. Several of the calibrated models showed close agreement between observed and predicted concentration distributions. The most complex models using the most information about the velocity distributions were no better (and in some cases worse) in prediction than the simplest models investigated. It would appear that our knowledge of the system, despite the quality of the experiments, is too uncertain to infer a precise model structure.  相似文献   

2.
Part 1 of this study assessed the effectiveness of dispersion models based around random particle tracking (RPT) applied an overbank flow in the Flood Channel Facility at HR Wallingford. Several of the RPT models that were used predicted the mixing behaviour of the tracer in the complex channel flow to within an acceptable accuracy. If there is no uniqueness in the model structure which can be used to represent the true system to within the limitations of the available observations, then this implies an inherent degree of ambiguity in our knowledge of the physically based model structure. This suggests that we should be less forthright in the optimisation of each individual model structure, and perhaps investigate more of the parameter combinations for each model which yield feasible simulations of the system. An alternative fuzzy calibration technique is introduced which avoids the optimisation process and takes account of uncertainties in the model structure, parameter sets and observed data in prediction.  相似文献   

3.
This paper illustrates how the acoustic Doppler current profiler (ADCP) and single-beam echo-sounder (SBES) recordings can be used for the calibration of existing software to assist in generalizing the morphodynamic processes in large rivers at key sites such as bifi.trcations and confluences. Calibration of the MIKE21C numerical model by the Danish Hydraulic Institute at the 25-km-long reach of Lower Paran~ near Rosario (Argentina) is presented. This reach includes two downstream confluences and two bifurcations. The model simulates a 2-D depth-averaged flow velocity and the related sediment fluxes to predict the bifurcation morphodynamics that affects the Paranh waterway. To investigate the river channel bathymetry, roughness, flow discharge allocation at bifurcations, suspended sediment concentration and grain size distributions, several instruments were used. These instruments included two ADCPs by Teledyne RDI working at frequencies of 600 and 1,200 kHz, a Sontek ADCP working at a frequency of 1,000 kHz and a SBES. The method to assess suspended sediment concentration and grain size distributions has been previously described. This paper focuses primarily on investigating dune morphology (by means of SBES depth measurements) and friction velocity (by means of ADCP profiling) to determine the river channel bed-roughness. The 2-D model results agree with observed values of bed-roughness, flow velocity and suspended sediment concentration distributions at the investigated sections, known data of water slope and total load of bed sediment are in good agreement with model results.  相似文献   

4.
Measurements of seismic anisotropy in fractured rock are used at present to deduce information about the fracture orientation and the spatial distribution of fracture intensity. Analysis of the data is based upon equivalent-medium theories that describe the elastic response of a rock containing cracks or fractures in the long-wavelength limit. Conventional models assume frequency independence and cannot distinguish between microcracks and macrofractures. The latter, however, control the fluid flow in many subsurface reservoirs. Therefore, the fracture size is essential information for reservoir engineers. In this study we apply a new equivalent-medium theory that models frequency-dependent anisotropy and is sensitive to the length scale of fractures. The model considers velocity dispersion and attenuation due to a squirt-flow mechanism at two different scales: the grain scale (microcracks and equant matrix porosity) and formation-scale fractures. The theory is first tested and calibrated against published laboratory data. Then we present the analysis and modelling of frequency-dependent shear-wave splitting in multicomponent VSP data from a tight gas reservoir. We invert for fracture density and fracture size from the frequency dependence of the time delay between split shear waves. The derived fracture length matches independent observations from borehole data.  相似文献   

5.
The movement of unconsolidated materials near the Earth's surface is often driven by disturbances that occur at a range of spatial and temporal scales. The nature of these disturbances ranges from highly variable, such as tree turnover, to periodic and predictable, such as frost heave or creep. To explore the effect of probabilistic disturbances on surface processes, we formulated a granular creep model with analogy to rate process theory (RPT) used for chemical reactions. According to the theory, individual particles must be energized to a height greater than adjacent particles in order for grain dilation and transport to occur. The height of neighbouring particles (which is akin to activation energy in chemical reactions) varies with slope angle such that energy barriers get smaller in the downslope direction as slopes steepen. When slopes approach the friction‐limited angle of repose, the height of energy barriers approaches zero and grains ?ow in the absence of disturbance. An exponential function is used to describe the probability distribution of particle excitation height although alternative distributions are possible. We tested model predictions of granular dynamics in an experimental sandpile. In the sandpile, acoustic energy serves as the disturbance agent such that grains dilate and shear in response. Particle velocities are controlled by the frequency of energy pulses that result in grain displacement. Using tracer particles, we observed a convex‐upward velocity pro?le near the surface of the sandpile, consistent with predictions of our RPT‐based velocity model. In addition, we depth‐integrated the velocity model to predict how ?ux rates vary with inclination of the sandpile and observed non‐linear ?ux–gradient curves consistent with model predictions. By varying the acoustic energy level in the experimental sandpile, we documented changes in the rate of grain movement; similar changes in modelled velocities were achieved by varying the exponent of the particle excitation probability distribution. The general agreement between observed and modelled granular behaviour in our simple laboratory sandpile supports the utility of RPT‐based methods for modelling transport processes (e.g. soil creep, frost heave, and till deformation), thus enabling us to account for the probabilistic nature of disturbances that liberate sediment in natural landscapes. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

6.
Influence of calibration methodology on ground water flow predictions   总被引:1,自引:0,他引:1  
We constructed a numerical model of transient ground water flow and solute transport for a portion of the Biscayne Aquifer in Florida, and calibrated the model with three different combinations of data from a 193-day period: head (h) data alone, data on h and ground water discharge to a canal (q), and data on h, q, and ground water chloride concentration (C). We used each of the three calibrated models to predict h and q during a 182-day test period separate from the calibration period. All three calibrated models predicted h equally well during the test period (r = 0.95, where r = 1 indicates perfect agreement between measured and simulated values), though the model calibrated on h alone had significantly different parameter values than the other two models. Predictions of q during the test period depended on calibration methodology; models calibrated with multiple targets simulated q more accurately than the model calibrated on h alone (r = 0.79 compared to r = 0.49). Based on the results of these simulations, we conclude: (1) Post-calibration prediction is important in assessing the value of different data types in automated calibration; (2) inverse-solution uniqueness is not a requirement for accurate h predictions; (3) relatively simple models can predict with reasonable accuracy transient ground water flow in a complex aquifer, and parameters governing this prediction can be estimated by nonlinear regression methods that incorporate both h and q data; (4) addition of C data to the calibration did not improve model predictive capacity because the information in the C data was similar to that in the q data, from the perspective of model calibration (the subsurface chemical signal in question was controlled mainly by seepage of high-chloride canal water into the low-chloride ground water system).  相似文献   

7.
Estimating overland flow erosion capacity using unit stream power   总被引:2,自引:0,他引:2  
Soil erosion caused by water flow is a complex problem. Both empirical and physically based approaches were used for the estimation of surface erosion rates. Their applications are mainly limited to experimental areas or laboratory studies. The maximum sediment concentration overland flow can carry is not considered in most of the existing surface erosion models. The lack of erosion capacity limitation may cause over estimations of sediment concentration. A correlation analysis is used in this study to determine significant factors that impact surface erosion capacity. The result shows that the unit stream power is the most dominant factor for overland flow erosion which is consistent with experimental data. A bounded regression formula is used to reflect the limits that sediment concentration cannot be less than zero nor greater than a maximum value. The coefficients used in the model are calibrated using published laboratory data. The computed results agree with laboratory data very well. A one dimensional overland flow diffusive wave model is used in conjunction with the developed soil erosion equation to simulate field experimental results. This study concludes that the non-linear regression method using unit stream power as the dominant factor performs well for estimating overland flow erosion capacity.  相似文献   

8.
Ground-motion models (GMMs) are widely used in probabilistic seismic hazard analysis (PSHA) to estimate the probability distributions of earthquake-induced ground-motion intensity measures (IMs) at a site, given an earthquake of a certain magnitude occurring at a nearby location. Accounting for spatial and cross-IM correlations in earthquake-induced ground motions has important implications on probabilistic seismic hazard and loss estimates. This study first develops a new Italian GMM with spatial correlation for 31 amplitude-related IMs, including peak ground acceleration (PGA), peak ground velocity (PGV), and 5%-damped elastic pseudo-spectral accelerations (PSAs) at 29 periods ranging from 0.01 to 4 seconds. The model estimation is performed through a recently developed one-stage nonlinear regression algorithm proposed by the authors, known as the Scoring estimation approach. In fact, current state-of-practice approaches estimate spatial correlation separately from the GMM estimation, resulting in inconsistent and statistically inefficient estimators of interevent and intraevent variances and parameters in the spatial correlation model. We test whether this affects the subsequent cross-IM correlation analysis. To this aim, based on the newly developed GMM, the empirical correlation coefficients from interevent and intraevent residuals are investigated. Finally, a set of analytical correlation models between the selected IMs are proposed. This is of special interest as several correlation models between different IMs have been calibrated and validated based on advanced GMMs and global datasets, lacking earthquakes in extensional regions; however, modeling the correlation between different IM types has not been adequately addressed by current, state-of-the-art GMMs and recent ground-motion records for Italy.  相似文献   

9.
Transport processes that lead to exchange of mass between surface water and groundwater play a significant role for the ecological functioning of aquatic systems, for hydrological processes and for biogeochemical transformations. In this study, we present a novel integral modeling approach for flow and transport at the sediment–water interface. The model allows us to simultaneously simulate turbulent surface and subsurface flow and transport with the same conceptual approach. For this purpose, a conservative transport equation was implemented to an existing approach that uses an extended version of the Navier–Stokes equations. Based on previous flume studies which investigated the spreading of a dye tracer under neutral, losing and gaining flow conditions the new solver is validated. Tracer distributions of the experiments are in close agreement with the simulations. The simulated flow paths are significantly affected by in- and outflowing groundwater flow. The highest velocities within the sediment are found for losing condition, which leads to shorter residence times compared to neutral and gaining conditions. The largest extent of the hyporheic exchange flow is observed under neutral condition. The new solver can be used for further examinations of cases that are not suitable for the conventional coupled models, for example, if Reynolds numbers are larger than 10. Moreover, results gained with the integral solver provide high-resolution information on pressure and velocity distributions at the rippled streambed, which can be used to improve flow predictions. This includes the extent of hyporheic exchange under varying ambient groundwater flow conditions.  相似文献   

10.
A novel, non-intrusive fluorescence imaging technique has been used to quantitatively measure the pore geometry, fluid velocity, and solute concentration within a saturated, three-dimensional porous medium. Discrete numerical averages of these quantities have been made over a representative volume of the medium and used to estimate macroscopic quantities that appear in conventional continuum models of flow and transport. The approach is meant to illustrate how microscopic information can be measured, averaged, and used to characterize medium-scale processes that are typically approximated constitutively. The experimental system consisted of a clear, cylindrical column packed with clear spherical beads and a refractive index-matched fluid seeded with fluorescent tracer particles and solute dye. By illuminating the fluid within the column with a scanning planar laser beam, details of flow and concentration within the pore spaces can be quantitatively observed, allowing for three-dimensional, dimensional, time dependent information to be obtained at good resolution. In time dependent information to be obtained at good resolution. In the current experiment, volumetrically averaged velocities and void-to-volume ratios are first compared with bulk measurements of fluid flux and medium porosity. Microscopic measurements of concentration are then used to construct cross-sectionally averaged profiles, mean breakthrough curves, and direct measurements of the dispersive flux, velocity variance, and concentration variance. In turn, the dispersive flux measurements are compared with mean concentration gradients to provide a basis for confirming the Fickian dispersion model and estimating dispersion coefficients for the medium. Coefficients determined in this manner are compared with others based upon traditional length-scale arguments, mean breakthrough analyses, and curve fits with numerical simulations.  相似文献   

11.
This paper reports the application of a two‐dimensional hydraulic model to a braided reach of the Avoca River, New Zealand. Field measurements of water surface elevation, depth and velocity obtained at low flow were used to validate the model and to optimize the parameterization of bed friction. The main systematic trends in the measured flow variables are reproduced by the model. However, field data are characterized by greater spatial variability than model output reflecting differences in the scale of processes measured in the field and represented by the model. Additional model runs were conducted to simulate flow patterns within the study reach at five higher discharges. The purpose of these simulations was to evaluate the potential for using two‐dimensional hydraulic models to quantify the reach‐scale hydraulic characteristics of braided rivers and their dependence on discharge. Changes in flow depth and velocity with increasing discharge exhibit trends that are consistent with the results of previous field investigations, although the tendency for the wetted area of the braidplain within particular depth and velocity categories to remain fixed as discharge rises, as has been noted for several braided rivers in New Zealand, was not observed. Modelled shear stress frequency distributions fit gamma functions that incorporate a distribution shape parameter, the value of which follows clear systematic trends with rising discharge. These results illustrate both the problems of, and potential for, using two‐dimensional hydraulic models in braided river applications. This leads to something of a paradox in that while such models provide a means of generating hydraulic information that would be difficult to obtain in the field at an equivalent spatial resolution, they are, due to the problems inherent to data collection, difficult to validate conclusively. Despite this limitation, the application of spatially distributed models to investigate relationships between discharge and reach‐scale form and process variables appears to have considerable potential. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

12.
This paper, based on a real world case study (Limmat aquifer, Switzerland), compares inverse groundwater flow models calibrated with specified numbers of monitoring head locations. These models are updated in real time with the ensemble Kalman filter (EnKF) and the prediction improvement is assessed in relation to the amount of monitoring locations used for calibration and updating. The prediction errors of the models calibrated in transient state are smaller if the amount of monitoring locations used for the calibration is larger. For highly dynamic groundwater flow systems a transient calibration is recommended as a model calibrated in steady state can lead to worse results than a noncalibrated model with a well-chosen uniform conductivity. The model predictions can be improved further with the assimilation of new measurement data from on-line sensors with the EnKF. Within all the studied models the reduction of 1-day hydraulic head prediction error (in terms of mean absolute error [MAE]) with EnKF lies between 31% (assimilation of head data from 5 locations) and 72% (assimilation of head data from 85 locations). The largest prediction improvements are expected for models that were calibrated with only a limited amount of historical information. It is worthwhile to update the model even with few monitoring locations as it seems that the error reduction with EnKF decreases exponentially with the amount of monitoring locations used. These results prove the feasibility of data assimilation with EnKF also for a real world case and show that improved predictions of groundwater levels can be obtained.  相似文献   

13.
A new parameter parsimonious rainfall–run‐off model, the Distance Distribution Dynamics (DDD) model, is used to simulate hydrological time series at ungauged sites in the Lygne basin in Norway. The model parameters were estimated as functions of catchment characteristics determined by geographical information system. The multiple regression equations relating catchment characteristics and model parameters were trained from 84 calibrated catchments located all over Norway, and all model parameters showed significant correlations with catchment characteristics. The significant correlation coefficients (with p‐value < 0.05) ranged from 0.22 to 0.55. The suitability of DDD for predictions in ungauged basins was tested for 17 catchments not used to estimate the multiple regression equations. For ten of the 17 catchments, deviations in Nash–Sutcliffe efficiency (NSE) criteria between the calibrated and regionalised model were less than 0.1, and for two calibrated catchments within the Lygne basin, the deviations were less than 0.08. The median NSE for the regionalized DDD for the 17 catchments for two time series was 0.66 and 0.72. Deviations in NSE between calibrated and regionalised models are well explained by the deviations between calibrated and regressed parameters describing spatial snow distribution and snowmelt respectively. The quality of the simulated run‐off series for the ungauged sites in the Lygne basin was assessed by comparing flow indices describing high, medium and low flow estimated from observed run‐off at the 17 catchments and for the simulated run‐off series. The indices estimated for the simulated series were generally well within the ranges defined by the 17 observed series. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

14.
We investigated the role of increasingly well‐constrained geologic structures in the subsurface (i.e., subsurface architecture) in predicting streambed flux and hyporheic residence time distribution (RTD) for a headwater stream. Five subsurface realizations with increasingly resolved lithological boundaries were simulated in which model geometries were based on increasing information about flow and transport using soil and geologic maps, surface observations, probing to depth to refusal, seismic refraction, electrical resistivity (ER) imaging of subsurface architecture, and time‐lapse ER imaging during a solute tracer study. Particle tracking was used to generate RTDs for each model run. We demonstrate how improved characterization of complex lithological boundaries and calibration of porosity and hydraulic conductivity affect model prediction of hyporheic flow and transport. Models using hydraulic conductivity calibrated using transient ER data yield estimates of streambed flux that are three orders of magnitude larger than uncalibrated models using estimated values for hydraulic conductivity based on values published for nearby hillslopes (10?4 vs. 10?7 m2/s, respectively). Median residence times for uncalibrated and calibrated models are 103 and 100 h, respectively. Increasingly well‐resolved subsurface architectures yield wider hyporheic RTDs, indicative of more complex hyporheic flowpath networks and potentially important to biogeochemical cycling. The use of ER imaging to monitor solute tracers informs subsurface structure not apparent from other techniques, and helps to define transport properties of the subsurface (i.e., hydraulic conductivity). Results of this study demonstrate the value of geophysical measurements to more realistically simulate flow and transport along hyporheic flowpaths.  相似文献   

15.
Many research tools for lahar hazard assessment have proved wholly unsuitable for practical application to an active volcanic system where field measurements are challenging to obtain. Two simple routing models, with minimal data demands and implemented in a geographical information system (GIS), were applied to dilute lahars originating from Soufrière Hills Volcano, Montserrat. Single-direction flow routing by path of steepest descent, commonly used for simulating normal stream-flow, was tested against LAHARZ, an established lahar model calibrated for debris flows, for ability to replicate the main flow routes. Comparing the ways in which these models capture observed changes, and how the different modelled paths deviate can also provide an indication of where dilute lahars, do not follow behaviour expected from single-phase flow models. Data were collected over two field seasons and provide (1) an overview of gross morphological change after one rainy season, (2) details of dominant channels at the time of measurement, and (3) order of magnitude estimates of individual flow volumes. Modelling results suggested both GIS-based predictive tools had associated benefits. Dominant flow routes observed in the field were generally well-predicted using the hydrological approach with a consideration of elevation error, while LAHARZ was comparatively more successful at mapping lahar dispersion and was better suited to long-term hazard assessment. This research suggests that end-member models can have utility for first-order dilute lahar hazard mapping.  相似文献   

16.
We investigate the importance of selecting two different methodologies for the determination of hydraulic conductivity from available grain-size distributions on the stochastic modeling of the depth-averaged breakthrough curve observed during a forced-gradient tracer test experiment. The latter was performed in the Lauswiesen alluvial aquifer, located near the city of Tübingen, Germany, by injecting NaBr into a well at a distance of about 50 m from a pumping well. We also examine the joint effect of the choice of the transport model adopted to describe solute transport at the site and the way the spatial distribution of porosity is assessed. In the absence of direct measurements of porosity, we consider: (a) the model used by Riva et al. (J Contam Hydrol 88:92–118, 2006; J Contam Hydrol 101:1–13, 2008), which relates the natural logarithms of effective porosity and conductivity through an empirical, experimentally-based, linear relationship derived for a nearby experimental site; and (b) a model based on a commonly used relationship linking the total porosity to the coefficient of uniformity of grain size distributions. Transport is described in terms of a purely advective process and/or by including mass exchange processes between mobile and immobile regions. Modeling of flow and transport is performed within a Monte Carlo framework, upon conceptualizing the aquifer as a random composite medium. Our results indicate that the model adopted to describe the correlation between conductivity and porosity and the way grain-sieve information are incorporated to depict the heterogeneous distribution of hydraulic conductivity can have relevant effects in the interpretation of the data at the site. All the conceptual models employed to describe the structural heterogeneity of the system and transport features can reasonably reproduce the global characteristics of the experimental depth-averaged breakthrough curve. Specific details, such as the peak concentration and the time of first arrival, can be better reproduced by a double porosity transport model when a correlation between conductivity and porosity based on grain size information at the site is considered. The best prediction of the late-time behavior of the measured breakthrough curves, in terms of the observed heavy tailing, is offered by directly linking porosity distribution to the spatial variability of particle size information.  相似文献   

17.
1 INTRODUCTION AND REVIEW Numerous investigations related to vertical velocity and suspension concentration distributions have been undertaken to the steady, uniform, open channel turbulent flow. The problems of sediment mixed water flows are of direct interest to the practical situations in the field of river sedimentation, coastalsediment transport and in the field of two-phase flow in particular. In hydraulic open channel flow, the vertical velocity profile is usually described by…  相似文献   

18.
Transport and losses of nitrate from sloped soils are closely linked to nitrogen fertilizer management. Previous studies have always focused on different types of fertilizer applications and rarely analysed various initial nitrate distributions as a result of nitrogen fertilizer applications. Under certain conditions, both subsurface lateral saturated flow and vertical leaching dominate nitrate losses. Soil tank experiments and HYDRUS‐2D modelling were used to better understand the subsurface nitrate transport and losses through lateral saturated flow and vertical leaching under various initial nitrate distributions. Low (L: 180 mg L?1), normal (N: 350 mg L?1), and high (H: 500 mg L?1) nitrate concentrations were used in five different distributions (NNNN, NLLN, LHHL, LNLN, and HNHN) along the slope of the tank. The first two treatments (NNNN and NLLN) were analysed both experimentally and numerically. Experiments were conducted under 12 rainfall events at intervals of 3 days. The HYDRUS‐2D model was calibrated and validated against the experimental data and demonstrated good model performance. The other three treatments (LHHL, LNLN, and HNHN) were investigated using the calibrated model. Nitrate concentrations in purple sloped soils declined exponentially with time under intermittent rainfalls, predominantly in the upper soil layers. Non‐uniform initial nitrate distributions contributed to larger differences between four locations along the slope in deeper soil layers. The non‐uniform nitrate distribution either enhanced or reduced decreases in nitrate concentrations in areas with higher or lower initial nitrate concentrations, respectively. Higher nitrate concentrations at the slope foot and along the slope were reduced mainly by lateral flow and vertical leaching, respectively. Increasing nitrogen application rates increased subsurface nitrate losses. Mean subsurface lateral nitrate fluxes were twice as large as mean vertical leaching nitrate fluxes. However, due to longer leaching durations, total nitrate losses due to vertical leaching were comparable with those due to lateral flow, which indicated comparable environmental risks to surface waters and groundwater.  相似文献   

19.
The occurrence of preferential flow in the subsurface has often been shown in field experiments. However, preferential flow is rarely included in models simulating the hydrological response at the catchment scale. If it is considered, preferential flow parameters are typically determined at the plot scale and then transferred to larger-scale simulations. Here, we successfully used the optimization algorithm DiffeRential Evolution Adaptive Metropolis (DREAM) to calibrate a 3D physics-based dual-permeability model directly at the catchment scale. In order to keep computational costs of the optimization routine at a reasonable level, we limited the number of parameters to be calibrated to the ones that had been shown before to be most influential for the simulation of discharge. We also calibrated parameters of the matrix domain and the macropore domain with a fixed parameter ratio between soil layers instead of calibrating every layer separately. These ratios reflected observed depth profiles of soil hydraulic properties at our study site. The dual-permeability parameter sets identified during calibration were able to simulate observed discharge time series satisfactorily but did not outperform a calibrated single-domain reference model scenario. Saturated hydraulic conductivities of the macropore domain were calibrated such that they became very similar to matrix saturated hydraulic conductivities, thereby effectively removing the effect of macropores. This suggests that the incorporation of vertical preferential flow as represented by the dual-permeability approach was not relevant for reproducing the hydrometric response reasonably well in the studied catchment. We also tested the scale-invariance of the calibrated dual-permeability parameter sets by using the parameter sets performing best at catchment scale to simulate plot-scale bromide depth profiles obtained from tracer irrigation experiments. This parameter transfer proved to be not successful, indicating that soil hydraulic parameters are scale-variant, independent of the direction of parameter transfer.  相似文献   

20.
电阻率和速度随机分布的MT与地震联合反演   总被引:10,自引:5,他引:5       下载免费PDF全文
在已有研究成果的基础上,为了适应物性参数剧烈变化的复杂模型并满足联合反演的要求,开发了速度和电阻率随机分布共网格单元模型的建模技术.基于这种统一的物性随机分布的网格介质模型,利用有限元方法和改进的射线追踪法分别正演计算大地电磁场和地震走时,结合改进的模拟退火算法,研究实现了电阻率和速度随机分布条件下的大地电磁与地震资料的同步联合反演.对物性界面不完全一致和物性变化剧烈的带地形复杂模型的试验,表明了该方法在精细反演复杂电阻率和速度结构方面的效果,克服了以往研究局限于简单模型的不足.对地震资料品质差的地区开展的实际资料联合反演,表明了方法的适用性,先验信息约束下的联合反演提高了反演精度.  相似文献   

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

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