一种土的非线性弹性本构模型参数的反演方法

旨在提出一种土的非线性弹性本构模型参数反演的方法。以现今普遍实行的地基载荷试验为基础,依据遗传算法的组合优化理论,采用正演计算和遗传算法优化相结合的方式,建立了土层非线性弹性本构模型参数反演的方法;并依据某黄土场地地基载荷试验数据,实施了黄土土层非线性弹性本构模型参数反演的全过程。计算结果表明,所建立的方法可以实现土层非线性弹性本构模型中相互关联的多个参数的组合优化,并在对初始值要求较低的情况下,可以获得良好的参数反演精度。从而为土的变形特性分析和土与其中及相邻结构的共同作用分析,提供了较好的土体本构模型参数的确定方法。     

Development and parameterisation of a complex hydrogeological model based on high-resolution direct-push data

Ongoing developments in geological and hydrogeological investigation techniques, especially direct-push methods, have led to an increase in the quality, density and spatial resolution of data available from such investigations. This has created new challenges in the development of numerical models in terms of accurately and efficiently translating detailed and complex conceptual models into effective numerical models. Suitable geometrical and numerical modelling tools are essential in order to meet these challenges. This paper describes the development of a three-dimensional hydrogeological flow model for a contaminated site near Berlin, Germany, based on high-resolution geological data obtained principally using direct-push methods. The available data were first interpreted to construct a detailed GIS-based geological model, which formed the basis of the conceptual site model. The conceptual model was then translated into a geometrical model, which was used to create a finite element numerical model. An innovative geometry object-based approach enabled the complex structural details of the conceptual model to be accurately reproduced in the numerical model domain. The resulting three-dimensional steady-state unconfined flow model was successfully calibrated using external automated calibration software, whereby parameter values for groundwater recharge and hydraulic conductivity were determined.     

地震孕育过程中若干地震学特征参量研究

从新疆实际出发提出一种简捷的预报程式,用以判定地震三要素。结果表明,(1)新疆3级地震累积频次低于均值6个月,1年内新疆发生6级地震的概率为90%;(2)以新疆地震活跃期响应地震分布区做为预备震源区,该区未来3年发生中强震的概率为86%;(3)预备震源区出现一年以上低6值,其250 km区域1年内发生6级地震的概率为75%;(4)当系统处于相对稳定状态时系统内部个体的涨落很难使系统出现非线性响应,而当系统处于亚失稳状态时,系统内部个体的涨落有可能触发系统中处于高应力状态下的地震活动。     

Assessment of radionuclide transport uncertainty in the unsaturated zone of Yucca Mountain

The present study assesses the uncertainty of flow and radionuclide transport in the unsaturated zone at Yucca Mountain using a Monte Carlo method. Matrix permeability, porosity, and sorption coefficient are considered random. Different from previous studies that assume distributions of the parameters, the distributions are determined in this study by applying comprehensive transformations and rigorous statistics to on-site measurements of the parameters. The distribution of permeability is further adjusted based on model calibration results. Correlation between matrix permeability and porosity is incorporated using the Latin Hypercube Sampling method. After conducting 200 Monte Carlo simulations of three-dimensional unsaturated flow and radionuclide transport for conservative and reactive tracers, the mean, variances, and 5th, 50th, and 95th percentiles for quantities of interest (e.g., matrix liquid saturation and water potential) are evaluated. The mean and 50th percentile are used as the mean predictions, and their associated predictive uncertainties are measured by the variances and the 5th and 95th percentiles (also known as uncertainty bounds). The mean predictions of matrix liquid saturation and water potential are in reasonable agreement with corresponding measurements. The uncertainty bounds include a large portion of the measurements, suggesting that the data variability can be partially explained by parameter uncertainty. The study illustrates propagation of predictive uncertainty of percolation flux, increasing downward from repository horizon to water table. Statistics from the breakthrough curves indicate that transport of the reactive tracer is delayed significantly by the sorption process, and prediction on the reactive tracer is of greater uncertainty than on the conservative tracer because randomness in the sorption coefficient increases the prediction uncertainty. Uncertainty in radionuclide transport is related to uncertainty in the percolation flux, suggesting that reducing the former entails reduction in the latter.     

CORS主要理论方法的比较分析

主要介绍了CORS的理论基础、基本概念以及国内外发展现状,对其主要理论方法进行了深入的分析比较,为我国各地的CORS建设提供了理论参考和技术借鉴.     

An approach to handling non-Gaussianity of parameters and state variables in ensemble Kalman filtering

The ensemble Kalman filter (EnKF) is a commonly used real-time data assimilation algorithm in various disciplines. Here, the EnKF is applied, in a hydrogeological context, to condition log-conductivity realizations on log-conductivity and transient piezometric head data. In this case, the state vector is made up of log-conductivities and piezometric heads over a discretized aquifer domain, the forecast model is a groundwater flow numerical model, and the transient piezometric head data are sequentially assimilated to update the state vector. It is well known that all Kalman filters perform optimally for linear forecast models and a multiGaussian-distributed state vector. Of the different Kalman filters, the EnKF provides a robust solution to address non-linearities; however, it does not handle well non-Gaussian state-vector distributions. In the standard EnKF, as time passes and more state observations are assimilated, the distributions become closer to Gaussian, even if the initial ones are clearly non-Gaussian. A new method is proposed that transforms the original state vector into a new vector that is univariate Gaussian at all times. Back transforming the vector after the filtering ensures that the initial non-Gaussian univariate distributions of the state-vector components are preserved throughout. The proposed method is based in normal-score transforming each variable for all locations and all time steps. This new method, termed the normal-score ensemble Kalman filter (NS-EnKF), is demonstrated in a synthetic bimodal aquifer resembling a fluvial deposit, and it is compared to the standard EnKF. The proposed method performs better than the standard EnKF in all aspects analyzed (log-conductivity characterization and flow and transport predictions).     

An Analytic Solution to Well-water Level Changes under Barometric Pressure

Under barometric pressure,groundwater flow in well-aquifer systems is a kind of hydro-mechanical coupling problem.Applying the flux boundary conditions on borehole wall and water pressure equilibrium conditions inside and outside the borehole wall under barometric pressure(BP),an analytic solution to well-water level changes has been proposed in this paper.The formulation shows that the BP coefficients increase with time and tend to BP constant.The Change of BP coefficients over time depends only on the ratio of transmissivity(T) to the well radius squared(r2w),and has nothing to do with the change in BP.The BP constant only relates to aquifer loading efficiency(B),and has nothing to do with the aquifer transmissivity and well radius.The BP coefficients change over time in the analytic formulation is consistent with the analysis of measured data from the Nanxi wells.Based on the BP coefficient changes over time,a parameter estimation method is suggested and discussed in its application to the estimation of the aquifer BP constant(or B) and transmissivity by using the Nanxi well data.     

Temporal and spatial variations of runoff and rainwash erosion on an agricultural field

This paper investigates the spatial and temporal variations of runoff, erosion and rate of sediment transport on an agricultural field submitted to natural rainfalls. The site, located in the Eastern Townships (Québec, Canada), is a corn field (10000 m²) where sheetwash erosion is active. Water (Q) and sediment (Qs) discharges were measured from June to October at eight locations on the field and for ten rainfall events. Analysis of the data was carried out on an aggregate data set and on the distributed measurements in time and space. The results showed that changes in vegetation, soil compaction and crusting are critical in determining temporal variations of runoff and erosion. Until August, the increase in soil compaction reduced infiltration capacity and depression storage and generated greater runoff for a given rainfall intensity (I). Sediment transport decreased as particle detachment is less likely to occur when vegetation breaks the drop impact and the soil surface is sealed. Later in the season, we observed an increase in sediment concentration associated with the presence of burrowing insects and harvest activity, providing loose sediments to the broken down surface. Intercepts and slopes of the relationship between Q and Qs also vary during the period of measurement. High sediment availability over the soil surface in June and October gives high intercept values. The slope of the relationship is more stable but difficult to estimate for extreme events (high values of I or low Q values) where the number of sampled points are small. During a rainfall, the response of the field is dominated by the topography and drainage area. The largest amount of runoff and erosion occurred on straight and steep slopes with small drainage areas, and on converging gentle slopes with large drainage areas. Although aggregate runoff and erosion values are decreasing with drainage area, parameters of the Qs-Q relationship for different locations on the field are not statistically different. These results bear important consequences for models of sheetwash erosion on agricultural fields.     

成都地面气温与风速年极大值的渐近分布及参数估计

利用１９３３～１９９２年期间成都地面最高气温和地面最大风速年极值的记录，通过统计推断，找出了成都地面最高气温年极值和地面最大风速年极值遵循的渐近分布——Ｗｅｉｂｕｌ分布和Ｇｕｍｂｅｌ分布，并讨论了它们的参数估计方法。     

A New Inversion Method of Sedimentary Strata For Deriving Double Parameters and Its Applications

The inverse problem of wave equation is the importance of study not only in seismic prospecting but also in applied mathematics. With the development of the research, the inverse methods of 1 - D wave equations have been trending towards the multiple parameters inversion . We have obtained an inverse method with double -parameter, in which medium density and wave velocity can be derived simultaneously. In this paper, to increase the inverse accuracy, the method is improved as follows. Firstly, the formula in which the Green Function is omitted are derived and used. Secondly, the regularizing method is reasonable used by choosing the stable function. With the new method, we may derive elastic parameter and medium density or medium density and wave velocity. Thus, lithology parameters for seismic prospecting may be obtained.After comparing the derived values from the new method with that from previous method, we obtain the new method through which substantially improve the derived accuracy . The new method has been applied to real depths inversion for sedimentary strata and volcanic rock strata in Chaoyanggou Terrace of Songliao Basin in eastern China. According to the inverse results,the gas - bearing beds are determlned.