CN111485850A - 基于大孔道数据的油井堵水方法及装置 - Google Patents

Abstract

本发明公开了一种基于连通的大孔道数据的油井堵水方法及装置，属于石油技术领域，该方法包括：根据取芯井岩芯样品的毛管压力数据，计算不同毛管压力下的孔隙半径，对孔隙半径及岩芯样品对应的渗透率、孔隙度进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型；根据大孔道行业标准、岩芯样品孔隙半径大小频率分布、多个孔隙半径和孔渗关系模型的相关性和二八定律，确定需要堵水的孔隙半径标准值及目标关系模型；采用目标关系模型，计算目标井各层位孔隙半径平均值，与孔隙半径标准值进行对比，确定目标井堵水层位，进而计算堵剂量。本发明能准确的定位堵水位置，提高了堵水的准确性；大大降低堵剂用量，节约成本的同时避免地层污染的风险。

Description

基于大孔道数据的油井堵水方法及装置

技术领域

本发明涉及石油技术领域，特别涉及一种基于大孔道数据的油井堵水方法及装置。

背景技术

在油藏开发过程中，注水开发作为最经济有效的开发方式被广泛应用。注水开发过程中，注入水在部分层位率先突破，在油水井之间形成水流优势通道，注入水沿水流优势通道直接进入油井，水驱波及面积减少，造成注水低效现象。堵水技术是改善水驱增油效果的重要手段，通过封堵水流优势通道，改变液流方向，增加水驱波及体积，提高驱油面积。堵水技术发展迅速，但如何准确确定堵水层位、堵剂用量是一项技术难题。

目前现场实施的堵水方法通常有两种，一种主要依据估算和经验制定，如依据油井生产数据和油层物性的对比，推测出水层位深度；利用孔隙度、有效厚度值计算目标层位地下孔隙体积，根据经验确定堵水剂用量。另一种采用油藏数值模拟技术，模拟整个油藏的生产过程，确定出水层位和目标层的孔隙体积。

上述第一种计算方法误差较大，易造成堵剂太少，堵水不成功，或者堵剂太多，造成地层污染；第二种方法需要全方位收集区块的基础资料，开展建模、数模工作，并对全区及单井的生产历史进行拟合，拟合精度达到规范后方可实施，工作量巨大。

发明内容

本发明实施例提供了一种油井抽油机的工作参数配置方法及装置，能够解决现有技术的问题，所述技术方案如下：

一方面，提供了一种基于连通的大孔道数据的油井堵水方法，所述方法包括：

对取芯井的岩芯样品进行压汞测试试验，获取所述岩芯样品的毛管压力数据，通过公式计算不同毛管压力下，汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径；

根据计算的最小孔隙半径，以及与其对应的所述岩芯样品的孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型，所述孔渗为孔隙度和渗透率的合称，所述孔渗测试数据包括所述岩芯样品不同层位的孔隙度和渗透率数据；

根据大孔道行业标准，确定油藏开发过程中需要堵水的孔隙半径标准值；

根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，从所述多个孔隙半径和孔渗关系模型中筛选出目标关系模型，所述目标关系模型用于计算层位前20％孔隙的孔隙半径值；

根据目标井的测井报告中不同层位的孔渗数据，采用所述目标关系模型，计算所述目标井不同层位的孔隙半径平均值；

基于所述目标井不同层位的孔隙半径平均值与所述需要堵水的孔隙半径标准值的对比，确定所述目标井的目标堵水层位；其中，所述目标堵水层位的孔隙半径平均值不小于所述需要堵水的孔隙半径标准值；

根据所述目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算需要的堵剂量。

在一种可能的实现方式中，对取芯井的岩芯样品进行压汞测试试验，获取所述岩芯样品的毛管压力数据，通过公式计算不同毛管压力下，汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径，包括：

在压汞测试中，针对所述岩芯样品的不同层位，获取所述层位下汞充满所述不同预设百分比时的毛管压力大小；

针对每个层位对应的不同毛管压力值，根据孔隙半径公式，计算所述毛管压力值对应的孔隙半径值，则所述孔隙半径值即为汞充满所述岩芯样品相应预设百分比时的最小孔隙半径。

在一种可能的实现方式中，根据计算的最小孔隙半径，以及与其对应的所述岩芯样品的孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型，包括:

根据计算的汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径，结合所述岩芯样品不同层位的孔隙度和渗透率，得到多组孔隙半径和孔渗数据；

针对每一组孔隙半径和孔渗数据进行多元线性回归处理，建立相应预设百分比下的孔隙半径和孔渗关系模型，得到多个孔隙半径和孔渗关系模型。

在一种可能的实现方式中，根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，从所述多个孔隙半径和孔渗关系模型中筛选出目标关系模型，包括：

根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，确定需要堵水的判别条件，所述判别条件为层位前20％孔隙的孔隙半径平均值大于或等于所述需要堵水的孔隙半径标准值；

基于所述判别条件，将所述多个孔隙半径和孔渗关系模型中孔隙半径占比在前20％的孔隙半径和孔渗关系模型确定为所述目标关系模型。

在一种可能的实现方式中，根据目标井的测井报告中不同层位的孔渗数据，采用所述目标关系模型，计算所述目标井不同层位的孔隙半径平均值，包括：

根据所述目标井的测井报告，获取所述目标井不同层位的孔隙度和渗透率值；

根据所述目标井每一个层位的孔隙度和渗透率，采用目标关系模型，计算所述层位相对于所述目标关系模型的孔隙半径值；

当目标关系模型数目为1个时，将所述目标关系模型的孔隙半径值确定为所述层位的孔隙半径平均值；当目标关系模型数目为多个时，计算多个目标关系模型对应的孔隙半径值的均值，将所述均值作为所述层位的孔隙半径平均值。

另一方面，提供了一种基于连通的大孔道数据的油井堵水装置，所述装置包括：

获取模块，用于对取芯井的岩芯样品进行压汞测试试验，获取所述岩芯样品的毛管压力数据，通过公式计算不同毛管压力下，汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径；

模型建立模块，用于根据计算的最小孔隙半径，以及与其对应的所述岩芯样品的孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型，所述孔渗为孔隙度和渗透率的合称，所述孔渗测试数据包括所述岩芯样品不同层位的孔隙度和渗透率数据；

孔隙半径标准值确定模块，用于根据大孔道行业标准，确定油藏开发过程中需要堵水的孔隙半径标准值；

筛选模块，用于根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，从所述多个孔隙半径和孔渗关系模型中筛选出目标关系模型，所述目标关系模型用于计算层位前20％孔隙的孔隙半径值；

孔隙半径计算模块，用于根据目标井的测井报告中不同层位的孔渗数据，采用所述目标关系模型，计算所述目标井不同层位的孔隙半径平均值；

目标堵水层位确定模块，用于基于所述目标井不同层位的孔隙半径平均值与所述需要堵水的孔隙半径标准值的对比，确定所述目标井的目标堵水层位；其中，所述目标堵水层位的孔隙半径平均值不小于所述需要堵水的孔隙半径标准值；

堵剂量计算模块，用于根据所述目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算需要的堵剂量。

在一种可能的实现方式中，获取模块，用于：

在压汞测试中，针对所述岩芯样品的不同层位，获取所述层位下汞充满所述不同预设百分比时的毛管压力大小；

针对每个层位对应的不同毛管压力值，根据孔隙半径公式，计算所述毛管压力值对应的孔隙半径值，则所述孔隙半径值即为汞充满所述岩芯样品相应预设百分比时的最小孔隙半径。

在一种可能的实现方式中，模型建立模块，用于:

根据计算的汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径，结合所述岩芯样品不同层位的孔隙度和渗透率，得到多组孔隙半径和孔渗数据；

针对每一组孔隙半径和孔渗数据进行多元线性回归处理，建立相应预设百分比下的孔隙半径和孔渗关系模型，得到多个孔隙半径和孔渗关系模型。

在一种可能的实现方式中，筛选模块，用于：

根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，确定需要堵水的判别条件，所述判别条件为层位前20％孔隙的孔隙半径平均值大于或等于所述需要堵水的孔隙半径标准值；

基于所述判别条件，将所述多个孔隙半径和孔渗关系模型中孔隙半径占比在前20％的孔隙半径和孔渗关系模型确定为所述目标关系模型。

在一种可能的实现方式中，孔隙半径计算模块，用于：

根据所述目标井的测井报告，获取所述目标井不同层位的孔隙度和渗透率值；

根据所述目标井每一个层位的孔隙度和渗透率，采用目标关系模型，计算所述层位相对于所述目标关系模型的孔隙半径值；

当目标关系模型数目为1个时，将所述目标关系模型的孔隙半径值确定为所述层位的孔隙半径平均值；当目标关系模型数目为多个时，计算多个目标关系模型对应的孔隙半径值的均值，将所述均值作为所述层位的孔隙半径平均值。

本发明实施例提供的技术方案带来的有益效果是：

对岩芯样品进行压汞测试试验，获取毛管压力数据，通过公式计算不同毛管压力下，汞充满岩芯样品不同预设百分比时分别对应的最小孔隙半径；将计算的最小孔隙半径，以及对应的岩芯样品孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型；根据大孔道行业标准，确定需要堵水的孔隙半径标准值；根据岩芯样品孔隙半径大小频率分布、多个孔隙半径和孔渗关系模型的相关性和二八定律，筛选出目标关系模型，目标关系模型用于计算层位前20％孔隙的孔隙半径值；根据目标井的测井报告中的孔渗数据，通过目标关系模型计算前20％孔隙的平均半径值；将目标井不同层位的孔隙半径平均值和需要堵水的孔隙半径标准值进行对比，确定目标堵水层位；根据该目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算所需堵剂量。本发明通过水流优势通道的刻画，可以定位参与流动的大孔道的位置，得到实施堵水的准确深度；提出只需要封堵连通的大孔隙的方法，大大降低了堵剂用量，节约了堵水成本。而且，过多的注入堵剂很容易造成泄漏，造成环境污染，本发明提供的方法准确计算堵剂用量，避免了污染的风险，为实际生产提供了简便、有效的技术手段。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例提供的一种基于大孔道数据的油井堵水方法的流程图；

图2是本发明实施例提供的一种基于大孔道数据的油井堵水方法的流程图；

图3是本发明实施例提供的一种对A井进行压汞测试的压汞曲线图；

图4是本发明实施例提供的一种多元线性回归处理结果图；

图5是本发明实施例提供的一种A井岩芯孔隙半径分布直方图；

图6是本发明实施例提供的一种目标井的测井报告图；

图7是本发明实施例提供的一种基于大孔道数据的油井堵水装置的结构示意图；

图8示出了本发明一个示例性实施例提供的计算机设备800的结构框图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。

图1是本发明实施例提供的一种基于连通的大孔道数据的油井堵水方法的流程图。参见图1，该方法包括：

101、对取芯井的岩芯样品进行压汞测试试验，获取岩芯样品的毛管压力数据，通过公式计算不同毛管压力下，汞充满岩芯样品不同预设百分比时分别对应的最小孔隙半径。

102、根据计算的最小孔隙半径，以及与其对应的岩芯样品的孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型。

其中，孔渗为孔隙度和渗透率的合称，孔渗测试数据包括岩芯样品不同层位的孔隙度和渗透率数据。

103、根据大孔道行业标准，确定油藏开发过程中需要堵水的孔隙半径标准值。

104、根据岩芯样品孔隙半径大小频率分布、多个孔隙半径和孔渗关系模型的相关性和二八定律，从多个孔隙半径和孔渗关系模型中筛选出目标关系模型，目标关系模型用于计算层位前20％孔隙的孔隙半径值。

105、根据目标井的测井报告中不同层位的孔渗数据，采用目标关系模型，计算目标井不同层位的孔隙半径平均值。

106、基于目标井不同层位的孔隙半径平均值与需要堵水的孔隙半径标准值的对比，确定目标井的目标堵水层位；

其中，目标堵水层位的孔隙半径平均值不小于需要堵水的孔隙半径标准值。

107、根据目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算需要的堵剂量。

本发明实施例提供的方法，对岩芯样品进行压汞测试试验，获取毛管压力数据，通过公式计算不同毛管压力下，汞充满岩芯样品不同预设百分比时分别对应的最小孔隙半径；将计算的最小孔隙半径，以及对应的岩芯样品孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型；根据大孔道行业标准，确定需要堵水的孔隙半径标准值；根据岩芯样品孔隙半径大小频率分布、多个孔隙半径和孔渗关系模型的相关性和二八定律，筛选出目标关系模型，目标关系模型用于计算层位前20％孔隙的孔隙半径值；根据目标井的测井报告中的孔渗数据，通过目标关系模型计算前20％孔隙的平均半径值；将目标井不同层位的孔隙半径平均值和需要堵水的孔隙半径标准值进行对比，确定目标堵水层位；根据该目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算所需堵剂量。本发明通过水流优势通道的刻画，可以定位参与流动的大孔道的位置，得到实施堵水的准确深度；提出只需要封堵连通的大孔隙的方法，大大降低了堵剂用量，节约了堵水成本。而且，过多的注入堵剂很容易造成泄漏，造成环境污染，本发明提供的方法准确计算堵剂用量，避免了污染的风险，为实际生产提供了简便、有效的技术手段。

在一种可能的实现方式中，对取芯井的岩芯样品进行压汞测试试验，获取该岩芯样品的毛管压力数据，通过公式计算不同毛管压力下，汞充满该岩芯样品不同预设百分比时分别对应的最小孔隙半径，包括：

在压汞测试中，针对该岩芯样品的不同层位，获取该层位下汞充满该不同预设百分比时的毛管压力大小；

针对每个层位对应的不同毛管压力值，根据孔隙半径公式，计算该毛管压力值对应的孔隙半径值，则该孔隙半径值即为汞充满该岩芯样品相应预设百分比时的最小孔隙半径。

在一种可能的实现方式中，根据计算的最小孔隙半径，以及与其对应的该岩芯样品的孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型，包括:

根据计算的汞充满该岩芯样品不同预设百分比时分别对应的最小孔隙半径，结合该岩芯样品不同层位的孔隙度和渗透率，得到多组孔隙半径和孔渗数据；

针对每一组孔隙半径和孔渗数据进行多元线性回归处理，建立相应预设百分比下的孔隙半径和孔渗关系模型，得到多个孔隙半径和孔渗关系模型。

在一种可能的实现方式中，根据该岩芯样品孔隙半径大小频率分布、该多个孔隙半径和孔渗关系模型的相关性和二八定律，从该多个孔隙半径和孔渗关系模型中筛选出目标关系模型，包括：

根据该岩芯样品孔隙半径大小频率分布、该多个孔隙半径和孔渗关系模型的相关性和二八定律，确定需要堵水的判别条件，该判别条件为层位前20％孔隙的孔隙半径平均值大于或等于需要堵水的孔隙半径标准值；

基于该判别条件，将该多个孔隙半径和孔渗关系模型中孔隙半径占比在前20％的孔隙半径和孔渗关系模型确定为目标关系模型。

在一种可能的实现方式中，根据目标井的测井报告中不同层位的孔渗数据，采用目标关系模型，计算该目标井不同层位的孔隙半径平均值，包括：

根据该目标井的测井报告，获取该目标井不同层位的孔隙度和渗透率值；

根据该目标井每一个层位的孔隙度和渗透率，采用目标关系模型，计算该层位相对于目标关系模型的孔隙半径值；

当目标关系模型数目为1个时，将目标关系模型的孔隙半径值确定为该层位的孔隙半径平均值；当目标关系模型数目为多个时，计算多个目标关系模型对应的孔隙半径值的均值，将该均值作为该层位的孔隙半径平均值。

上述所有可选技术方案，可以采用任意结合形成本公开的可选实施例，在此不再一一赘述。

图2是本发明实施例提供的一种基于连通的大孔道数据的油井堵水方法的流程图。参见图2，该方法包括：

201、对取芯井的岩芯样品进行压汞测试，获取该岩芯样品的毛管压力数据，通过公式计算出不同毛管压力下，汞充满岩芯样品不同预设百分比时分别对应的最小孔隙半径。

该步骤201包括以下流程：在压汞测试中，针对岩芯样品的不同层位，获取层位下汞充满不同预设百分比时的毛管压力大小；针对每个层位对应的不同毛管压力值，根据孔隙半径公式，计算毛管压力值对应的孔隙半径值，则孔隙半径值即为汞充满岩芯样品相应预设百分比时的最小孔隙半径。

具体地，对取芯井进行压汞测试，针对岩芯样品的不同层位，根据压汞测试数据可以得到以下关系：(1)岩芯压汞测试中，汞会率先进入岩芯中的大孔隙，随着压力的增加，再充满岩芯中的中小孔隙。压汞测试过程中，可以获取到汞率先充满5％、10％、15％、20％、25％、30％、35％、40％、45％、50％岩芯时毛管压力的大小；(2)根据毛管压力与孔隙半径的公式，计算不同毛管压力下，汞充满5％、10％、15％、20％、25％、30％、35％、40％、45％、50％岩芯时的最小孔隙半径r5、r10、r15、r20、r25、r30、r35、r40、r45、r50。

其中孔隙半径是通过压汞过程中的毛管压力计算出来；压汞测试过程中，压力不断增加，汞是逐渐按顺序充满岩芯中的大、中、小孔隙，当汞充满5％岩芯时，显示的毛管压力值是充满5％岩芯的最大压力，所以计算出来的数据r5是指岩芯孔隙由大到小排列，前5％孔隙中的最小孔隙半径。

其中，压汞测试过程中，根据毛管力计算孔隙半径的孔隙半径公式如(1)所示。

r(d)＝2δcosθ/Pc 公式(1)

其中，Pc表示毛管压力，单位Pa(帕斯卡)；r(d)表示汞充满d％岩芯时的最小孔隙半径，单位mm(毫米)；θ表示润湿角，单位(度)；δ表示表面张力，单位mN/m(毫牛/米)。其中，压汞测试中还会得到润湿角、表面张力的值，这两种数据一般都是常规的数据，另外每次压汞测试中汞和空气的表面张力也基本上都是恒定的数据，也可以在压汞测试过程中获取。

需要说明的是，一般情况下需要计算孔隙半径的平均值，本发明采用最小孔隙半径值代替孔隙半径的平均值，当然，具体实施过程中，在求解汞充满岩芯样品不同预设百分比的孔隙半径时，也可以采用多个孔隙半径的平均值，例如，在求解汞充满5％的孔隙半径值时，计算汞充满1％、3％和5％的孔隙半径的平均值，本发明对孔隙半径平均值的求解方式不作限定。

压汞测试过程中，汞率先充满大孔隙然后再充满中小孔隙；如对取芯井A井D层位71块目的层岩芯进行压汞测试，所测压汞曲线如图3所示。基于该压汞曲线，通过上述公式(1)可计算出A井不同层位的岩芯，汞充满5％、10％、15％、20％、25％、30％、35％、40％、45％、50％岩芯时的最小孔隙半径r5、r10、r15、r20、r25、r30、r35、r40、r45、r50，如表1所示。

表1中记录了取芯井A井不同样品深度的渗透率、孔隙度、以及汞分别进入前5％、10％、15％、20％、25％、30％、35％、40％、45％和50％孔隙体积时的最小孔隙半径值。例如，表1中第1行(除去表头)、第6列r5um对应的单元格表示孔隙按照从大到小排列，压汞进入前5％大孔隙体积时的最小孔隙半径为5.65um(微米)。需要说明的是，岩芯样品的渗透率、孔隙度的值都是实验测得的，其中，孔隙度是通过液体饱和法测得，能够得到岩芯样品不同层位对应的孔隙度。渗透率是通过达西定律的方法测得，能够得到岩芯样品不同层位对应的渗透率。

表1

序号	层位	样品深度m	渗透率md	孔隙度％	r5μm	r10μm	r15μm	r20μm	r25μm	r30μm	r35μm	r40μm	r45μm	r50μm
															1	Ⅳ_D	3040.2	20.6	19.7	5.65	4.60	3.74	3.04	2.48	2.01	1.64	1.33	1.08	0.59
2	Ⅳ_D	3041.6	311.3	23.7	16.44	13.70	11.42	9.52	7.94	6.62	5.51	4.60	3.83	1.89
															3	Ⅳ_D	3042.6	492.6	24.7	17.44	15.26	13.35	11.68	10.22	8.95	7.83	6.85	5.99	4.04
4	Ⅳ_D	3043.5	621.8	23.3	23.03	19.60	16.67	14.19	12.07	10.27	8.74	7.44	6.33	3.84
															5	Ⅳ_D	3044.5	590.7	22.7	17.63	15.97	14.47	13.12	11.89	10.77	9.76	8.85	8.02	7.07
6	Ⅳ_D	3045.3	273.8	21.5	13.17	11.67	10.34	9.16	8.12	7.19	6.37	5.65	5.00	3.89
															7	Ⅳ_D	3046.2	530.1	21.4	18.69	16.67	14.86	13.26	11.82	10.54	9.40	8.39	7.48	6.43
8	Ⅳ_D	3047.4	334.1	21.5	14.14	12.80	11.58	10.48	9.48	8.58	7.76	7.02	6.36	5.33
															9	Ⅳ_D	3048.1	442.2	22.0	15.33	13.97	12.73	11.60	10.57	9.63	8.77	7.99	7.28	5.66
10	Ⅳ_D	3048.8	395.1	23.0	14.63	12.78	11.16	9.75	8.51	7.43	6.49	5.67	4.95	2.67
															11	Ⅳ_D	3049.9	783.1	24.4	18.18	16.41	14.81	13.37	12.07	10.89	9.83	8.87	8.01	6.59
12	Ⅳ_D	3050.6	28.4	18.3	5.55	4.71	4.00	3.39	2.88	2.44	2.07	1.76	1.49	1.16
															13	Ⅳ_D	3051.2	13.4	17.6	4.57	3.71	3.02	2.46	2.00	1.63	1.32	1.08	0.88	0.63
14	Ⅳ_D	3052.1	58.0	19.4	6.77	5.86	5.08	4.40	3.81	3.30	2.86	2.48	2.15	1.52
															15	Ⅳ_D	3053.6	149.9	21.9	9.88	8.74	7.74	6.85	6.06	5.37	4.75	4.20	3.72	3.00
16	Ⅳ_D	3054.4	87.3	20.8	9.20	7.80	6.61	5.61	4.75	4.03	3.42	2.90	2.46	1.54
															17	Ⅳ_D	3055.4	119.6	21.6	12.08	9.96	8.22	6.78	5.59	4.61	3.81	3.14	2.59	1.54
18	Ⅳ_D	3056.3	102.4	21.0	9.94	8.42	7.13	6.04	5.12	4.34	3.67	3.11	2.63	2.30
															19	Ⅳ_E	3094.1	170.6	20.0	8.62	7.28	6.15	5.19	4.38	3.70	3.13	2.64	2.23	1.31
20	Ⅳ_E	3096.8	526.1	22.3	20.27	17.84	15.71	13.83	12.17	10.71	9.43	8.30	7.31	6.65
															21	Ⅳ_E	3099.7	1420.0	24.6	24.97	22.56	20.38	18.42	16.64	15.03	13.58	12.27	11.09	10.68
22	Ⅳ_E	3100.4	721.6	23.4	15.28	14.15	13.10	12.13	11.23	10.40	9.63	8.91	8.25	7.78
															23	Ⅴ_F	3128.7	671.2	23.6	19.97	18.03	16.29	14.71	13.28	11.99	10.83	9.78	8.83	6.20
24	Ⅴ_F	3130.1	591.4	22.8	18.96	16.87	15.00	13.35	11.87	10.56	9.40	8.36	7.44	4.85
															25	Ⅴ_F	3130.6	2980.0	26.8	76.15	64.93	55.35	47.19	40.24	34.30	29.25	24.93	21.26	18.04
26	Ⅴ_F	3131.0	2220.0	24.2	41.14	36.13	31.72	27.86	24.46	21.48	18.86	16.56	14.54	12.79
															27	Ⅴ_F	3132.6	2220.0	22.2	41.65	33.61	27.12	21.89	17.66	14.25	11.50	9.28	7.49	5.69
28	Ⅴ_F	3133.6	93.1	17.3	19.87	13.61	9.32	6.38	4.37	2.99	2.05	1.40	0.96	0.62
															29	Ⅴ_F	3134.3	8.4	17.5	3.89	3.17	2.59	2.11	1.72	1.41	1.15	0.94	0.76	0.55
30	Ⅴ_F	3135.7	428.9	20.2	17.27	14.91	12.87	11.11	9.59	8.28	7.15	6.17	5.33	3.61
															31	Ⅴ_F	3136.8	858.5	23.6	23.28	20.61	18.24	16.14	14.29	12.65	11.20	9.91	8.77	7.82
32	Ⅴ_F	3137.2	1450.0	25.3	25.35	23.14	21.11	19.27	17.58	16.05	14.64	13.36	12.19	10.92
															33	Ⅴ_F	3138.4	810.9	23.4	22.93	20.30	17.97	15.90	14.08	12.46	11.03	9.76	8.64	6.29
34	Ⅴ_F	3140.3	234.6	21.4	15.54	13.34	11.44	9.82	8.43	7.23	6.21	5.33	4.57	2.63
															…	…	…										…	…	…

202、根据计算的最小孔隙半径，以及与其对应的岩芯样品的孔渗测试数据，对孔隙半径、渗透率和孔隙度进行多元线性回归处理，建立多个孔隙半径和孔渗的关系模型。

其中，孔渗为孔隙度和渗透率的合称，孔渗测试数据包括岩芯样品不同层位的孔隙度和渗透率数据。

对不同岩芯的r5、r10、r15、r20、r25、r30、r35、r40、r45、r50，与对应岩芯的渗透率、孔隙度进行多元线性回归，建立取芯井孔隙半径和孔隙度、渗透率的多个关系模型。

取芯井是对油田中最具有代表性的井进行取芯，所以取芯井的岩芯样品孔隙半径和孔隙度、渗透率的关系模型，可以代表整个油田的孔隙半径和孔隙度、渗透率的相互关系。

具体地，步骤202的过程包括：根据计算的汞充满岩芯样品不同预设百分比时分别对应的最小孔隙半径，结合岩芯样品不同层位的孔隙度和渗透率，得到多组孔隙半径和孔渗数据；针对每一组孔隙半径和孔渗数据进行多元线性回归处理，建立相应预设百分比下的孔隙半径和孔渗关系模型，得到多个孔隙半径和孔渗关系模型。

以A井为例，依据表1的数据，对A井不同层位的r5、r10、r15、r20、r25、r30、r35、r40、r45、r50与对应的渗透率和孔隙度进行多元线性回归处理，得到如图4所示的多元线性回归处理结果，包括10个孔隙半径和孔渗的关系模型，以及每个模型对应的回归相关系数。其中，每个关系模型中k表示渗透率，单位md、

表示孔隙度，％、R表示相关系数，r5表示孔隙半径从大到小排列，汞充满前5％孔隙体积时的最小孔隙半径，r10、r15等含义以此类推，在此不做赘述。

由图4可知，回归相关系数R在0.9678至0.9870之间，结果表明孔隙半径与渗透率和孔隙度相关程度高，所建立出的关系模型贴近实际情况，完全可用。

203、根据大孔道行业标准，确定油藏开发过程中需要堵水的孔隙半径标准值。

行业标准SY/T 6285-1997《油气储层评价方法》中描述砂岩和碳酸盐岩含油储层划分，依据孔隙半径中值r50，r50≥15μm的为大孔道储层。

油藏水驱开发过程中，注入水率先通过大孔道，完成驱油过程以后，形成水流高速通道，注入水直接进入油井，造成油井水淹；中小孔隙中原油未实现水驱过程，无需堵水。本发明意在实施高精度堵水，仅准确的堵住见水的大孔道，使注入水顺利驱替到中小孔隙，降低油井产水量的同时，提高中小孔隙中的原油动用程度。因此，本发明确定15μm为需要堵水的孔隙半径标准值，将大于或等于15μm做为筛选堵水孔隙的条件之一。

204、根据岩芯样品孔隙半径大小频率分布、多个孔隙半径和孔渗关系模型的相关性和二八定律，确定需要堵水的判别条件，该判别条件为层位前20％孔隙的孔隙半径平均值大于或等于需要堵水的孔隙半径标准值。

根据测试油藏水驱生产过程中，孤立的大孔隙不会造成水淹，只有参与流动的、连通的孔隙才能造成水淹，本发明通过筛选出此类大孔隙，实施准确堵水。行业标准中认为R50≥15μm为大孔道，本发明通过研究孔隙半径分布规律、孔隙半径和孔渗关系模型的相关性和二八定律，并结合需要堵水的孔隙半径标准值和堵水实践经验发现实际连通的、需要堵水的大孔道为r20≥15μm，也即需要堵水的判别条件为r20≥15μm，下面为具体论证过程：

虽然行业标准SY/T 6285-1997《油气储层评价方法》中描述砂岩和碳酸盐岩含油储层划分，依据孔隙半径中值R50，R50≥15μm的为大孔道储层。而从堵水实践中认识到，只有部分连通的大孔道才会造成水淹，远远到不了50％。从A井岩芯孔喉半径大小频率分布图来看，如图5，孔喉半径大于15μm的连通的大孔道，只占该岩芯的17％，即对A井来说r17≥15μm，为大孔道储层，考虑实际操作中的误差，界定研究区内连通的大孔道标准为r20≥15μm。

其中，二八定律是19世纪末20世纪初意大利经济学家帕累托发现的。他认为，在任何一组东西中，最重要的只占其中一小部分，约20％，其余80％尽管是多数，却是次要的。应用此理论，油井的水驱开发过程中，油井的产油、出水主要来自于前20％的孔隙；尤其水驱过程中的高速水流通道，主要来自于孔隙中前20％的大孔隙，因此计算前20％的孔隙大小，即可判断是否需要堵水。

另外，如图4孔隙半径和孔渗关系模型的回归相关系数也显示，在r20附近回归相关系数最高，说明r20时，最贴近油藏的驱替情况，因此通过计算前20％孔隙大小来判断是否需要堵水的方式是合理的，因此本发明将层位前20％孔隙的孔隙半径平均值大于或等于需要堵水的孔隙半径标准值作为需要堵水的判别条件。

205、基于上述判别条件，将所述多个孔隙半径和孔渗关系模型中孔隙半径占比在前20％的孔隙半径和孔渗关系模型确定为所述目标关系模型。

以上论证认为，计算前20％孔隙的孔隙半径大小即可用来判断是否需要堵水，为使发明更具操作性并降低实施风险，本发明采用前20％孔隙半径的平均值与堵水孔隙标准值进行对比，来判断是否需要堵水。因此，本发明将前20％孔隙半径对应的孔隙半径和孔渗关系模型确定为目标关系模型，也即从图4所示的多个孔隙半径和孔渗关系模型中筛选出r5、r10、r15、r20关系模型作为本发明堵水方案所采用的目标关系模型，以A井为例，最终筛选出的目标关系模型包括r5、r10、r15、r20所对应的关系模型。

其中，目标关系模型用于计算目标井不同层位前20％孔隙的孔隙半径值。后续，根据目标关系模型计算目标井前20％孔隙半径平均值，将孔隙半径平均值与需要堵水的孔隙半径标准值进行对比，以确定目标堵水层位，详细过程参见步骤206和步骤207。

206、根据目标井的测井报告中不同层位的孔渗数据，采用目标关系模型，计算目标井不同层位的孔隙半径平均值。

其中，目标井是指实际油藏中的生产井。目标井的测井报告包括该目标井不同层位的深度、孔隙度和渗透率的对应数据，如图6所示为目标井的测井报告。本发明通过目标关系模型，计算目标井不同深度(层位)下的孔隙半径平均值。

其中，目标关系模型是通过对取芯井的岩芯样品压汞测试实验得到的，基于取芯井的岩芯样品代表性，此目标关系模型可以应用到该油田的所有目标井。

具体过程包括：根据目标井的测井报告，获取目标井不同层位的孔隙度和渗透率值；根据每一个层位的孔隙度和渗透率，采用目标关系模型，计算该层位相对于目标关系模型的孔隙半径值；当目标关系模型数目为1个时，将所述目标关系模型的孔隙半径值确定为所述层位的孔隙半径平均值；当目标关系模型数目为多个时，计算多个目标关系模型对应的孔隙半径值的均值，将该均值作为该层位的孔隙半径平均值。

其中，当目标关系模型个数为多个时，计算每一个层位的孔隙半径的过程为：根据目标井每一个层位的孔隙度和渗透率，采用该多个目标关系模型，分别计算该层位相对于每一个目标关系模型的孔隙半径，得到该层位的多个孔隙半径，如r5、r10、r15、r20，然后通过算术平均值，计算前20％的孔隙半径的均值，从而得到该层位的孔隙半径平均值。

下面对上述过程进行举例说明，例如，步骤204中筛选出了r5、r10、r15、r20对应的关系模型为目标关系模型，针对目标井A1的每一个层位，将该层位的孔隙度和渗透率分别代入r5、r10、r15、r20对应的目标关系模型中，得到该层位相对于r5、r10、r15、r20关系模型的孔隙半径，将r5、r10、r15、r20对应的孔隙半径取算术平均值，确定为该层位的前20％的孔隙半径平均值，记为R20(平均)。如表2所示为目标井A不同层位R20(平均)信息表。表2中不同井段代表不同层位，第一行对应的编号为砂体为D1,第二行对应的砂体编号为D2，第三行对应的砂体编号为D3，以此类推，在此不做赘述。

表2

207、基于目标井不同层位的孔隙半径平均值与需要堵水的孔隙半径标准值的对比，确定所述目标井的目标堵水层位。

其中，目标堵水层位的孔隙半径平均值不小于需要堵水的孔隙半径标准值，需要堵水的孔隙半径标准值为孔隙半径大于等于15μm。

在本发明实施例中，将目标井不同层位的孔隙半径平均值与需要堵水的孔隙半径标准值进行比较，针对每一个层位，当该层位的孔隙半径平均值大于或等于需要堵水的孔隙半径标准值时，将该层位确定为目标堵水层位；当该层位的孔隙半径平均值小于需要堵水的孔隙半径标准值时，无需堵水，即该层位不是目标堵水层位。

仍以目标井A井为例，步骤203中确定需要堵水的孔隙半径标准值为15μm，由表2可知，R20(平均)满足不小于15μm条件的井段所在层位为砂体D1、D2、D6、D7，如表2中1、2、6、7行，因此，确定砂体D1、D2、D6、D7所指示层位为目标堵水层位，需要实施堵水。

208、根据目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算需要的堵剂量。

在本发明实施例中，根据目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，采用预设体积算法确定堵剂体积。其中，预设体积算法如公式(2)所示。

其中，V表示堵剂体积，单位m³(立方米)；π为圆周率；γ表示堵剂推进有效范围；h表示目标堵水层位的有效厚度，单位m(米)；

表示目标堵水层位的孔隙度，％。其中需要说明的是，该目标井的井距取平均井距，堵剂推进有效范围γ取1/4平均井距，单位m(米)；堵剂量采用堵剂体积表示。

具体地，针对目标井的每一个目标堵水层位，根据该目标堵水层位的厚度、孔隙度和堵剂推进有效范围，利用上述公式(2)计算该目标堵水层位的堵剂量。

仍以A井为例，采用公式(2)计算目标堵水层位D1、D2、D6、D7的堵剂体积，其中，A井的井距为800m，如表3所示为A1井目标堵水层位的堵剂体积结果表。

表3

下面结合上述步骤201至步骤208，对本发明的研究思路进行解释说明:本发明采用反演孔隙结构法，利用油藏的孔隙度、渗透率反推油藏孔隙半径。具体方法如下：利用实验室岩芯样品压汞测试报告中的毛管压力数据，计算汞充满岩芯样品不同预设百分比时分别对应的最小孔隙半径，结合对应的岩芯孔渗测试数据，通过多元线性回归处理建立孔隙半径和孔渗关系模型；进而根据大孔道行业标准、岩芯样品孔隙半径频率分布、孔隙半径和孔渗关系模型相关性和二八定律，筛选出目标关系模型；然后将该目标关系模型应用到目标井，目标井是指实际油藏中的生产井；根据油田内每口油井测井报告中目标层位的孔隙度和渗透率，利用目标关系模型，计算目标井各层位前20％孔隙半径平均值，当某层位前20％孔隙半径平均值大于或等于需要堵水的孔隙半径标准值时，确定该层位需要堵水；进而根据需要堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算所需堵剂量。本发明实施例提供的堵水方法论证了，只有连通的大孔道才会造成水淹，堵水只需针对连通的大孔道即可，不仅可以提高中小孔隙中的原油动用程度，而且大大降低了堵剂用量，降低了堵水成本。且过多的注入堵剂易发生泄露造成环境污染，本发明提供的方法，可简便、准确的计算堵剂量，提高经济效益，降低环境污染的风险。

本发明实施例提供的方法，对岩芯样品进行压汞测试试验，获取毛管压力数据，通过公式计算不同毛管压力下，汞充满岩芯样品不同预设百分比时分别对应的最小孔隙半径；将计算的最小孔隙半径，以及对应的岩芯样品孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型；根据大孔道行业标准，确定需要堵水的孔隙半径标准值；根据岩芯样品孔隙半径大小频率分布、多个孔隙半径和孔渗关系模型的相关性和二八定律，筛选出目标关系模型，目标关系模型用于计算层位前20％孔隙的孔隙半径值；根据目标井的测井报告中的孔渗数据，通过目标关系模型计算前20％孔隙的平均半径值；将目标井不同层位的孔隙半径平均值和需要堵水的孔隙半径标准值进行对比，确定目标堵水层位；根据该目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算所需堵剂量。本发明通过水流优势通道的刻画，可以定位参与流动的大孔道的位置，得到实施堵水的准确深度；提出只需要封堵连通的大孔隙的方法，大大降低了堵剂用量，节约了堵水成本。而且，过多的注入堵剂很容易造成泄漏，造成环境污染，本发明提供的方法准确计算堵剂用量，避免了污染的风险，为实际生产提供了简便、有效的技术手段。

图7是本发明实施例提供的一种基于大孔道数据的油井堵水装置的结构示意图，参见图7，所述装置包括：

获取模块701，用于对取芯井的岩芯样品进行压汞测试试验，获取岩芯样品的毛管压力数据，通过公式计算不同毛管压力下，汞充满岩芯样品不同预设百分比时分别对应的最小孔隙半径；

模型建立模块702，用于根据计算的最小孔隙半径，以及与其对应的岩芯样品的孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型，孔渗为孔隙度和渗透率的合称，孔渗测试数据包括岩芯样品不同层位的孔隙度和渗透率数据；

孔隙半径标准值确定模块703，用于根据大孔道行业标准，确定油藏开发过程中需要堵水的孔隙半径标准值；

筛选模块704，用于根据岩芯样品孔隙半径大小频率分布、多个孔隙半径和孔渗关系模型的相关性和二八定律，从多个孔隙半径和孔渗关系模型中筛选出目标关系模型，目标关系模型用于计算层位前20％孔隙的孔隙半径值；

孔隙半径计算模块705，用于根据目标井的测井报告中不同层位的孔渗数据，采用目标关系模型，计算目标井不同层位的孔隙半径平均值；

目标堵水层位确定模块706，用于基于目标井不同层位的孔隙半径平均值与需要堵水的孔隙半径标准值的对比，确定目标井的目标堵水层位；其中，目标堵水层位的孔隙半径平均值不小于需要堵水的孔隙半径标准值；

堵剂量计算模块707，用于根据目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算需要的堵剂量。

在一种可能的实现方式中，获取模块，用于：

在压汞测试中，针对该岩芯样品的不同层位，获取该层位下汞充满该不同预设百分比时的毛管压力大小；

针对每个层位对应的不同毛管压力值，根据孔隙半径公式，计算该毛管压力值对应的孔隙半径值，则该孔隙半径值即为汞充满该岩芯样品相应预设百分比时的最小孔隙半径。

在一种可能的实现方式中，模型建立模块，用于:

根据计算的汞充满该岩芯样品不同预设百分比时分别对应的最小孔隙半径，结合该岩芯样品不同层位的孔隙度和渗透率，得到多组孔隙半径和孔渗数据；

针对每一组孔隙半径和孔渗数据进行多元线性回归处理，建立相应预设百分比下的孔隙半径和孔渗关系模型，得到多个孔隙半径和孔渗关系模型。

在一种可能的实现方式中，筛选模块，用于：

根据该岩芯样品孔隙半径大小频率分布、该多个孔隙半径和孔渗关系模型的相关性和二八定律，确定需要堵水的判别条件，该判别条件为层位前20％孔隙的孔隙半径平均值大于或等于需要堵水的孔隙半径标准值；

基于该判别条件，将该多个孔隙半径和孔渗关系模型中孔隙半径占比在前20％的孔隙半径和孔渗关系模型确定为目标关系模型。

在一种可能的实现方式中，孔隙半径计算模块，用于：

根据该目标井的测井报告，获取该目标井不同层位的孔隙度和渗透率值；

根据该目标井每一个层位的孔隙度和渗透率，采用目标关系模型，计算该层位相对于目标关系模型的孔隙半径值；

当目标关系模型数目为1个时，将所述目标关系模型的孔隙半径值确定为所述层位的孔隙半径平均值；当目标关系模型数目为多个时，计算多个目标关系模型对应的孔隙半径值的均值，将所述均值作为所述层位的孔隙半径平均值。

本发明实施例提供的装置，对岩芯样品进行压汞测试试验，获取毛管压力数据，通过公式计算不同毛管压力下，汞充满岩芯样品不同预设百分比时分别对应的最小孔隙半径；将计算的最小孔隙半径，以及对应的岩芯样品孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型；根据大孔道行业标准，确定需要堵水的孔隙半径标准值；根据岩芯样品孔隙半径大小频率分布、多个孔隙半径和孔渗关系模型的相关性和二八定律，筛选出目标关系模型，目标关系模型用于计算层位前20％孔隙的孔隙半径值；根据目标井的测井报告中的孔渗数据，通过目标关系模型计算前20％孔隙的平均半径值；将目标井不同层位的孔隙半径平均值和需要堵水的孔隙半径标准值进行对比，确定目标堵水层位；根据该目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算所需堵剂量。本发明通过水流优势通道的刻画，可以定位参与流动的大孔道的位置，得到实施堵水的准确深度；提出只需要封堵连通的大孔隙的方法，大大降低了堵剂用量，节约了堵水成本。而且，过多的注入堵剂很容易造成泄漏，造成环境污染，本发明提供的方法准确计算堵剂用量，避免了污染的风险，为实际生产提供了简便、有效的技术手段。

需要说明的是：上述实施例提供的基于大孔道数据的油井堵水装置在水驱开发制定堵水方案过程中，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的基于大孔道数据的油井堵水装置与基于大孔道数据的油井堵水方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。

图8示出了本发明一个示例性实施例提供的计算机设备800的结构框图。该计算机设备800可以是：智能手机、平板电脑、MP3播放器(Moving Picture Experts Group AudioLayer III，动态影像专家压缩标准音频层面3)、MP4(Moving Picture Experts GroupAudio Layer IV，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。计算机设备800还可能被称为用户设备、便携式计算机设备、膝上型计算机设备、台式计算机设备等其他名称。

通常，计算机设备800包括有：处理器801和存储器802。

处理器801可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器801可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器801也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器801可以在集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器801还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。

存储器802可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器802还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器802中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器801所执行以实现本申请中方法实施例提供的油井抽油机的工作参数配置方法。

在一些实施例中，计算机设备800还可选包括有：外围设备接口803和至少一个外围设备。处理器801、存储器802和外围设备接口803之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口803相连。具体地，外围设备包括：射频电路804、触摸显示屏805、摄像头806、音频电路807、定位组件808和电源809中的至少一种。

外围设备接口803可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器801和存储器802。在一些实施例中，处理器801、存储器802和外围设备接口803被集成在同一芯片或电路板上；在一些其他实施例中，处理器801、存储器802和外围设备接口803中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。

射频电路804用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。射频电路804通过电磁信号与通信网络以及其他通信设备进行通信。射频电路804将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路804包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路804可以通过至少一种无线通信协议来与其它计算机设备进行通信。该无线通信协议包括但不限于：城域网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例中，射频电路804还可以包括NFC(Near Field Communication，近距离无线通信)有关的电路，本申请对此不加以限定。

显示屏805用于显示UI(User Interface，用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏805是触摸显示屏时，显示屏805还具有采集在显示屏805的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器801进行处理。此时，显示屏805还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏805可以为一个，设置计算机设备800的前面板；在另一些实施例中，显示屏805可以为至少两个，分别设置在计算机设备800的不同表面或呈折叠设计；在再一些实施例中，显示屏805可以是柔性显示屏，设置在计算机设备800的弯曲表面上或折叠面上。甚至，显示屏805还可以设置成非矩形的不规则图形，也即异形屏。显示屏805可以采用LCD(Liquid Crystal Display，液晶显示屏)、OLED(Organic Light-EmittingDiode,有机发光二极管)等材质制备。

摄像头组件806用于采集图像或视频。可选地，摄像头组件806包括前置摄像头和后置摄像头。通常，前置摄像头设置在计算机设备的前面板，后置摄像头设置在计算机设备的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件806还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。

音频电路807可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器801进行处理，或者输入至射频电路804以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在计算机设备800的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器801或射频电路804的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路807还可以包括耳机插孔。

定位组件808用于定位计算机设备800的当前地理位置，以实现导航或LBS(Location Based Service，基于位置的服务)。定位组件808可以是基于美国的GPS(GlobalPositioning System，全球定位系统)、中国的北斗系统、俄罗斯的格雷纳斯系统或欧盟的伽利略系统的定位组件。

电源809用于为计算机设备800中的各个组件进行供电。电源809可以是交流电、直流电、一次性电池或可充电电池。当电源809包括可充电电池时，该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。

在一些实施例中，计算机设备800还包括有一个或多个传感器810。该一个或多个传感器810包括但不限于：加速度传感器811、陀螺仪传感器812、压力传感器813、指纹传感器814、光学传感器815以及接近传感器816。

加速度传感器811可以检测以计算机设备800建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器811可以用于检测重力加速度在三个坐标轴上的分量。处理器801可以根据加速度传感器811采集的重力加速度信号，控制触摸显示屏805以横向视图或纵向视图进行用户界面的显示。加速度传感器811还可以用于游戏或者用户的运动数据的采集。

陀螺仪传感器812可以检测计算机设备800的机体方向及转动角度，陀螺仪传感器812可以与加速度传感器811协同采集用户对计算机设备800的3D动作。处理器801根据陀螺仪传感器812采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。

压力传感器813可以设置在计算机设备800的侧边框和/或触摸显示屏805的下层。当压力传感器813设置在计算机设备800的侧边框时，可以检测用户对计算机设备800的握持信号，由处理器801根据压力传感器813采集的握持信号进行左右手识别或快捷操作。当压力传感器813设置在触摸显示屏805的下层时，由处理器801根据用户对触摸显示屏805的压力操作，实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。

指纹传感器814用于采集用户的指纹，由处理器801根据指纹传感器814采集到的指纹识别用户的身份，或者，由指纹传感器814根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器801授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器814可以被设置计算机设备800的正面、背面或侧面。当计算机设备800上设置有物理按键或厂商Logo时，指纹传感器814可以与物理按键或厂商Logo集成在一起。

光学传感器815用于采集环境光强度。在一个实施例中，处理器801可以根据光学传感器815采集的环境光强度，控制触摸显示屏805的显示亮度。具体地，当环境光强度较高时，调高触摸显示屏805的显示亮度；当环境光强度较低时，调低触摸显示屏805的显示亮度。在另一个实施例中，处理器801还可以根据光学传感器815采集的环境光强度，动态调整摄像头组件806的拍摄参数。

接近传感器816，也称距离传感器，通常设置在计算机设备800的前面板。接近传感器816用于采集用户与计算机设备800的正面之间的距离。在一个实施例中，当接近传感器816检测到用户与计算机设备800的正面之间的距离逐渐变小时，由处理器801控制触摸显示屏805从亮屏状态切换为息屏状态；当接近传感器816检测到用户与计算机设备800的正面之间的距离逐渐变大时，由处理器801控制触摸显示屏805从息屏状态切换为亮屏状态。

本领域技术人员可以理解，图8中示出的结构并不构成对计算机设备800的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。

在示例性实施例中，还提供了一种计算机可读存储介质，例如包括指令的存储器，上述指令可由终端中的处理器执行以完成下述实施例中的基于大孔道数据的油井堵水方法。例如，所述计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等。

本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，本发明具体实施方式只是结合附表、附图对本发明进行了示例性描述，显然本发明具体实现并不受上述方式的限制，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于连通的大孔道数据的油井堵水方法，其特征在于，所述方法包括：

对取芯井的岩芯样品进行压汞测试试验，获取所述岩芯样品的毛管压力数据，通过公式计算不同毛管压力下，汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径；

根据计算的最小孔隙半径，以及与其对应的所述岩芯样品的孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型，所述孔渗为孔隙度和渗透率的合称，所述孔渗测试数据包括所述岩芯样品不同层位的孔隙度和渗透率数据；

根据大孔道行业标准，确定油藏开发过程中需要堵水的孔隙半径标准值；

根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，从所述多个孔隙半径和孔渗关系模型中筛选出目标关系模型，所述目标关系模型用于计算层位前20％孔隙的孔隙半径值；

根据目标井的测井报告中不同层位的孔渗数据，采用所述目标关系模型，计算所述目标井不同层位的孔隙半径平均值；

基于所述目标井不同层位的孔隙半径平均值与所述需要堵水的孔隙半径标准值的对比，确定所述目标井的目标堵水层位，所述目标堵水层位的孔隙半径平均值不小于所述需要堵水的孔隙半径标准值；

根据所述目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算需要的堵剂量。

2.根据权利要求1所述的方法，其特征在于，对取芯井的岩芯样品进行压汞测试试验，获取所述岩芯样品的毛管压力数据，通过公式计算不同毛管压力下，汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径，包括：

在压汞测试中，针对所述岩芯样品的不同层位，获取所述层位下汞充满所述不同预设百分比时的毛管压力大小；

针对每个层位对应的不同毛管压力值，根据孔隙半径公式，计算所述毛管压力值对应的孔隙半径值，则所述孔隙半径值即为汞充满所述岩芯样品相应预设百分比时的最小孔隙半径。

3.根据权利要求1所述的方法，其特征在于，根据计算的最小孔隙半径，以及与其对应的所述岩芯样品的孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型，包括:

根据计算的汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径，结合所述岩芯样品不同层位的孔隙度和渗透率，得到多组孔隙半径和孔渗数据；

针对每一组孔隙半径和孔渗数据进行多元线性回归处理，建立相应预设百分比下的孔隙半径和孔渗关系模型，得到多个孔隙半径和孔渗关系模型。

4.根据权利要求1所述的方法，其特征在于，根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，从所述多个孔隙半径和孔渗关系模型中筛选出目标关系模型，包括：

根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，确定需要堵水的判别条件，所述判别条件为层位前20％孔隙的孔隙半径平均值大于或等于所述需要堵水的孔隙半径标准值；

基于所述判别条件，将所述多个孔隙半径和孔渗关系模型中孔隙半径占比在前20％的孔隙半径和孔渗关系模型确定为所述目标关系模型。

5.根据权利要求1所述的方法，其特征在于，根据目标井的测井报告中不同层位的孔渗数据，采用所述目标关系模型，计算所述目标井不同层位的孔隙半径平均值，包括：

根据所述目标井的测井报告，获取所述目标井不同层位的孔隙度和渗透率值；

根据所述目标井每一个层位的孔隙度和渗透率，采用目标关系模型，计算所述层位相对于所述目标关系模型的孔隙半径值；

当目标关系模型数目为1个时，将所述目标关系模型的孔隙半径值确定为所述层位的孔隙半径平均值；当目标关系模型数目为多个时，计算多个目标关系模型对应的孔隙半径值的均值，将所述均值作为所述层位的孔隙半径平均值。

6.一种基于连通的大孔道数据的油井堵水装置，其特征在于，所述装置包括：

获取模块，用于对取芯井的岩芯样品进行压汞测试试验，获取所述岩芯样品的毛管压力数据，通过公式计算不同毛管压力下，汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径；

模型建立模块，用于根据计算的最小孔隙半径，以及与其对应的所述岩芯样品的孔渗测试数据，对孔隙半径、孔隙度和渗透率进行多元线性回归处理，建立多个孔隙半径和孔渗关系模型，所述孔渗为孔隙度和渗透率的合称，所述孔渗测试数据包括所述岩芯样品不同层位的孔隙度和渗透率数据；

孔隙半径标准值确定模块，用于根据大孔道行业标准，确定油藏开发过程中需要堵水的孔隙半径标准值；

筛选模块，用于根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，从所述多个孔隙半径和孔渗关系模型中筛选出目标关系模型，所述目标关系模型用于计算层位前20％孔隙的孔隙半径值；

孔隙半径计算模块，用于根据目标井的测井报告中不同层位的孔渗数据，采用所述目标关系模型，计算所述目标井不同层位的孔隙半径平均值；

目标堵水层位确定模块，用于基于所述目标井不同层位的孔隙半径平均值与所述需要堵水的孔隙半径标准值的对比，确定所述目标井的目标堵水层位，所述目标堵水层位的孔隙半径平均值不小于所述需要堵水的孔隙半径标准值；

堵剂量计算模块，用于根据所述目标堵水层位的孔隙度、有效厚度和堵剂推进有效范围，计算需要的堵剂量。

7.根据权利要求6所述的装置，其特征在于，获取模块，用于：

在压汞测试中，针对所述岩芯样品的不同层位，获取所述层位下汞充满所述不同预设百分比时的毛管压力大小；

针对每个层位对应的不同毛管压力值，根据孔隙半径公式，计算所述毛管压力值对应的孔隙半径值，则所述孔隙半径值即为汞充满所述岩芯样品相应预设百分比时的最小孔隙半径。

8.根据权利要求6所述的装置，其特征在于，模型建立模块，用于:

根据计算的汞充满所述岩芯样品不同预设百分比时分别对应的最小孔隙半径，结合所述岩芯样品不同层位的孔隙度和渗透率，得到多组孔隙半径和孔渗数据；

针对每一组孔隙半径和孔渗数据进行多元线性回归处理，建立相应预设百分比下的孔隙半径和孔渗关系模型，得到多个孔隙半径和孔渗关系模型。

9.根据权利要求6所述的装置，其特征在于，筛选模块，用于：

根据所述岩芯样品孔隙半径大小频率分布、所述多个孔隙半径和孔渗关系模型的相关性和二八定律，确定需要堵水的判别条件，所述判别条件为层位前20％孔隙的孔隙半径平均值大于或等于所述需要堵水的孔隙半径标准值；

基于所述判别条件，将所述多个孔隙半径和孔渗关系模型中孔隙半径占比在前20％的孔隙半径和孔渗关系模型确定为所述目标关系模型。

10.根据权利要求6所述的装置，其特征在于，孔隙半径计算模块，用于：

根据所述目标井的测井报告，获取所述目标井不同层位的孔隙度和渗透率值；

根据所述目标井每一个层位的孔隙度和渗透率，采用目标关系模型，计算所述层位相对于所述目标关系模型的孔隙半径值；

当目标关系模型数目为1个时，将所述目标关系模型的孔隙半径值确定为所述层位的孔隙半径平均值；当目标关系模型数目为多个时，计算多个目标关系模型对应的孔隙半径值的均值，将所述均值作为所述层位的孔隙半径平均值。

Images