裂缝发育程度对低孔隙度岩石渗流特性的影响

低孔隙度岩石中的裂缝对储层渗透率具有重要影响，但裂缝的存在导致岩心代表性样品选取和高精度岩石物理参数测量困难。为研究裂缝对低孔隙度岩石渗透率的影响，本文通过高精度CT扫描实验构建了低孔隙度岩石的三维数字岩心模型，采用添加平板裂缝的方法构建了不同裂缝参数的低孔隙度岩石数字岩心，并利用格子玻尔兹曼（LBM）方法计算了不同裂缝参数数字岩心模型的渗流场分布和绝对渗透率。结果表明，尽管低孔隙度岩石的数字岩心模型基质渗透率低，但裂缝的存在对岩石渗透率有一定程度的提高。然而，裂缝发育程度对渗透率影响规律不同：当单条裂缝孔隙度在0-0.4 %时，裂缝对模型渗透率影响不明显；当单条裂缝孔隙度大于0.4 %时，裂缝对模型渗透率具有显著影响；模型渗透率随裂缝开度增大而增大，随裂缝倾角增大而减小，随裂缝数量增加而增大。另外，裂缝与基质存在耦合作用，与裂缝相连的孔隙中流体流速明显提高，显示裂缝对基质孔隙的强连通作用。本研究结果对含裂缝的低孔隙度储层渗透率精确计算及储层压裂后的油气产能评价有指导意义。     

低渗透油藏动态渗吸机理实验研究及数字岩心模拟

在低渗油藏的开发过程中,渗吸作用对于产量提升起到重要作用.传统的岩心渗流实验难以从微观的角度较为清晰地描述岩石物性、孔喉分布、岩石组分构成等特征.利用数字岩心技术建立孔隙网络模型,对真实岩心孔隙结构进行渗吸模拟分析具有重要意义.对松辽盆地某区域岩心样品通过综合微计算机断层扫描技术(micro computed tomography,Micro-CT)、数字岩心分析技术,实现低渗储层数字岩心渗吸过程的模拟,与实验室岩心自发渗吸实验对照,结果表明,利用数字岩心技术模拟得到的孔隙度、渗透率参数与实验室测量结果相差较小,提取的孔隙网络模型与真实岩心误差较小,从而可以保证数字岩心得到的两相渗流模拟结果具有较高的可信度,为油藏渗吸作用的机理研究提供了新的思路,对指导实际油藏开发具有重要意义.     

应用核磁T2谱与数字岩心技术计算粒度分布方法

提出一种利用核磁T2谱结合数字岩心技术计算岩石粒度分布的方法。岩石粒度分布近似符合Weibull分布,改变分布参数可得到不同粒度分布。给出初始Weibull分布后采用过程法构建数字岩心,并用随机行走法对构建的数字岩心进行核磁数值模拟,遍历分布参数区间使模拟得到T2谱与试验T2谱逼近,误差最小时获得的Weibull分布即为岩石粒度分布。最后根据实际核磁测井资料对井下岩石粒度进行连续计算,粒度计算结果与实验室粒度分析结果对比良好,验证了方法准确性。     

基于砂岩数字岩心的地层参数模拟与粗化方法研究

基于micro-CT成像的数字岩心技术的发展,拓展了岩石物性参数的测量和分析方法。基于三维数字岩心的物性参数模拟,可较为准确的反应岩石的物性,针对砂岩三维数字岩心图像的地层参数模拟与粗化,对研究岩心的流体分布等具有十分重要的意义。首先,基于三维数字岩心图像,采用闵可夫斯基泛函和K-means聚类法对岩心进行精细分层;其次,采用驱替模拟法,模拟岩心在不同驱替状态下湿相和非湿相的分布,针对不同的湿相分布状态,对各个小层在不同方向上的地层参数进行模拟;最后,采用数学平均法、曲线回归法和去除薄层后的曲线回归法对岩心的地层参数进行粗化并对比分析各种粗化方式结果的差异。     

砂岩数字岩心的地层参数模拟与粗化方法研究

基于micro-CT成像的数字岩心技术的发展,拓展了岩石物性参数的测量和分析方法。基于三维数字岩心的物性参数模拟,可较为准确地反应岩石的物性,针对砂岩三维数字岩心图像的地层参数模拟与粗化,对研究岩心的流体分布等具有十分重要的意义。首先,基于三维数字岩心图像,采用闵可夫斯基泛函和K-means聚类法对岩心进行精细分层;其次,采用驱替模拟法,模拟岩心在不同驱替状态下湿相和非湿相的分布,针对不同的湿相分布状态,对各个小层在不同方向上的地层参数进行模拟;最后,采用数学平均法、曲线回归法和去除薄层后的曲线回归法对岩心的地层参数进行粗化并对比分析各种粗化方式结果的差异。     

加压时间对储层岩心渗透率的影响

目前储层应力敏感性研究仅限于有效覆压对渗透率的影响,忽略了加压时间对岩石变形的作用.在高压孔渗仪上用拟三轴岩心夹持器和天然岩心,考察了加压时间对岩石渗透率的影响,提出并计算了单位覆压下渗透率绝对变化速率和相对变化速率.结果表明,岩石变形或渗透率变化与时间有关,渗透率随加压时间的延长而降低.在加压初期渗透率变化幅度较大,随着时间的延长,渗透率变化逐渐变缓;幂函数在总体趋势上能拟合渗透率随时间的变化,其拟合误差在工程应用的合理范围内;基础渗透率大于4.6×10-3 μm2的岩心,渗透率下降主要发生在4 h内;渗透率绝对变化速率随时间的延长而降低;渗透率绝对变化速率和相对变化速率与岩石物性有关,岩心基础渗透率越高,渗透率绝对变化速率越大,相对变化速率越小.渗透率相对变化速率的建立,为储层应力敏感性评价提供了一个新的定量参数.     

基于Kozeny-Carman方程的水相相对渗透率计算方法

岩石孔隙中流体的相对渗透率是储层评价、油气田开发的关键参数之一,与流体性质、岩石的润湿性、孔隙结构等有关。一般通过稳态法与非稳态法实验获得。由于相渗实验的周期长、成本高,且条件较为苛刻,难以在现场中进行大量的应用。以Kozeny–Carman方程和阿尔奇公式为基础,将复杂的孔隙结构及流体分布进行简化,得到了毛管渗流中外油内水和外水内油状态下水相的相对渗透率模型。并通过G油田4块岩心的实验数据对模型进行了分析。研究表明:外油内水状态下的水相相对渗透率模型与孔隙中油水渗流规律相符,用该模型计算的水相相对渗透率与实测数据较为吻合。     

考虑渗透率张量的各向异性油藏流线模拟方法

针对各向异性油藏,考虑全张量形式渗透率,建立三维两相不可压缩流体流线模拟方法。采用十九点差分格式对压力方程差分离散并进行求解得到压力场,对全张量渗透率下运动方程求解得到速度场,采用Pollock方法确定并追踪流线,最后根据一维水驱油理论得到任一流线内的含油饱和度分布。不同渗透率张量油藏模型计算结果表明:油藏压力沿最大主渗透率方向传播较快;由于流速计算时交叉项的引入,流线与压力梯度方向一般不一致;注采井方向与最大主渗透率方向垂直布置可有效减轻注入水沿最大主渗透率方向的突进,使流线在整个油藏中分布更加均匀,改善水驱开发效果。     

数字岩心物性数值模拟

针对岩心物性实验周期长且实验结果难以复现等缺点,利用数字岩心对储层岩石物性进行数值模拟。使用计算机断层扫描(CT)技术,获得CT序列图;并由序列图构建三维岩心图像。利用最大内切球方法建立岩心图像的孔隙网络模型,计算岩心渗透性以及岩心中孔喉参数;利用分形维数描述岩心图像的均质性。将数字岩心物性模拟结果和毛管压力曲线数据进行对比,分析了当前物性数值模拟和真实物性实验之间的差距。利用多组岩心样本进行了实验。结果表明,由于成像设备的原因,岩心图像只能表征岩石中的大尺寸孔隙特征;此外,CT图像得到的均质性和物性实验具有很好的一致性,渗透性和物性实验具有一定的出入。     

利用泥质含量与孔隙度计算砂岩渗透率

砂岩储层的渗透率受多个因素的影响,单相关模型往往有其局限性,为保证渗透率模型具有优良的解释能力和预测效果,本文提出一种综合利用岩心粒度分析与物性分析资料,根据泥质含量与孔隙度计算渗透率的方法,并且通过实际井资料的试处理验证,说明其计算方法是可靠的。     

An Index mechanism for multi-scale view in spatial databases

There are numerous geometric objects stored in the spatial databases. An importance function in a spatial database is that users can browse the geometric objects as a map efficiently. Thus the spatial database should display the geometric objects users concern about swiftly onto the display window. This process includes two operations:retrieve data from database and then draw them onto screen. Accordingly, to improve the efficiency, we should try to reduce time of both retrieving object and displaying them. The former can be achieved with the aid of spatial index such as R-tree, the latter require to simplify the objects. Simplification means that objects are shown with sufficient but not with unnecessary detail which depend on the scale of browse. So the major problem is how to retrieve data at different detail level efficiently. This paper introduces the implementation of a multi-scale index in the spatial database SISP (Spatial Information Shared Platform) which is generalized from R-tree. The difference between the generalization and the R-tree lies on two facets: One is that every node and geometric object in the generalization is assigned with a importance value which denote the importance of them, and every vertex in the objects are assigned with a importance value,too. The importance value can be use to decide which data should be retrieve from disk in a query. The other difference is that geometric objects in the generalization are divided into one or more sub-blocks, and vertexes are total ordered by their importance value. With the help of the generalized R-tree, one can easily retrieve data at different detail levels.Some experiments are performed on real-life data to evaluate the performance of solutions that separately use normal spatial index and multi-scale spatial index. The results show that the solution using multi-scale index in SISP is satisfying.     

Characterization of Ordovician carbonate reservoirs, southeastern Saskatchewan, Canada

The discovery of the prolific Ordovician Red River reservoirs in 1995 in southeastern Saskatchewan was the catalyst for extensive exploration activity which resulted in the discovery of more than 15 new Red River pools. The best yields of Red River production to date have been from dolomite reservoirs. Understanding the processes of dolomitization is, therefore, crucial for the prediction of the connectivity, spatial distribution and heterogeneity of dolomite reservoirs.The Red River reservoirs in the Midale area consist of 3～4 thin dolomitized zones, with a total thickness of about 20 m, which occur at the top of the Yeoman Formation. Two types of replacement dolomite were recognized in the Red River reservoir: dolomitized burrow infills and dolomitized host matrix. The spatial distribution of dolomite suggests that burrowing organisms played an important role in facilitating the fluid flow in the backfilled sediments. This resulted in penecontemporaneous dolomitization of burrow infills by normal seawater. The dolomite in the host matrix is interpreted as having occurred at shallow burial by evaporitic seawater during precipitation of Lake Almar anhydrite that immediately overlies the Yeoman Formation. However, the low δ18O values of dolomited burrow infills (-5.9‰～ -7.8‰, PDB) and matrix dolomites (-6.6‰～ -8.1‰, avg. -7.4‰ PDB) compared to the estimated values for the late Ordovician marine dolomite could be attributed to modification and alteration of dolomite at higher temperatures during deeper burial, which could also be responsible for its 87Sr/86Sr ratios (0.7084～0.7088) that are higher than suggested for the late Ordovician seawaters (0.7078～0.7080). The trace amounts of saddle dolomite cement in the Red River carbonates are probably related to "cannibalization" of earlier replacement dolomite during the chemical compaction.     

A computer generator for randomly layered structures

AcomputergeneratorforrandomlylayeredstructuresYUJia shun1,2,HEZhen hua2(1.TheInstituteofGeologicalandNuclearSciences,NewZealand;2.StateKeyLaboratoryofOilandGasReservoirGeologyandExploitation,ChengduUniversityofTechnology,China)Abstract:Analgorithmisintrod…     

Perspherical Approaches for Electrons in Atoms and Charged Excitons in Quantum Wells

Instead of following Fock’s expansion,we solve the Schrodinger equation for some quantum mechanical manybody systems such as electrons in atoms and charged excitons in quantum wells in a similar way in hyperspherical coordinates by expanding the wave functions into orthonormal complete basis sets of the hyperspherical hannonics(HHs)of hyperangles and generalized Laguerre polynomials(GLPs)of the hyperradius.This leads the equation to     

Advances in broadband wireless communications under high-mobility scenarios

Future mobile communication systems aim at providing very high data transmission rates, even in high-mobility scenarios such as high-speed wheel-track trains, maglev trains, highway vehicles, airplanes, guided missiles or spacecraft. A particularly important commercial application is the strong and increasing worldwide demand for high- speed broadband wireless communications （up to 574.8 km/ h test speeds or 380 km/h commercial speeds） in railways, providing data, voice and video services for applications such as onboard entertainment services to passengers, train control, train dispatch, train sensor status handling and sur- veillance. In such high-mobility scenarios, there are a number of communication challenges, including fast hand- over, location updating, high-speed channel modeling, estimation and equalization, anti-Doppler spreading tech- niques, fast power control, and dedicated network architec- ture. Because signal transmission in very high-speed scenarios will inevitably experience serious deterioration, it is imperative to develop key broadband mobile communi- cation techniques for such high-speed vehicles.     

Wuhan University Journal of Natural Sciences TOTAL CONTENTS

<正>~~     

海南岛地体及其毗邻陆缘晚中生代—新生代古地磁研究和构造演化

本文叙述了对海南岛及其毗邻大陆边缘白垩纪到第四纪地层岩石进行古地磁研究的全部工作过程。通过分析岩石中剩余磁矢量的磁偏角及磁倾角的变化,提出海南岛白垩纪以来经历的构造演化模式如下:早期伴随顺时针旋转而向南迁移,后期伴随逆时针转动并向北运移。联系该地区及邻区的地质、地球物理资料,对海南岛上述的构造地体运动提出以下认识:北部湾内早期有一拉张作用,主要是该作用使湾内地壳显著伸长减薄,形成北部湾盆地。从而导致了海南岛的早期构造运动,而海南岛后期的构造运动则主要是受南海海底扩张的影响。海南地体运动规律的阐明对于了解北部湾油气盆地的形成演化有重要的理论和实际意义。     

Scale-based spatial data model for GIS

Being the primary media of geographical information and the elementary objects manipulated, almost all of maps adopt the layer-based model to represent geographic information in the existent GIS. However, it is difficult to extend the map represented in layer-based model. Furthermore, in Web-Based GIS, It is slow to transmit the spatial data for map viewing. In this paper, for solving the questions above, we have proposed a new method for representing the spatial data. That is scale-based model. In this model we represent maps in three levels: scale-view, block, and spatial object, and organize the maps in a set of map layers, named Scale-View, which associates some given scales.Lastly, a prototype Web-Based GIS using the proposed spatial data representation is described briefly.     

Exciton Seebeck in Molecular Systems

Various applications relevant to the exciton dynamics,such as the organic solar cell,the large-area organic light-emitting diodes and the thermoelectricity,are operating under temperature gradient.The potential abnormal behavior of the exicton dynamics driven by the temperature difference may affect the efficiency and performance of the corresponding devices.In the above situations,the exciton dynamics under temperature difference is mixed with     

Achieving Unusual States of Matter Under High Pressure

The elongation method,originally proposed by Imamura was further developed for many years in our group.As a method towards O(N)with high efficiency and high accuracy for any dimensional systems.This treatment designed for one-dimensional(ID)polymers is now available for three-dimensional(3D)systems,but geometry optimization is now possible only for 1D-systems.As an approach toward post-Hartree-Fock,it was also extended to