多孔紊动射流的数值模拟与实验研究进展

多孔射流作为流体运动的一种重要类型,在工程中具有非常广泛的用途。由于多孔射流各股水流之间的卷吸和掺混的存在,增加了流场的复杂性。长期以来研究人员从理论分析、实验测量和数值模拟方面对多孔紊动射流进行了大量的研究工作,目前已经认识了流场中的许多流动特性和流动机理。从数值模拟和实验研究两个方面,比较分析了国内外多孔紊动射流的研究现状和研究结果,评述了不同湍流模型以及不同实验测量方法对紊动射流的预测能力,讨论了存在的问题并对该领域的研究方向进行了展望。     

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

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

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.     

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.     

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

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

顾及DEM误差自相关的坡度计算模型精度分析

基于DEM的坡度计算,其误差来源于DEM误差、DEM结构和坡度计算模型。在顾及DEM误差自相关的前提下,对四种DEM坡度计算模型进行了分析和评价。研究表明,三阶不带权差分能给出较高的坡度计算精度;在局部窗口中,格网点数量越多,坡度计算越准确;等权比不等权的坡度计算模型更准确;DEM误差自相关结构形式对坡度计算无影响。进一步的理论分析和试验分析还表明：DEM误差自相关性的存在,不仅能够改善地形分析的精度,也能改善DEM自身精度。     

人工边坡潜在滑动面研究——以广州科学城某人工高边坡为例

运用强度参数的改变对边坡破坏面形迹影响不明显这一特点,在数值模拟过程中通过改变岩体强度参数,有效地获取潜在滑动面的位置和形态,较好地解决了滑动面搜索的难题。将该法应用于广州科学城某人工高边坡稳定性的研究,在三维数值模拟过程中,将强度参数大幅度折减,计算后获得各剖面的剪应变增量图,从这些图中可获得潜在滑动面。这与人们通常将此类边坡的中风化面作为滑动面存在较大差别。将该滑动面运用极限平衡法进行计算,计算结果显示各剖面的安全系数基本都大于1.2,边坡稳定但仍需要加固处理,与三维数值模拟结果相一致。由此认为用这种分析法确定出的潜在滑动面合理、计算结果可靠,可作为搜索边坡潜在滑动面并计算安全系数的方法之一。     

非均匀月壤介质的被动微波辐射传输模拟

基于非均匀月壤物理模型和辐射传输方程,模拟月壤介质中的微波辐射传输特性,探讨频率、月壤厚度等与月表亮温的关系。结果表明:在低频段,月壤微波辐射亮温的动态变化范围较大,可探测的月壤厚度大,3 GHz时的最大可探测月壤厚度达12.4 m;在高频段对应的可探测月壤厚度较小,特别是从50GHz往后的频率段内,最大可探测月壤厚度均小于2 m。不同频率的亮温-厚度变化曲线没有交叉点,且频率越高,所能探测的月壤厚度越小。根据模拟结果,建立了月壤厚度与亮温的查找表。基于查找表,利用单个波段的亮温数据即可得到月壤厚度信息。