[发明专利]利用纵横波速度的奇异性识别天然气水合物赋存类型方法

权利要求书

1.一种利用纵横波速度的奇异性识别天然气水合物赋存类型方法，其特征在于，包括以下步骤：

S1、针对纵波速度，逐点计算奇异性值；针对横波速度，逐点计算奇异性值；

S2、制作纵波速度奇异性和横波速度奇异性交会图；

S3、在所述纵波速度奇异性和横波速度奇异性交会图基础上做统计分析，利用阈值方法识别天然气水合物赋存类型；

步骤S3的具体步骤如下所述：

a.对不含天然气水合物地区测井曲线的纵波速度的奇异性值α(Vp)和横波速度的奇异性值α(Vs)，做交会图分析，并得到一个最小范围值a，使得不含天然气水合物区域的散点分布在|α(Vp)-α(Vs)|≤a区间内；

b.对要分析的含天然气水合物地区测井段内的纵波速度奇异性值α(Vp)和横波速度的奇异性值α(Vs)交会图做统计分析：设散点落在|α(Vp)-α(Vs)|≤a区间内的个数为N₁，设散点落在|α(Vp)-α(Vs)|≤a区间外的个数为N₂；

c.给定一个阈值b，当N₁/N₂＝b时，判定为孔隙型水合物，当N₁/N₂b时，判定为裂隙型水合物；

所述阈值b的大小根据理论实验结果得出。

2.如权利要求1所述的利用纵横波速度的奇异性识别天然气水合物赋存类型方法，其特征在于，步骤S1的具体步骤如下所述：

a.针对实际测井的纵波速度曲线，将纵波速度信号中的每一点都看作奇点，利用小波变换计算得到不同尺度上的模极大值，直接通过分析纵波速度数据中每一点上的模极大值随尺度的变化情况，进而进行奇异性分析；

从奇点之上的模极大值出现的最小尺度处开始，根据Mallat的奇异性值计算公式，利用最小二乘法拟合，得到以下奇异性计算公式中的奇异性值α，如式(2)所示，

其中，x为进行奇异性分析的点，m为用来计算奇异性值α的小波变换模极大值的个数，s_i为第i个尺度，|Wf(x,s_i)|为在尺度s_i情况下x点之上的小波变换模极大值，α为要计算的奇异性值；

b.同理，针对实际测井的横波速度曲线，将横波速度信号中的每一点都看作奇点，利用小波变换计算得到不同尺度上的模极大值，直接通过分析横波速度数据中每一点上的模极大值随尺度的变化情况，进而进行奇异性分析；

从奇点之上的模极大值出现的最小尺度处开始，根据Mallat的奇异性值计算公式，利用最小二乘法拟合，利用公式(2)得到奇异性值α。

3.如权利要求2所述的利用纵横波速度的奇异性识别天然气水合物赋存类型方法，其特征在于：步骤S2的具体步骤如下所述：

a.假设纵波速度为Vp、横波速度为Vs，则纵波速度Vp对应的奇异性值为α(Vp)，横波速度Vs对应的奇异性值为α(Vs)；

b.以纵波速度奇异性值α(Vp)为横坐标，以横波速度奇异性值α(Vs)为纵坐标，得到α(Vp)和α(Vs)的交会图；

当α(Vp)与α(Vs)相等时，在所述交会图上，形成直线α(Vp)＝α(Vs)；

通过天然气水合物的数据计算的α(Vp)和α(Vs)形成交会散点；

当所述天然气水合物为裂隙型水合物时，所述交会散点散乱分布在直线α(Vp)＝α(Vs)之外；当所述天然气水合物为孔隙型水合物时，所述交会散点集中分布在直线α(Vp)＝α(Vs)附近。