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。不同频率的亮温-厚度变化曲线没有交叉点,且频率越高,所能探测的月壤厚度越小。根据模拟结果,建立了月壤厚度与亮温的查找表。基于查找表,利用单个波段的亮温数据即可得到月壤厚度信息。     

贵州糯东井田水文地质条件分析

通过对井田内泉和暗河流量的调查研究,认为永宁镇组一、三段碳酸盐岩裂隙溶洞水为井田的主要含水层,含水层主要接受大气降水及河流的渗透补给;煤层上覆地层飞仙关组为相对隔水层,厚度大于500m,下覆无强含水层,其地下水对煤层的开采影响不大,井田水文地质条件较简单。采用比拟法对井田的涌水量进行了预算,推算该矿井开采设计涌水量为13 238m^3/d。     

塔河油田2区三叠系油藏类型初步探讨

塔河油田2区三叠系油藏是西北分公司开发较早的碎屑岩油藏之一,开发层系有上油组、中油组,油柱高度在15～20m,储层为辫状三角洲相辫状三角洲前缘亚相水下分流河道、河口坝、分流间湾微相沉积,非均质性严重。本区上油组砂体厚度40～0m,中油组25～35m,分布较稳定,物性为中孔、中高渗储层,原油属低粘度轻质油。但2区三叠系油藏实钻变化大,见水快、产量递减快、含水上升快、采出程度低,开发效果较差,经过多次调整虽有一定程度改善,仍存在一些问题。本文认为2区三叠系油藏开发效果较差的原因之一是对本区油藏类型的认识上存在问题即2区三叠系油藏是以断裂控制的断块油藏而非背斜油藏,并延伸分析正在投入开发评价的S96-T901井区、TK232-T453井区及T206-T208井区等都为断块（断背斜）油藏。应用断块油藏的观点对塔河地区三叠系油藏勘探开发有重要指导意义。