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.     

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

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

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

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

有限元法模拟开挖引起的基坑水平位移和邻近建筑物沉降

以徐州某采用水泥土墙支护的基坑为例,采用有限元数值模拟的方法,对基坑开挖引起的水平位移量和邻近建筑物沉降量进行了计算,并与实际监测值进行了对比分析,得到了一些关于如何控制基坑水平位移和减少基坑周边建筑沉降的结论。     

数值模拟技术在海坨子油田开发调整研究中的应用

讨论分析数值模拟技术在海坨子油田开发调整研究中的思路和方法,采用美国Landmark公司VIP2003．4版本数值模拟软件包进行模拟研究,通过参数准备,静、动态拟合,模拟出海坨子油田海S9—11井组1—9砂组各小层剩余储量大小及不同阶段流体饱和度分布,结合可视化技术形象描述了海坨子油田剩余油变化规律及潜力区,在此基础上预测不同井网加密方式开发指标,并优选出最佳方案,为海坨子油田下步开发调整提供了重要依据,取得了较好的效果。     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.     

Effects of onshore and offshore topography on sea breeze circulation: An observational study at eastern Attica,Greece

Field measurements of wind, air temperature and humidity were taken at the eastern part of the Attika district in June 1991, to examine the topographic influences exerted on the local sea breezes. These influences are due either to the elongated Evia island, faced by the northern half of Attica coastline some tens of kilometers offshore, or to the coast-parallel range of Hymettos mountain, rising steeply 12 km onshore. The instrumentation consisted mainly of three tethered meteorological balloons released at characteristic sires (i.e., the coast, a location between shoreline and mountain foot and the mountain top) and three autographic ground-based anemometers operating at selected locations. Data from the ground-based and upper air stations of the Hellenic National Meteorological Service, as well as the diurnal weather maps were also obtained and analyzed. Observations were made under different synoptic wind and the latter was found to determine remarkably the significance of the topographic effects. A preliminary two-dimensional numerical approach was also made concerning the sea breeze capability to reach the Hymettos mountain top in the case of a weak opposing geostrophic flow.     

A stability study of asteroid families near the 3:1 and 5:2 resonance with Jupiter

By using theD-criterion Lindblad (1992) has identified 14 asteroid families from a sample of 4100 numbered asteroids with proper elements from Milani and Kneevi (1990). Taxonomic types and other physical properties for a significant number of objects in five of the families show strong homogeneity within each family, further strengthening their internal relationship.To test the hypothesis of a common origin in, e.g., a catastrophic collision event, we have set out to integrate the orbits of the members of the Maria, Dora and Oppavia-Gefion families over some 10⁶ years. The mean distance for the Maria family is close to the 3:1 mean-motion resonance with Jupiter, while the other two families lie close to the 5:2 resonance.We used a simplified solar system model which included the perturbations by Jupiter and Saturn only and implemented Everhart's variable stepsize integrator RA15. All close encounters between the family members (within 0.1 AU) were recorded as well. Preliminary results from integrations over 4×10⁵ years are presented here.The statistics of close encounters show pronounced peaks for several members within each family, while for others no significant levels above the background of random encounters or even very low frequencies were found. This indicates a subclustering within the families. Quite a lot of very close (<0.005 AU) mutual encounters are found, which suggest that, at least for the larger members in a family, the mutual gravitational interactions could be of some importance for the real orbital evolutions.The encounter statistics between the Dora and Oppavia family members suggest a possible interrelationship between this two groups.     

新疆地区夏季平流层低层风场特征分析及数值模拟

利用2017年6—8月的FNL再分析资料分析了新疆地区夏季平流层低层风场的时空演变特征，开展了平流层低层风场的高分辨率数值模拟和检验。结果表明：新疆地区夏季平流层纬向风随时间的变化存在经向差异，同时准零风层开始和结束的时间也存在经向差异，准零风层开始时间南部地区早于北部地区，而准零风层结束时间则相反。整个夏季新疆地区上空的准零风层处于70～40 hPa之间，其高度随时间呈先降低而后升高的变化趋势。准零风层数值模拟结果表明，模拟的准零风层参数与探空资料分析结果相比二者存在一致的变化趋势，准零风层起始高度的平均绝对误差为467 m，该高度对应的风速均方根误差为1.75 m/s。