非线性波消波及波场分布研究   总被引：1，自引：1，他引：0  

王巧莎  李明海 《海洋科学》2018,42(6):30-38

本文基于雷诺平均N-S方程,并结合RNG k-e方程建立了粘性数值波浪水槽,对不同波陡、不同相对水深、不同相对波高的非线性规则波的阻尼消波问题和波场分布进行研究。文中提出了两种描述消波区内部阻尼变化的阻尼函数,分别适用于小波陡情形和高波陡情形。研究结果表明,小波陡组消波区可设为一个波长,阻尼系数取10~4～10~5即可满足消波要求,计算结果与实验结果及造波理论吻合良好;高波陡组消波区可设为两个波长,阻尼系数取10~4～10~5亦可满足消波要求,计算结果与实验结果吻合良好。此外,当波陡较小时,波场内反射情况的小幅改变即可对整个波场造成影响,特别是当水深较浅时这种影响极为明显,需谨慎考虑。当波陡较大时,水波能量较高,整个波场沿水波传播方向可观测到明显的衰减现象,在具体试验中需进行考虑。  相似文献   

Three dimensional simulation of fluid flow in X‐ray CT images of porous media     

A. Al‐Omari  E. Masad 《国际地质力学数值与分析法杂志》2004,28(13):1327-1360

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

Numerical manifold method for dynamic nonlinear analysis of saturated porous media     

H. W. Zhang  L. Zhou 《国际地质力学数值与分析法杂志》2006,30(9):927-951

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

About Darcy's law in non-Galilean frame     

C. Geindreau  E. Sawicki  J.-L. Auriault  P. Royer 《国际地质力学数值与分析法杂志》2004,28(3):229-249

This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek⁻¹ is 𝒪(1). These variations are estimated to be less than 5% if Ek⁻¹<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

Sedimentary characteristics of tide-dominated estuary in Donghetang Formation（Upper Devonian）, central Tarim Basin     

WANG Li-ping LI Chuang TANG Bing-yong 《东北亚地学研究》2005,8(1):78-93

The Donghetang Formation （Upper Devonian） in central Tarim Basin has been thought an important oil and gas reservoir since the abundant oil and gas resources were found in the wells W16, W20, W34, and other fields. However, the sedimentary environment of the Donghetang Formation has been disputed because it suffered from both tidal and fluvial actions and there were not rich fossils in the sandstone. After the authors analyzed sedimentary features by means of drill cores, well logging data, paleosols, and with SEM obseruations, three kinds of sedimentary environments were distinguished： alluvial fan, tide-dominated estuary, and shelf. Particularly, the sedimentary features of tide-dominated estuary were studied in detail. Besides, the authors discussed sedimentary characteristics of the Donghetang Formation which was divided into two fourth-order sequences and five system tracts. At the same time, according to the forming process of five system tracts, the whole vertical evolution and lateral transition of tide-dominated estuary were illustrated clearly. Finally, the reservoir quality was evaluated based on porosity and permeability.  相似文献   

Upscaling of elastic properties for large scale geomechanical simulations     

F. Chalon  M. Mainguy  P. Longuemare  P. Lemonnier 《国际地质力学数值与分析法杂志》2004,28(11):1105-1119

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

鄂尔多斯盆地西峰地区长8砂岩微观非均质性的实验分析     

刘林玉  王震亮  柳益群 《岩矿测试》2008,27(1):29-32

通过真实砂岩微观模型实验,客观地分析了鄂尔多斯盆地西峰地区长8砂岩的微观非均质性特征,发现该区长8储层为致密型砂岩,直观地显示出长8砂岩储层具有很强的微观非均质性,砂岩的成岩作用和孔隙结构是影响鄂尔多斯盆地西峰地区长8砂岩储层微观非均质性的主要原因。在低渗透砂岩中,引起砂岩产生强烈的非均质性,而这些砂岩则形成长8砂岩储层的高渗带。  相似文献   

鄂尔多斯盆地东缘煤储层微孔隙结构特征及其影响因素   总被引：2，自引：0，他引：2  

张松航  汤达祯  唐书恒  许浩  林文姬  张彪 《地质学报》2008,82(10):1341-1349

在对鄂尔多斯盆地东缘58套采自不同层位不同地域煤样显微组分测试和低温液氮比表面测试的基础上,从比表面积、吸附曲线形态、阶段孔径结构分布、阶段孔径比表面积贡献、FHH分形特征几个方面刻画了煤储层微孔隙特征。同时,就“孔比表面积—孔隙结构—煤物质组成与化学结构”这一因果链,分析了鄂尔多斯盆地东缘煤微孔隙特征与煤变质程度、煤岩显微组分以及矿物组分与类型的关系。研究发现,Ro在0.60%~1.91%范围内,比表面积及FHH分形分维数先明显降低,Ro至1.3%后又略有回升,同时吸附回线“滞后环”及分形尺度对应的压力范围明显减小,孔径3~4 nm“墨水瓶”型孔明显减少乃至消失,并成为比表面积减少的主要贡献者。在相似煤级条件下,惰质组分较镜质组分有更多的比表面积贡献。粘土矿物组分含量>10%时,孔径3~4 nm“墨水瓶”型孔对比表面贡献明显,且随煤级增高影响作用相对加大。  相似文献   

高、低煤阶煤层气藏地质特征及控气作用差异性研究   总被引：4，自引：0，他引：4  

王勃  巢海燕  郑贵强  李景明  郭志斌  王红岩  刘洪林  李贵中 《地质学报》2008,82(10):1396-1401

高、低煤阶煤层气地质特征及控气作用差异性是研究煤层气富集成藏的重要组成部分,是煤层气勘探开发理论研究过程中重要的基础性研究领域之一。本文以中国沁水、阜新盆地和美国粉河盆地等典型的含气盆地为例,探讨了高、低煤阶煤层气的储层物性差异,分析了构造控气和水文地质控气作用的差异性。研究表明,高煤阶气藏含气量高,CH4百分含量高,δ13C1值大于-38.75‰,储层渗透率变化小,储层改造难,构造热事件对煤层气的生成、富集贡献大,持续的水动力使气藏遭到破坏,且破坏幅度大,现今地下水格局对气藏的形成具有一定的影响;低煤阶气藏含气量低,CH4百分含量低,δ13C1值大于-49.11‰,储层渗透率变化大,储层易改造,煤热演化史及煤阶影响着煤层气的生成、富集,在煤层气生成过程中活跃的水动力是甲烷生成的主要的水文地质条件之一,但持续的水动力使气藏遭到破坏,且破坏幅度小,而合适的地层水矿化度则是低煤阶煤层气生成的重要条件,地下水格局对气藏的调整和改造起到决定性的影响。  相似文献   

松辽盆地尚家—太平川地区石油地质条件及勘探方向分析   总被引：17，自引：9，他引：8       下载免费PDF全文

施立志  林铁锋  吴海波  王卓卓  祁越 《地质科学》2008,43(3):558-568

通过对尚家—太平川地区构造演化、烃源岩特征、沉积特征、储盖组合、油水关系及油藏类型和油气成藏特征的分析,研究了本区的石油地质条件,认为烃源岩具有较好的生、排烃能力,分流河道砂体为油气聚集提供较好的沉积环境,良好的储盖组合有利于油气的聚集与保存。研究区没有统一油水界面,尚家地区以岩性—构造油藏为主,太平川地区以构造—岩性油藏为主。区内扶余—杨大城子油层油气成藏受有效性圈闭、优势运移通道、过剩压力、储集条件和构造诱导作用控制。据上述认识确定了研究区下一步的勘探方向,并已在指导勘探实践中取得了较好的效果。  相似文献