基于频变AVO反演的页岩储层含气性地震识别技术

结合地震岩石物理技术,研究了叠前频变AVO反演在四川盆地龙马溪组页岩储层含气性识别中的应用.首先,应用Backus平均理论将测井数据粗化为地震尺度储层模型,应用传播矩阵理论进行高精度地震正演及井震标定,分析页岩气储层地震响应特征.其次,基于岩心观测结果,应用Chapman多尺度裂缝理论设计页岩气储层理论模型,研究储层衰减、频散以及对应的地震反射特征.应用该理论模型测试频变AVO反演方法,计算结果表明：对于研究区地层结构和地震数据,区分流体类型的优势频率不是地震子波的主频,还受层间调谐干涉等储层结构因素控制,也进一步说明理论模型测试和标定的重要性.最后,将频变AVO反演技术应用到四川盆地龙马溪组页岩地层,计算得到的频散属性为页岩气储层含气性识别提供依据.

     

Effects of petrophysical parameters on attenuation and dispersion of seismic waves in the simplified poroelastic theory: Comparison between the Biot's theory and viscosity‐extended theory

The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly.     

Rock-physics forward modelling to predict seismic behaviour: A case study for exploration target in Mahanadi basin,East Coast of India

One of the major aspects of rock-physics forward modelling is to predict seismic behaviour at an undrilled location using drilled well data. It is important to model the rock and fluid properties away from drilled wells to characterize the reservoir and investigate the root causes of different seismic responses. Using the forward modelling technique, it is possible to explain the amplitude responses of present seismic data in terms of probable rock and reservoir properties. In this context, rock-physics modelling adds significant values in the prospect maturation process by reducing the risk of reservoir presence in exploration and appraisal phases. The synthetic amplitude variation with offset gathers from the forward model is compared with real seismic gathers to ensure the fidelity of the existing geological model. ‘Prospect A’ in the study area has been identified from seismic interpretation, which was deposited as slope fan sediments in Mahanadi basin, East Coast of India. The mapped prospect has shown class-I amplitude variation with offset response in seismic without any direct hydrocarbon indicator support. The existing geological model suggests the presence of an excellent gas reservoir with proven charge access from the fetch area, moderate porosity and type of lithology within this fan prospect. But, whether the seismic response from this geological model will exhibit a class-I amplitude variation with offset behaviour or ‘dim spot’ will be visible; the objective of the present study is to investigate these queries. A rock-physics depth trend analysis has been done to envisage the possibilities of class-I reservoir in ‘Prospect A’. Forward modelling, using a combination of mechanical and chemical compaction, shows the synthetic gas gathers at ‘Prospect A’, which are class I in nature. The study has also depicted 2D forward modelling using lithology and fluid properties of discovery well within similar stratigraphy to predict whether ‘dim spot’ will be seen in seismic. The estimated change in synthetic amplitude response has been observed as ∼5% at contact, which suggests that the changes will not be visible in seismic. The study connects the existing geological model with a top-down seismic interpretation using rock-physics forward modelling technique to mature a deep-water exploratory prospect.     

Generalized anisotropy parameters and approximations of attenuations and velocities in viscoelastic media of arbitrary anisotropy type – theoretical and experimental aspects*

Seismic anisotropy in geological media is now widely accepted. Parametrizations and explicit approximations for the velocities in such media, considered as purely elastic and moderately anisotropic, are now standards and have even been extended to arbitrary types of anisotropy. In the case of attenuating media, some authors have also recently published different parametrizations and velocity and attenuation approximations in viscoelastic anisotropic media of particular symmetry type (e.g., transversely isotropic or orthorhombic). This paper extends such work to media of arbitrary anisotropy type, that is to say to triclinic media. In the case of homogeneous waves and using the so‐called ‘correspondence principle’, it is shown that the viscoelastic equations (for the phase velocities, phase slownesses, moduli, wavenumbers, etc.) are formally identical to the corresponding purely elastic equations available in the literature provided that all the corresponding quantities are complex (except the unit vector in the propagation direction that remains real). In contrast to previous work, the new parametrization uses complex anisotropy parameters and constitutes a simple extension to viscoelastic media of previous work dealing with non‐attenuating elastic media of arbitrary anisotropy type. We make the link between these new complex anisotropy parameters and measurable parameters, as well as with previously published anisotropy parameters, demonstrating the usefulness of the new parametrization. We compute the explicit complete directional dependence of the exact and of the approximate (first and higher‐order perturbation) complex phase velocities of the three body waves (qP, qS1 and qS2). The exact equations are successfully compared with the ultrasonic phase velocities and phase attenuations of the three body waves measured in a strongly attenuating water‐saturated sample of Vosges sandstone exhibiting moderate velocity anisotropy but very strong attenuation anisotropy. The approximate formulas are checked on experimental data. Compared to the exact solutions, the errors observed on the first‐order approximate velocities are small (<1%) for qP‐waves and moderate (<10%) for qS‐waves. The corresponding errors on the quality factor Q are moderate (<6%) for qP‐waves but critically large (up to 160%) for the qS‐waves. The use of higher‐order approximations substantially improves the accuracy, for instance typical maximum relative errors do not exceed 0.06% on all the velocities and 0.6% on all the quality factors Q, for third‐order approximations. All the results obtained on other rock samples confirm the results obtained on this rock. The simplicity of the derivations and the generality of the results are striking and particularly convenient for practical applications.     

Effect of fracture fill on frequency‐dependent anisotropy of fractured porous rocks

In fractured reservoirs, seismic wave velocity and amplitude depend on frequency and incidence angle. Frequency dependence is believed to be principally caused by the wave‐induced flow of pore fluid at the mesoscopic scale. In recent years, two particular phenomena, i.e., patchy saturation and flow between fractures and pores, have been identified as significant mechanisms of wave‐induced flow. However, these two phenomena are studied separately. Recently, a unified model has been proposed for a porous rock with a set of aligned fractures, with pores and fractures filled with two different fluids. Existing models treat waves propagating perpendicular to the fractures. In this paper, we extend the model to all propagation angles by assuming that the flow direction is perpendicular to the layering plane and is independent of the loading direction. We first consider the limiting cases through poroelastic Backus averaging, and then we obtain the five complex and frequency‐dependent stiffness values of the equivalent transversely isotropic medium as a function of the frequency. The numerical results show that, when the bulk modulus of the fracture‐filling fluid is relatively large, the dispersion and attenuation of P‐waves are mainly caused by fractures, and the values decrease as angles increase, almost vanishing when the incidence angle is 90° (propagation parallel to the fracture plane). While the bulk modulus of fluid in fractures is much smaller than that of matrix pores, the attenuation due to the “partial saturation” mechanism makes the fluid flow from pores into fractures, which is almost independent of the incidence angle.     

Sensitivity of the elastic anisotropy and seismic reflection amplitude of the Eagle Ford Shale to the presence of kerogen

The Eagle Ford Shale of Central and South Texas is currently of great interest for oil and gas exploration and production. Laboratory studies show that the Eagle Ford Shale is anisotropic, with a correlation between anisotropy and total organic carbon. Organic materials are usually more compliant than other minerals present in organic‐rich shales, and their shapes and distribution are usually anisotropic. This makes organic materials an important source of anisotropy in organic‐rich shales. Neglecting shale anisotropy may lead to incorrect estimates of rock and fluid properties derived from inversion of amplitude versus offset seismic data. Organic materials have a significant effect on the PP and PS reflection amplitudes from the Austin Chalk/Upper Eagle Ford interface, the Upper Eagle Ford/Lower Eagle Ford interface, and the Lower Eagle Ford/Buda Limestone interface. The higher kerogen content of the Lower Eagle Ford compared with that of the Upper Eagle Ford leads to a negative PP reflection amplitude that dims with offset, whereas the PS reflection coefficient increases in magnitude with increasing offset. The PP and PS reflection coefficients at the Austin Chalk/Upper Eagle Ford interface, the Upper Eagle Ford/Lower Eagle Ford interface, and the Lower Eagle Ford/Buda Limestone interface all increase in magnitude with increasing volume fraction of kerogen.     

Estimating azimuthal stress‐induced P‐wave anisotropy from S‐wave anisotropy using sonic log or vertical seismic profile data

Most sedimentary rocks are anisotropic, yet it is often difficult to accurately incorporate anisotropy into seismic workflows because analysis of anisotropy requires knowledge of a number of parameters that are difficult to estimate from standard seismic data. In this study, we provide a methodology to infer azimuthal P‐wave anisotropy from S‐wave anisotropy calculated from log or vertical seismic profile data. This methodology involves a number of steps. First, we compute the azimuthal P‐wave anisotropy in the dry medium as a function of the azimuthal S‐wave anisotropy using a rock physics model, which accounts for the stress dependency of seismic wave velocities in dry isotropic elastic media subjected to triaxial compression. Once the P‐wave anisotropy in the dry medium is known, we use the anisotropic Gassmann equations to estimate the anisotropy of the saturated medium. We test this workflow on the log data acquired in the North West Shelf of Australia, where azimuthal anisotropy is likely caused by large differences between minimum and maximum horizontal stresses. The obtained results are compared to azimuthal P‐wave anisotropy obtained via orthorhombic tomography in the same area. In the clean sandstone layers, anisotropy parameters obtained by both methods are fairly consistent. In the shale and shaly sandstone layers, however, there is a significant discrepancy between results since the stress‐induced anisotropy model we use is not applicable to rocks exhibiting intrinsic anisotropy. This methodology could be useful for building the initial anisotropic velocity model for imaging, which is to be refined through migration velocity analysis.     

地震频段弹性模量测试系统改进与升级

获取地震频段弹性模量对于地震数据定量解释、研究地震波传播特征及油气勘探开发具有重要意义.前期建立的测试系统基于应力应变法可以获得地震频段内岩心的杨氏模量、泊松比和衰减.该系统仅适用于ϕ38 mm岩心, 而测量物性参数(孔隙度、渗透率等)、测量超声纵横波速度时一般使用ϕ25 mm岩心.为将ϕ25 mm岩心应用于地震频段弹性模量的测试, 统一各项岩石物理测试的数据, 本文结合有限元数值模拟和实验测试对地震频段弹性模量测试系统的机械结构与应变采集系统(惠斯通电桥)做了改进与升级.通过数值模拟与标定实验(铝和有机玻璃)论证了在ϕ25 mm岩心上进行低频实验的可靠性和准确性.利用改进后的测试系统对常规砂岩在干燥(空气)、充气(N₂)与完全饱和流体(白油)条件开展了高频(MHz)超声和低频(1~3 kHz)应力应变测量实验, 其结果表明: 干燥岩样和充气岩样的弹性参数在低频和超声频段一致, 以此验证了改进后测量系统的可靠性; 而饱和白油的砂岩弹性参数在1~3000 Hz频段范围内表现出明显的频散特性.改进后的地震频段弹性模量测试系统可以为定量研究含流体岩石的弹性频散特征提供有效的实验测量工具, 其实验结果不仅可应用于生产实践, 也可为检验、修正岩石物理模型的提供重要的实验依据.

     

流体饱和度对裂缝-孔隙岩石频变纵波速度及其各向异性的影响

流体饱和度会改变裂缝性储层的纵波速度,从而影响地层速度的频散特性及各向异性程度,导致储层地震响应特征复杂,储层预测多解性强,流体识别难度大.本文根据多相流体饱和裂缝性储层的特点,借助于Norris和KG模型,建立部分饱和裂缝-孔隙等效介质模型,给出频变地震波速度随流体饱和度变化的精确关系式.数值模拟结果表明,当气、水两相共存时,随着含水饱和度的增加,高频段纵波相速度逐渐增大,各向异性程度逐渐减小;低频段纵波相速度逐渐减小,各向异性程度不变;相速度频散及其各向异性程度逐渐增强.组合已有的孔隙弹性理论模型,对实验室人工裂缝-孔隙砂岩岩样的纵波速度进行拟合,计算得到的曲线与实验室测量散点值吻合度较高,表明组合模型在给定参数下的有效性.该研究能够为多相流体饱和裂缝性储层的地震响应特征分析奠定扎实的理论基础,为提高储层预测的确定性和流体识别的准确性提供可靠的理论依据.

     

Influence of subsurface injection on time‐lapse seismic: laboratory studies at seismic and ultrasonic frequencies

Seismic monitoring of reservoir and overburden performance during subsurface CO₂ storage plays a key role in ensuring efficiency and safety. Proper interpretation of monitoring data requires knowledge about the rock physical phenomena occurring in the subsurface formations. This work focuses on rock stiffness and elastic velocity changes of a shale overburden formation caused by both reservoir inflation induced stress changes and leakage of CO₂ into the overburden. In laboratory experiments, Pierre shale I core plugs were loaded along the stress path representative for the in situ stress changes experienced by caprock during reservoir inflation. Tests were carried out in a triaxial compaction cell combining three measurement techniques and permitting for determination of (i) ultrasonic velocities, (ii) quasistatic rock deformations, and (iii) dynamic elastic stiffness at seismic frequencies within a single test, which allowed to quantify effects of seismic dispersion. In addition, fluid substitution effects connected with possible CO₂ leakage into the caprock formation were modelled by the modified anisotropic Gassmann model. Results of this work indicate that (i) stress sensitivity of Pierre shale I is frequency dependent; (ii) reservoir inflation leads to the increase of the overburden Young's modulus and Poisson's ratio; (iii) in situ stress changes mostly affect the P‐wave velocities; (iv) small leakage of the CO₂ into the overburden may lead to the velocity changes, which are comparable with one associated with geomechanical influence; (v) non‐elastic effects increase stress sensitivity of an acoustic waves; (iv) and both geomechanical and fluid substitution effects would create significant time shifts, which should be detectable by time‐lapse seismic.     

Effects of fracture intersections on seismic dispersion: theoretical predictions versus numerical simulations

The detection and characterisation of domains of intersecting fractures are important goals in several disciplines of current interest, including exploration and production of unconventional reservoirs, nuclear waste storage, CO₂ sequestration, and groundwater hydrology, among others. The objective of this study is to propose a theoretical framework for quantifying the effects of fracture intersections on the frequency‐dependent elastic properties of fluid‐saturated porous and fractured rocks. Three characteristic frequency regimes for fluid pressure communication are identified. In the low‐frequency limit, fractures are in full pressure communication with the embedding porous matrix and with other fractures. Conversely, in the high‐frequency limit, fractures are hydraulically isolated from the matrix and from other fractures. At intermediate frequencies, fractures are hydraulically isolated from the matrix porosity but can be in hydraulic communication with each other, depending on whether fracture sets are intersecting. For each frequency regime, the effective stiffness coefficients are derived using the linear‐slip theory and anisotropic Gassmann equations. Explicit mathematical expressions for the two characteristic frequencies that separate the three frequency regimes are also determined. Theoretical predictions are then applied to two synthetic 2D samples, each containing two orthogonal fracture sets: one with and another without intersections. The resulting stiffness coefficients, Thomsen‐style anisotropy parameters, and the transition frequencies show good agreement with corresponding numerical simulations. The theoretical results are applicable not only to 2D but also to 3D fracture systems and are amenable to being employed in inversion schemes designed to characterise fracture systems.     

Integrating basin modeling with seismic technology and rock physics

We explore the link between basin modelling and seismic inversion by applying different rock physics models. This study uses the E‐Dragon II data in the Gulf of Mexico. To investigate the impact of different rock physics models on the link between basin modelling and seismic inversion, we first model relationships between seismic velocities and both (1) porosity and (2) effective stress for well‐log data using published rock physics models. Then, we build 1D basin models to predict seismic velocities derived from basin modelling with different rock physics models, in a comparison with average sonic velocities measured in the wells. Finally, we examine how basin modelling outputs can be used to aid seismic inversion by providing constraints for the background low‐frequency model. For this, we run different scenarios of inverting near angle partial stack seismic data into elastic impedances to test the impact of the background model on the quality of the inversion results. The results of the study suggest that the link between basin modelling and seismic technology is a two‐way interaction in terms of potential applications, and the key to refine it is establishing a rock physics models that properly describes changes in seismic signatures reflecting changes in rock properties.     

多波时移地震AVO反演研究

数值模拟了油藏含油饱和度与有效压力变化时移地震AVO的响应,确定利用时移地震AVO区分油藏参数的变化、实现油藏定量解释的可行性.从Aki等 AVO近似方程出发,详细推导了P_P波和P_S转换波时移地震AVO计算公式.结合岩石物理近似关系和本文推导的时移地震AVO计算公式,推导了利用多波时移地震AVO反演油藏含油饱和度和压力变化的方程.数据试验表明,文中推导的多波时移地震AVO方程能较好地反演油藏含油饱和度变化和有效压力变化,实现油藏定量解释.     

Compaction trend versus seismic anisotropy in shaly formations

Shales comprise more than 60% of sedimentary rocks and form natural seals above hydrocarbon reservoirs. Their sealing capacity is also used for storage of nuclear wastes. The world's most important conventional oil and gas reservoirs have their corresponding source rocks in shale. Furthermore, shale oil and shale gas are the most rapidly expanding trends in unconventional oil and gas. Shales are notorious for their strong elastic anisotropy, i.e., so‐called vertical transverse isotropy. This vertical transverse isotropy, characterised by a vertical axis of invariance, is of practical importance as it is required for correct surface seismic data interpretation, seismic to well tie, and amplitude versus offset analysis. A rather classical paradigm makes a clear link between compaction in shales and the alignment of the clay platelets (main constituent of shales). This would imply increasing anisotropy strength with increasing compaction. Our main purpose is to check this prediction on two large databases in shaly formations (more than 800 samples from depths of 0–6 km) by extracting the major trends in the relation between seismic anisotropy and compaction. The statistical analysis of the database shows that the simultaneous increase in density and velocity, a classical compaction signature, is quite weakly correlated with the anisotropy strength. As a consequence, compaction can be excluded as a major cause of seismic anisotropy, at least in shaly formations. Also, the alignment of the clay platelets can explain most of the anisotropy measurements of both databases. Finally, a method for estimating the orientation distribution function of the clay platelets from the measurement of the anisotropy parameters is suggested.     

A strategy for automated analysis of passive microseismic data to image seismic anisotropy and fracture characteristics

Monitoring of induced seismicity is gaining importance in a broad range of industrial operations from hydrocarbon reservoirs to mining to geothermal fields. Such passive seismic monitoring mainly aims at identifying fractures, which is of special interest for safety and productivity reasons. By analysing shear‐wave splitting it is possible to determine the anisotropy of the rock, which may be caused by sedimentary layering and/or aligned fractures, which in turn offers insight into the state of stress in the reservoir. We present a workflow strategy for automatic and effective processing of passive microseismic data sets, which are ever increasing in size. The automation provides an objective quality control of the shear‐wave splitting measurements and is based on characteristic differences between the two independent eigenvalue and cross‐correlation splitting techniques. These differences are summarized in a quality index for each measurement, allowing identification of an appropriate quality threshold. Measurements above this threshold are considered to be of good quality and are used in further interpretation. We suggest an automated inversion scheme using rock physics theory to test for best correlation of the data with various combinations of fracture density, its strike and the background anisotropy. This fully automatic workflow is then tested on a synthetic and a real microseismic data set.     

中国大陆地壳与上地幔地震各向异性研究

地壳与上地幔各向异性研究在地球动力学领域有广泛的应用和重要的科学意义,其研究有助于许多地质和地球物理基本问题的解释.中国大陆的各向异性研究,在地壳与上地幔两个尺度上都有许多重要进展,并在基本理论和岩石实验研究中取得进步.文章认为,各向异性研究揭示的科学问题将会是广泛的和深远的.     

Research progress on layered seismic anisotropy - A review

Seismic anisotropy is an effective feature to study the inner structure of the Earth. In complex tectonic area, the assumption of single-layer anisotropy is sometimes not well consistent with the observed data; thus, the assumption of multi-layered (i.e. stratified) anisotropy should be considered. At present, the main methods to study anisotropy include receiver functions, shear wave splitting from local and teleseismic events (SKS, SKKS, and PKS, hereafter collectively called XKS), P- and Pn wave travel time inversion, surface wave inversion from far-field earthquakes and ambient noise. Each of the above method has its own advantages and limitations. Thus, one or more of the above methods are often combined to characterize multi-layered anisotropy, of which the depth range of anisotropic layers are different. This paper reviews the research progress of multi-layered anisotropy for the purpose of providing a basis for future seismic anisotropy investigations.     

汶川地震余震序列的地震各向异性

利用2008年5月12日汶川地震震源区及周边地震台站记录的余震序列资料,使用剪切波分裂系统分析法,分析了汶川地震发震构造龙门山断裂带及周边地区的地壳各向异性特征,推断了地壳最大主压应力方向及空间分布特征.研究结果表明：大致以安县为界,位于龙门山北东段的台站快剪切波的偏振方向为北东向,与断裂带走向一致;而位于龙门山西南段的台站快剪切波的偏振方向为北西向,与断裂带走向垂直;这个特征同样揭示出龙门山断裂带西南段逆冲、北东段带有明显走滑性质的特征.研究还显示,靠近龙门山与鲜水河、安宁河小江断裂交汇区附近的台站快剪切波的偏振方向表现比较离散,这可能是由震源区局部的复杂地质构造引起,与该地区复杂的主压应力方向特点一致.     

鄂尔多斯块体北缘与西缘地区地壳各向异性特征

本研究使用内蒙古自治区数字测震台网2010年1月至2017年10月区域小地震的波形记录资料,采用SAM方法,进行了地壳剪切波分裂的分析,得到鄂尔多斯块体北缘与西缘地区地壳介质地震各向异性的初步研究结果.根据15个台站161个有效地震记录的分析,鄂尔多斯块体北缘与西缘地区的快剪切波平均偏振方向为NE44.4°±38.4°,慢剪切波平均时间延迟为1.7±1.6 ms·km^-1.研究区域的快剪切波偏振显示出两个优势方向,一个是NE方向,另一个是近NS方向.区内的逆冲凸起与走滑正倾断层构造对剪切波分裂产生了直接的影响,造成了剪切波分裂参数的复杂分布,反映了剪切波分裂参数受到区域应力和构造共同作用的影响.鄂尔多斯块体北缘的快波偏振特征有NE和近NS两个优势偏振方向,其东区与西区的快剪切波偏振表现出明显不同的特征.东区的第一快剪切波优势偏振方向为NE,第二快剪切波优势偏振方向为近NS;西区的第一快剪切波优势偏振方向为近EW,第二快剪切波优势偏振方向为近NS.鄂尔多斯块体北缘的区域背景主压应力方向可能总体上为近NS方向,但空间分布有差异,东区NE方向的优势偏振与西区近EW方向的优势偏振更可能反映了断裂与构造的影响.鄂尔多斯块体西缘的快剪切波偏振特征显示出非常清楚的NE向的优势偏振方向,近NS向的优势偏振方向则不太明显,反映出该地区复杂构造对各向异性分布的影响.慢波时间延迟呈现出西低东高的特点,时间延迟的高值出现在鄂尔多斯块体北缘的东部,时间延迟的这种西低东高的各向异性强度变化,可能反映了区域构造活动西强东弱的特性.

     

基于两级尺度粗化的裂缝-多孔隙介质地震响应分析方法

实际地震勘探中,储层物性参数的差异是导致地震波响应特征发生变化的根本原因,而建立储层物性参数与地震响应特征之间的联系,需要跨越微观孔隙尺度、介观测井尺度以及宏观地震尺度等三个不同尺度空间.本文基于已知井的岩石物理实验数据和测井数据,利用复杂多孔隙介质理论将微观尺度孔隙岩石粗化到介观测井尺度,利用Backus平均理论将介观测井尺度的模型进一步粗化到宏观地震尺度,最终,得到地震尺度裂缝-多孔隙介质模型.其数值计算结果与测井数据和地震数据的对比表明：基于两级尺度粗化算法的裂缝多孔隙介质模型在给定参数下是有效的,且基于该模型的地震响应特征分析方法能够对储层的地震响应特征随物性参数的变化进行分析.