用地震资料估算珊溪水库岩层孔隙度和饱和度

本文首先从室内岩石力学实验数据出发,利用Biot-Gassmann方程和Pride模型建立的岩石骨架与岩石基质体积模量函数关系式求出珊溪水库区岩石基质模量和固结系数,然后分别利用联合地震波走时方程组和多台和达曲线求得P波速度和波速比,最后以流纹斑岩为例,计算了岩石饱和度、孔隙度,并分析了地震中造成的饱和度和孔隙度的变化。结果表明:(1)地震过程中岩石饱和度均大于99%,且大多为100%,即地震中岩石处于水饱和或接近水饱和状态,波速和波速比的变化主要由岩石孔隙度的变化引起。(2)珊溪水库震中区流纹斑岩孔隙度的上限值为6.6%~15.5%,下限值为3.6%~4.2%。下限值与地震前由实验室测得的孔隙度相近。如果实验室测得的流纹斑岩孔隙度是其初始值,则地震中岩石孔隙度增大了3.2%~12.09%。(3)2008年1月地震以后震中区岩石孔隙度逐渐闭合,直到2014年8月孔隙度开始迅速增大,9月15日起原震区再次发生一组新的地震活动,震群活动前岩石出现了扩容现象。     

岩体裂隙充水后对拐角频率的影响

2002年在浙江温州的珊溪水库发生了一系列水库诱发地震，其拐角频率明显小于同震级的构造地震。本从二相介质理论出发，讨论了水库岩体在扩容、流体渗入饱和和强度弱化过程中波速和波速比的变化以及对拐角频率的影响。根据实验结果，当孔隙度为0．15，岩体弹性模量减小23％时，拐角频率可以减小约29％。由此解释了水库诱发地震拐角频率偏小主要是由裂隙发育和库水渗透扩散使地壳浅部岩体强度弱化引起的。因此拐角频率可以作为区别两的一个依据。     

温州珊溪水库地震重新定位与速度结构联合反演

通过震源位置和速度结构联合反演方法,利用浙江和福建区域地震台网和珊溪水库台网给出的P波走时资料,得到了珊溪水库地区的三维速度结构,重新确定了珊溪水库地震的震源参数。结果表明:①震中总体呈现NW向分布,NW走向的双溪-焦溪垟断裂可能为珊溪水库地震序列的发震断层。②珊溪水库地震震源深度最大为9.5km,平均为5.4km,小于华南地震区10km的平均震源深度。③水库北、南两岸的地震较浅,震源深度均小于5km,水库淹没区地震较深。水库诱发地震之初的几年中震源深度有一个逐渐变大的过程,这可能是由于库水逐渐往深部下渗,从而进一步诱发了更深处的地震所致。④研究区存在一个低P波速度异常区,低速区位于水库淹没区内的多组断裂交汇部位,地震大多发生在低速异常区内。这可能与水库蓄水后库水下渗有关。     

双相介质的AVO正演模拟

岩石的孔隙度、流体饱和度等信息是影响地震波振幅随炮检距变化（AVO）的重要因素.本文在实验给定了岩石的物性参数（孔隙度及孔隙流体的不同相态）,利用Gassmann方程计算储层条件下的纵、横波速度,通过模拟不同类型的孔隙流体的地震响应,研究双相介质中流体成分的变化对地震反射波AVO的影响.     

水库诱发地震中水库水体下方的地震波速比

以大桥水库MS4.6诱发地震和紫坪铺水库ML2.5诱发震群为例,研究了穿过水体下方的地震波的波速比在地震前后的变化.由彝海子台记录的穿过大桥水库水体下方的地震波所计算的波速比,在蓄水后4.6级主震前存在高值异常,反映库水对波速比有影响;由八角台记录的穿过紫坪铺水库水体下方的地震波所计算的波速比显示,在2006年8月28日后一时段库区东北部出现微小的波速比高值异常,结果于2006年10月至2007年2月在八角台附近的汶川水磨发生了2.5级水库诱发震群.因此认为,根据穿过水库水体下方的地震波射线计算的单台波速比,能够反映地震孕育过程中震源及周围地区介质特性的改变,可作为监测水库地震的一种手段.     

松原前郭地震区孕震构造的地震层析成像研究

本文选取2013年10月到2018年12月期间松原市周边17个地震台记录到的515个地震事件的2926条P波和2665条S波到时数据,采用双差成像方法联合反演松原地震区震源参数和中上地壳(0～25 km)高分辨率三维体波速度和波速比结构,进而利用O'Connell-Budiansky理论估计了该区岩石介质的裂缝密度和饱和度参数.重定位后的震源参数精度有了显著提高.松原震区P波和S波的速度、波速比结构均表现出明显的不均匀特性,主要表现为低P波、低S波、低波速比结构,推测松原地震区发生更大地震的可能性较低.松原震群区浅地表区域(＜5 km)整体表现为低P波、低S波和高波速比结构,而岩石介质的裂缝密度和饱和度表现为高值.浅层高波速比、高裂缝密度和高饱和度结构可能与流体有关;从区域构造出发表明松原地震发震构造与东北地区的区域构造应力场密切相关.松原震区的三维精细速度和高精度震源位置信息为深入了解该地区的地壳结构、孕震构造和地震活动性提供重要参考.     

含流体砂岩地震波频散实验研究

为了研究孔隙流体对不同渗透率岩石地震波速度的影响,在实验室利用跨频带岩石弹性参数测试系统得到了应变幅值10-6的2~2000Hz频段下的地震波速度和1 MHz频率下的超声波速度,利用差分共振声谱法得到了频率600Hz岩石干燥和完全饱水情况下岩石声学参数.实验表明,在低饱和度下,致密砂岩在地震和超声频段下没有明显的频散;在高饱和度下纵波速度的频散变得明显.从干燥到完全水饱和条件,不同频率测量的致密砂岩的体积模量随岩石孔隙度增高而降低,且体积模量的变化量受岩石微观孔隙结构的影响较大.高孔、高渗砂岩无论在低含水度下还是在高含水饱和度下频散微弱,并且在地震频段下围压对于岩石纵横波速度的影响要大于频率的影响.高孔、高渗砂岩和致密砂岩不同含水饱和度下的频散差异可应用于储层预测,油气检测等方面,同时该研究可以更好地帮助理解岩石的黏弹性行为,促进岩石物理频散理论的发展,提高地震解释的精度.     

基于等效Thomsen参数的P-SV波AVO属性研究

在等效Thomsen各向异性参数的P-SV波反射系数近似公式基础上研究了反射系数的多种AVO(振幅随炮检距变化)属性特征,针对不同的属性特征构建了多属性AVO交绘图;并利用反射系数公式对三类含气砂岩AVO的特征进行分析.结果表明P-SV波反射系数公式可以有效的区分第三类含气砂岩;岩石的孔隙度、流体饱和度等信息是影响地震波AVO的重要因素.利用Gassmann方程进行了对上层为HTI介质,下层为孔隙储层的介质模型进行了流体替换计算,分析了孔隙度、含气饱和度和各向异性参数变化对地层AVO的影响.     

连续胶结声速理论胶结半径表达式的一般性推广及应用

针对Dvorkin 和Nur （1996） 提出的连续胶结声速理论（CCT）中的胶结半径表达式仅适用于孔隙式胶结疏松砂岩的局限性,通过胶结物与颗粒之间的空间几何关系,在胶结物中心厚度不为零的前提下,推导了胶结物胶结半径的一般表达式,并考察了胶结物中心厚度对疏松砂岩声波速度的影响,探讨了该表达式在基底式胶结疏松砂岩纵横波速度预测中的应用.结果表明：随胶结物中心厚度变大,纵、横波速度均减小,且孔隙度越大,胶结物含量越少,速度下降越快;孔隙度或胶结物含量不变时,纵横波速比随胶结物中心厚度的增大而增大;纵横波速比随孔隙度的变化规律与胶结物中心厚度有关,当胶结物中心厚度较小时,孔隙度增大,纵横波速比略微减小,而当胶结物中心厚度较大时,孔隙度的增大会引起纵横波速比略微增大;对于基底式胶结疏松砂岩,使用推广后表达式修正的CCT模型相对原CCT模型能够更好地预测其声波速度,同时发现,利用其纵横波速比无法预测由孔隙度或胶结物含量变化引起的疏松砂岩软硬的变化,但是在胶结物含量或孔隙度不变的情况下,纵横波速比可指示由胶结物中心厚度变化引起的疏松砂岩软硬的变化.     

珊溪水库M_L4.6震群精确定位与发震构造研究

利用浙江、福建区域地震台网和珊溪水库台网给出的P波走时资料,通过震源位置和速度结构联合反演方法,重新确定了珊溪水库地震的震源参数。并通过假设发震断层可以用一个平面进行拟合的思路,利用重新确定的震源参数拟合得到了发震断层面参数。收集了2267条P波到时数据,求解了ML4.0以上地震震源机制解,结合水库区的断层活动性特征和宏观调查资料,讨论了珊溪水库地震序列的发震构造。     

基于气体包裹体模型的低孔低渗储层声波速度研究

随着石油勘探工业的持续发展和技术水平的日益提高,低孔低渗油气藏已成为我国油气勘探开发的重要领域之一,但是该类油气藏的储层岩石物理关系复杂,对其评价也相对较难。本文针对中深层气藏低孔低渗储层评价存在的困难,采用修正的White气体包裹体模型开展岩石物理研究。首先对气体包裹体模型进行分析,得到纵横波速度的计算公式,进而理论计算并分析纵横波速度与孔隙度、饱和度、压力及温度参数的变化关系,最后结合岩心实验对理论计算结果进行了验证。研究结果表明基于气体包裹体模型的理论计算结果与实测数据吻合较好,可以较好地为低孔低渗复杂储层声波速度测量及解释提供技术支持。      

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

A physical model for shear-wave velocity prediction1

The clay-sand mixture model of Xu and White is shown to simulate observed relationships between S-wave velocity (or transit time), porosity and clay content. In general, neither S-wave velocity nor S-wave transit time is a linear function of porosity and clay content. For practical purposes, clay content is approximated by shale volume in well-log applications. In principle, the model can predict S-wave velocity from lithology and any pair of P-wave velocity, porosity and shale volume. Although the predictions should be the same if all measurements are error free, comparison of predictions with laboratory and logging measurements show that predictions using P-wave velocity are the most reliable. The robust relationship between S- and P-wave velocities is due to the fact that both are similarly affected by porosity, clay content and lithology. Moreover, errors in the measured P-wave velocity are normally smaller than those in porosity and shale volume, both of which are subject to errors introduced by imperfect models and imperfect parameters when estimated from logs. Because the model evaluates the bulk and shear moduli of the dry rock frame by a combination of Kuster and Toksöz’ theory and differential effective medium theory, using pore aspect ratios to characterize the compliances of the sand and clay components, the relationship between P- and S-wave velocities is explicit and consistent. Consequently the model sidesteps problems and assumptions that arise from the lack of knowledge of these moduli when applying Gassmann's theory to this relationship, making it a very flexible tool for investigating how the v_P-v_s relationship is affected by lithology, porosity, clay content and water saturation. Numerical results from the model are confirmed by laboratory and logging data and demonstrate, for example, how the presence of gas has a more pronounced effect on P-wave velocity in shaly sands than in less compliant cleaner sandstones.     

岩石的等效孔隙纵横比反演及其应用

通过融合Gassmann方程和由微分等效介质理论建立的干岩石骨架模型--DEM解析模型,本文提出根据纵波(和横波)速度反演岩石等效孔隙纵横比进行储层孔隙结构评价和横波速度预测的方法.首先,利用Gassmann方程和DEM解析模型建立岩石的纵、横波速度与密度、孔隙度、饱和度和矿物组分等各参数之间的关系;其次,将岩石孔隙等效为具有单一纵横比的理想椭球孔,应用非线性全局寻优算法来寻找最佳的等效孔隙纵横比使得理论预测与实际测量的弹性模量之间的误差最小;最后,将反演得到的等效孔隙纵横比代入到Gassmann方程和DEM解析模型中构建横波速度.实验室和井孔测量数据应用表明,反演得到的等效孔隙纵横比可准确反映储层的孔隙结构,对于裂缝型储层如花岗岩,其孔隙纵横比通常小于0.025,而对于孔隙型储层如砂岩,其孔隙纵横比通常大于0.08.只利用纵波与同时利用纵、横波反演得到的孔隙纵横比结果几乎完全一致,而且由纵波构建的横波与实测横波吻合良好,说明本文提出的等效孔隙纵横比反演及其横波速度预测方法是有效的.     

南海北部神狐海域天然气水合物分解的测井异常

南海北部神狐海域GMGS-1钻探揭示SH3井天然气水合物位于稳定带上部,厚度约为10 m.氯离子异常计算的水合物饱和度最高达26%,高水合物饱和度层出现高电阻率和低纵波速度.为分析该低纵波速度异常,本文基于简化的三相介质理论计算了饱和水纵波速度,在深度195 m附近,测量的纵波速度小于饱和水纵波速度.利用阿尔奇公式,基于原位温度、盐度、密度孔隙度和测量的电阻率,利用交会分析确定了该井的阿尔奇常数为a=1.1和m=2.3.基于该参数,利用阿尔奇方程计算的水合物饱和度占孔隙空间5%~20%,局部地层水合物饱和度达26.8%,在垂向上分布不均匀.由于钻探可能导致水合物发生分解而产生游离气,原位游离气和水合物分解产生的气体都能造成低纵波速度异常.由于地震资料采集在测井之前完成,利用不同速度制作合成地震记录并与地震资料进行对比,能够确定水合物稳定带上部的低速异常形成原因.     

Compressional- and shear-wave velocities as a function of confining stress in dry sandstones

A laboratory study was carried out to investigate the influence of confining stress on compressional- and shear-wave velocities for a set of rock samples from gas-producing sandstone reservoirs in the Cooper Basin, South Australia. The suite of samples consists of 22 consolidated sublitharenites with helium porosity ranging from 2.6% to 16.6%. We used a pulse-echo technique to measure compressional- and shear-wave velocities on dry samples (cylindrical 4.6 × 2 cm) at room temperature and at elevated confining stress (≤ 60 MPa). Compressional- and shear-wave velocities in samples increase non-linearly with confining stress. A regression equation of the form V = A ? Be^?DP gives a good fit to the measured velocities with improved prediction of velocities at high confining stresses compared with equations suggested by other studies. The predicted microcrack-closure stresses of the samples show values ranging from 70 MPa to 95 MPa and insignificant correlation with porosity, permeability or clay content. There is a positive correlation between change in velocity with core porosity and permeability, but this association is weak and diminishes with increasing confining stress. Experimental results show that pore geometry, grain-contact type, and distribution and location of clay particles may be more significant than total porosity and clay content in describing the stress sensitivity of sandstones at in situ reservoir effective stress. The stress dependence of Cooper Basin sandstones is very large compared with data from other studies. The implication of our study for hydrocarbon exploration is that where the in situ reservoir effective stress is much less than the microcrack-closure stress of the reservoir rocks, the variation of reservoir effective stress could cause significant changes in velocity of the reservoir rocks. The velocity changes induced by effective stress in highly stress-sensitive rocks can be detected at sonic-log and probably surface-seismic frequencies.     

基于常规测井数据的纵波本征衰减估计

地震波本征衰减反映了地层及其所含流体的一些特性,对油气勘探开发有重要意义.已有的理论研究与实验发现,地震频带内的衰减主要与中观尺度(波长与颗粒尺度之间)的斑状部分饱和、完全饱和岩石弹性非均匀性情况下波诱导的局部流体流有关.这种衰减与岩石骨架、孔隙度及充填流体的性质密切相关.本文着重讨论均匀流体分布、斑状或非均匀流体分布两种情况下部分饱和岩石的纵波模量差异.以经典岩石物理理论和衰减机制认识为基础,通过分析低频松弛状态、高频非松弛状态岩石的弹性模量,讨论储层参数(如孔隙度、泥质含量以及含水饱和度等)与纵波衰减之间的确定性关系.上述方法与模型在陆相砂泥岩地层与海相碳酸盐岩地层中的适用性通过常规测井资料得到了初步验证.     

Geophysical constraints on deep weathering and water storage potential in the Southern Sierra Critical Zone Observatory

The conversion of bedrock to regolith marks the inception of critical zone processes, but the factors that regulate it remain poorly understood. Although the thickness and degree of weathering of regolith are widely thought to be important regulators of the development of regolith and its water‐storage potential, the functional relationships between regolith properties and the processes that generate it remain poorly documented. This is due in part to the fact that regolith is difficult to characterize by direct observations over the broad scales needed for process‐based understanding of the critical zone. Here we use seismic refraction and resistivity imaging techniques to estimate variations in regolith thickness and porosity across a forested slope and swampy meadow in the Southern Sierra Critical Zone Observatory (SSCZO). Inferred seismic velocities and electrical resistivities image a weathering zone ranging in thickness from 10 to 35 m (average = 23 m) along one intensively studied transect. The inferred weathering zone consists of roughly equal thicknesses of saprolite (P‐velocity < 2 km s^?1) and moderately weathered bedrock (P‐velocity = 2–4 km s^?1). A minimum‐porosity model assuming dry pore space shows porosities as high as 50% near the surface, decreasing to near zero at the base of weathered rock. Physical properties of saprolite samples from hand augering and push cores are consistent with our rock physics model when variations in pore saturation are taken into account. Our results indicate that saprolite is a crucial reservoir of water, potentially storing an average of 3 m³ m^?2 of water along a forested slope in the headwaters of the SSCZO. When coupled with published erosion rates from cosmogenic nuclides, our geophysical estimates of weathering zone thickness imply regolith residence times on the order of 10⁵ years. Thus, soils at the surface today may integrate weathering over glacial–interglacial fluctuations in climate. Copyright © 2013 John Wiley & Sons, Ltd.     

τ-p SEISMIC DATA FOR VISCOELASTIC MEDIA – PART 1: MODELLING1

Viscoelastic modelling reveals that the interaction of compressional-wave velocity C_p, compressional-wave quality factor Q_p, shear-wave velocity C_s, shear-wave quality factor Q_s and Poisson's ratio as a function of time intercept τ and ray parameter p, is complicated; however, distinct, potentially diagnostic behaviours are seen for different combinations of viscoelastic parameters. Synthetic seismograms for three viscoelastic reservoir models show that variations in the Poisson's ratio produce visible differences when compared to the corresponding elastic synthetic seismograms; these differences are attributable to interaction of the elastic parameters with Q_p and Q_s. When the P-wave acoustic impedance contrast is small, viscoelastic effects become more apparent and more useful for interpretation purposes. The corresponding amplitude and net phase spectra reveal significant differences between the elastic and the viscoelastic responses. When P-wave reflectivities are large, they tend to dominate the total response and to mask the Q reflectivity effects. The attenuation effects are manifested as an amplitude decay that increases with both time and ray parameter. The sensitivity of the computed seismic responses for various combinations of viscoelastic parameters suggests the opportunity for diagnostic interpretation of τ-p seismic data. The interpretation of the viscoelastic parameters can permit a better understanding of the rock types and pore fluid distribution existing in the subsurface.