页岩气勘探开发中的地震技术

1、页岩气勘探开发中的地震技术页岩气作为非常规油气资源，与常规油气资源相比，一般情况下，页岩气的富集并不形成类似于常规油气的圈闭，存在局部富集的差异，没有明显的边界，属于“自生自储”型的低丰度连续的油气藏，在勘探开发过程中，有着不同的评价标准、工作流程和技术流程。因此，需要结合页岩气勘探开发的实际情况，认真审视常规油气勘探开发中发挥巨大作用的地球物理技术(尤其是地震技术)的定位，分析地震技术在页岩气勘探开发中的作用。岩气储集在厚层的泥页岩中(中间可以有薄夹层),解释、识别这样的地层相对比较容易,从勘探到开发都离不开地震勘探技术。以地震技术为主体的气藏描述是页岩气储层识别与评价的核心，具体地震勘探任

2、务包括查明页岩层的深度、厚度、分布范围、产状形态，寻找页岩层内有机质丰度高、裂缝发育、渗透性好、脆性大的部位，即页岩气“甜点区”。微地震监测是一种主要用于油气田开发的新的地震方法，主要包括井中微地震监测和地面微地震监测，主要用于在水力压裂作业过程中，了解裂缝的走向和评价压裂的效果，对诱导裂缝的方位、几何形态进行监测。1 页岩构造解释与优选1.1 页岩层位解释利用地震、地质以及测井资料，准确标定页岩层的顶底界面，分析目的层页岩层段的地震波响应特征在地震剖面上识别和追踪页岩储层，确定页岩层的深度、厚度以及分布范围，获得页岩层构造形态、断层展布以及沉积厚度特征等。结合区域构造及沉积环境的研究，单井地

3、质剖面测井解释以及地震属性分析，划分页岩层段沉积微相类型和沉积微相变化规律，确定有利的页岩层沉积微相发育特征以及展布范围。1.2 优质页岩预测图1 龙马溪组页岩层自然伽马反演剖面图（据李志荣，2011）在查明了页岩层的构造与沉积特征后，需要在页岩的构造-层序格架内寻找优质页岩发育区域，即圈定页岩有机质丰度、成熟度高的位置。通过页岩实验测试数据与测井分析，建立测井-岩电关系模板，在分析优质页岩储层与地震反射波响应特征的基础上，通过井震联合反演确定优质页岩发育位置。例如，优质页岩与普通页岩相比具有有机质含量高、自然伽马数值高的特点，据此进行井震联合反演，能够对优质页岩层进行很好的预测（图1）。2

4、页岩含气性检测页岩储层的含气量决定了页岩气开发是否具有商业价值，所以对页岩含气性检测直接决定了页岩气的勘探开发价值。目前对页岩储层进行含气性检测的地震技术主要有叠后波阻抗反演、叠前AVO反演、地震属性分析技术和叠前弹性阻抗反演以及频谱分解技术。2.1 叠后波阻抗反演叠后反演的基础是褶积模型，即地震数据可以看作地震子波与反射系数的褶积。通过压缩子波的反褶积处理，将地震数据转换为近似的反射系数序列，然后再由反射系数序列得到波阻抗剖面。随着页岩层含气量的增大，储层体积密度和层速度会降低，从而导致波阻抗值减小，所以在页岩层的地质模型约束下拾取页岩层波阻抗数据，其波阻抗低值区代表低密度或低速区（图2），

5、也是预测的储图2 龙马溪组页岩气层波阻抗反演剖面层含气区。2.2 叠前AVO反演AVO反演所依据岩石物理学理论和振幅随偏移距变化理论,借助于Zoeppritz方程或近似式,对CDP道集反射振幅的变化作最小平方拟合,直至理论值与观测值很好的拟合,利用振幅随炮检距变化关系曲线计算出截距P和梯度G两个参数,再通过这两个参数反演出所需的弹性参数,进而进行岩性、流体识别.由于页岩气富集，导致储层体积密度减小。弹性波速度降低，对含气性检测参数具有明显的影响。美国Eagle Ford)页岩气层通过叠前AVO技术,反演出与储层孔隙流体性质相关的参数,由于含气储层具有低值,据此预测页岩储层含气有利区(图3).2

6、.3 叠前弹性阻抗反演Patrick Connolly将纵波反射系数随炮检距变化的思想引入叠后地震道的正-反演问题中，率先提出了弹性阻抗(Ze)为纵波速度、横波速度、密度和入射角的函数。图3 Eagle Ford 页岩气层属性图（据Galen Treadgold,2010）弹性阻抗函数是对声波阻抗概念的推广，它是入射角的函数,声波阻抗是弹性阻抗入射角为00时的特例,其不仅具有叠后波阻抗反演的优点，而且还弥补了叠前AVO反演技术稳定性和分辨率较低的不足,同时弹性阻抗较波阻抗包含更多的岩性和物性信息,增强了反演技术预测和描述储层的能力.2.4 频谱分解技术频谱分解技术是一项基于频率谱分解的储层特色

7、解释技术，主要依据是含油气储层的高频吸收特性，即当地震波经过含油气储层时，其高频成分能量衰减较地震波通过不含油气储层时严重。目前该技术在常规储层流体检测中已初见端倪，但在页岩气检测中应用较少。在频谱分解技术与常规AVO反演技术基础上，综合两种方法各自的优势，产生了分频AVO技术与频变AVO技术。分频AVO技术，是采用分频技术，先对叠前地震进行分频处理，然后在分频数据的基础上实现AVO油气检测。频变AVO技术是基于AVO的Zoeppritz方程，建立反射系数与频率之间的数学关系，推导出截距、梯度、碳氢检测因子等属性与频率之间的数学关系，综合地质、地震、测井等数据，反演出高精度的频变AVO属性，检

8、测页岩气储层。程冰洁（2012）等人在川西新场地区3D3C（三维三分量）叠前地震资料进行频变AVO处理并获得了流体检测因子和拟泊松比两种频变AVO属性，并对含气性进行了预测解释。图4中，流体检测因子的值在彩色区域较大，在空白区域较小，这与钻井吻合较好，这说明流体检测因子越大，含气的可能性越大；图5中，高产井都具有较强的频变拟泊松比异常区（红色区域），低产井区并未出现出现异常（程冰洁，2012）。图4 新场Tx24层段流体10HZ检测因子含气性识别平面图(据陈冰洁等，2012)图5 新场Tx24层段流体10HZ才与15HZ频变拟泊松比含气性识别平面图2.5 地震属性分析技术地震属性分析就是利用属

9、性提取手段，从地震资料中提取隐藏的信息，并将其转换成与岩性、物性相关的参数，可以为地质解释或油藏工程直接服务，从而达到充分发挥地震资料潜力，提高地震资料进行储层预测和表征的能力。地震属性优化及其综合解释技术，全方位地利用地震资料，将多属性提取与分析融为一体，这在一定程度上克服了地震解释的多解性，能更好地为地质综合研究服务。3 页岩总有机碳含量（TOC）预测页岩的总有机碳(TOC)含量(以下称TOC质量分数)是评价页岩生烃能力最重要的指标参数。页岩是一种有机质丰度很高的烃源岩，因而对页岩TOC质量分数的预测可参照烃源岩的评价方法来进行。对页岩的总有机碳评价对油气的勘探开发有很大意义，可以预测页岩

10、的生烃潜力。除实验室检测方法外，对烃源岩TOC质量分数的评价主要以测井资料为主，且评价方法很多，其中采用电阻率-孔隙度重叠法进行单井TOC质量分数评价1,2最为广泛。对于页岩TOC质量分数的横向分布预测，研究发现其可通过地震多属性反演方法来实现。地震多属性反演就是在三维地震资料和钻井、录井、测井等资料解释的储层物性、含气性和岩石地化指标等目标参数的基础上，经过构造解释拾取层位，提取和优化地震属性，建立地震属性或组合与已知目标参数之间的统计关系来反演目标参数分布的技术。金吉能（2012）等通过地震多属性反演建立起了多属性组合和总有机碳含量的相关关系，且相关系数很高，并最终得到了页岩TOC质量分数

11、预测平面图（图6），通过井震结合对比分析验证了这个方法的可行性，相关系数很高（图7）。图6 某区块的TOC质量分数平面分布预测图7 各井点处的TOC质量分数实测值和预4 页岩裂缝预测相对于常规油气藏, 页岩气的成因类型、富集机理及生产机制等都具有一定的特殊性, 页岩气产量高低直接与泥页岩内部天然微裂缝发育程度有关,这说明微裂缝的存在某种程度上提高了水力压裂效应的有效性, 从而极大改善了泥页岩的渗流能力,为页岩气从基岩孔隙进入井孔提供了必要的运移通道 3-6 。同时, 泥页岩中的小孔洞、微裂缝和纳米级微孔隙也是页岩气的重要聚集空间, 除了以吸附状态赋存于岩石颗粒和有机质表面外,还有大部分页岩气以

12、游离状态赋存于微小的孔隙之中3,7-8.而页岩气产量高低直接与泥页岩内部天然微裂缝发育程度有关,裂缝的存在某种程度上提高了水力压裂效应的有效性,从而极大改善了泥页岩的渗流能力，故裂缝预测对页岩气的勘探开发具有很重要的指导意义。目前应用的地震裂缝预测技术方法主要有三大类，即转换波裂缝检测、纵波方位各向异性检测和叠后地震属性分析9。4.1 多分量转换波检测裂缝20世纪80年代以来，随着各向异性研究的深入，特别是横波分裂现象的发现和多分量地震检波器的研制成功，使得应用多分量转换波研究裂缝成为可能。转换波勘探克服了纯横波勘探激发难、成本高、静校正量大等缺陷，同时具有所得信息多、兼有纵波和横波的优势、资

13、料信噪比较高、频带较宽、勘探深度较大等优点，使得多分量转换波勘探成为油气储层探测的有力工具。从多分量转换波资料中确定裂缝的方向和密度也成为地震勘探的重要目标之一21。Ata等22给出了在委内瑞拉通过三分量地震进行裂缝检测的应用实例；Van Dok等23、Gaiser等24-25给出了利用转换波进行裂缝表征的实例；Grechka等26-27探讨了在利用多方位、多分量地震资料进行裂缝参数预测时不同裂缝等效介质模型的适用性；Vasconcelos等28、Engelder等29、Bale等30、Roxana等31、Tang等32和蔡希源等33也都给出了利用转换波数据进行裂缝预测的实例(图8)。图8 横

14、波分裂时差和纵波相干体预测的裂缝分布(据蔡希源等，2010)在处理多分量资料时，判别快、慢横波是否完全分离是多分量裂缝预测的关键。其中比较有代表性的方法包括最小熵旋转法34、正交基旋转法35及全局寻优法36。Dumitru等34提出的最小熵旋转法，无需互相关法中快、慢波波型相似这一假设条件，可以对快、慢横波的方向作出评估。用最小熵旋转法计算的裂缝走向较为准确，尤其是快、慢波时差较大的情况下。然而，如果快、慢波时差较小，该方法受干涉影响，精度降低。黄中玉等35提出正交基旋转方法基于快、慢横波为正交偏振这一假设条件，推导出野外采集系统坐标与自然坐标之间夹角的解析关系式，从而实现快、慢横波分离，估算

15、地层裂缝发育方位，实现裂缝检测的目的。理论模型测试结果表明，正交基旋转法具有较高的准确性与可靠性，但受资料信噪比的影响较大。如果地下有多层各向异性介质，且方位角各不相同，则快、慢波分离过程要复杂得多。Dariu36把全局寻优法中的模拟退火技术引入到裂缝探测工作当中，开发了一种自动逐层估算极化方向和快、慢波时间延迟技术。转换波裂缝检测常用的方法有相对时差梯度法和层剥离法32-33，37-39。相对时差梯度法是一种在数据体上计算裂缝发育方位和密度的方法，其步骤为：首先分时窗扫描裂缝发育方向；然后将各方位的径向分量和横向分量数据旋转到裂缝方向上得到快、慢波数据体；最后计算快、慢波的时差数据体，并计算

16、其梯度以获得反映裂缝发育密度的相对时差梯度数据体，用以检测裂缝在空间上的发育情况。层剥离法是一种沿层位进行检测的方法，其步骤为：首先沿目的层顶面计算裂缝发育方位和密度，并通过时间补偿和旋转分析，消除上覆地层各向异性的影响；然后在目的层底界面分析裂缝发育方位和密度，以获得某一特定勘探目的层的裂缝发育方位和密度。4.2 纵波方位各向异性检测裂缝尽管裂缝诱导的方位各向异性对纵波、横波和转换波的传播特征有着重要的影响，但早期研究关注的重点是近垂直入射情况下分裂横波的时间延迟及反射振幅的变化。虽然此类研究能给出垂直裂缝系统的方位和强度（裂隙密度）信息，但横波类的检测方法对裂缝充填物的响应并不敏感。随着研

17、究的深入，研究者逐渐认识到纵波信号所携带的与方位相关的变化特征不仅可用于解决裂缝的方位、密度问题，而且对了解裂缝充填状况有所帮助15。纵波检测裂缝方法由于费用和资料品质等方面较之于横波、转换波更具优势而日益受到重视。应该说，纵波裂缝检测方法具有非常好的发展前景，但对由地震得到的各向异性参数与裂缝系统参数之间关系的充分研究是相关方法成功应用的先决条件。目前利用纵波各向异性进行裂缝检测的方法有动校正（NMO）速度方位变化裂缝检测、正交地震测线纵波时差裂缝检测、纵波方位AVO（AVOZ/AVAZ）和纵波阻抗随方位角变化（IPVA）裂缝检测方法。在各向同性介质中，反射波的时距曲线是抛物线形式，在进行叠

18、加之前，需要进行动校正(NMO)，这时的叠加速度(NMO速度)是一个常量。在裂缝诱导的EDA介质中，由于各向异性的存在，不同方向的NMO速度随EDA介质对称轴和地震测线之间的夹角不同而改变40-43.。此时，NMO速度在水平面内的轨迹是一个椭圆（HTI介质情况下），椭圆的长轴平行于裂缝的走向，短轴平行于HTI介质的对称轴（即垂直于裂缝走向）。据此对三维地震数据沿不同方向的CMP道集进行精细速度分析，求取NMO速度（即叠加速度），并在平面上模拟NMO速度的椭圆轨迹，椭圆的长轴指向即为裂缝发育方向。正交地震测线纵波时差裂缝检测的原理和动校正(NMO)速度方位变化裂缝检测方法类似，也是依据速度的各向

19、异性变化特征，但着眼点是因方位速度差异引起的纵波方位时差响应。对于固定的炮检距而言，垂直排列裂缝发育目的层（HTI）的层间纵波方位时差(Azimuthal Moveout Rseponse,AMR)是地震测线和裂缝走向之间倍夹角、偏移距以及各向异性系数的函数43-45。通过对几组正交测线的层间AMR进行计算分析，可以推断裂缝走向，进而反演出各向异性参数43-47。纵波方位AVO（AVOZ/AVAZ）裂缝检测方法相对于基于速度方位异常的纵波检测方法来说，在一定程度上减少了上覆地层的影响，能够提供关于裂缝及其充填流体的更为准确、详尽的信息16,19,48（图9）。在各向同性介质中，由于泊松比的变化

20、造成纵波反射系数随入射角变化13。在裂缝诱导的各向异性介质中，弹性参数如泊松比的变化只是引起方位AVO响应特征的一个因素。尽管纵波方位AVO响应在裂缝检测中是非常重要的地震响应特征，但在实际应用中，由于岩石弹性参数及流体存在的影响，有时由于裂缝各向异性造成的方位AVO特征并不是特别明显19。即便观测到了较为清晰的方位AVO响应特征，也需要结合NMO速度的椭圆轨迹分析和相应的流体替换模型来进行进一步的研究49-50.图9 AVO斜率G方位性差异值T6反射层层位切片(刘朋波等，2008)红色反映北东向裂缝，蓝色反映北西向裂缝由于AVOZ方法受噪声影响较大而稳定性差51，速度随方位角变化分析方法（V

21、VA）只能识别大套储层，对于薄层分辨率不够。在AVOZ方法基础上出现了波阻抗随方位角变化来预测裂缝型储层的方法（IPVA）52-53。将该方法应用于罗家裂缝发育地区的全方位三维P波地震资料分析中，分析结果表明，与AVOZ和VVA方法相比，IPVA方法具有更好的稳定性和分辨率。该方法可以作为纵波方位各向异性裂缝检测方法的重要补充。4.3 叠后地震属性分析裂缝叠后三维地震属性分析是裂缝识别与预测的方法之一。在已有研究中，振幅类、频率类、相位类地震属性以及一些分析技术（如地震波形分类、时频分析、沿层切片等）已被广泛应用于识别和预测裂缝发育带54。由于裂缝形态的特殊性，以上地震属性及分析技术更适用于推

22、断裂缝发育区的概貌，而较为精细的裂缝地震属性分析主要围绕地震反射波形的突变（不连续性）来开展，主要有曲率分析、相干分析和频谱分解等技术(图10)。图10 沿着不同属性提取的Woodford层切片（据郭艳霞，2011）（a）相干属性，（b）波阻抗，（c）显示的是（a）和（b）的叠合图，注意到断层处一般表现出低波阻抗特征（黄色、绿色）；（d）构造负曲率，（e）波阻抗能量负曲率。（f）显示的是（d）和（e）的叠合图从几何地震学的角度看，反射点集合可以视为一个时间标量场，该标量场某一反射面的梯度反映的是该反射面的起伏变化率，即单位反射时间内反射面沿不同方向的变化增量，表示的是反射曲面沿方向矢量所在法截

23、面截取曲线的一阶导数视倾角的大小；而该方向上的曲率定义为该曲线上密切圆半径的倒数，亦即为该方向上该曲线的二阶导数。由此可见，看似复杂的地震几何属性系列不过是沿不同方向计算的一阶、二阶导数体。很多研究者，如Roberts55，Hart56和Sigismondi等57发表了运用地震反射曲率属性预测裂缝的研究实例。在诸多曲率属性之中，最大正曲率和最小负曲率对裂缝识别最有价值。近几年较为突出的进展是Marfurt58，Chopra59和Al-Dossary等60在三维曲率体计算、构造倾角滤波、多尺度曲率分析等方面的研究成果。相干分析技术通过对地震波形纵向和横向相似性的判别，得到地震相干性的估计值，相似

24、地震道具有较高的相干系数，而不连续性强的地方具有较低的相干系数。Bahorich61等于1995年介绍了第一代相干算法，该算法建立在归一化互相关基础上，效率高，但抗噪性能差。Marfurt等62在1998年给出了基于多道相似性的第二代相干算法，该算法提高了抗噪能力和计算结果的垂向分辨率，但是计算道数的增加降低了侧向分辨率60。1999年，Gersztenkorn等63提出了基于本征结构的第三代相干算法。与第一代、第二代算法相比，第三代相干算法在抗噪性能、纵横向分辨率等方面有了明显改善，但计算耗时较大。此后又有一些新的改进算法，如Randen等64的几何结构张量方法Cohen等65局部结构熵估算

25、等。目前，相干分析在断裂解释及裂缝预测中的应用已经非常普遍，在裂缝预测中，相干分析结果多和地震几何属性倾角/方位角叠合显示(图11)。图11 YC凹陷某亚段相干属性裂缝预测图(据陆明华等，2013)频谱分解提供了一种在地震带宽范围内从不同离散频率尺度对地下地质特征进行解释的手段。该技术被广泛应用于沉积旋回性分析66、薄层厚度预测67、地震古地貌分析68、烃类检测69等方面。在应用过程中，出现了短时傅里叶变换、小波变换等多种频谱分解算法。频谱分解应用于裂缝检测的基础是地震反射不连续性特征相应于地质异常具有多尺度性，通过生成不同频率数据体，利用纵横向上时频点或时频段上的频谱差、频谱比、频谱下降率等

26、描述不同尺度的地震波衰减特征，揭示裂缝发育带，乃至对其含油气性进行检测(图12)。除此之外，频谱分解也常与曲率分析、相干分析配合使用，作为研究特定频率成分不连续性特征的一种数据预处理手段54.与裂缝识别相关的叠后地震属性分析是近年来业界的一个研究热点。按照Bob Peebler（2008）的说法，适应于石油工业的需要，地震勘探可划分为1970年之前的二维地震、2000年之前的三维地震、2000年之后特别是2007年之后的全波地震3个时代70。今天，相对于全波数据，三维纵波勘探已有的庞大的数据量以及较为完备的技术基础，使得三维叠后属性分析手段在裂缝预测方面仍然具有较为广阔的发展空间，这也是200

27、0年之后，倾角/方位角/曲率、相干分析、地震结构/纹理属性（能量、熵、对比度、均质性）等地震几何类属性不断呈现的原因所在71-72。图12 川南威远构造上的W112井裂缝及油气检测(据李煜伟等，2013)5 页岩气开发井中地震技术井中地震技术是在地面地震技术基础上向“高分辨率、高信噪比、高保真”发展的一种地球物理手段,在油气勘探开发中,可将钻、测井和地震技术很好的结合起来,成为有机联系钻测井资料和地面地震资料对储层进行综合解释的有效途径。5.1 微地震监测技术微地震监测是一种用于油气田开发的新地震方法，该方法优于利用测井方法监测压裂裂缝效果,在压裂施工中,可在邻井(或在增产压裂措施井中)布置井

28、下地震检波器,也可在地面布设常规地震检波器,监测压裂过程中地下岩石破裂所产生的微地震事件,记录在压裂期间由岩石剪切造成的微地震或声波传播情况,通过处理微地震数据确定压裂效果,实时提供压裂施工过程中所产生的裂缝位置、裂缝图13 地面接收信号微地震监测原理图（据贾利春等，2012）方位、裂缝大小(长度、宽度和高度)、裂缝复杂程度,评价增产方案的有效性,并优化页岩气藏多级改造的方案。图13是地面接收信号的井下微地震裂缝监测工作原理图，压裂井和监测井位于同一井区，压裂井压裂施工过程中产生的微地震信号通过地层向周围传播，位于邻井中的接收器接收这些信号并传至地面数据采集器，处理后可得到微地震监测图。图14

29、是在以邻井布置井下地震检波器为列展示了微地震监测压裂裂缝的微地震事件。从图15可以看出微地震活动性表征的复杂裂缝系统显示,裂缝模式随时间推移而扩展。图14 邻井接收信号微地震监测原理图（据罗蓉等，2011）图15 微地震监测压裂裂缝的微地震事件图（据weatherford公司）5.2 其他井中地震技术由Zero-Offset VSP、Offset VSP、Walkaway VSP、Walkaround VSP逐步发展到3D VSP技术，都是较为成熟的井中地震技术。其中,3D VSP技术和微地震采集配套施工配合监测储层改造人工裂缝发育分布状况是国外石油大公司的通常做法。3D VSP观测是一种可靠

30、的识别裂隙方向和裂隙密度分布的方法,3D VSP P-P和P-S成像用于陆上构造解释,可大大改善纵、横向分辨率和断裂系统分辨率。3D VSP测井与地面地震结合体现了综合地震勘探能力。此外,四维地震可用于检测在生产过程中,随着温度压力变化页岩气(游离气及吸附气)的变化情况,以助页岩气开发优化开采。井驱动地震数据处理是一种提高地震数据处理水平和质量的手段,也是发展趋势,使用这种技术配套,需要提高地震资料处理技术人员的整体水平。6 结语全国页岩气资源调查显示我国页岩气资源潜力巨大，但其勘探开发是一个全新课题，目前研究程度低#勘探开发技术较薄弱，特别是页岩气地震技术仍处在探索阶段。针对页岩气地质特点，

31、综合运用地震属性、地震反演、频谱分解、各向异性等技术，对页岩储层分布、厚度、含气性、裂缝发育、力学特征以及压裂缝监测等方面进行研究，进而寻找页岩气富集“甜点区”将是页岩气地震勘探技术的主要研究方向。参考文献1 许晓宏，黄海平，卢松年测井资料与烃源岩有机碳含量的定量关系研究J江汉石油学院学报,1998,20(3):8122 Passey Q R,许春艳，吴红文用测井孔隙度和电阻率求取有机质丰度的模型J国外测井技术,2002,17(5):43463 金吉能，潘仁芳，王鹏，黎建.地震多属性反演预测页岩总有机碳含量J.石油天然气学报,2012,34(11):68-72.4 李志荣,邓小江,杨晓,巫芙蓉

32、,刘定锦,张红,谭荣彪,周跃宗.四川盆地南部页岩气地震勘探新进展J.天然气工业,2011:40-43.5 王延光，杜启振.泥岩裂缝性储层地震勘探方法初探J.地球物理学进展,2006,21(2)：494-501.6 程冰洁,徐天吉,李曙光.频变AVO含气性识别技术研究与应用J.地球物理学报,2012,55(2):608-613.7 Curtis J B. Fractur ed shale-gas systems J . AAPG Bulletin , 2002, 86( 11) : 1921-1938.8 Mont gomery S L, Jarvie D M, Bowker K A, etal

33、 . Miss issipp-ian Barnett shale, Fort Worth basin, north-central T ex as: Gasshaleplay with multi trill ion cubic foot potential J . AAPGBullet in, 2005, 89( 2) : 155-175.9 Bow ker K A. Barnett shale gas production, Fort Worth basin :Issues and discussionJ . AAPG Bulletin, 2007, 91( 4) :523-533.10

34、Wang Xiang, Liu Yuhu a, Zhang Min, et al . Conditions of formation and accumulation for shale gas J . Natu ral Gas Geoscience, 2010, 21(2) : 350-356. 王祥, 刘玉华, 张敏, 等. 页岩气形成条件及成藏影响因素研究 J . 天然气地球科学, 2010, 21( 2) : 350-356. 11 Zh ang Jin chuan, Xu e Hui, Zhang Deming, etal . Shale gas an dits accu mulat

35、ionm echanism J . Geoscience, 2003, 17( 4) : 466. 张金川, 薛会, 张德明, 等. 页岩气及其成藏机理 J . 现代地质, 2003, 17( 4) : 466. 12 Zhang Jin chuan, Jin Zhijun, Yuan M ingsheng. Reser voirin gmechanism of shale gas and its distribut ion J . Natural GasIndustry, 2004, 24( 7) : 15- 18. 张金川, 金之钧, 袁明生. 页岩气成藏机理和分布 J . 天然气工业,

36、2004, 24( 7) : 15-18. 13 刘振峰,曲寿利,孙建国,董宁.地震裂缝预测技术研究进展J.石油物探,2012,51(2):193-198.14张公社等. 利用全三维纵波资料进行裂缝检测. 石油地球物理勘探,2004 ,39 (1) :414415 Schoenberg M A ,Dean S ,Sayers. Azimuth_dependent turning of seismic waves reflected f rom f ractured res2 ervoirs. Geop hysics , 1999 ,64 (4) :1160117116 Gray D , Hea

37、d Kim. Fracture detection in Manderson Field :A 32D AVAZ case history. The L eading Edge ,2000 ,19 (11) :1214122117 Li X Y. Fracture detection using azimuthal variation of P2wave move out f rom orthogonal seismic survey lines. Geop hysics ,1999 ,64 (4) :1193120118 Thomsen L. Weak elastic anisot ropy

38、. Geophysics ,1986 ,51 (10) :1954196619 杜世通. 地震波动力学. 山东东营:石油大学出版社,1996 ,35750020 Shuey R T. A Simplification of the Zoeppritz equa2tions. Geophysics ,1985 ,50 (4) : 609634.21 Stewart R R. Gaiser J E. Brown R J. et al. Converted-wave seismic explorations applicationsJ. Geophysics.2003.68(1), 40-5722

39、Ata E.Michelena R J. Mapping distribution of fractures in a reservoir with P-S converted waves J.The Leading Edge. 1995.14 (6). 664-67623 Van Dok R R. Gaiser J E. Jackson A R. et al. 3-D converted-wave processing: Wind River Basin case history J. Expanded Abstracts of the 670 Annual Internat SEG Mtg

40、.1997.1206-120924 Gaiser J E. Applicatijons for vector coordinate systems of 3-D converted-wave data J. The Leading Edge. 999.1801), 1290-130025 Gaiser J E. Van Dok R R. Borehole calibration of PS-waves for fracture characterizations Pinedale field.Wyoming J. Expanded Abstract of 750 Annual Ice tern

41、at SEG Mtg.2005.873-87626 Grechka V.Kachanov M. Seismic characterization of multiple fracture sets, does orthotropy sufficeJ.Geophysics. 2006.71(3), D93-D10527 Grechka V.Kachanov M. Effective elasticity of fractured rocks, a snapshot of the work in progress J.Geophysics. 2006.71(6), W45-W5828 Vascon

42、celos I. Grechka V. Seismic characterization of multiple fracture sets at Rulison Field. ColoradoJ. Geophysics. 2007.72(2), B19-B3029 Engelder T.Lash GG.Uzcdtegui R S. Joint sets that enhance production from middle and upper Devonian gas shales of the Appalachian basin J. AAPG Bulletin.2009.93(7) ,8

43、57-88930 Bale R.Gratakos B. Mattocks B. et al. Shear wave splitting applications for fracture analysis and improved imaging, some onshore examplesJ). First Break.2009.27(9), 73-8331 Roxana M V. Using multrcomponent seismic data to delineate hydrocarbon reservoirs: 2D-3C Willesden Green. Alberta and

44、3D-3C Manitou Lake. Saskatchewan D. Alberta.Canada:Urivemity of Calgary.200932 Tang J M.Huang Y.Xu X R.et al. Application of converted-wave 3D/3C data for fracture detection in a deep tight-gas reservoirJ. The Leading Edge.2009.28(7), 826-83733 蔡希源.启建明.陈本池.三维多分里地君技术在川西新场地区深层致密砂岩裂缝检测中的应用J.石油学报.2010.3

45、1(5):737-743Cai X Y. Tang J M.Chen B C. A 3D multi-component seismic exploration technique applied to the detection of deep-seated tight sandstone fractures inthe Xinchang area of the western Sichuan BasinJ.Acts Petrolei Sinica.2010.31(5),737-74334 Dumitru G. Bale R. Minimum entropy Rotation: a new

46、shear-wave splitting technique for converted wave data J. Expanded Abstract of 700 Annual internat SEG Mtg.2000,1229-123235 黄中玉.赵金洲.正交基旋转的横波分裂检测技术J.石油地珠物理勘探.2004.39(2).143-152.Huang Z Y. Zhao J Z. Technique of S-wave splitting detection by orthonormal basis rotation J. Oil Geophysical Prospecting. 2

47、004.39(2),149 15236 Dariu H P.Granger Y. Bireftingenee analysis using simulated annealing J. Expanded Abstract of 67EAGE Annual Conference. 2005 ,Z9937 Dewangan P.Tsvankin 1. Velocity independent layer stripping of PP and PS reflection traveltimes J.Geophysics.2006.71(4), U59-U6538 Haacke R R. Westb

48、rook G K. Peacock S. Layer stripping of shearwave splitting in marine PSwavesJ3. Geophysical Journal International. 2009. 176(3), 782-80439 Wang X.Tsvankin I. Estimation of interval anisotropy parameters using velocity independent layer strippingJ.Geophysics.2009.74(5), WBI17-WB12740 Grechka V.Tsvan

49、kin I. NMO-velocity surfaces and Dix-type formulas in anisotropic heterogeneous media J" Geophysics. 2002.67 (3). 939-95141 Tsankin L. Reflection moveut and parameter estimation for horizontal tranverse isotropy J. Geophysics.1996.62(2), 479-49242 sconcelos I. Tsvankin I. Nonhyperbolc moveoutin

50、version of wide-azimuth P-wave data for orthorhombic media J" Geophysical Prospecting. 2006.54(5), 535-552.43 桂先志.段天友.肠远元.林.裂缝性偏层纵波地震检测方法研究JJ.石油天然气学报，2007.29(4):75-79. Gui X Z. Duan T Y.Yi Y Y.et aL On Pwave Seismic Detection Methods for Fractured ReservoirsJ.Journal of Oil and Gas Technology.

51、2007.29(4),75-79.44 Li X Y. Fracture reservoir delineation using mult-component seismic data J. Geophysical Prosperting.1997.45(1). 39-64.45 Li X Y. Fracture detection using azimuthal variation of P-wave moveout from orthogonal seismic survey lines J" Geophysics. 1999.64(4): 1193-1201.46 Zheng

52、Y. Seismic azimuthal anisotropy and fracture analysis from PP reflection data D. Alberta.Canada:University of Calgary.200647 Kozlov E A.Varivoda D Y. Dense 3D residual move out analysis as a tool for HTI parameter estimationJ ). Geophysical Prospecting. 2005. 53 (1 ):131-148.48 Perez M A. Detection

53、of fracture orientation using azimuthal variation of Pwave AVO responses D.Massachusetts. USA: Massachusetts institute of technology.1997.49 Xu X. Tsvankin I. A case study of azimuthal AVO analysis with anisotropic spreading correctionJ. The Leading Edge. 2007.26(12): 1552-1561.50 杨劫勇.赵群.王世且.等.纵波方位各

54、向异性及其在裂缝检侧中的应用J.石油物探.2006.45(2).177-181. Yang Q Y. Zhao Q. Wang S X. et al. P- wave azimuthal anisotropy and its application in detection of fractures J. Geophysical Prospecting for Petroleum.2006.45(2) ,177-181.51Zheng Y.Todorovic-Marinic D.Larson G. Fracturte detection:Ambiguity and practical solu

55、tion J. Expanded Abstract of 74th Annual Internat SEG Mtg.2004,1575-1578.52 曲寿利.季玉新.王鑫.等.全方位P波属性裂缝检测方法【J.石油地球物理助探.2001.36(4),390-397. Qu S L.Ji Y X.Wang X.et al. Seismic method for using full azimuth P-wave attribution to detect frayture J. Oil Geophysical Prospecting. 2001.36(4) ,390-397.53 Qu S L.

56、Ji Y X.Wang X.et al. Fracture detection by using full azimuth P wave attributes J. Applied Geophysics. 2007.4 (3) , 238-24354 Salamoff S I. the use of complex seimic reflection attributes to delineate subsurface fracture networks,an example from Teapot Dome. Wyoming D. Colorado. USA ,Colorado State

57、University.200655 Roberts A. Curvature attributes and their application to 3D interpreted horizons J. First Hreak,2001.19(2), 85-10056 Hart B S.Pearson R. Rawling G C.3-D seismic horizon based approaches to fracture-swarm sweet spot definition in tight-gas reservoirs J. The Leading Edge.2002.21(1), 28-3557 Sigismondi M E.Soldo J C. Curvature attributes and seismic interpretation, Case studies from Argentina basins J. The Leading Edge. 2003. 22 (11 ),1122-112658 Marfurt K J. Robust estimates of reflector dip and azimuth J