含水率对页岩动态弹性力学性质的影响研究

声波测试是研究岩石动态力学性质的一种重要的无损检测方法。在有机质丰富的页岩中,页岩的力学性质受含水量、岩石层理、总有机碳含量(TOC)和微观结构变化的影响。本研究在声波测试前,首先使用CT扫描技术对鹰滩页岩的每组测试样品进行扫描,对实验样品的完整性进行筛选,利用核磁共振法(NMR)测定样品的含水率。然后利用超声波脉冲传输技术在不同围压下对每组测试岩心进行声波测试。测试后,利用X射线衍射(XRD)分析了样品的矿物组分,利用扫描电子显微镜(SEM)对微结构进行了表征。研究结果表明:围压和含水率对纵波(P波)和横波(S波)速度均有显著影响,两种速度都随着围压的增加而增加。含水率的增加会降低横波与纵波的波速,含水率对横波速度比纵波速度的影响更大。含水率降低了杨氏模量和剪切模量,增加了体积模量、拉梅参数和泊松比。分析表明,引起岩石力学性质含水率变化的机制可能是粘土-水相互作用、化学反应和毛细压力变化的综合结果。     

致密岩石纵横波波速各向异性的比较研究

通过对板岩、千枚岩、糜棱岩和变质砂岩等4种岩石在平行和垂直层理、板理的3个正交方向上的纵、横波波速试验,提出了致密岩石纵波波速的各向异性与横波波速的各向异性普遍存在着一致性和差异性两种特征。一致性特征主要表现为致密岩石的横波各向异性随纵波各向异性的增强而增强;差异性特征主要表现为在平行和垂直层理、板理的两个方向上,纵波波速的各向异性指数大于横波,亦即纵波波速的各向异性比横波更为显著。从横观各向同性弹性介质的波速方程出发,通过具体的算例对这两种特征进行了较好的理论分析。     

基于微形态理论的各向异性散体波动分析

基于散粒体微观力学理论,忽略颗粒转动引起的相对位移,考虑颗粒接触的组构各向异性,根据宏微观能量守恒推导得到了散体材料各向异性微形态本构关系,进而通过单位接触方向积分的递推公式推导出了各向异性本构张量表达式;在此基础上,根据哈密顿原理得到了各向异性散体材料的运动平衡方程和边界条件,从而求得了平面波在各向异性散粒体中的传播规律和频散关系,最后对波的频散关系和频率带隙进行了详细的参数分析。研究表明：该模型预测了散体中包含3类12种位移波：3种纵波、6种横波和3种平面内横向剪切波;横观各向同性条件下,接触各向异性参数a20越大,纵波和横波的频率越大,而平面内横向剪切波的频率越小;正交各向异性条件下,随着接触各向异性参数a22的增大,与2方向运动相关的横波频率增大,而与3方向运动相关的横波频率则减小;但a22的变化对纵波频率影响很小。材料各向异性程度对横波带宽影响不大,但对纵波带宽影响较大：a20的增大使得声?光学波间的带宽减小,而光学波间的带宽增大,当a20>0.84时,声?光学波间的带隙消失;但是a22的增大则使得声?光学波间的带宽增大,而光学波间的带宽减小。退化为各向同性模型后,预测3类波的频散曲线与其他各向同性模型的结果基本一致。     

焦石坝地区页岩各向异性及其影响的分析研究

黑色富含有机质页岩由于沉积方式的特殊性,一般具有较强的各向异性。但现有大部分研究是基于各向同性理论开展页岩储层预测工作,这与页岩的实际地质特点不符。文中联合利用实验室岩石力学测试数据以及各向异性参数经验公式计算焦石坝地区页岩的各向异性参数,并进一步探讨了各向异性对页岩岩石速度、储层AVO响应特征以及岩石力学参数计算的影响。结果表明,各向异性对页岩速度、脆性指数预测结果影响较大,在地层倾角较大时进行速度分析及页岩可压性预测时应考虑各向异性参数的影响;同时,AVO正演结果表明地震反射入射角小于40°时,各向异性算法与各向同性算法计算的页岩AVO响应差异较小,利用各向同性叠前反演方法进行页岩储层预测具有可行性。     

利用高阶交错网格有限差分法数值模拟VTI介质井孔声场

横向各向同性（TI）介质是岩石地球物理中常见的一种现象,研究其井孔声场传播特征对声波测井理论以及为声波测井解释提供依据具有重要意义。针对具有垂直对称轴的横向各向同性（VTI）介质,根据柱坐标系条件下的弹性波波动方程,推导了速度-应力交错有限差分公式,采用时间二阶、空间十阶的交错有限差分算法对VTI介质中的井孔声场进行数值模拟。给出了在均匀介质中井孔声场不同时刻的波场快照,以及不同各向异性系数的VTI介质中的波场快照,计算了井轴上声源激发出的声波全波列波形。结果表明,在其他条件不变的条件下,VTI地层的各向异性系数的增大对横波的传播影响不大,但会使得纵波在纵向上的传播速度相对变小,径向上变化不大。各向异性系数的增大会使声波测井全波列首波信号时差变大,声波幅度略变小。      

围压下裂隙对地震波速度影响的物理实验研究

构建了等直径不同裂隙密度和等裂隙密度不同裂隙直径两组物理模型,进行不同围压条件下多方向的超声波速度测试,并运用Hudson理论进行了理论模型计算。结果显示,计算与实测结果吻合较好。随围压的增大,纵、横波速度均近线性增加,纵、横波各向异性基本保持不变;裂隙密度从2%增大到6%,纵波速度不同程度降低,其中慢纵波降低幅度相对较大,快横波变化不明显,而慢横波则大幅降低。随着裂隙密度的增大,纵、横波各向异性均增大,且横波各向异性增加速率大于纵波;裂隙直径从2 mm增大到3 mm,快纵波速度增加很小,慢纵波增加明显,横波速度均不发生改变。随着裂隙直径的增大,纵波各向异性逐渐降低,横波各向异性保持不变。最后结合试验结果分析了Hudson理论在不同深度进行参数预测的必要条件。研究结果有助于油气生产、地下水的开采与控制、污染处理等。      

泥页岩储层岩石力学特性及脆性评价

泥页岩储层的岩石力学特性对油气开发影响极大,进行泥页岩力学特性和脆性评价方面的研究,可以为泥页岩油钻井和压裂设计工作提供技术支撑。实验研究表明,泥页岩抗压强度与围压、杨氏模量成正相关;体积应变量随杨氏模量减小而增大,随泊松比增加而增加;泥页岩破坏在低围压下以劈裂式破坏为主,高围压时多出现剪切式破坏。泥页岩的脆性与其弹性参数和矿物组成关系密切,通过数值模拟和实验测量,综合弹性参数和矿物组分两种方法提出了一种新的脆性评价方法-弹性参数与矿物成分组合法(EP&MC Method),并实现了单井脆性评价,效果较好。脆性评价既是储层岩石力学特性分析的重要内容,也是压裂选层的重要依据。     

不同温压下岩石弹性波速度、衰减及各向异性与组构的关系

结合岩石组构分析 ,阐述了岩石弹性波传播速度和衰减以及它们的各向异性与岩石组构之间的关系。在不同温压条件下对具有很强晶格优选方位的岩石样品的研究表明 ,随着围压的增加 ,波速和Q值均增大 ,但是在相互正交的 3个方向上 (垂直或平行于层理面及线理方向 )增大的速度并不相同 ,这与微裂隙的逐渐闭合密切相关。观测到的波速和Q值的各向异性具有不同的形成机理 ,波速各向异性主要与定向分布的微裂隙和主要矿物的晶格优选方位等构造因素有关 ;高围压下Q值各向异性与速度各向异性正好相反 ,可能是由于定向排列的矿物晶体沿不同方向其边界之间接触程度不同造成的。对岩石组构的研究不仅可以揭示岩体的变形机制、变形的动力学过程及其有关的热力学信息 ,还可以对宏观岩石的各种物理性质 ,尤其是力学特性 ,从微观机理上加以解释。文中特别强调了岩石组构分析对研究岩石物理性质的各向异性具有十分重要的意义。     

岩石孔隙流体对纵横波速度影响的实验研究及意义

我们用ＭＴＳ的超声波测试系统研究了地表岩样饱和水及孔隙压力对纵横波速度的影响。在实验研究中,分别在常温常压条件下测量了干燥岩石和完全饱水岩石的纵横波速度,并在模拟地层条件下,测量了岩石纵横波速度随孔隙压力的变化。实验结果显示：（１）岩石完全饱和水后,与干燥时相比,纵波速度将明显增加,最大可增加３０％,且岩石在干燥时的纵波速度越低速度增大的比例越大,但纵波速度增大的比例与孔隙度无明显关系。饱水前后横波速度几乎保持不变。（２）在模拟地层的温压条件下,当孔隙压力由低到高,纵横波速度均将线性减小,且横波速度比纵波速度下降的比例更大。据此,得到了识别地层异常高压和地层含气的地震速度标志。     

湘西桑页1井页岩弹性波速实验研究

随着页岩气勘探开发工作的不断深入,地震相关勘探技术在页岩气钻井部署和页岩气富集带或"甜点"区评价中发挥重要作用。不同应力环境下岩石的弹性波速及其各向异性实验数据是地震资料叠前反演的重要约束性基础资料,对深部页岩气地球物理探测、储层评价及甜点预测具有重要意义。本文以湘西桑页1井龙马溪组页岩、五峰组页岩及底板泥灰岩钻孔岩芯为研究对象,实验测量了不同应力条件下岩芯的弹性波速特征,计算了应力环境下的波速差异,并对比了页岩与泥灰岩的波速各向异性。结果表明,在5~100 MPa应力环境下,粉砂质、硅质和含粉砂页岩的纵波速度分别为4.8~5.1 km/s、5.1~5.3 km/s和5.2~5.5 km/s,泥灰岩的纵波速度在5.9~6.4 km/s之间;原位应力环境下,页岩与泥灰岩的纵波速度差异介于14.9%~21.7%之间;受微裂缝和层理面的影响,页岩垂直层理方向纵波速度明显低于平行层理方向,而泥灰岩的纵波速度为平行层理方向明显高于垂直层理方向。     

Generalized solution to the anisotropic Mandel's problem

Existing solutions to Mandel's problem focus on isotropic, transversely isotropic, and orthotropic materials, the last two of which have one of the material symmetry axes coincide with the vertical loading direction. The classical plane strain condition holds for all these cases. In this work, analytical solution to Mandel's problem with the most general matrix anisotropy is presented. This newly derived analytical solution for fully anisotropic materials has all the three nonzero shear strains. Warping occurs in the cross sections, and a generalized plane strain condition is fulfilled. This solution can be applied to transversely isotropic and orthotropic materials whose material symmetry axes are not aligned with the vertical loading direction. It is the first analytical poroelastic solution considering mechanical general anisotropy of elasticity. The solution captures the effects of material anisotropy and the deviation of the material symmetry axes from the vertical loading direction on the responses of pore pressure, stress, strain, and displacement. It can be used to match, calibrate, and simulate experimental results to estimate anisotropic poromechanical parameters. This generalized solution is capable of reproducing the existing solutions as special cases. As an application, the solution is used to study the responses of transversely isotropic and orthotropic materials whose symmetry axes are not aligned with the vertical loading direction. Examples on anisotropic shale rocks show that the effects of material anisotropy are significant. Mandel-Cryer's effects are highly impacted by the degree of material anisotropy and the deviation of the material symmetry axes from the vertical loading direction.     

A modi?ed failure criterion for transversely isotropic rocks

A modified failure criterion is proposed to determine the strength of transversely isotropic rocks. Me-chanical properties of some metamorphic and sedimentary rocks including gneiss, slate, marble, schist, shale, sandstone and limestone, which show transversely isotropic behavior, were taken into consider-ation. Afterward, introduced triaxial rock strength criterion was modified for transversely isotropic rocks. Through modification process an index was obtained that can be considered as a strength reduction parameter due to rock strength anisotropy. Comparison of the parameter with previous anisotropy in-dexes in literature showed reasonable results for the studied rock samples. The modified criterion was compared to modified Hoek-Brown and Ramamurthy criteria for different transversely isotropic rocks. It can be concluded that the modified failure criterion proposed in this study can be used for predicting the strength of transversely isotropic rocks.     

页岩各向异性特征的试验研究

为研究彭水页岩气区块储层的各向异性特征,开展了石柱县龙马溪组页岩的单轴和三轴压缩试验,分析了其力学特性、强度特征和破裂模式的各向异性,并揭示了其破坏机制的各向异性。结果表明：(1) 龙马溪组页岩具有明显的各向异性特征,弹性模量在平行层理方向最大,垂直层理方向最小,且围压的增加使其增加速率不断减小;0°、30°和60°、90°页岩的泊松比随围压的增加呈现出了相反的变化规律,这可能与页岩层理间孔隙和微裂缝的良好发育有关。(2) 相同围压下,0°试样强度最高,90°次之,30°最低,总体上呈现出两边高、中间低的U型变化规律,而不同角度的Hoek-Brown强度准则能较好地反映其强度的各向异性特征。(3) 页岩破裂模式的各向异性是由破坏机制的各向异性引起的,而强度的各向异性是由破坏机制的各向异性控制的。单轴压缩时,0°页岩为沿层理的张拉劈裂破坏,30°为沿层理的剪切滑移破坏,60°为贯穿层理和沿层理的复合剪切破坏,90°为贯穿层理的张拉破坏。三轴压缩时,0°为贯穿层理的共轭剪切破坏,30°为沿层理的剪切滑移破坏,60°和90°为贯穿层理的剪切破坏;页岩地层的层状沉积结构和层理间的弱胶结作用是破坏机制各向异性的根源。研究结果为水平井井壁的稳定性分析和水力压裂施工设计等提供了技术参考。     

Porothermoelastic analyses of anisotropic hollow cylinders with applications

It has been known that material anisotropy and thermal stresses affect borehole stability significantly. Aiming at the experimental studies associated with borehole stability in anisotropic (transversely isotropic) poroelastic materials subject to non‐isothermal conditions, this paper details and applies an anisotropic porothermoelastic solution to an unjacketed hollow cylinder in a triaxial set‐up. Numerical analyses are presented to demonstrate thermal and material anisotropy effects on the pore pressure and the stress concentrations in and around the geometry of a hollow cylinder subjected to thermal and stress perturbations. Copyright © 2004 John Wiley & Sons, Ltd.     

不同固结条件下尾矿动孔压演化规律

地震作用下尾矿坝有液化失稳的危险,尾矿动孔压的演化规律可以间接体现其液化过程。为研究尾矿动孔压的演化规律,开展了一系列的动三轴试验。结果表明：尾矿动孔压演化具有明显的阶段性,在等压固结时可用S形反正弦函数曲线描述,在偏压固结时可用J形指数函数曲线描述;尾矿材料在等压固结条件下临界孔压接近于围压,在偏压固结条件下临界孔压小于围压,且随围压及固结应力比的增大而减小,随尾矿平均粒径的增大而增大。基于极限平衡理论推导了等压固结和偏压固结条件下的临界孔隙水压力方程,从理论上阐释了试验结果所揭示的动孔压演化规律,可以为地震区上游式尾矿坝的抗震设计提供参考。     

Experimental study of the seismic properties of the Eastern Alps (Italy) along the Aurina–Tures–Badia Valleys transect

In order to better constrain the interpretation and the nature of the seismic reflectors, experimental measurements at high confining pressure (up to 300 MPa) and room temperature of the compressional wave velocity (Vp) on 10 samples representative of the most common lithologies along the Aurina (Ahrntal), Tures (Tauferer Tal), and Badia (Abtei Tal) Valleys profile (Eastern Alps, Italy) have been performed. For each sample, the speed of ultrasonic waves was measured in three mutually perpendicular directions, parallel and normal to the rock foliation and lineation.The main results are:(a) Good agreement between the calculated vs. measured modal compositions of the considered rocks, indicating that they were presumably equilibrated at the estimated P–T conditions; therefore, the seismic properties are representative of the crustal level indicated by the thermobarometry.(b) Measured and calculated average Vp are in good agreement, and are typical of mid-crustal level (6.0–6.5 km/s). Only the amphibolites show Vp typical of the lower crust (7.2 km/s).(c) The seismic anisotropy of metapelites is very high (12–27%), both with orthorhombic and transverse isotropy symmetry; amphibolites are transversely isotropic with an anisotropy of 8%; orthogneisses and granitoids are isotropic or weakly anisotropic.(d) The contacts between amphibolites and all other rock types may generate good reflections, provided they are not steeply inclined. Although the metamorphic foliation remains steeply inclined, discordant buried sub-horizontal igneous contacts may be detected.     

横观各向同性水敏性地层斜井眼坍塌压力确定

在以泥页岩为目标储层的非常规油气藏水平井开发中, 钻遇井壁的稳定性问题突出, 一方面是由于在设计钻井液密度时简单采用各向同性介质模型和Mohr-Coulomb强度准则, 另一方面则是由于选用的钻井液性能欠佳, 难以保证建井全过程的井壁稳定。为此, 将地层视为横观各向同性介质, 采用Mogi-Coulomb强度准则判别岩石基体和层理面的稳定性, 同时考虑泥页岩水化作用对岩石力学强度参数的影响, 建立了水敏性地层井壁稳定分析模型, 分析了坍塌压力随井斜、方位、钻井液性能、钻井时间等因素的变化规律。研究结果表明:对于具有显著层理面结构的泥页岩地层, 采用横观各向同性模型比各向同性介质模型更能描述强度的各向异性;采用Mogi-Coulomb准则比Mohr-Coulomb强度准则得到的坍塌压力值与实际吻合程度更高, 在保证井壁稳定的前提下有利于降低钻井液密度, 实现钻井提速;提高钻井液封堵能力、减小泥页岩吸水扩散能力, 有利于延长坍塌周期。利用东濮、威远、焦石坝等工区多口页岩气水平井的实钻资料对本文方法进行了验证, 证实本文方法确定的钻井液密度窗口下值更能满足安全钻井需求。     

Bayesian inversion of synthetic AVO data to assess fluid and shale content in sand-shale media

Reservoir characterization of sand-shale sequences has always challenged geoscientists due to the presence of anisotropy in the form of shale lenses or shale layers. Water saturation and volume of shale are among the fundamental reservoir properties of interest for sand-shale intervals, and relate to the amount of fluid content and accumulating potentials of such media. This paper suggests an integrated workflow using synthetic data for the characterization of shaley-sand media based on anisotropic rock physics (T-matrix approximation) and seismic reflectivity modelling. A Bayesian inversion scheme for estimating reservoir parameters from amplitude vs. offset (AVO) data was used to obtain the information about uncertainties as well as their most likely values. The results from our workflow give reliable estimates of water saturation from AVO data at small uncertainties, provided background sand porosity values and isotropic overburden properties are known. For volume of shale, the proposed workflow provides reasonable estimates even when larger uncertainties are present in AVO data.     

Anisotropic mechanical properties and brittleness evaluation of layered phylliteEI北大核心

为了研究层状千枚岩的力学特性与各向异性特征,开展了不同层理倾角与不同围压下的千枚岩力学试验。对比分析了千枚岩试样强度、变形、脆性与破坏模式等各向异性特征。结果表明:(1)随着层理倾角β增加,岩样的强度、变形特征曲线形状呈U型;随围压增加,岩样强度及塑性增强,各向异性度逐渐减弱稳定。(2)采用多种强度准则描述岩样强度各向异性,其中Saeidi准则和改进Ramamurthy准则能很好地预测岩样在不同层理倾角下的抗压强度。(3)基于岩样峰前应力-应变曲线与能量特征提出了综合脆性评价指标,在层理倾角β=45°左右时,岩样脆性指标较低,更易发生剪切滑移破坏,得出脆性下降顺序为:沿层理面拉伸劈裂>穿层理面拉伸劈裂>沿层理面剪切>穿层理面剪切。(4)千枚岩的破坏模式与层理倾角和围压有关,单轴条件下,岩样劈裂破坏后易形成复杂裂纹网络;高围压下,岩样破裂后多形成单一的沿层理面或贯穿多层理面的剪切破坏。