专利名称：一种获取煤岩的速度特征的控制方法
技术领域：
本发明涉及煤田勘探领域，尤其涉及一种获取煤岩的速度特征的控制方法。
背景技术：
随着煤炭资源开采深度的加大，深部高温热害和深层煤岩与瓦斯突出等动力学灾害的危险程度逐渐提高，对煤田地震勘探精度的要求也逐渐提高。目前，煤田勘探不仅要求地震勘探预测小断层，还要求对深层煤顶底板岩性、构造煤分布、瓦斯富集程度等情况进行定性或定量的预测。随着煤层气资源勘探开发力度的加大，水平井技术得以广泛的应用，煤层气开发方案的确定需要提前精确预测煤厚和煤岩中裂缝的发育程度与发育方向。因此，煤田地震与测井约束反演方法逐渐在煤田勘探中得到应用与推广。正如油气田地震反演技术应用一样，煤田地震的反演也需要煤岩密度、速度及其相互换算的定量关系支持，以降低地震反演中的不确定性与多解性。
由于煤岩是典型的软岩，因此取芯和制样难度较大。而且由于以前煤田地震勘探主要以构造勘探为主，对于煤岩的声学弹性模量并无要求，使得国内外对于煤岩的声学特征的研究相对较少。在煤岩弹性参数测量的两种方法中，以静态法较多，主要是利用三轴应力一应变测试技术进行煤岩变形测量与弹性模量计算，目的是通过煤岩的力学参数的测量服务于井下巷道开采工程。由于动态法可以方便地模拟地层条件，从而使得动态法测量的弹性参数对实际地球物理勘探的指导意义更为明显。由于煤炭并非国外的主要能源，加之国外的煤田构造简单，国际上对煤岩的弹性测量研究较少，而主要通过与测井数据的对比和合成地震记录重点分析煤系地层的地震反射特征及其对油气储层反射特征的影响。对于国内煤田勘探而言，随着阻抗反演的推广应用，对于煤岩弹性模量特征的研究需求越来越多。但目前通过实验室进行超声测试才能得到煤岩的速度特征，由于涉及设备过多且过程复杂，这样对于国内在野外的实际煤田勘探来说是不现实的，而煤岩的速度特征对于煤田勘探来说是十分重要的参数。因此当前需要一种获取煤岩的速度特征的技术方案来解决当前国内不能在野外的实际煤田勘探中，获取国内不同变质程度煤岩的速度特征的问题。
发明内容
本发明所要解决的技术问题是提供一种获取煤岩的速度特征的控制方法，解决了当前国内不能在野外的实际煤田勘探中，获取国内不同变质程度煤岩的速度特征的问题。为了解决上述问题，本发明提供了一种获取煤岩的速度特征的控制方法，包括获取需要勘探的煤岩的密度P并进行存储；根据获取的煤岩的密度P，通过下述方式计算得到该煤岩的垂直层理方向纵波速度Vp:VP=a P Vp+b,
其中P为密度，单位g/cm3 ；VP为垂直层理方向纵波速度，单位m/s ;a与b为统计系数，对于中国国内煤岩，a=0. 416，b=791。进一步地，上述方法还可包括所述获取需要勘探的煤岩的密度是通过钻孔取芯进行质量和体积测量，获得煤岩的密度；或者是通过野外现场的地球物理密度测井直接获得煤岩的密度。进一步地，上述方法还可包括根据得到的垂直层理方向的纵波速度，通过下述方式计算得到垂直层理方向的横波速度VP=a+bVs,其中Vp为垂直层理方向的纵波速度，Vs为垂直层理方向的横波速度，单位m/s ；a与b为统计系数，对于中国国内煤岩，a=-490, b=2. 23。进一步地，上述方法还可包括根据获取的煤岩的密度P，通过下述方式计算得 到该煤岩的垂直层理方向横波速度Vs=a P Vs+b,其中P为密度，单位g/cm3 ；VS为垂直层理方向横波速度，单位m/s ;a与b为统计系数，对于中国国内煤岩，a=0. 342，b=575。与现有技术相比，应用本发明，通过测得的煤岩的密度，可以较精确地得到煤岩的纵横波速度，不再需要利用煤岩样品进行实验室的超声测量，本发明适合野外的实际煤田勘探的需求；从实际测试中可以看出，对于国内不同变质程度煤岩弹性参数的预测精度较闻。
图I是本发明实际煤岩测试中煤岩测速示意图；图2是基于本发明实际煤岩测试的纵波速度与密度关系的示意图；图3是基于本发明实际煤岩测试的横波速度与密度关系的示意图；图4是基于本发明实际煤岩测试的纵横波速度间的关系的示意图；图5是本发明的获取煤岩的速度特征的控制方法的流程图。
具体实施例方式本发明的基于超声测试获取煤岩参数的控制方法在具体实现中要通过系统中各设备之间信息交互来进行信息和/或数据的收集，并通过其内的控制器(可以是CPU等进行控制处理信息和/或数据，本发明对此不作任何限定)，其间还可以通过各种存储器(可以是内存、硬盘或其他存储设备)进行信息和/或数据的储存和传送，本发明对此不作任何限定。真密度(即真相对密度，True Relative Density，TRD)，指材料在绝对密实状态下的体积内固体物质的实际体积,不包括内部空隙。多数材料为多孔物质,具有与外部相通的开口孔和不通的闭孔，将含有闭孔材料的密度称为视密度(即视相对密度，ApparentRelative Density,ARD)。煤的视密度是表示煤物理特性的一项指标,用于煤矿及地勘部门计算煤的埋藏量；在贮煤仓的设计及煤的运输、磨碎和燃烧等过程的计算中，也都需要这项测值。此外视密度还用来计算煤的孔隙率，作为煤层气计量基准。本发明中所指的密度为视密度，又称视相对密度，即地质勘探行业所使用的密度。
下面通过实际煤岩测试对本发明作进一步说明。本发明针对采自国内不同地区的不同变质程度煤岩的弹性速度进行了实验室测量，探索煤岩速度一密度、纵波速度一横波速度之间的关系，并通过与经典Gardner公式、Lindseth公式和Castagna公式的对比，尝试总结具有普遍适用性的煤岩速度一密度、纵波速度一横波速度的经验公式，以方便于煤田地震勘探的应用；并着重说明煤岩弹性参数的特殊性及其应用于煤田地震反演需要注意的问题。在本发明中采取下列相对完整数据的8块煤样的超声测试结果与密度数据进行分析。
样品编号采样地点煤样变质程度加工规格
A义马煤矿褐煤6cm见方 ~平顶山八矿戊组IeI6cm见方~
~~Cl平顶山八矿己组β6cm见方~
~~C2平顶山八矿己组β6cm见方~
D新疆气煤6cm见方
~鹤壁六矿2143工作面贫瘦煤6cm见方~
~~Fl焦作方庄矿无烟煤6cm见方~
~~F2焦作九里山矿无烟煤6cm见方~表I测试样品信息在实际中，超声测量采用常温常压行波传播一脉冲透射的方法进行测试。整套仪器由脉冲信号发生器、超声换能器、放大器、计数器和示波器组成，实验使用的是压电陶瓷柱状纵横波换能器。为保证样品与换能器耦合良好，测试纵波时采用凡士林进行耦合；测试横波时采用蜂蜜耦合。由于测试煤样为边长6cm的立方体，选用超声的低频段，主频为IOOKHz ;整个测量系统误差小于1% ;考虑到煤岩的特殊性，最大误差不超过3%。本发明中超声测量按两个平行截面分别测量纵波与横波速度，即分别测量了煤样沿煤层走向、倾向和垂向(垂直层理)3个方向的纵横波速度。为与煤田人工地震的观测方式相统一，分别以X、I、z分别代表煤层的走向、倾向和垂直层理的方向。如图3所示，Vx,Vy>Vz分别表示沿煤层走向、倾向和垂直层理的纵波速度。横波振动方向与波前方向垂直，由于煤样中裂隙的存在，横波通过煤样传播会分裂成两个相互垂直的横波，所以横波沿煤样某个方向传播时会有两个速度值。以沿X方向传播为例，沿X方向传播的横波有沿I、z两个方向的偏振，对应两个速度，记作Vxy与Vxz,下标的第一个字母X代表横波传播的方向，第二个字母代表与传播方向垂直的方向(即横波振动方向)。即Vxy表不横波沿X传播，振动方向与I平行；Vxz表示沿X传播，振动方向与z平行。纵横波速度计算采用T, L ,, LVp =-=--,( I )
tp - to ts-to其中VP为纵波速度，单位m/s ；VS为横波速度，单位m/s ；L为发射、接收换能器中心间的距离，单位m ；tp为纵波在样品中的旅行时单位s ；ts为横波在样品中的旅行时，单位s ；t0为仪器系统的零延时，单位S。在地震数据处理、解释和反演中，速度一密度的关系，尤其纵波速度一密度关系是经常要用到的公式。为此，针对煤岩分析它的纵横波速度与密度的关系，以认识煤岩的特殊性。如图2和图3所示为三方向纵波速度、横波速度与密度之间的试验散点图和线性回归方程。由图2、图3可以看出各个方向上的纵波速度与横波速度都随着密度增大而增大；即纵波速度、横波速度与密度均存在正线性相关性，平均相关系数O. 8。但三方向纵波速度-密 度相关性与横波速度-密度相关性规律不同纵波速度与密度的相关性按走向、垂向、倾向三方向划分依次降低，且三方向相关性差异不大，均在80%附近；而横波速度一密度的相关性按垂向、走向、倾向依次降低，且倾向方向的横波速度-密度的相关性比其它2个方向要差得较大。而且从图2和图3还可以看出，三方向纵波速度与横波速度存在以下的统一性规律走向速度大于倾向速度，垂向速度最�。蛔呦蛩俣扔氪瓜蛩俣炔钜旖洗螅阆蛩俣冉咏骄俣�。在煤田地震数据处理及与测井数据的联合反演中，经常需要用到密度一速度的关系式，以减少反演参数的数量。此外，由于很多煤田测井缺少声波曲线，还经常需要根据密度曲线进行纵波速度的换算。但到目前为止，按照油气领域井震联合反演的方法和商业软件，工程技术人员还习惯于使用Gardner公式进行两者关系的替换。而且在常规煤田地震勘探中一般是使用纵波勘探地下结构，其中纵波速度是地震勘探需要预测的主要参数之
OGardner 公式如下p=cl( 2 )其中P为密度,单位g/cm3 ;VP为垂向纵波速度,单位m/s ;a与b为统计系数,无量纲，一般取 a=0.31，b=0. 25。虽然Gardner公式在油气勘探领域的应用得到了广泛的认可，但对于煤岩会产生巨大的误差，平均误差高达78. 5%。该公式是由国外专家根据国外的大量沉积岩(主要是砂岩和泥岩)的超声测量结果统计得到的，由于国内外沉积环境和构造条件的差异，该公式并不适合于中国地区的沉积岩使用；更不适合煤这种软岩的速度-密度特征。在Gardner公式的基础上,针对软岩的沉积特征，Lindseth提出了一个改进的经验公式，对煤岩应用效果较为理想。Lindseth公式如下Vp = a P Vp+b (3)其中P为密度,单位g/cm3 ;VP为垂向纵波速度,单位m/s ;a与b为统计系数,无量纲，一般取 a=0. 308，b=1054。由于我国煤岩沉积的复杂性，Lindseth公式预测的纵波速度与实际煤样也会产生平均10. 5%的相对误差,但该精度已远远高于Gardner公式。综上所述，根据本发明实测的煤岩垂向纵波速度与密度的关系数据，本发明得到如下的回归经验公式VP=aP Vp+b (4)其中P为密度,单位g/cm3 ;VP为垂向纵波速度,单位m/s ;a与b为统计系数,无量纲；针对中国6种典型变质程度煤取a=0. 416，b=791。从表2可以看到，显然该公式的统计系数不同于Lindseth公式，平均相对误差不到10%,远高于Gardner公式的精度,也略高于Lindseth公式的精度。
权利要求
1.一种获取煤岩的速度特征的控制方法，其特征在于，包括 获取需要勘探的煤岩的密度P并进行存储； 根据获取的煤岩的密度P，通过下述方式计算得到该煤岩的垂直层理方向纵波速度Vp: VP=a P Vp+b， 其中P为密度，单位g/cm3 ;VP为垂直层理方向纵波速度,单位m/s ;a与b为统计系数，对于中国国内煤岩，a=0. 416，b=791。
2.如权利要求I所述的方法，其特征在于， 所述获取需要勘探的煤岩的密度是通过钻孔取芯进行质量和体积测量，获得煤岩的密度；或者是通过野外现场的地球物理密度测井直接获得煤岩的密度。
3.如权利要求I所述的方法，其特征在于， 还包括根据得到的垂直层理方向的纵波速度，通过下述方式计算得到垂直层理方向的横波速度VP=a+bVs, 其中Vp为垂直层理方向的纵波速度，Vs为垂直层理方向的横波速度，单位m/s ;a与b为统计系数，对于中国国内煤岩，a=-490, b=2. 23。
4.如权利要求3所述的方法，其特征在于， 进一步包括根据获取的煤岩的密度P，通过下述方式计算得到该煤岩的垂直层理方向横波速度 Vs=a P Vs+b, 其中P为密度，单位g/cm3 ;VS为垂直层理方向横波速度,单位m/s ;a与b为统计系数，对于中国国内煤岩，a=0. 342，b=575。
全文摘要
本发明公开了一种获取煤岩的速度特征的控制方法，包括获取需要勘探的煤岩的密度ρ并进行存储，根据获取的煤岩的密度ρ，通过VP=aρVP+b的方式计算得到该煤岩的垂直层理方向纵波速度VP，其中a=0.416，b=791。本发明解决了当前国内不能在野外的实际煤田勘探中，获取国内不同变质程度煤岩的速度特征的问题。
文档编号G01V1/30GK102809761SQ20121025684
公开日2012年12月5日 申请日期2012年7月23日 优先权日2012年7月23日
发明者王赟 申请人:中国科学院地质与地球物理研究所
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
