土壤吸湿-湿过程中水力传导度的拟合研究

非饱和含水介质在吸湿和脱湿的不同过程中,水力传导度随介质含水率的变化规律是不同的。本文选取大沽河流域两种代表性砂土,对其吸湿-脱湿两个过程分别进行试验观测,并使用VanGenuchten—Mualem模型对其吸湿-脱湿过程进行水力传导度的拟合。研究结果表明：砂土在吸湿过程的非饱和水力传导度与脱湿过程的非饱和水力传导度大小并不相同,在含水率相同的情况下,吸湿过程的水力传导度要比脱湿过程大,在含水率较低时（小于15％）相差可达1个数量级。     

鲁北平原区井灌回归系数分析

利用井灌回归系数计算回归补给水量,是地下水资源评价计算井灌回归水量的基本方法。井灌回归系数确定的合理与否,直接影响到地下水资源评价成果的合理性。以试验站的观测资料为基础,结合引水灌溉调查成果,对影响井灌回归系数大小的因素和区域分布规律进行了系统分析,依地下水埋深和灌溉方式划分引黄灌区、井河结合灌区、井灌区三种类型,分析确定鲁北平原区适用的井灌回归系数。     

爆破震动对某边坡稳定性影响的数值模拟

爆破震动是影响边坡稳定性的一个重要因素之一。结合实际工程文章利用了拉格朗日有限差分法FLAC程序对爆破震动作用下某边坡进行了数值模拟。讨论了利用FLAC进行边坡动力分析时设置边界条件设置、动力时程的合成、输入以及转化,阻尼的合理选取。从位移场和速度场进行了详细的动力分析。获得该边坡爆破震动特征及位移场、振动速度场的变化规律。研究成果对保证边坡安全具有较好的指导作用。     

Modelling Positional Uncertainty of Line Features by Accounting for Stochastic Deviations from Straight Line Segments

The assessment of positional uncertainty in line and area features is often based on uncertainty in the coordinates of their elementary vertices which are assumed to be connected by straight lines. Such an approach disregards uncertainty caused by sampling and approximation of a curvilinear feature by a sequence of straight line segments. In this article, a method is proposed that also allows for the latter type of uncertainty by modelling random rectangular deviations from the conventional straight line segments. Using the model on a dense network of sub‐vertices, the contribution of uncertainty due to approximation is emphasised; the sampling effect can be assessed by applying it on a small set of randomly inserted sub‐vertices. A case study demonstrates a feasible way of parameterisation based on assumptions of joint normal distributions for positional errors of the vertices and the rectangular deviations and a uniform distribution of missed sub‐vertices along line segments. Depending on the magnitudes of the different sources of uncertainty, not accounting for potential deviations from straight line segments may drastically underestimate the positional uncertainty of line features.     

软弱夹层几何参数对试样力学行为影响颗粒元模拟研究

基于离散元理论的颗粒元程序,采用颗粒接触连接模型和滑动模型,建立了含软弱夹层试样的颗粒元模型。通过数值模型实验,对不同软弱夹层几何参数试样进行加载模拟,探讨了不同夹层厚度及不同倾角条件下试样的破坏过程和破坏形式,并分析其对试样强度的影响。模拟结果表明,裂隙首先产生于软弱夹层处,随着荷载的增加,再逐渐扩展,呈现渐进性破坏特征;夹层越厚,倾角越大,试样沿软弱夹层呈现较强的滑动特征,如果夹层薄且缓,软弱夹层成为试样破坏的非控制因素;软弱夹层厚度和夹角对试样的峰值强度都有影响,随软弱夹层厚度和夹角的增大,试样的峰值强度降低。     

基于GIS和VRML技术的东圳水库3维仿真建模

以福建省莆田市东圳水库为例,讨论了基于VRML的水库及其周围地理环境、流域状况仿真建模过程,以及所采用的技术和方法。通过仿真系统可以直观分析水库的环境情况,水库管理者可在场景中漫游,从不同位置进行观察分析,进而提出相应解决措施。同时该系统作为对水库进行研究的人机交互平台,可进行实时完善和扩展。     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Numerical simulation of migration–accumulation of oil resources

Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value to the evaluation of oil resources and to the determination of the location and amount of oil deposits. This thesis discusses the characteristics of petroleum geology and permeation fluid mechanics. For the three-dimensional problems of Dongying hollow of Shengli Petroleum Oil Field, it puts forward a new model and a kind of modified method of upwind finite difference fractional steps implicit interactive scheme. For the famous hydraulic experiment of secondary migration–accumulation, the numerical simulation test has been done, and both the computational and experimental results are basically identical. For the actual problem of Dongying hollow, the numerical simulation test and the actual conditions are basically coincident. Thus, the well-known problem has been solved.     

Accurate local upscaling with variable compact multipoint transmissibility calculations

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.     

西藏札达盆地控盆断裂有限元数值模拟

札达盆地是喜马拉雅构造带中的一个山间断陷盆地,其演化过程与盆地两侧的控盆构造密切相关。对控盆断裂的构造应力场进行模拟计算,将有助于进一步深化对本区构造控盆的认识。因此,在对盆地构造地质进行详细调查的基础上,结合本区的深部地质与地球物理资料,对札达盆地控盆断裂的构造应力场进行了模拟。计算结果表明,札达盆地的演化明显受盆地两侧边界断裂的控制,札达盆地是在整体南北向挤压应力的作用下,不同块体差异隆升作用的结果。其南侧的控盆断裂为北倾的正断层,北侧的控盆断裂为南倾的逆断层,二者共同形成了南降北升的翘板式断陷盆地运动过程,是喜马拉雅地块在陆内汇聚挤压构造环境中构造应力场调整的一种方式。