阅读说明：本技术 基于流体替换的井中流体因子敏感度计算方法及系统 (Well fluid factor sensitivity calculation method and system based on fluid replacement ) 是由 毕进娜 陈科 张瑶 林庆富 张杨 朱凌燕 王佳 于 2018-09-10 设计创作，主要内容包括：公开了一种基于流体替换的井中流体因子敏感度计算方法及系统。该方法包括：根据包含初始纵波速度、初始横波速度、初始密度的测井曲线,确定储层位置,建立初始模型；对初始模型进行Gassmann流体替换,获得替换模型,计算替换模型的纵波速度、横波速度、密度；根据初始纵波速度、初始横波速度、初始密度与替换模型的纵波速度、横波速度、密度,分别计算初始模型与替换模型的多个流体识别因子；根据流体因子敏感性公式,计算每个流体识别因子对应的敏感性。本发明通过Gassmann流体替换,确定流体识别因子及其敏感性,优选敏感性高的流体识别因子,有助于精确反应储层特征,寻找含油或者含气的储层,为常规油气、煤层气及页岩油气的勘探、开发服务。(A method and system for calculating fluid factor sensitivity in a well based on fluid replacement is disclosed. The method comprises the following steps: determining the position of a reservoir according to a logging curve comprising an initial longitudinal wave velocity, an initial transverse wave velocity and an initial density, and establishing an initial model; carrying out Gassmann fluid replacement on the initial model to obtain a replacement model, and calculating the longitudinal wave velocity, the transverse wave velocity and the density of the replacement model; respectively calculating a plurality of fluid identification factors of the initial model and the replacement model according to the initial longitudinal wave velocity, the initial transverse wave velocity, the initial density and the longitudinal wave velocity, the transverse wave velocity and the density of the replacement model; and calculating the sensitivity corresponding to each fluid identification factor according to a fluid factor sensitivity formula. The method determines the fluid identification factor and the sensitivity thereof through Gassmann fluid replacement, preferably selects the fluid identification factor with high sensitivity, is favorable for accurately reflecting the characteristics of the reservoir stratum, finds the reservoir stratum containing oil or gas and provides service for the exploration and development of conventional oil gas, coal bed gas and shale oil gas.)

基于流体替换的井中流体因子敏感度计算方法及系统

技术领域

本发明涉及油气、页岩油气及煤层气地震勘探与开发领域，更具体地，涉及一种基于流体替换的井中流体因子敏感度计算方法及系统。

背景技术

地下岩石实际上是由固体矿物和孔隙流体组成的多孔双相介质，孔隙流体的存在必然会影响多孔岩石的力学特征，从而使介质岩石表现出一定的弹性异常，在地震岩石物理模型构建及储层流体地震波响应分析基础上，可将与孔隙流体有关的异常特性表征为流体因子，因此，如何从构建的流体因子中优选出敏感的流体因子是关键。

在油气储层的识别研究中，大都通过获得储层AVO信息或地震资料信息中三个基本弹性参数(纵波速度、横波速度、密度)进行流体识别，而在实际应用中，对于不同流体类型的储层，这三个参数的应用有限，有时不是很敏感等，不能直接有效地进行流体识别分析。为此，近年来国内外的专家为了有效突出不同流体的物性差异，对上述的纵波速度、横波速度、密度进行了各种组合，如纵波阻抗、横波阻抗、泊松比等，这些组合的流体因子在一定程度上提高了流体识别的敏感程度。但是，面对数量日益增多的流体识别因子，在实际应用中，如何选择最灵敏的流体识别因子，是目前急需解决的重要问题。因此，有必要开发一种基于流体替换的井中流体因子敏感度计算方法及系统。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种基于流体替换的井中流体因子敏感度计算方法及系统，其能够通过Gassmann流体替换，确定流体识别因子及其敏感性，优选敏感性高的流体识别因子，有助于精确反应储层特征，寻找含油或者含气的储层，为常规油气、煤层气及页岩油气的勘探、开发服务。

根据本发明的一方面，提出了一种基于流体替换的井中流体因子敏感度计算方法。所述方法可以包括：根据包含初始纵波速度、初始横波速度、初始密度的测井曲线，删除所述测井曲线中的数据中的异常值，确定储层位置，建立初始模型；对所述初始模型进行Gassmann流体替换，获得替换模型，计算所述替换模型的纵波速度、横波速度、密度；根据所述初始纵波速度、所述初始横波速度、所述初始密度与所述替换模型的纵波速度、横波速度、密度，分别计算初始模型与替换模型的多个流体识别因子；根据流体因子敏感性公式，计算每个流体识别因子对应的敏感性。

优选地，所述替换模型的纵波速度为：

其中，V_p为替换模型的纵波速度，K_sat为替换模型的体积模量，μ_sat为替换模型的剪切模量，ρ_sat为替换模型的密度。

优选地，所述替换模型的横波速度为：

其中，V_s为替换模型的横波速度，μ_sat为替换模型的剪切模量，ρ_sat为替换模型的密度。

优选地，所述替换模型的密度为：

ρ_sat＝(1-φ)ρ_ma+φρ_fl (3)

其中，ρ_sat为替换模型的密度，ρ_fl为流体密度，ρ_ma为骨架密度，φ为孔隙度。

优选地，通过公式(4)计算每个流体识别因子对应的敏感性：

其中，mean_water为水层的流体因子的均值，mean_gasoil为油气层的流体因子的均值，std_gasoil为油气层流体因子的方差。

根据本发明的另一方面，提出了一种基于流体替换的井中流体因子敏感度计算系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：根据包含初始纵波速度、初始横波速度、初始密度的测井曲线，删除所述测井曲线中的数据中的异常值，确定储层位置，建立初始模型；对所述初始模型进行Gassmann流体替换，获得替换模型，计算所述替换模型的纵波速度、横波速度、密度；根据所述初始纵波速度、所述初始横波速度、所述初始密度与所述替换模型的纵波速度、横波速度、密度，分别计算初始模型与替换模型的多个流体识别因子；根据流体因子敏感性公式，计算每个流体识别因子对应的敏感性。

优选地，所述替换模型的纵波速度为：

其中，V_p为替换模型的纵波速度，K_sat为替换模型的体积模量，μ_sat为替换模型的剪切模量，ρ_sat为替换模型的密度。

优选地，所述替换模型的横波速度为：

其中，V_s为替换模型的横波速度，μ_sat为替换模型的剪切模量，ρ_sat为替换模型的密度。

优选地，所述替换模型的密度为：

ρ_sat＝(1-φ)ρ_ma+φρ_fl (3)

其中，ρ_sat为替换模型的密度，ρ_fl为流体密度，ρ_ma为骨架密度，φ为孔隙度。

优选地，通过公式(4)计算每个流体识别因子对应的敏感性：

其中，mean_water为水层的流体因子的均值，mean_gasoil为油气层的流体因子的均值，std_gasoil为油气层流体因子的方差。

本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的

具体实施方式

中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。

附图说明

通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。

图1示出了根据本发明的基于流体替换的井中流体因子敏感度计算方法的步骤的流程图。

图2示出了根据本发明的一个实施例的18种流体识别因子的敏感性的直方图。

具体实施方式

下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。

图1示出了根据本发明的基于流体替换的井中流体因子敏感度计算方法的步骤的流程图。

在该实施例中，根据本发明的基于流体替换的井中流体因子敏感度计算方法可以包括：步骤101，根据包含初始纵波速度、初始横波速度、初始密度的测井曲线，删除测井曲线中的数据中的异常值，确定储层位置，建立初始模型；步骤102，对初始模型进行Gassmann流体替换，获得替换模型，计算替换模型的纵波速度、横波速度、密度；步骤103，根据初始纵波速度、初始横波速度、初始密度与替换模型的纵波速度、横波速度、密度，分别计算初始模型与替换模型的多个流体识别因子；以及步骤104，根据流体因子敏感性公式，计算每个流体识别因子对应的敏感性。

在一个示例中，替换模型的纵波速度为：

其中，V_p为替换模型的纵波速度，K_sat为替换模型的体积模量，μ_sat为替换模型的剪切模量，ρ_sat为替换模型的密度。

在一个示例中，替换模型的横波速度为：

其中，V_s为替换模型的横波速度，μ_sat为替换模型的剪切模量，ρ_sat为替换模型的密度。

在一个示例中，替换模型的密度为：

ρ_sat＝(1-φ)ρ_ma+φρ_fl (3)

其中，ρ_sat为替换模型的密度，ρ_fl为流体密度，ρ_ma为骨架密度，φ为孔隙度。

在一个示例中，通过公式(4)计算每个流体识别因子对应的敏感性：

其中，mean_water为水层的流体因子的均值，mean_gasoil为油气层的流体因子的均值，std_gasoil为油气层流体因子的方差。

具体地，根据本发明的基于流体替换的井中流体因子敏感度计算方法可以包括：

根据包含初始纵波速度、初始横波速度、初始密度的测井曲线，删除测井曲线中的数据中的异常值，确定储层位置，输入储层的油气水含量，三者相加总含量为100％，油气水的弹性参数和模量参数(体积模量、剪切模量)，设置岩石的骨架和泥质特性，即岩石的弹性参数(纵波速度、横波速度、密度)和模量参数(体积模量、剪切模量)，获得初始模型。

设置Gassmann流体替换后储层孔隙中的油气水的含量，三者相加总含量为100％吗，根据Gassmann流体替换公式和Wood模型对初始模型进行Gassmann流体替换，获得替换模型，通过公式(1)计算替换模型的纵波速度，通过公式(2)计算替换模型的横波速度，通过公式(3)计算替换模型的密度。

根据初始纵波速度、初始横波速度、初始密度与替换模型的纵波速度、横波速度、密度，通过流体识别因子计算公式，分别计算初始模型与替换模型的多个流体识别因子。

根据流体因子敏感性公式，通过公式(4)计算每个流体识别因子对应的敏感性，敏感度的数值越大，说明流体识别因子的敏感性越好，区分这两种流体的能力越强，可以将优选出来流体识别因子，应用于同工区中其它未知的测井。

本方法通过Gassmann流体替换，确定流体识别因子及其敏感性，优选敏感性高的流体识别因子，有助于精确反应储层特征，寻找含油或者含气的储层，为常规油气、煤层气及页岩油气的勘探、开发服务。

应用示例

为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。

根据本发明的基于流体替换的井中流体因子敏感度计算方法可以包括：

根据包含初始纵波速度、初始横波速度、初始密度的测井曲线，删除测井曲线中的数据中的异常值，确定储层位置，输入储层的油气水含量，三者相加总含量为100％，油气水的弹性参数和模量参数(体积模量、剪切模量)，设置岩石的骨架和泥质特性，即岩石的弹性参数(纵波速度、横波速度、密度)和模量参数(体积模量、剪切模量)，获得初始模型。

设置Gassmann流体替换后储层孔隙中的油气水的含量，三者相加总含量为100％吗，根据Gassmann流体替换公式和Wood模型对初始模型进行Gassmann流体替换，获得替换模型，通过公式(1)计算替换模型的纵波速度，通过公式(2)计算替换模型的横波速度，通过公式(3)计算替换模型的密度。

根据初始纵波速度、初始横波速度、初始密度与替换模型的纵波速度、横波速度、密度，通过流体识别因子计算公式，分别计算初始模型与替换模型的18个流体识别因子，如表1所示。

表1

序号	流体因子	序号	流体因子
				1	纵波速度	10	拉梅系数泊松比
2	纵波阻抗	11	Gassman流体项泊松比
				3	拉梅系数密度	12	阻抗比
4	剪切模量密度	13	Gassman流体项孔隙度
				5	Gassmann流体项	14	密度
6	泊松比	15	阻抗差
				7	横波阻抗	16	Russell因子
8	横波速度	17	HSFIF流体因子
				9	流体压缩系数	18	泊松阻抗

其中，纵波速度：V_p；纵波阻抗：I_p＝ρ×V_p；第一拉梅系数密度λρ：

第二拉梅系数密度μρ：

Gassman流体项：

c为干岩石纵横波速度平方比；泊松比：横波阻抗：I_s＝ρ×V_s；横波速度：V_s；流体压缩系数：第一拉梅系数泊松比λσ：Gassman流体项泊松比fσ：

为干岩石纵横波速度平方比；阻抗比(纵横波速度比)：

Gassman流体项孔隙度K_f：

φ_c临界孔隙度；密度：ρ；阻抗差：I_p-I_s；Russell因子ρf：

c为调节参数；HSFIF流体因子σ_HSFIF：

为调节参数；泊松阻抗：PI＝ρ(V_p-cV_s)。

图2示出了根据本发明的一个实施例的18种流体识别因子的敏感性的直方图。

根据流体因子敏感性公式，通过公式(4)计算每个流体识别因子对应的敏感性，如图2所示，敏感度的数值越大，说明流体识别因子的敏感性越好，区分这两种流体的能力越强，可以将优选出来流体识别因子，应用于同工区中其它未知的测井。

综上所述，本发明通过Gassmann流体替换，确定流体识别因子及其敏感性，优选敏感性高的流体识别因子，有助于精确反应储层特征，寻找含油或者含气的储层，为常规油气、煤层气及页岩油气的勘探、开发服务。

本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。

根据本发明的实施例，提供了一种基于流体替换的井中流体因子敏感度计算系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：根据包含初始纵波速度、初始横波速度、初始密度的测井曲线，删除测井曲线中的数据中的异常值，确定储层位置，建立初始模型；对初始模型进行Gassmann流体替换，获得替换模型，计算替换模型的纵波速度、横波速度、密度；根据初始纵波速度、初始横波速度、初始密度与替换模型的纵波速度、横波速度、密度，计算多个流体识别因子；根据流体因子敏感性公式，计算每个流体识别因子对应的敏感性。

在一个示例中，替换模型的纵波速度为：

其中，V_p为替换模型的纵波速度，K_sat为替换模型的体积模量，μ_sat为替换模型的剪切模量，ρ_sat为替换模型的密度。

在一个示例中，替换模型的横波速度为：

其中，V_s为替换模型的横波速度，μ_sat为替换模型的剪切模量，ρ_sat为替换模型的密度。

在一个示例中，替换模型的密度为：

ρ_sat＝(1-φ)ρ_ma+φρ_fl (3)

其中，ρ_sat为替换模型的密度，ρ_fl为流体密度，ρ_ma为骨架密度，φ为孔隙度。

在一个示例中，通过公式(4)计算每个流体识别因子对应的敏感性：

其中，mean_water为水层的流体因子的均值，mean_gasoil为油气层的流体因子的均值，std_gasoil为油气层流体因子的方差。

本系统通过Gassmann流体替换，确定流体识别因子及其敏感性，优选敏感性高的流体识别因子，有助于精确反应储层特征，寻找含油或者含气的储层，为常规油气、煤层气及页岩油气的勘探、开发服务。

本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。

以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。