一种利用工程参数计算泥岩地层纵波速度的方法

本发明涉及一种利用工程参数来计算泥岩地层纵波速度的方法，可以为随钻三压力评价提供基础参数，属于钻井工程领域。

背景技术：

1、在油气钻井过程中，及时准确的获取地层三压力剖面是钻井成功的重要基础。地层三压力指的是地层孔隙压力、地层破裂压力和地层坍塌压力，当钻进过程中钻柱与裸眼井壁之间的钻井液柱压力低于地层孔隙压力时，可能导致溢流甚至井喷；当钻井液柱压力低于地层坍塌压力时，可能导致泥岩地层失稳；当钻井液柱压力高于地层破裂压力时，可能导致地层破裂形成井漏。目前的地层三压力评价工作可分为钻前评估、钻中监测和钻后评价三个方面，其中钻前和钻中阶段缺乏实时的地层声波速度数据，监测阶段一般只能计算地层孔隙压力和地层破裂压力、无法计算地层坍塌压力，由于钻井液密度低于地层坍塌压力导致的钻井复杂事故层出不穷，严重影响了钻进效率和勘探进程。

2、目前获取泥岩地层纵波速度的方式有以下几类：

3、1、在电缆测井作业后获得。在下入下级套管或完钻之前，为了取全取准地质资料，作业者们的普遍做法是针对裸眼地层进行电缆测井作业，声波测井一般是常测或必测项目之一。由于地层刚刚打开，且电缆测井工具在测量过程中紧贴裸眼井壁，电缆测井获得的地层声波数据被认为最接近地层真实状态，是技术人员常用的基础数据类型。尽管电缆声波测井具有测量密度高、连续性好和一致性强的优点，但易受井眼不规则、泥浆滤液侵入等因素的困扰，且只能用于钻后或中途完钻阶段。

4、2、通过随钻测井作业获得。为了在随钻过程中获得泥岩地层的速度数据，测井工作者开发了随钻声波测井工具，相对于电缆测井技术，随钻测井不受井斜和泥浆滤液侵入地层的影响，更真实的反映了原状地层的地质特征。但受限于钻井噪声、仪器偏心和声波频散等因素的制约，并且该技术受井斜和井眼环境影响较大，其仪器工具、数据传输和处理解释技术仍不成熟，随钻声波测井技术仅应用于海域中浅部钻井等部分钻井市场，在高温高压井地层三压力评价中的案例极少。

5、3、通过岩屑声波测量获得。为了及时获取钻进过程中的地层声波数据，研究者开发了针对岩屑进行超声波测量的技术，该技术测量成本低、实时性较好，但由于目前普遍采用pdc钻头造成岩屑细小，实际测量效果不佳，未普遍推广应用。

6、4、提取井旁地震道获得。利用地震层速度资料进行地层压力预测是新探区或者新层系探井的唯一选择，但由于地震资料受岩性和孔隙度等非地层压力因素影响较大，且地震层速度分辨率较井筒资料低，因此仅应用于钻前预测阶段。

7、5、通过元素录井等随钻资料进行回归分析获得。利用岩屑的元素录井资料和实测的电缆声波测井数据进行多元线性回归，可以得到二者之间的转换模型。由于元素录井资料在实质上是地层岩石矿物组分的反映，极易受地层在纵向和横向上的非均质性影响，因而计算的声波速度结果精度不高。

8、6、使用区域经验公式进行计算获得。这种方法是利用区域已钻井的声波速度资料进行深度回归，进而推测正钻井的地层声波速度分布情况。该方法在本质上是通过邻井对比得到本井的地层属性信息，在构造稳定、区域异常地层压力分布较为规律的地区和层系中可以使用，但无法应用于地层非均质性强、深层和地层高压的探井地层压力监测。

9、由于泥岩地层纵波速度是计算地层三压力的关键参数，而目前获取泥岩地层纵波速度的测量工具或计算方法无法在随钻过程中满足资料需求，严重制约了随钻过程中的地层三压力计算和井壁稳定性评价工作。工程参数是钻进过程中自然产生和实施记录的基础数据，提出利用工程参数来计算可靠精度的泥岩纵波速度的方法，对于保障钻井安全和提速增效具有重要意义。本发明所述的有益效果：

10、1、本发明依赖的计算数据来自于钻机钻进过程中的伴生参数，不需要增加其他的测试项目即可完成计算，具有明显的经济性；

11、2、本发明可在钻井过程中实时评价变质岩地层裂缝的发育情况，相较于其他方法具有明显的实时性。

12、3、本发明提供一种在随钻过程中，利用实时的工程参数来识别和评价地层裂缝的方法，通过对钻井过程中实时采集的工程参数进行分析，建立校正的修正钻压指数与裂缝发育特征之间的响应关系，从而根据实时的钻井工程参数对钻遇的地层裂缝进行识别和评价。

技术实现思路

1、本发明设计开发了一种利用工程参数计算泥岩地层纵波速度的方法，根据已钻井资料，建立地层孔隙压力和纵波速度的数学转换模型，结合正钻井的岩石可钻性和岩屑录井数据计算得到的地层孔隙压力数据，通过转换模型得到正钻地层泥岩层的纵波速度，可以为随钻地层三压力评价提供实时且可靠的基础数据。

2、本发明提供的技术方案为：

3、一种利用工程参数计算泥岩地层纵波速度的方法，包括如下具体步骤：收集已钻井的地层孔隙压力、声波测井、井径测井和岩屑录井数据，提取已钻井未扩径段泥岩层的地层孔隙压力和纵波速度数据的步骤；根据常压且未扩径段泥岩层的纵波速度值与深度数据，建立泥岩层纵波速度趋势线的步骤；参照提出的回归模型，建立已钻井泥岩层地层孔隙压力和纵波速度量化关系的步骤；获取正钻井的随钻工程参数和岩屑录井资料，提取正钻井泥岩层岩石可钻性参数的步骤；基于eaton法，以正钻井泥岩层的岩石可钻性参数计算泥岩层地层孔隙压力的步骤；根据已钻井资料建立的回归模型和正钻井资料计算的地层孔隙压力数据，计算正钻井泥岩层纵波速度的步骤。

4、其中，收集已钻井的地层孔隙压力、声波测井、井径测井和岩屑录井数据，提取已钻井未扩径段泥岩层的地层孔隙压力和纵波速度数据的步骤具体包括：搜集已钻井资料的步骤，包括钻头尺寸、地层孔隙压力、声波测井、井径测井和岩屑录井数据；提取已钻井未扩径段深度序列的步骤，根据实测井径测井数据与钻头尺寸之差除以钻头尺寸小于10％为通过标准，确定未扩径段的深度序列；根据岩屑录井数据确定未扩径段泥岩点深度序列的步骤，以单层厚度不小于2m为判别标准，确定未扩径段泥岩点深度序列；根据地层孔隙压力数据，提取常压且未扩径段泥岩点声波速度的步骤，地层孔隙压力系数大于1.05的最小深度数据为常压段下限，该下限以上地层为常压段，然后根据常压段深度范围提取常压且未扩径段泥岩点纵波速度数据；

5、根据常压且未扩径段泥岩层的纵波速度值与深度数据，建立泥岩层纵波速度趋势线的步骤，利用如下公式，计算得到泥岩点纵波速度趋势线

6、vn＝e^(a×h+b)

7、式中，vn为纵波速度趋势值，单位为m/s，h为深度，单位为m，e为自然对数，a、b为常量。

8、参照提出的回归模型，建立已钻井泥岩层地层孔隙压力和纵波速度量化关系的步骤，根据未扩径段泥岩点深度序列，从地层孔隙压力数据和声波测井数据中提取未扩径段泥岩点的相应数据，按照以下公式回归得到模型相关系数。

9、v＝vn×[(po-pp)/(po-ph)]^c

10、式中：v为纵波速度测量值，单位为m/s，po为上覆地层压力，单位为g/cm3；pp为静水压力，单位为g/cm3；ph为地层孔隙压力，单位为g/cm3；c为回归系数，无量纲。

11、获取正钻井的随钻工程参数和岩屑录井资料，提取正钻井泥岩层岩石可钻性参数的步骤，具体包括：搜集正钻井的dc指数和岩屑录井数据的步骤；提取泥岩点深度序列数据的步骤，以单层厚度不小于2m为判别标准，确定泥岩点深度序列；提取泥岩点深度序列的dc指数数据的步骤，参照泥岩点深度序列，提取对应的dc指数数据；

12、基于eaton法，以正钻井泥岩层的岩石可钻性参数计算泥岩层地层孔隙压力的步骤参照sy/t5623—2009，具体包括：根据中浅层或常压段泥岩点dc指数与深度的变化规律，建立dc指数趋势线；根据dc指数和dc指数趋势线计算地层孔隙压力的步骤。

13、根据已钻井资料建立的回归模型和正钻井资料计算的地层孔隙压力数据，计算正钻井泥岩层纵波速度的步骤包括：根据泥岩点深度序列提取泥岩点地层孔隙压力数据的步骤；根据泥岩点地层孔隙压力数据和已钻井回归模型，计算得到正钻井泥岩层纵波速度。