岩石储层中值半径预测模型建立方法、装置及系统与流程

本发明涉及地质勘测领域，具体而言，涉及一种岩石储层中值半径预测模型建立方法、装置及系统。

背景技术：

测井，也叫地球物理测井或矿场地球物理，简称测井，是利用岩层的电化学特性、导电特性、声学特性、放射性等地球物理特性，测量地球物理参数的方法，属于应用地球物理方法(包括重、磁、电、震、核)之一。根据地质和地球物理条件，合理地选用综合测井方法，可以详细研究钻孔地质剖面、探测有用矿产、详细提供计算储量所必需的数据，如油层的有效厚度、孔隙度、含油气饱和度和渗透率等，以及研究钻孔技术情况等任务。此外，井中磁测、井中激发激化、井中无线电波透视和重力测井等方法还可以发现和研究钻孔附近的盲矿体。测井方法在石油、煤、金属与非金属矿产及水文地质、工程地质的钻孔中，都得到广泛的应用。特别在油气田、煤田及水文地质勘探工作中，可以测量岩石的中值半径已成为不可缺少的勘探方法之一。

岩石中值半径是反映储层岩石渗透性能的重要指标，但是要获得此指标首先要用专用设备获取岩心，再用取样设备在岩心上取样，然后将取得的样品放入专门的分析化验设备—压汞仪进行测试，费工费时，成本很高，且对人体及环境均有污染。

技术实现要素：

有鉴于此，本发明实施例的目的在于提供一种岩石储层中值半径预测模型建立方法、装置及系统。

第一方面，本发明实施例提供了一种岩石储层中值半径预测模型建立方法，所述岩石储层中值半径预测模型建立方法包括：

接收测井仪器发送的不同深度的岩石储层样本对应的多个不同类别的测井参数；

接收孔隙分析仪发送的不同深度的岩石储层样本的中值半径；

将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前n个测井参数样本；

对找出的各个深度的对应的不同测井参数样本以及中值半径进行二元回归分析，建立岩石储层中值半径预测模型。

第二方面，本发明实施例还提供了一种岩石储层中值半径预测模型建立装置，所述岩石储层中值半径预测模型建立装置包括：

信息收发单元，用于接收测井仪器发送的不同深度的岩石储层样本对应的多个不同类别的测井参数；

所述信息收发单元还用于接收孔隙分析仪发送的不同深度的岩石储层样本的中值半径；

线性回归分析单元，用于将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前n个测井参数样本；

预测模型建立单元，用于对找出的各个深度的对应的不同测井参数样本以及中值半径进行二元回归分析，建立岩石储层中值半径预测模型。

第三方面，本发明实施例还提供了一种岩石储层中值半径预测模型建立系统，所述岩石储层中值半径预测模型建立系统包括岩石储层中值半径预测模型建立装置、测井仪器以及孔隙分析仪，所述岩石储层中值半径预测模型建立装置分别与所述测井仪器、所述孔隙分析仪建立通信连接，所述岩石储层中值半径预测模型建立装置包括信息收发单元、线性回归分析单元、预测模型建立单元，

所述测井仪器用于发送不同深度的岩石储层样本对应的多个不同类别的测井参数至所述信息收发单元；

所述孔隙分析仪用于发送不同深度的岩石储层样本的中值半径至所述信息收发单元；

所述信息收发单元用于接收测井仪器发送的不同深度的岩石储层样本对应的多个不同类别的测井参数；

所述信息收发单元还用于接收孔隙分析仪发送的不同深度的岩石储层样本的中值半径；

所述线性回归分析单元用于将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前n个测井参数样本；

所述预测模型建立单元用于对找出的各个深度的对应的不同测井参数样本以及中值半径进行线性回归分析，建立岩石储层中值半径预测模型。

与现有技术相比，本发明提供的岩石储层中值半径预测模型建立方法、装置及系统，通过将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前n个测井参数样本，对找出的各个深度的对应的不同测井参数样本以及中值半径进行线性回归分析，建立岩石储层中值半径预测模型。当再次对同一区域的岩石储层采集中值半径时，无需再利用专用设备获取岩心，及取样设备在岩心上取样，再用将取得的样品放入压汞仪进行测试，省时省力，成本很低，并且对人体及环境无污染，仅需将测井参数输入预建立好的岩石储层中值半径预测模型即可将岩石中值半径计算得出，省时省力。

为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。

附图说明

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

图1为本发明实施例提供的岩石储层中值半径预测模型建立系统的交互示意图；

图2为本发明实施例提供的服务器的结构框图；

图3为本发明实施例提供的岩石储层中值半径预测模型建立方法的流程图；

图4为本发明实施例提供的岩石储层中值半径预测模型建立装置的功能单元示意图。

图标：100-测井仪器；200-服务器；300-孔隙分析仪；400-岩石储层中值半径预测模型建立装置；101-处理器；102-存储器；103-存储控制器；104-外设接口；401-信息收发单元；402-线性回归分析单元；403-预测模型建立单元。

具体实施方式

下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明较佳实施例所提供的岩石储层中值半径预测模型建立方法、装置及系统可应用于如图1所示的应用环境中。如图1所示，测井仪器100、服务器200、孔隙分析仪300位于网络中，通过该网络，测井仪器100与服务器200进行数据交互。该服务器200可以是，但不限于，网络服务器、数据库服务器、云端服务器等等。测井仪，是用于凿井领域，可对井壁进行连续扫描，也可对任意水平进行横向扫描，给出井筒竖直剖面、水平断面、井筒有效断面、井筒偏斜距离，并且可以测出各类测井参数，例如，声波参数、伽马参数、den密度参数、rt电阻率参数等等。

图2示出了一种可应用于本发明实施例中的服务器200的结构框图。所述服务器包括岩石储层中值半径预测模型建立装置400、存储器102、存储控制器103及处理器101。

所述存储器102、存储控制器103及处理器101，各元件相互之间直接或间接地电性连接，以实现数据的传输或交互。例如，这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。所述岩石储层中值半径预测模型建立装置400包括至少一个可以软件或固件(firmware)的形式存储于所述存储器102中。所述处理器101用于执行存储器102中存储的可执行模块，例如，所述岩石储层中值半径预测模型建立装置400包括的软件功能模块或计算机程序。

其中，存储器102可以是，但不限于，随机存取存储器(randomaccessmemory，ram)，只读存储器102readonlymemory，rom)，可编程只读存储器(programmableread-onlymemory，prom)，可擦除只读存储器(erasableprogrammableread-onlymemory，eprom)，电可擦除只读存储器(electricerasableprogrammableread-onlymemory，eeprom)等。其中，存储器102用于存储程序，所述处理器101在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的流过程定义的服务器所执行的方法可以应用于处理器101中，或者由处理器101实现。

处理器101可能是一种集成电路芯片，具有信号的处理能力。上述的处理器101可以是通用处理器，包括中央处理器(centralprocessingunit，简称cpu)、网络处理器(networkprocessor，简称np)等；还可以是数字信号处理器(dsp)、专用集成电路(asic)、现成可编程门阵列(fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器101也可以是任何常规的处理器101等。

外设接口104将各种输入/输出装置耦合至处理器101以及存储器102。在一些实施例中，外设接口104、处理器101以及存储控制器103可以在单个芯片中实现。在其他一些实例中，他们可以分别由独立的芯片实现。

可以理解，图2所示的结构仅为示意，服务器200还可包括比图2中所示更多或者更少的组件，或者具有与图2所示不同的配置。图2中所示的各组件可以采用硬件、软件或其组合实现。

请参阅图3，本发明实施例提供了一种岩石储层中值半径预测模型建立方法，应用于服务器200，所述岩石储层中值半径预测模型建立方法包括：

步骤s301：接收测井仪器100发送的不同深度的岩石储层样本对应的多个不同类别的测井参数。

本实施例中，测井参数电阻率参数、伽马参数、密度参数、中子参数、井径参数以及声波参数等等。

例如，位于地下1800m的岩石储层的测井参数分别为电阻参数10.231，伽马参数90.179，密度参数2.489，中子参数15.38，井径参数10.367，声波参数71.02。

具体实施例中，可将测井仪器100安装于钻井机的钻头内，当钻井机进行钻井时，测井仪器100即可对钻出的岩心的各种测井参数进行采集，测井仪器100安装有无线通信模块，通过无线通信模块即可将采集到的各种测井参数发送至服务器200。

步骤s302：接收孔隙分析仪300发送的不同深度的岩石储层样本的中值半径。

本实施例中，孔隙分析仪300为压汞仪，压汞仪内置强大的数据处理和报告程序包，快速升压、灵活、可控的真空系统，和高性能的低/高压系统，可使用汞侵入法来测定总孔体积、孔径分布、中值半径、孔隙率、密度和传输性。

步骤s303：将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前n个测井参数样本。

较佳地，由于对测井参数选择的种类越多，建立的岩石储层中值半径预测模型越具备价值，而测井参数与岩石中值半径相关性越好，建立的岩石储层中值半径预测模型越具备价值，因而将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前两个测井参数样本，最后得出的岩石储层中值半径预测模型越具备预测价值。

具体地，依据算式计算出各个深度对应的不同类别的测井参数与对应的中值半径之间的相关性系数，并找出最接近于1的两个相关性系数分别对应的两类测井参数样本，其中，rij为相关性系数，x为测井参数，y为中值半径，

其中，根据发明人试验采集到的测井参数，各个深度对应的不同类别的测井参数以及中值半径，如表1所示：

表1

具体地，回归分析是研究变量之间的相关关系的一种多元统计分析方法，当假定某个因变量与p个自变量之间有线性相关关系时，可以将这种关系用下式表示：

上式称为线性回归方程，式中的b0,b1,...,bp称为回归系数，它们是待定的，如果已经获得了关于这些变量的n组数据(x1i,x2i,...,xpi,yii＝1,2,...,n)，将第i组数据代入上式就是

式中称为回归值，也叫计算值，它与观测值yi之间是有误差的，回归系数的确定就是要使误差平方和达到最小，也就是利用最小二乘原理确定回归系数，归结为如下方程组：

任意两个变量xi与xj之间的相关系数由下式确定：

其中，其中，相关系数的分布区间在(-1，+1)之间，当相关系数小于0时称为负相关，相关系数大于0时称为正相关。将表1中的每种类别的测井参数的不同深度的测井参数代入上式中，当相关系数越趋近于1时，说明该测井参数与中值半径的相关性越好，反之，则说明相关性越差。本实施例中，相关性较好的两类测井参数为声波参数和伽马参数。

步骤s304：对找出的各个深度的对应的不同测井参数样本以及中值半径进行二元回归分析，建立岩石储层中值半径预测模型。

具体地，所述岩石储层中值半径预测模型为其中，x1为声波参数、x2为伽马参数，a、b、m为常数，y为中值半径。

本实施例中，将上述的不同深度的声波参数、伽马参数代入上式中进行二元回归分析得到岩石储层中值半径预测模型当然地，声波参数及伽马参数的数据组越多，岩石储层中值半径预测模型越具备参考价值。当需要再次测量中值半径时，仅需将测井参数的声波参数、伽马参数代入即可。

请参阅图4，本发明实施例还提供了一种岩石储层中值半径预测模型建立装置400，所述岩石储层中值半径预测模型建立装置400包括信息收发单元401、线性回归分析单元402以及预测模型建立单元403。

所述信息收发单元401用于接收测井仪器100发送的不同深度的岩石储层样本对应的多个不同类别的测井参数。

可以理解地，步骤s301可以由信息收发单元401执行。

所述信息收发单元401还用于接收孔隙分析仪300发送的不同深度的岩石储层样本的中值半径。

可以理解地，步骤s302可以由信息收发单元401执行。

所述线性回归分析单元402用于将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前n个测井参数样本。

可以理解地，步骤s303可以由线性回归分析单元402执行。

具体地，所述线性回归分析单元402用于将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前两个测井参数样本。

具体地，所述线性回归分析单元402用于依据算式计算出各个深度对应的不同类别的测井参数与对应的中值半径之间的相关性系数，并找出最接近于1的两个相关性系数分别对应的两类测井参数样本，其中，rij为相关性系数，x为测井参数，y为中值半径，

预测模型建立单元403用于对找出的各个深度的对应的不同测井参数样本以及中值半径进行二元回归分析，建立岩石储层中值半径预测模型。

可以理解地，步骤s304可以由预测模型建立单元403执行。

本实施例中，所述两类测井参数样本分别为声波参数和伽马参数，所述岩石储层中值半径预测模型为其中，x1为声波参数、x2为伽马参数，a、b、m为常数，y为中值半径。

请参阅图1，本发明实施例提供了一种岩石储层中值半径预测模型建立系统，所述岩石储层中值半径预测模型建立系统包括岩石储层中值半径预测模型建立装置400、测井仪器100以及孔隙分析仪300。所述岩石储层中值半径预测模型建立装置400分别与所述测井仪器100、所述孔隙分析仪300建立通信连接，所述岩石储层中值半径预测模型建立装置400包括信息收发单元401、线性回归分析单元402、预测模型建立单元403。

所述测井仪器100用于发送不同深度的岩石储层样本对应的多个不同类别的测井参数至所述信息收发单元401，所述孔隙分析仪300用于发送不同深度的岩石储层样本的中值半径至所述信息收发单元401；所述信息收发单元401用于接收测井仪器100发送的不同深度的岩石储层样本对应的多个不同类别的测井参数；所述信息收发单元401还用于接收孔隙分析仪300发送的不同深度的岩石储层样本的中值半径；所述线性回归分析单元402用于将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前n个测井参数样本；所述预测模型建立单元403用于对找出的各个深度的对应的不同测井参数样本以及中值半径进行二元回归分析，建立岩石储层中值半径预测模型。

所述线性回归分析单元402用于依据算式计算出各个深度对应的不同类别的测井参数与对应的中值半径之间的相关性系数，并找出最接近于1的两个相关性系数分别对应的两类测井参数样本，其中，rij为相关性系数，x为测井参数，y为中值半径，

综上所述，本发明提供的岩石储层中值半径预测模型建立方法、装置及系统，通过将位于各个深度对应的不同类别的测井参数以及中值半径进行多元线性回归分析，找出与中值半径相关性最好的前n个测井参数样本，对找出的各个深度的对应的不同测井参数样本以及中值半径进行二元回归分析，建立岩石储层中值半径预测模型。当再次对同一区域的岩石储层采集中值半径时，无需再利用专用设备获取岩心，及取样设备在岩心上取样，再用将取得的样品放入压汞仪进行测试，省时省力，成本很低，并且对人体及环境无污染，仅需将测井参数输入预建立好的岩石储层中值半径预测模型即可将岩石中值半径计算得出，省时省力。

在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。

所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-onlymemory)、随机存取存储器(ram，randomaccessmemory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。