[发明专利]页岩气储层非均质性的确定方法、装置和设备

权利要求书

1.一种页岩气储层非均质性的确定方法，其特征在于，包括：

获取目标页岩气储层的水力压裂施工曲线在预设时间段内的多个压裂段的套压数据；

基于各个压裂段的套压数据进行经验模态分解，得到各个压裂段分解得到的多个本征模态函数并确定各个压裂段分解得到的本征模态函数的数量；

对所述各个压裂段的多个本征模态函数进行希尔伯特变换，得到各个压裂段的至少一个希尔伯特谱图；

基于所述各个压裂段分解得到的本征模态函数的数量和所述各个压裂段的至少一个希尔伯特谱图，确定所述目标页岩气储层的非均质性；

其中，在得到各个压裂段分解得到的多个本征模态函数之后，还包括：根据所述各个压裂段的套压数据和分解得到多个本征模态函数，确定所述各个压裂段分解得到的各个本征模态函数的相关性系数；将所述各个压裂段中相关性系数大于预设阈值至少一个本征模态函数作为各个压裂段的特征本征模态函数；

根据对所述各个压裂段的多个本征模态函数进行希尔伯特变换，获取各个压裂段的至少一个希尔伯特谱图，包括：对所述各个压裂段的特征本征模态函数进行希尔伯特变换，得到所述各个压裂段的特征本征模态函数的希尔伯特瞬时频率和归一化瞬时振幅；基于所述各个压裂段的特征本征模态函数的希尔伯特瞬时频率和瞬时振幅，得到各个压裂段的特征希尔伯特谱图；

基于所述各个压裂段分解得到的本征模态函数的数量和所述各个压裂段的至少一个希尔伯特谱图，确定所述目标页岩气储层的非均质性，包括：获取任意两个压裂段分解得到的本征模态函数的数量；按照如下规则确定所述目标页岩气储层的非均质性：在所述任意两个压裂段分解得到的本征模态函数的数量不相等的情况下，分解得到的本征模态函数的数量越多的压裂段的非均质性程度越高。

2.根据权利要求1所述的方法，其特征在于，基于所述各个压裂段分解得到的本征模态函数的数量和所述各个压裂段的至少一个希尔伯特谱图，确定所述目标页岩气储层的非均质性，还包括：

按照如下规则确定所述目标页岩气储层的非均质性：在所述任意两个压裂段分解得到的本征模态函数的数量相等的情况下，获取所述任意两个压裂段的特征希尔伯特谱图；比较所述任意两个压裂段的特征希尔伯特谱图，其中，特征希尔伯特谱图中瞬时频率和归一化瞬时振幅越高的压裂段的非均质性程度越高。

3.根据权利要求1所述的方法，其特征在于，所述预设时间段的起始点为所述水力压裂施工曲线中套压数据从最高点下降百分之十对应的时间点，所述预设时间段的终点为所述水力压裂施工曲线中套压数据下降至储层压力对应的时间点。

4.根据权利要求1所述的方法，其特征在于，基于各个压裂段的套压数据进行经验模态分解，得到各个压裂段分解得到的多个本征模态函数，包括：

获取目标压裂段的套压数据随时间变化的所有极大值点和极小值点；

基于所述目标压裂段的套压数据随时间变化的所有极大值点，利用三次样条插值函数拟合形成上包络线；

基于所述目标压裂段的套压数据随时间变化的所有极小值点，利用三次样条插值函数拟合形成下包络线；

基于所述上包络线和所述下包络线，确定所述上包络线和所述下包络线的均值；

基于所述上包络线和所述下包络线的均值，确定所述目标压裂段的套压数据残差；

确定所述目标压裂段的套压数据残差局部零极点与过零点的数目是否相等或相差一个，和所述上包络线和所述下包络线的均值是否为0；

在确定所述目标压裂段的套压数据残差局部零极点与过零点的数目相等或相差一个，和所述上包络线和所述下包络线的均值为0的情况下，将所述目标压裂段的套压数据残差作为分解得到的一个本征模态函数。