基于反演基质矿物模量和多约束条件的双相介质AVO正演方法（英文）

基于双相介质理论的AVO正演技术是储层性质描述和流体预测的有效技术手段之一,但是输入参数中基质矿物模量的准确性和双相介质模型的的合理性极大地影响双相介质AVO正演效果。因此,本文采用基于流体因子的基质矿物模量反演方法,自适应反演基质矿物体积模量。引入具有岩石物理意义的多约束条件,使得流体替换技术制作的双相介质模型具有岩石物理意义。保证获得的双相介质AVO特征反映实际地层响应,真实可靠。通过不同岩性岩样的对比分析,说明反演方法的优越性和准确性。同时LH地区实际资料应用,获得孔隙度和流体饱和度等重要岩性参数变化时双相介质AVO特征,特别是不同储层孔隙度在同一入射角对应快纵波和横波反射系数幅值的大小差异和突变角差异是分辨储层孔隙度大小的依据。     

地层模量分解及在流体识别中的应用

储层流体识别是确定油气水分布,合理布设井位,提高钻井成功率的关键之一.本文基于流体饱和孔隙介质岩石物理模型,对地震反演的地层体积模量进行分解,获得孔隙流体体积模量,并依据油、气、水(尤其是气-油、气-水)模量的显著差异进行识别.文中简要分析了Gassmann模型和Kuster-Toksz模型的特征,详细讨论了孔隙形态和饱和度对弹性模量的影响,提出了联合Kuster-Toksz方程和Gassmann方程的体积模量分解方法.该方法通过Kuster-Toksz方程从测井数据中反演地层骨架固体和干骨架的弹性模量,再利用Gassmann方程对地层体积模量进行分解,既考虑了孔隙形态,又充分利用了Gassmann方程的易用性.理论模型结果表明方法是可行的.方法应用于西部地区某气田,流体识别与地层含气性预测结果与钻井基本一致,进一步证实了方法的有效性.     

流体体积模量计算方法研究

Wood模型、Patchy模型、Domenico模型及Brie的经验公式是常用的流体体积模量计算模型,目前低孔低渗或致密储层一般采用Brie的经验公式来计算流体体积模量.通过深入研究这几种模型,计算出流体体积模量的上下界,将上下界分别带入Gassmann方程反推出Brie干岩石剪切模量模型指数值范围,从指数范围内寻找一个最优值,使得纵横波预测误差最小,这个最优值即为剪切模型中的指数值.Brie剪切模型中采样点的指数值为固定值,将该固定值表示为随深度变化的变量,优化了Brie干岩石模量的计算方法.将优化后的Brie干岩石模型与Gassmann方程相结合反推出流体的体积模量.本文对Weyburn油田常规储层、胜利油田低孔低渗储层及苏里格气田致密储层进行研究,得出如下结论:(1)流体体积模量除了受各相流体的体积模量、含水饱和度、压力的影响外,还与孔隙的连通程度有关,即在有效压力不大的情况下,流体体积模量随含水饱和度的变化规律主要是连通性决定的;(2)低孔低渗、致密储层流体体积模量岩石物理模型与常规储层有很大的区别,Wood模型适用于常规储层流体体积模量的计算,而Wood模型和Domenico模型相结合的形式适用于低孔低渗和致密储层流体体积模量的计算.     

求取岩石基质体积模量的线形拟合方法(英文)

岩石基质的体积模量或其倒数--压缩系数,在进行油气预测的流体替换和孔隙度反演时,是重要的输入参数,但是利用现有方法很难准确求得。文中提出了一种求取该参数的线形拟合方法,该方法通过对Gassmann方程的合理简化并引入Eshelby-Walsh干燥岩石椭球包体近似公式,获得了计算岩石基质压缩系数的拟合公式,可方便地利用公式计算该参数。实际碳酸盐岩岩样的岩石物理测试分析显示:利用饱和岩样和干燥岩样测得的基质压缩系数的差异小于1%,说明所求参数是正确的,可靠的。     

基于约束条件横波速度反演和流体替代

依据辽东湾地区一组在地层条件和不同流体相态下（气饱和与水饱和）实验室测定的岩石纵波速度、横波速度及密度等数据,讨论Gassmann方程在横波速度反演和流体替代方面的应用效果,指出岩石骨架剪切模量和体积模量估算精度是影响纵波、横波速度预测精度的关键因素之一,提出基于约束条件横波速度反演和流体替代方法,速度预测精度的改善十分明显.     

砂泥岩干岩石模量估计

流体替换是了解和预测地震波速度和波阻抗如何依赖孔隙流体变化的有效工具,参数多、不确定性大是流体替代的特点,其中,干岩石模量是链接流体和饱和岩石的关键,也是Gassmann方程的基础,因此干岩石模量值的确定是流体替代的难点.前人对干岩石模量也进行了大量的研究,提出了许多经典模型:疏松砂岩模型、Kuster-Toksoz模型、自相容模型(selfconsistent模型)、微分有效介质模型(DEM模型)等,但这些模型都具有一定的局限性,本文通过对Gassmann方程图形分析方法(Mavko and Mukerji,1995)的研究,提出了一种计算干岩石模量的新方法.通过实际的岩样数据分析这几种算法的应用效果,研究结果认为求取干岩石模量时,本文提出的新方法具有很好的应用效果.     

测井资料Xu-White模型预测横波速度的一些新观点

目前横波预测的方法大致可以分为两种:经验公式预测和理论岩石物理模型。由于经验公式预测一般具有区域性,研究者更重视岩石物理模型预测。目前大多数岩石物理模型预测横波的方法假定地下流体的物性参数(速度和密度)不受地层深度的影响,且孔隙扁率是恒定的,实际上这并不科学。因为矿物的体积模量和剪切模量随所处地层深度发生改变,而对于孔隙扁率则随颗粒形状、孔隙度等的变化有较大变化。针对这些情况,提出一种新的改进的Xu-White横波预测方法,并可取得较好的效果。     

岩石基质模量与临界孔隙度的联合预测方法（英文）

岩石的基质模量和临界孔隙度是岩石物理模型的两项重要输入参数,但临界孔隙度通常很难获得。本文基于有效体积模量与孔隙度的线性模型,提出一种快速计算基质模量和临界孔隙度的方法。通过对实验室或现场测量的有效体积模量与孔隙度数据进行最小二乘拟合,利用拟合系数即可计算基质模量和临界孔隙度。本文方法适用的样品孔隙度范围广,并且计算结果能够准确反映样品的粘土含量、压力大小和饱和状态的差异。粘土含量高的样品的基质模量低,临界孔隙度小。压力对基质模量影响不大,但是临界孔隙度会随着压力的升高而略微增大。饱和样品计算得到的基质模量要比干燥的更大。     

基于变形P-L模型的矩阵方程迭代精细横波预测

横波速度预测问题的关键有两个,一是如何建立合理的岩石物理模型,二是针对建立的横波预测目标函数,如何准确高效地求解.针对第一个问题,对Pride模型和Lee模型(P-L模型)进行变形,提出拟固结指数的概念,将干岩石模量和岩石基质模量相联系,变形后的P-L模型在没有降低P-L模型准确度的情况下简化了问题的复杂度,建立起了饱和流体岩石弹性模量与干岩石模量、岩石基质模量、混合流体模量之间的关系,进而计算理论上的纵波速度,并通过比较实测纵波速度与计算的理论纵波速度大小,最终建立了横波预测的目标函数.针对第二个问题,借鉴了地震反演的思路,将该目标函数的最优化问题转化为线性矩阵方程组迭代求解问题,通过几步迭代就可以求解出合适的拟固结指数,进而得到预测横波速度.实际验证和应用表明,该横波预测方法具有很好的稳定性和准确性,并且岩石物理模型的构建和目标函数的求解思路可用于其他储集类型地层的横波预测.     

基于碳酸盐岩孔隙结构预测孔隙度方法研究

许多研究都已经证实在碳酸盐岩储层中孔隙结构对声波速度影响很大,因此在孔隙度反演时必须考虑孔隙结构的影响.本文通过对Gassmann方程的合理简化并引入Eshelby-Walsh干燥岩石椭球包体近似公式,推导出包含岩石孔隙结构参数、饱和岩石压缩系数与岩石基质压缩系数三个参数的计算孔隙度的新公式,该式岩石基质压缩系数是通过Russell流体因子和Gassmann-Boit-Geer tsma方程计算式建立目标函数反演获得的,然后通过测井交汇图技术把岩石孔隙结构参数与岩石基质压缩系数优选转换成纵横波速度与密度关系式,进而导出具体地区考虑碳酸盐岩孔隙结构孔隙度具体计算公式,测井与地震资料应用表明,基于碳酸盐岩孔隙结构的孔隙度预测方法的精度高于常规方法.     

Effective bulk modulus of the pore fluid at patchy saturation

We explore a package of parallel porous layers, each filled with a different fluid. Assume that this package is sampled by an elastic wave with the wavelength much larger than the thickness of an individual layer. Also assume that the layers are hydraulically isolated from each other, meaning that the diffusion length is smaller than that of the individual layer. This assumption is relevant to a patchy saturation scenario. Suppose that we wish to conduct the fluid substitution operation on this package treated as a single porous elastic body. What is the effective bulk modulus of the pore fluid to be used in this operation that will result in the same elastic modulus as computed by Backus averaging the individual moduli of the layers? We address this question analytically by assuming that the porosity, dry frame, and the mineral matrix properties of the individual layers are the same for all layers. The only difference between the layers is the pore fluid. We find that the resulting effective bulk modulus of the fluid thus derived falls between the arithmetic and harmonic averages of the fluid bulk moduli in the layers. It can be approximated by a linear combination of these two bounds where the weights are 0.50 and 0.50 or 0.75 for the arithmetic average and 0.25 for the harmonic average, depending on the elastic moduli of the dry frame, the mineral, and the pore fluids. This solution also provides a relation between the effective bulk modulus of the pore fluid in the system under examination and water saturation to be used in the fluid substitution operation at a coarse spatial scale.     

Recycling and “re-hydration” of degassed magma inducing transient dissolution/crystallization events at Stromboli (Italy)

Intrusive degassing and recycling of degassed and dense magma at depth have been proposed for a long time at Stromboli. The brief explosive event that occurred at the summit craters on 9 January 2005 threw out bombs and lapilli that could be good candidates to illustrate recycling of shallow degassed magma at depth. We present an extensive data set on both the textures and the mineral, bulk rock and glassy matrix chemistry of the “9 Jan” products. The latter have the common shoshonitic–basaltic bulk composition of lavas and scoriae issued from typical strombolian activity. In contrast they differ by the heterogeneous chemistry of their matrix glasses and their crystal textures that testify to crystal dissolution event(s) just prior magma crystallization upon ascent and eruption. Comparison between mineral paragenesis of the natural products and experimental phase equilibria suggest water-induced magma re-equilibration. We propose that mineral dissolution is related to water enrichment of the recycled degassed magma, via differential gas bubble transfer and to some extents its physical mixing with volatile-rich magma blobs. However, all these features illustrate transient processes. Even though evidence of mineral dissolution is ubiquitous at Stromboli, its effect on the bulk magma chemistry is minor because of the subtle interplay between mineral dissolution and crystallization in magmas having comparable bulk chemistry.     

中国南海北部陆坡孔隙度的求取

孔隙度在岩石物理和地震勘探领域都是一个非常重要的参数,水饱和地层的孔隙度可以用于求海底地层的弹性模量进而求取纵横波速度、识别水合物以及估计水合物的储量等.然而在没有打钻的情况下要相对准确地了解海底未固结成岩地层的孔隙度变化趋势是困难的.历史上有学者提出了几个模型,都是基于孔隙度随深度的增加呈现指数下降的趋势,并且需要知道指数模型中的若干参数,因此就产生了两个问题,一是孔隙度是否随深度增加而呈现指数下降,二是指数模型中参数如何选取.本文从另外一个角度,即利用海底介质的压实性从表层开始往下逐层求取水饱和地层的孔隙度,此方法被称之为"体积模量法",它仅需要知道表层孔隙度和地层矿物组分,从而节省了勘探成本.本文将该方法应用于我国南海地区和美国ODP井数据,结果与测井数据以及地震反演数据吻合得较好.     

北京塔院井水位微动态资料的分析和应用

对塔院井地下水水位微动态资料和气压观测资料,采用三种方法计算了气压效率,作了谱分析和调和分析,并利用所得结果求出了该井所在含水层的水文地质参数     

Rock-physics models of hydrate-bearing sediments in permafrost,Qilian Mountains,China

Rock-physics models are constructed for hydrate-bearing sediments in the Qilian Mountains permafrost region using the K–T equation model, and modes I and II of the effective medium model. The K–T equation models the seismic wave propagation in a twophase medium to determine the elastic moduli of the composite medium. In the effective medium model, mode I, the hydrate is a component of the pore inclusions in mode I and in mode II it is a component of the matrix. First, the P-wave velocity, S-wave velocity, density, bulk modulus, and shear modulus of the sediment matrix are extracted from logging data.. Second, based on the physical properties of the main components of the sediments, rockphysics model is established using the K–T equation, and two additional rock-physics models are established assuming different hydrate-filling modes for the effective medium. The model and actual velocity data for the hydrate-bearing sediments are compared and it is found that the rock-physics model for the hydrate-filling mode II well reproduces the actual data.     

Rock physics modelling of porous rocks with multiple pore types: a multiple-porosity variable critical porosity model

A critical porosity model establishes the empirical relationship between a grain matrix and a dry rock by the concept of critical porosity. The model is simple and practical and widely used. But the critical porosity in the model is a fixed value that cannot relate to pore structure. The aim of this paper is to establish the theoretical relationship between critical porosity and pore structure by combining Kuster–Toksöz theory with the critical porosity model. Different from the traditional critical porosity model, critical porosity is not an empirical value but varied with pore shape and the ratio of bulk modulus versus shear modulus of the grain matrix. The substitution of the theoretical relationship into Kuster–Toksöz theory will generate the formulae for the bulk and shear moduli of multiple-porosity dry rocks, which is named the multiple-porosity variable critical porosity model. The new model has been used to predict elastic moduli for sandstone and carbonate rock. We compare the modelling results for P- and S-wave velocities and elastic moduli with the experimental data. The comparison shows that the new model can be used to describe the elastic properties for the rocks with multiple pore types.     

Estimation of Dry Fracture Weakness,Porosity, and Fluid Modulus Using Observable Seismic Reflection Data in a Gas-Bearing Reservoir

Fracture detection and fluid identification are important tasks for a fractured reservoir characterization. Our goal is to demonstrate a direct approach to utilize azimuthal seismic data to estimate fluid bulk modulus, porosity, and dry fracture weaknesses, which decreases the uncertainty of fluid identification. Combining Gassmann’s (Vier. der Natur. Gesellschaft Zürich 96:1–23, 1951) equations and linear-slip model, we first establish new simplified expressions of stiffness parameters for a gas-bearing saturated fractured rock with low porosity and small fracture density, and then we derive a novel PP-wave reflection coefficient in terms of dry background rock properties (P-wave and S-wave moduli, and density), fracture (dry fracture weaknesses), porosity, and fluid (fluid bulk modulus). A Bayesian Markov chain Monte Carlo nonlinear inversion method is proposed to estimate fluid bulk modulus, porosity, and fracture weaknesses directly from azimuthal seismic data. The inversion method yields reasonable estimates in the case of synthetic data containing a moderate noise and stable results on real data.