一种页岩水化损伤系数的预测模型建立方法及预测方法与流程

1.本发明涉及页岩气勘探开发领域，尤其涉及一种页岩水化损伤系数的预测模型建立方法及预测方法。


背景技术：

2.我国页岩气资源丰富，为应对我国对于天然气的强劲需求，加大页岩气勘探开发是我国石油工业的重要任务。其中，页岩地层井壁失稳问题是阻碍页岩气高效开发的重大技术难题之一。页岩地层井壁失稳的主要原因是页岩水化损伤特征。由于页岩的强水化特征，钻井液与井壁岩石接触，形成强水化损伤，造成井壁稳定性降低，从而导致井壁垮塌失稳。
3.因此，为建立合理的控制井壁稳定技术，需要对页岩地层水化损伤进行精准评价，从而形成对应的技术措施。现有水化损伤评价方法的考虑因素较为单一，主要针对水化时间、流体类型进行分析。但是，实际水化过程中，水化损伤程度与外部条件、流体类型、岩石自身特征都具有相关性。外部温度、压差、流体特征、岩石密度、孔隙度、渗透率等参数都会对水化损伤造成影响。目前，尚未存在一套综合方法对水化损伤进行定理评价。


技术实现要素：

4.本技术为了解决上述技术问题提供一种页岩水化损伤系数的预测模型建立方法及预测方法。
5.本技术通过下述技术方案实现：
6.一种页岩水化损伤系数的预测模型建立方法，包括以下步骤：
7.制备页岩标准圆柱岩心样品。
8.对页岩岩样开展超声波透射实验，获取原状页岩的频谱特征。
9.基于室内实验手段，获取页岩水化参数，主要分为页岩岩石参数(页岩活度、页岩密度、孔隙度、渗透率)、流体参数(流体活度)及外部环境参数(浸泡时间、浸泡温度、浸泡压差)。
10.基于不同页岩岩石参数、流体参数及外部环境参数，开展浸泡实验。
11.对浸泡后的页岩岩心，开展超声波透射实验，获取水化后的页岩频谱特征。
12.根据原状页岩频谱特征、水化后页岩频谱特征，基于水化损伤系数计算模型，获得水化损伤系数实验值；
13.基于页岩岩石参数、流体参数、外部环境参数、水化损伤系数，采用神经网络算法，建立基于页岩水化参数(岩石参数、流体参数、外部环境参数)的页岩水化损伤系数预测模型。
14.将待预测岩样的页岩水化参数，输入建立的页岩水化损伤系数预测模型，即可获得页岩水化损伤系数，实现页岩水化损伤效应评价。
15.与现有技术相比，本技术具有以下有益效果：
16.1，本技术考虑了现有页岩地层水化效应评价方法的缺陷，利用水化前后的声波频谱特征，借助神经网络算法，构建了页岩水化损伤系数的预测模型；使用该预测模型预测时，所有输入参数均为原状页岩参数，即是说，使用本技术的预测方法获取页岩水化损伤系数，不会对页岩岩样造成损坏，从而增加了页岩岩样利用率；
17.2，本技术的预测方法，输入参数包括了岩石特征、流体特征、外部环境特征，能够较为全面系统地实现水化损伤评价；
18.3，本技术通过水化前后的声波频谱特征，借助bp神经网络算法，以不同页岩岩石参数、流体参数及外部环境参数为输入参数，建立了页岩水化损伤系数的预测模型，实现了页岩水化损伤的定量评价，对形成控制页岩地层井壁稳定技术具有重要意义。
附图说明
19.此处所说明的附图用来提供对本技术实施方式的进一步理解，构成本技术的一部分，并不构成对本发明实施方式的限定。
20.图1为本发明流程框图；
21.图2是页岩频谱示意图；
22.图3中a是基于页岩频谱的水化损伤系数示意图，b是12块岩样的水化损伤系数实验值分布图；
23.图4是bp神经网络计算与学习的示意图；
24.图5是页岩水化损伤系数预测结果与实测结果对比图。
具体实施方式
25.为使本技术的目的、技术方案和优点更加清楚，下面将结合实施方式中的附图，对本发明实施方式中的技术方案进行清楚、完整地描述。显然，所描述的实施方式是本发明一部分实施方式，而不是全部的实施方式。通常在此处附图中描述和示出的本发明实施方式的组件可以以各种不同的配置来布置和设计。
26.因此，以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。
27.需要说明的是，在不冲突的情况下，本发明中的实施方式及实施方式中的特征可以相互组合。需要说明的是，本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。
28.实施例1
29.本实施例以12块页岩岩样与6种流体类型为基础，获取输入层参数，进而开展水化损伤系数的预测。具体包括以下步骤：
30.步骤一、制备页岩标准圆柱岩心。采用页岩野外露头，通过钻机钻取标准圆柱样岩心12块。
31.步骤二、对原状页岩岩样开展超声波透射实验，获取原状页岩的频谱特征。
32.其中，超声波透射实验可以直接得到页岩波形图，借助傅里叶转换，可以得到页岩频谱图。基于傅里叶变化的页岩频谱图如图2所示。
33.步骤三、基于室内实验，获取页岩水化参数。页岩水化参数主要分为页岩岩石参数、流体参数和外部环境参数。
34.页岩岩石参数具体包括页岩活度、页岩密度、孔隙度、渗透率；
35.流体参数包括流体活度；
36.外部环境参数具体包括浸泡时间、浸泡温度、浸泡压差。
37.(a)针对页岩岩石参数：
38.其中，页岩密度采用体积法测量，通过天平测量页岩岩样质量。采用游标卡尺分别测量页岩岩样直径与长度，进而计算岩样体积。基于岩样质量与体积，计算得到页岩岩样密度。
39.孔隙度采用气测法获取：首先通过游标卡尺测量得到岩样体积v
f
，然后将页岩放入样品室开始气测，当样品室两端压力平衡时，基于平衡压力查阅图版，获得岩样骨架体积v
s
。
40.根据岩样体积v
f
和骨架体积v
s
，采用以下公式计算岩样孔隙度
[0041][0042]
渗透率采用气测法获取：基于达西定理，气测渗透率公式为：
[0043][0044]
上式中，k为渗透率；p
o
为大气压力，固定为0.1mpa；q为大气压下的气体流量，cm3/s；μ为气体粘度，本实施例采用氮气，粘度为0.0178mpa.s；l为岩样长度，a横截面积，l与a通过游标卡尺测试计算；p1为气测时岩样室里的进气端压力；p2为气测时岩样室里的出气端压力。
[0045]
页岩活度测试方法基于饱和湿度法原理，根据湿度与活度的线性相关性获取。
[0046]
湿度与活度的线性相关性如下所示：
[0047]
a
w
＝ar
h
b
[0048]
上式中，a
w
为活度，r
h
为湿度；a、b为线性拟合系数，分别为0.0099和
‑
0.0038。基于湿度与活度相关性，利用湿度仪首先测试页岩湿度，进而得到页岩活度。
[0049]
(b)针对流体参数：
[0050]
流体活度测试方法与页岩岩样活度测试方法一致，基于流体湿度与流体活度的线性相关性，首先测试流体湿度，进而得到流体活度。
[0051]
(c)针对外部环境参数：
[0052]
选取不同浸泡时间、浸泡温度及浸泡压差，作为外部环境参数。其中，浸泡温度通过烘箱加热，浸泡压差通过向浸泡容器内部充入氮气，形成压差。
[0053]
综上，页岩水化参数表1所示。
[0054]
表1：页岩水化参数表
[0055][0056][0057]
步骤四、基于不同页岩岩石参数、流体参数及外部环境参数，开展浸泡实验。
[0058]
步骤五、对浸泡后的页岩岩心，开展超声波透射实验，获取水化后的页岩频谱特征，根据步骤二获得的原状页岩频谱与水化后页岩频谱，基于水化损伤系数计算模型，获取水化损伤系数实验值。
[0059]
考虑水化损伤导致页岩内部孔隙、裂缝增多，页岩内部结构变化必然导致声波频谱变化。采用声波频谱特征中，水化前(原状)页岩与水化后页岩主频与最大振幅为依据，建立了水化损伤系数计算模型，如下所示：
[0060][0061]
式中，s为水化损伤系数，m
a
为原状页岩主频，m
c
为水化页岩主频，f
a
为原状页岩的最大振幅，f
b
为水化页岩的最大振幅。
[0062]
水化损伤系数的构建示意与分布如图3所示。
[0063]
步骤六、基于页岩岩石参数、流体参数、外部环境参数、水化损伤系数实验值，采用神经网络算法，建立基于页岩水化参数(岩石参数、流体参数、外部环境参数)的页岩水化损伤效应评价方法。
[0064]
其中，bp神经网络算法主要分为输入层、隐含层、输出层。基于实验测试数据，通过对不同页岩岩石参数、流体参数及外部环境参数进行归一化处理，归一化后参数为输入层参数(8层输入参数)。水化损伤系数实验值作为输出层参数，即水化损伤系数实验值为理想输出参数。隐含层是用于神经网络学习，基于输入层参数与输出层的理想输出参数，开展神经网络学习。隐含层的层数为：
[0065]
[0066]
上式中，n
y
为隐含层层数；n
in
为输入层层数；m
out
为输出层层数；η为经验常数。实例分析中，输入层为8层，输出层1层，经验常数取6，选定隐含层为9层。
[0067]
确定bp神经网络的输入层、输出层和隐含层后，如图4所示，采用200次最大学习次数，精度为0.001，开展bp神经网络计算，建立基于页岩水化参数的页岩水化损伤预测模型。
[0068]
结果验证：基于本实施例的水化损伤预测方法，选用7块岩样进行验证分析。具体如下：
[0069]
获取7块岩样的页岩水化参数，页岩水化参数包括页岩活度、页岩密度、孔隙度、渗透率、流体活度、浸泡时间、浸泡温度、浸泡压差；
[0070]
将7块岩样的页岩水化参数输入本发明方法建立的页岩水化损伤预测模型，获得水化损伤系数的预测值；
[0071]
同时，按照前述方法，开展超声波透射实验，获取7块岩样的原状页岩频谱与水化后页岩频谱，根据水化损伤系数计算模型，获得7块岩样的水化损伤系数实验值。
[0072]
对比水化损伤系数的预测值与水化损伤系数实验值，如图5所示。由图5可知：预测结果与实验测试结果的整体吻合度较好，证明了本发明的水化损伤系数预测方法的准确性。
[0073]
以上所述的具体实施方式，对本技术的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。