1.一种各向异性地层含水饱和度计算方法,其特征在于,包括如下步骤:步骤1:获取待检测地层的电阻率张量矩阵R,其中, ρij表示待检测地层的i方向对于j方向上叠加的电流密度所能产生的电场强度的能力,i,j∈{x,y,z};
步骤2:根据步骤1中获得的电阻率张量矩阵R计算待检测地层的一阶不变量I1、二阶不变量I2和三阶不变量I3;
步骤3:根据待检测地层的一阶不变量I1、二阶不变量I2和三阶不变量I3计算待检测地层的含水饱和度Sw,如式Ⅰ所示:其中,a为绕曲程度,m为胶结指数,φ为孔隙度,ρw为水电阻率。
2.如权利要求1所述的各向异性地层含水饱和度计算方法,其特征在于,步骤2中一阶不变量I1、二阶不变量I2和三阶不变量I3通过式Ⅱ获得:I1=tr(R) 式ⅡI3=det(R)
其中,tr(R)表示R的迹,det(R)表示R的行列式。
3.如权利要求1所述的各向异性地层含水饱和度计算方法,其特征在于,式Ⅰ中的参数a和m通过水饱和样本进行标定。
4.如权利要求3所述的各向异性地层含水饱和度计算方法,其特征在于,a=100,m=
1.08。
5.一种各向异性地层含水饱和度计算装置,其特征在于,包括采集及预处理模块和含水饱和度计算模块;
所述的采集及预处理模块用于获取待检测地层的电阻率张量矩阵R,其中,ρij表示待检测地层的i方向对于j方向上叠加的电流密度所能产生的电场强度的能力,i,j∈{x,y,z};还用于根据电阻率张量矩阵R计算待检测地层的一阶不变量I1、二阶不变量I2和三阶不变量I3;
所述的含水饱和度计算模块用于根据待检测地层的一阶不变量I1、二阶不变量I2和三阶不变量I3计算待检测地层的含水饱和度Sw,如式Ⅰ所示:其中,a为绕曲程度,m为胶结指数,φ为孔隙度,ρw为水电阻率。
6.如权利要求1所述的各向异性地层含水饱和度计算装置,其特征在于,采集及预处理模块中一阶不变量I1、二阶不变量I2和三阶不变量I3通过式Ⅱ获得:I1=tr(R) 式ⅡI3=det(R)
其中,tr(R)表示R的迹,det(R)表示R的行列式。
7.如权利要求1所述的各向异性地层含水饱和度计算装置,其特征在于,式Ⅰ中的参数a和m通过水饱和样本进行标定。
8.如权利要求3所述的各向异性地层含水饱和度计算装置,其特征在于,a=100,m=
1.08。