CN116816336A - 地层渗透率测井计算方法、装置、电子设备及介质 - Google Patents

Abstract

本申请公开了一种地层渗透率测井计算方法、装置、电子设备及介质，包括：针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线；通过校正后曲线，计算岩心地面孔隙度；建立岩心地面孔隙度和覆压孔隙度的第一转换公式，计算岩心地面孔隙度对应的覆压孔隙度；建立覆压孔隙度和覆压渗透率的第二转换公式；将岩心地面孔隙度对应的覆压孔隙度带入第二转换公式，获得岩心地面孔隙度对应的覆压渗透率，即为地层渗透率。本发明基于有限的岩心实验，利用环境校正后的测井曲线开展大斜度井条件下弱成岩地层渗透率的计算，为储层特征和产能评价提供精确的地层渗透率参数，具有资料覆盖广、成本低和精度高等优势，从而指导油气田的勘探开发工作。

Description

地层渗透率测井计算方法、装置、电子设备及介质

技术领域

本发明涉及石油勘探开发领域，更具体地，涉及一种地层渗透率测井计算方法、装置、电子设备及介质。

背景技术

在油气田勘探开发中，渗透率是评价储层特征、确定储层产能的关键参数之一。通过岩心实验、试井分析、测井解释等技术手段均可求取地层的渗透率值。实际应用中，由于岩心和试井资料较少，且获取成本高、时效慢。利用测井曲线计算地层渗透率是油田现场更为常用和便捷的手段，具有资料覆盖广、成本较低和快速准确等优势。

目前，测井技术计算地层渗透率的一般计算方法，是以地面岩心数据为刻度，利用测井曲线计算地层孔隙度，通过孔隙度-渗透率数学转换公式，求取地层渗透率。但在大斜度井、弱成岩储层条件下，利用测井资料求取地层渗透率面临两个问题：一是大斜度井的测井曲线需要环境校正。由于受重力或井轨迹的影响，中子、密度测井仪器在斜井段推靠井壁不紧密，造成测井曲线录取质量不高，和直井的测井响应特征存在较大差异，进而影响了大斜度井测井计算地层渗透率的精度。二是弱成岩储层的地面岩心渗透率不能代表地层实际渗透率。由于弱成岩储层的岩石胶结差、未充分压实固结成岩，岩心取到地面后由于泄压作用，地面条件下测得的渗透率、孔隙度远大于地层实际渗透率、孔隙度。因此，受压力变化的影响，弱成岩储层的地面岩心渗透率不能等同于地层实际渗透率。如果以岩心地面渗透率刻度测井计算地层渗透率，将会产生较大误差。

标准层法是目前进行中子、密度等常规测井曲线校正的常用方法。该方法多用于直井常规测井曲线校正，前提是认为在工区范围内标准层的中子、密度测井响应值是相同的，不会随着不同井标准层深度的差异而变化，而忽略了地层压实作用对标准层地层性质的影响，即不同深度的标准层应该具有不同的测井响应值，特别是在岩石成岩较弱的地层，如果不考虑地层深度的变化对测井响应值的影响，将极大影响测井曲线校正的精度，给包括渗透率在内的储层参数计算带来很大困难。

试井分析技术可以得到包括地层渗透率、表皮系数、地层系数、地层压力在内的地层参数，可以有效表征地层的动态特征。但试井成本高、测试时间长，油田不会大量录取，且试井所得的地层参数为多层合试层段的参数平均值，精度有限，不能用于小层的储层及产能的精细评价。

岩心覆压实验是能够得到较为精确地层渗透率的主要方法之一。通过对取到地面的岩心不断加压，测量出不同压力条件下渗透率和孔隙度值。当有效压力等于油气藏地层压力时，可认为该条件下测得的覆压渗透率近似等于地层渗透率。但受覆压渗透率实验成本高、花费时间长等因素的影响，覆压实验数据非常少，油田现场一般做地面孔渗实验，用地面岩心孔渗作为地层条件下的岩心孔渗来拟合计算测井渗透率，从而造成较大误差。特别是对于弱成岩地层，该方法计算的地层渗透率误差更大。

因此，有必要开发一种地层渗透率测井计算方法、装置、电子设备及介质。

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

发明内容

本发明提出了一种地层渗透率测井计算方法、装置、电子设备及介质，其能够基于有限的岩心实验，利用环境校正后的测井曲线开展大斜度井条件下弱成岩地层渗透率的计算，为储层特征和产能评价提供较为精确的地层渗透率参数，具有资料覆盖广、成本低和精度高等优势，从而指导油气田的勘探开发工作。

第一方面，本公开实施例提供了一种地层渗透率测井计算方法，包括：

针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线；

通过所述校正后曲线，计算岩心地面孔隙度；

建立所述岩心地面孔隙度和覆压孔隙度的第一转换公式，计算所述岩心地面孔隙度对应的覆压孔隙度；

建立所述覆压孔隙度和覆压渗透率的第二转换公式；

将所述岩心地面孔隙度对应的覆压孔隙度带入所述第二转换公式，获得所述岩心地面孔隙度对应的覆压渗透率，即为所述地层渗透率。

优选地，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线包括：

确定区域标准层，获取所述区域标准层内直井标准层的中子值、密度值和深度值；

建立所述直井标准层的中子值随深度变化的关系式以及密度值随深度变化的关系式；

将测量深度校正到垂直深度，获取大斜度井标准层的中子理论值与密度理论值；

获取大斜度井标准层的中子实测值与密度实测值；

计算大斜度井标准层的中子校正量与密度校正量；

根据所述中子校正量与所述密度校正量对大斜度井斜井段的中子曲线和密度曲线进行校正，获得所述校正后曲线。

优选地，通过公式(1)计算大斜度井标准层的中子校正量：

ΔCNL＝CNL1–CNL2 (1)

其中，ΔCNL为中子校正量，CNL1为中子理论值，CNL2为中子实测值。

优选地，通过公式(2)计算大斜度井标准层的密度校正量：

ΔDEN＝DEN1–DEN2 (2)

其中，ΔDEN为密度校正量，DEN1为密度理论值，DEN2为密度实测值。

优选地，通过所述校正后曲线，采用中子-密度交会法进行储层孔隙度计算：

其中，φ_t为岩心地面孔隙度，CNL为环境校正后的中子值，DEN为环境校正后的密度值。

优选地，所述第一转换公式为：

其中，φ_t为岩心地面孔隙度，φ为覆压孔隙度，E和F为拟合系数。

优选地，所述第二转换公式为：

K＝G×e^H×φ (5)

其中，φ为覆压孔隙度，K为覆压渗透率，G和H为拟合系数。

作为本公开实施例的一种具体实现方式，

第二方面，本公开实施例还提供了一种地层渗透率测井计算装置，包括：

环境校正模块，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线；

岩心地面孔隙度计算模块，通过所述校正后曲线，计算岩心地面孔隙度；

第一转换模块，建立所述岩心地面孔隙度和覆压孔隙度的第一转换公式，计算所述岩心地面孔隙度对应的覆压孔隙度；

第二转换模块，建立所述覆压孔隙度和覆压渗透率的第二转换公式；

预测模块，将所述岩心地面孔隙度对应的覆压孔隙度带入所述第二转换公式，获得所述岩心地面孔隙度对应的覆压渗透率，即为所述地层渗透率。

优选地，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线包括：

确定区域标准层，获取所述区域标准层内直井标准层的中子值、密度值和深度值；

建立所述直井标准层的中子值随深度变化的关系式以及密度值随深度变化的关系式；

将测量深度校正到垂直深度，获取大斜度井标准层的中子理论值与密度理论值；

获取大斜度井标准层的中子实测值与密度实测值；

计算大斜度井标准层的中子校正量与密度校正量；

根据所述中子校正量与所述密度校正量对大斜度井斜井段的中子曲线和密度曲线进行校正，获得所述校正后曲线。

优选地，通过公式(1)计算大斜度井标准层的中子校正量：

ΔCNL＝CNL1–CNL2 (1)

其中，ΔCNL为中子校正量，CNL1为中子理论值，CNL2为中子实测值。

优选地，通过公式(2)计算大斜度井标准层的密度校正量：

ΔDEN＝DEN1–DEN2 (2)

其中，ΔDEN为密度校正量，DEN1为密度理论值，DEN2为密度实测值。

优选地，通过所述校正后曲线，采用中子-密度交会法进行储层孔隙度计算：

其中，φ_t为岩心地面孔隙度，CNL为环境校正后的中子值，DEN为环境校正后的密度值。

优选地，所述第一转换公式为：

其中，φ_t为岩心地面孔隙度，φ为覆压孔隙度，E和F为拟合系数。

优选地，所述第二转换公式为：

K＝G×e^H×φ (5)

其中，φ为覆压孔隙度，K为覆压渗透率，G和H为拟合系数。

第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：

存储器，存储有可执行指令；

处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的地层渗透率测井计算方法。

第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的地层渗透率测井计算方法。

其有益效果在于：本发明能够基于有限的岩心实验，利用环境校正后的测井曲线开展大斜度井条件下弱成岩地层渗透率的计算。一是考虑了地层压实作用对不同深度标准层中子、密度测井响应的影响，避免了以往标准层法忽略了压实作用对标准层性质变化影响的做法，能够准确、快速实现中子和密度曲线的环境校正，提供了正确的测井曲线；二是考虑了利用覆压孔隙度、覆压渗透率对测井基于地面孔隙度计算的储层孔隙度、渗透率值进行校正，计算得到的地层渗透率与覆压渗透率的误差小，为储层特征和产能评价提供较为精确的地层渗透率参数，具有资料覆盖广、成本低和精度高等优势，从而指导油气田的勘探开发工作。

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

附图说明

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

图1示出了根据本发明的一个实施例的地层渗透率测井计算方法的步骤的流程图。

图2示出了根据本发明的一个实施例的标准层的中子值与深度变化关系示意图。

图3示出了根据本发明的一个实施例的标准层的密度值与深度变化关系示意图。

图4示出了根据本发明的一个实施例的环境校正后的中子、密度曲线的示意图。

图5示出了根据本发明的一个实施例的岩心地面孔隙度和覆压孔隙度的关系示意图。

图6示出了根据本发明的一个实施例的覆压孔隙度和覆压渗透率的关系示意图。

图7示出了根据本发明的一个实施例的测井计算的地层孔隙度、渗透率与覆压孔隙度、覆压渗透率的对比示意图。

图8示出了根据本发明的一个实施例的一种地层渗透率测井计算装置的框图。

附图标记说明：

201、环境校正模块；202、岩心地面孔隙度计算模块；203、第一转换模块；204、第二转换模块；205、预测模块。

具体实施方式

下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。

本发明提供一种地层渗透率测井计算方法，包括：

针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线；

通过校正后曲线，计算岩心地面孔隙度；

建立岩心地面孔隙度和覆压孔隙度的第一转换公式，计算岩心地面孔隙度对应的覆压孔隙度；

建立覆压孔隙度和覆压渗透率的第二转换公式；

将岩心地面孔隙度对应的覆压孔隙度带入第二转换公式，获得岩心地面孔隙度对应的覆压渗透率，即为地层渗透率。

在一个示例中，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线包括：

确定区域标准层，获取区域标准层内直井标准层的中子值、密度值和深度值；

建立直井标准层的中子值随深度变化的关系式以及密度值随深度变化的关系式；

将测量深度校正到垂直深度，获取大斜度井标准层的中子理论值与密度理论值；

获取大斜度井标准层的中子实测值与密度实测值；

计算大斜度井标准层的中子校正量与密度校正量；

根据中子校正量与密度校正量对大斜度井斜井段的中子曲线和密度曲线进行校正，获得校正后曲线。

在一个示例中，通过公式(1)计算大斜度井标准层的中子校正量：

ΔCNL＝CNL1–CNL2 (1)

其中，ΔCNL为中子校正量，CNL1为中子理论值，CNL2为中子实测值。

在一个示例中，通过公式(2)计算大斜度井标准层的密度校正量：

ΔDEN＝DEN1–DEN2 (2)

其中，ΔDEN为密度校正量，DEN1为密度理论值，DEN2为密度实测值。

在一个示例中，通过校正后曲线，采用中子-密度交会法进行储层孔隙度计算：

其中，φ_t为岩心地面孔隙度，CNL为环境校正后的中子值，DEN为环境校正后的密度值。

在一个示例中，第一转换公式为：

其中，φ_t为岩心地面孔隙度，φ为覆压孔隙度，E和F为拟合系数。

在一个示例中，第二转换公式为：

K＝G×e^H×φ (5)

其中，φ为覆压孔隙度，K为覆压渗透率，G和H为拟合系数。

具体地，对大斜度井的中子和密度测井曲线进行环境校正：

选取区域标准层，并获取区域标准层内直井标准层的中子值、密度值和深度值。选取岩性纯、厚度大、分布稳定的泥岩层作为区域标准层。为保证每口直井标准层读值的准确性，每口井的读值层段不少于5个。按照深度顺序，读取每个层段内中子和密度曲线的平均值，深度值取每个层段的中部深度值。

建立直井标准层的中子值随深度变化的关系式以及密度值随深度变化的关系式。

将中子值与深度值做线性回归，建立标准层中子值随深度变化的关系式：

CNL＝A*DEPTH+B

式中：CNL为中子值，DEPTH为深度值，A、B为常数。

将密度值与深度值做线性回归，建立标准层密度值随深度变化的关系式：

DEN＝C*DEPTH+D

式中：DEN为密度值，DEPTH为深度值，C、D为常数。

将大斜度井标准层内测量深度校正到垂直深度，获取大斜度井标准层的中子理论值与密度理论值。

将大斜度井标准层内任一深度点的测量深度校正到垂直深度，深度值的采样间隔为0.125m，与补偿密度、密度曲线实际测量值的采样间隔一致。如果大斜度井的标准层在垂直井段内，则不需要进行垂深校正。

将大斜度井标准层内任一深度点的垂直深度值代入中子和深度关系式，计算得到标准层每个垂直深度点的中子理论值CNL1；将大斜度井标准层内任一深度点的垂直深度值代入密度与深度关系式，计算得到标准层每个垂直深度点的密度理论值DEN1。

分别读取大斜度井标准层中子实测值CNL2和密度实测值DEN2。

根据大斜度井标准层中子的理论值和实测值获取大斜度井斜井段中子的校正量，即为公式(1)，并根据大斜度井标准层密度的理论值和实际测量值获取大斜度井斜井段密度的校正量，即为公式(2)。

对大斜度井斜井段的中子曲线和密度曲线进行校正，将中子校正量ΔCNL加到大斜度井斜井段中子实际测量值，即完成中子曲线的环境校正；将密度校正量ΔDEN加到大斜度井斜井段密度实际测量值，即完成密度曲线的环境校正。

通过足够多的取心资料进行岩心归位，利用环境校正后中子和密度测井曲线，采用中子-密度交会法进行储层孔隙度计算，并利用岩心分析地面孔隙度数据校正孔隙度的计算结果，此时得到的孔隙度相当于未做覆压校正的岩心地面孔隙度。

选取与上覆岩压和地层压力一致或者接近的覆压岩石实验数据，构建地面岩心孔隙度和覆压孔隙度的交会图版，得到岩心地面孔隙度和覆压孔隙度的第一转换公式为公式(4)，计算岩心地面孔隙度对应的覆压孔隙度。

选择地层条件下的岩心覆压孔隙度和覆压渗透率实验数据，构建覆压孔隙度和覆压渗透率交会图版，得到覆压孔隙度和覆压渗透率的第二转换公式为公式(5)。

将岩心地面孔隙度对应的覆压孔隙度带入第二转换公式，获得岩心地面孔隙度对应的覆压渗透率，计算得到的覆压渗透率即为弱成岩储层的地层渗透率。

本发明还提供一种地层渗透率测井计算装置，包括：

环境校正模块，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线；

岩心地面孔隙度计算模块，通过校正后曲线，计算岩心地面孔隙度；

第一转换模块，建立岩心地面孔隙度和覆压孔隙度的第一转换公式，计算岩心地面孔隙度对应的覆压孔隙度；

第二转换模块，建立覆压孔隙度和覆压渗透率的第二转换公式；

预测模块，将岩心地面孔隙度对应的覆压孔隙度带入第二转换公式，获得岩心地面孔隙度对应的覆压渗透率，即为地层渗透率。

在一个示例中，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线包括：

确定区域标准层，获取区域标准层内直井标准层的中子值、密度值和深度值；

建立直井标准层的中子值随深度变化的关系式以及密度值随深度变化的关系式；

将测量深度校正到垂直深度，获取大斜度井标准层的中子理论值与密度理论值；

获取大斜度井标准层的中子实测值与密度实测值；

计算大斜度井标准层的中子校正量与密度校正量；

根据中子校正量与密度校正量对大斜度井斜井段的中子曲线和密度曲线进行校正，获得校正后曲线。

在一个示例中，通过公式(1)计算大斜度井标准层的中子校正量：

ΔCNL＝CNL1–CNL2 (1)

其中，ΔCNL为中子校正量，CNL1为中子理论值，CNL2为中子实测值。

在一个示例中，通过公式(2)计算大斜度井标准层的密度校正量：

ΔDEN＝DEN1–DEN2 (2)

其中，ΔDEN为密度校正量，DEN1为密度理论值，DEN2为密度实测值。

在一个示例中，通过校正后曲线，采用中子-密度交会法进行储层孔隙度计算：

其中，φ_t为岩心地面孔隙度，CNL为环境校正后的中子值，DEN为环境校正后的密度值。

在一个示例中，第一转换公式为：

其中，φ_t为岩心地面孔隙度，φ为覆压孔隙度，E和F为拟合系数。

在一个示例中，第二转换公式为：

K＝G×e^H×φ (5)

其中，φ为覆压孔隙度，K为覆压渗透率，G和H为拟合系数。

具体地，对大斜度井的中子和密度测井曲线进行环境校正：

选取区域标准层，并获取区域标准层内直井标准层的中子值、密度值和深度值。选取岩性纯、厚度大、分布稳定的泥岩层作为区域标准层。为保证每口直井标准层读值的准确性，每口井的读值层段不少于5个。按照深度顺序，读取每个层段内中子和密度曲线的平均值，深度值取每个层段的中部深度值。

建立直井标准层的中子值随深度变化的关系式以及密度值随深度变化的关系式。

将中子值与深度值做线性回归，建立标准层中子值随深度变化的关系式：

CNL＝A*DEPTH+B

式中：CNL为中子值，DEPTH为深度值，A、B为常数。

将密度值与深度值做线性回归，建立标准层密度值随深度变化的关系式：

DEN＝C*DEPTH+D

式中：DEN为密度值，DEPTH为深度值，C、D为常数。

将大斜度井标准层内测量深度校正到垂直深度，获取大斜度井标准层的中子理论值与密度理论值。

将大斜度井标准层内任一深度点的测量深度校正到垂直深度，深度值的采样间隔为0.125m，与补偿密度、密度曲线实际测量值的采样间隔一致。如果大斜度井的标准层在垂直井段内，则不需要进行垂深校正。

将大斜度井标准层内任一深度点的垂直深度值代入中子和深度关系式，计算得到标准层每个垂直深度点的中子理论值CNL1；将大斜度井标准层内任一深度点的垂直深度值代入密度与深度关系式，计算得到标准层每个垂直深度点的密度理论值DEN1。

分别读取大斜度井标准层中子实测值CNL2和密度实测值DEN2。

根据大斜度井标准层中子的理论值和实测值获取大斜度井斜井段中子的校正量，即为公式(1)，并根据大斜度井标准层密度的理论值和实际测量值获取大斜度井斜井段密度的校正量，即为公式(2)。

对大斜度井斜井段的中子曲线和密度曲线进行校正，将中子校正量ΔCNL加到大斜度井斜井段中子实际测量值，即完成中子曲线的环境校正；将密度校正量ΔDEN加到大斜度井斜井段密度实际测量值，即完成密度曲线的环境校正。

通过足够多的取心资料进行岩心归位，利用环境校正后中子和密度测井曲线，采用中子-密度交会法进行储层孔隙度计算，并利用岩心分析地面孔隙度数据校正孔隙度的计算结果，此时得到的孔隙度相当于未做覆压校正的岩心地面孔隙度。

选取与上覆岩压和地层压力一致或者接近的覆压岩石实验数据，构建地面岩心孔隙度和覆压孔隙度的交会图版，得到岩心地面孔隙度和覆压孔隙度的第一转换公式为公式(4)，计算岩心地面孔隙度对应的覆压孔隙度。

选择地层条件下的岩心覆压孔隙度和覆压渗透率实验数据，构建覆压孔隙度和覆压渗透率交会图版，得到覆压孔隙度和覆压渗透率的第二转换公式为公式(5)。

将岩心地面孔隙度对应的覆压孔隙度带入第二转换公式，获得岩心地面孔隙度对应的覆压渗透率，计算得到的覆压渗透率即为弱成岩储层的地层渗透率。

本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述的地层渗透率测井计算方法。

本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述的地层渗透率测井计算方法。

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

实施例1

图1示出了根据本发明的地层渗透率测井计算方法的步骤的流程图。

如图1所示，该地层渗透率测井计算方法包括：步骤101，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线；步骤102，通过校正后曲线，计算岩心地面孔隙度；步骤103，建立岩心地面孔隙度和覆压孔隙度的第一转换公式，计算岩心地面孔隙度对应的覆压孔隙度；步骤104，建立覆压孔隙度和覆压渗透率的第二转换公式；步骤105，将岩心地面孔隙度对应的覆压孔隙度带入第二转换公式，获得岩心地面孔隙度对应的覆压渗透率，即为地层渗透率。

以某区某井为实施例，对大斜度井的中子和密度测井曲线进行环境校正：

选取区域标准层，并获取区域标准层内直井标准层的中子值、密度值和深度值。

选取岩性纯、厚度大、分布稳定的流沙港组流二段泥页岩层作为区域标准层。读取了研究区内涠A、涠B、涠C、涠D、涠E等5口直井流二段的中子值、密度值和深度值，每口井的读值层段均不少于5个。按照深度顺序，读取每个层段内中子和密度曲线的平均值，深度值取每个层段的中部深度值。

建立直井标准层的中子值随深度变化的关系式以及密度值随深度变化的关系式。

图2示出了根据本发明的一个实施例的标准层的中子值与深度变化关系示意图。

将中子值与深度值做线性回归，如图2所示，建立标准层中子值随深度变化的关系式：

CNL＝-0.0109*DEPTH+59.475

图3示出了根据本发明的一个实施例的标准层的密度值与深度变化关系示意图。

将密度值与深度值做线性回归，如图3所示，建立标准层密度值随深度变化的关系式：

DEN＝0.0001*DEPTH+2.1372

将大斜度井标准层内测量深度校正到垂直深度，获取大斜度井标准层的中子理论值与密度理论值。

将标准层内任一深度点的测量深度校正到垂直深度，深度值的采样间隔为0.125m，与补偿密度、密度曲线实际测量值的采样间隔一致。将标准层内任一深度点的垂直深度值代入中子和深度关系式，计算得到标准层每个垂直深度点的中子理论值；将标准层内任一深度点的垂直深度值代入密度与深度关系式，计算得到标准层每个垂直深度点的密度理论值。

获取大斜度井标准层中子实测值以及密度实测值。根据大斜度井标准层中子的理论值和实测值获取大斜度井斜井段中子的校正量ΔCNL，并根据大斜度井标准层密度的理论值和实测值获取大斜度井斜井段密度的校正量ΔDEN。

求取标准层中子理论值的平均值CNL1和密度理论值的平均值DEN1，计算得到CNL1等于29.5％，DEN1等于2.49g/cm³。

求取标准层中子实测值的平均值CNL2和密度实测值的平均值DEN2，计算得到CNL2等于24.1％，DEN1等于2.60g/cm³。

计算斜井段中子值的校正量ΔCNL、密度值的校正量ΔDEN

ΔCNL＝CNL1–CNL2＝29.5-24.1＝5.4(％)

ΔDEN＝DEN1–DEN2＝2.49–2.60＝-0.11(g/cm³)

图4示出了根据本发明的一个实施例的环境校正后的中子、密度曲线的示意图。

对大斜度井斜井段的中子和密度曲线进行环境校正，如图4所示。

将得到的中子校正量5.4％加到斜井段中子实测值，即完成中子曲线的环境校正。将得到的密度校正量-0.11g/cm³加到斜井段密度实测值，即完成密度曲线的环境校正。

通过足够多的取心资料进行岩心归位，利用环境校正后中子和密度测井曲线，采用中子-密度交会法进行1824-1860m储层孔隙度计算，并利用岩心分析地面孔隙度数据校正孔隙度的计算结果。此时得到的孔隙度相当于未做覆压校正的岩心地面孔隙度，计算得到1824-1860m平均孔隙度27.6％，与岩心地面孔隙度平均值27.4％基本一致。

图5示出了根据本发明的一个实施例的岩心地面孔隙度和覆压孔隙度的关系示意图。

涠西区块弱成岩储层的地层压力平均为20Mpa，选取上覆岩压为20Mpa的覆压岩石实验数据，构建地面岩心孔隙度和覆压孔隙度的交会图版，如图5所示，得到覆压孔隙度φ计算公式，将地面孔隙度曲线转换为覆压孔隙度曲线：

转换得到涠*井1824-1860m的测井计算覆压孔隙度平均为23.87％，与岩心覆压孔隙度平均值23.94％相差误差很小，绝对误差小于1％，表明精度较高。

图6示出了根据本发明的一个实施例的覆压孔隙度和覆压渗透率的关系示意图。

选择地层条件下的岩心覆压孔隙度和覆压渗透率实验数据，构建覆压孔隙度和覆压渗透率交会图版，如图6所示，得到弱成岩储层的覆压渗透率K计算公式：

K＝0.0847×e^0.3331×φ

图7示出了根据本发明的一个实施例的测井计算的地层孔隙度、渗透率与覆压孔隙度、覆压渗透率的对比示意图。

将测井计算的储层孔隙度转换为地层条件下的覆压孔隙度，将覆压孔隙度曲线转换为覆压渗透率曲线，计算得到的覆压渗透率即为弱成岩储层的地层渗透率，如图7所示。计算结果表明，岩心地面渗透率平均为389mD，覆压渗透率平均为289mD，测井计算的平均地层渗透率301mD，与覆压渗透率的绝对误差12mD，相对误差小于5％，满足此类弱成岩油气藏的储层特征和产能评价的精度要求。

实施例2

图8示出了根据本发明的一个实施例的一种地层渗透率测井计算装置的框图。

如图8所示，该地层渗透率测井计算装置，包括：

环境校正模块201，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线；

岩心地面孔隙度计算模块202，通过校正后曲线，计算岩心地面孔隙度；

第一转换模块203，建立岩心地面孔隙度和覆压孔隙度的第一转换公式，计算岩心地面孔隙度对应的覆压孔隙度；

第二转换模块204，建立覆压孔隙度和覆压渗透率的第二转换公式；

预测模块205，将岩心地面孔隙度对应的覆压孔隙度带入第二转换公式，获得岩心地面孔隙度对应的覆压渗透率，即为地层渗透率。

作为可选方案，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线包括：

确定区域标准层，获取区域标准层内直井标准层的中子值、密度值和深度值；

建立直井标准层的中子值随深度变化的关系式以及密度值随深度变化的关系式；

将测量深度校正到垂直深度，获取大斜度井标准层的中子理论值与密度理论值；

获取大斜度井标准层的中子实测值与密度实测值；

计算大斜度井标准层的中子校正量与密度校正量；

根据中子校正量与密度校正量对大斜度井斜井段的中子曲线和密度曲线进行校正，获得校正后曲线。

作为可选方案，通过公式(1)计算大斜度井标准层的中子校正量：

ΔCNL＝CNL1–CNL2 (1)

其中，ΔCNL为中子校正量，CNL1为中子理论值，CNL2为中子实测值。

作为可选方案，通过公式(2)计算大斜度井标准层的密度校正量：

ΔDEN＝DEN1–DEN2 (2)

其中，△DEN为密度校正量，DEN1为密度理论值，DEN2为密度实测值。

作为可选方案，通过校正后曲线，采用中子-密度交会法进行储层孔隙度计算：

其中，φ_t为岩心地面孔隙度，CNL为环境校正后的中子值，DEN为环境校正后的密度值。

作为可选方案，第一转换公式为：

其中，φ_t为岩心地面孔隙度，φ为覆压孔隙度，E和F为拟合系数。

作为可选方案，第二转换公式为：

K＝G×e^H×φ (5)

其中，φ为覆压孔隙度，K为覆压渗透率，G和H为拟合系数。

实施例3

本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述地层渗透率测井计算方法。

根据本公开实施例的电子设备包括存储器和处理器。

该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。

该处理器可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。

本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。

有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。

实施例4

本公开实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的地层渗透率测井计算方法。

根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。

上述计算机可读存储介质包括但不限于：光存储介质(例如：CD－ROM和DVD)、磁光存储介质(例如：MO)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置ROM的媒体(例如：ROM盒)。

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

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

Claims (10)

1.一种地层渗透率测井计算方法，其特征在于，包括：

针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线；

通过所述校正后曲线，计算岩心地面孔隙度；

建立所述岩心地面孔隙度和覆压孔隙度的第一转换公式，计算所述岩心地面孔隙度对应的覆压孔隙度；

建立所述覆压孔隙度和覆压渗透率的第二转换公式；

将所述岩心地面孔隙度对应的覆压孔隙度带入所述第二转换公式，获得所述岩心地面孔隙度对应的覆压渗透率，即为所述地层渗透率。

2.根据权利要求1所述的地层渗透率测井计算方法，其中，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线包括：

确定区域标准层，获取所述区域标准层内直井标准层的中子值、密度值和深度值；

建立所述直井标准层的中子值随深度变化的关系式以及密度值随深度变化的关系式；

将测量深度校正到垂直深度，获取大斜度井标准层的中子理论值与密度理论值；

获取大斜度井标准层的中子实测值与密度实测值；

计算大斜度井标准层的中子校正量与密度校正量；

根据所述中子校正量与所述密度校正量对大斜度井斜井段的中子曲线和密度曲线进行校正，获得所述校正后曲线。

3.根据权利要求2所述的地层渗透率测井计算方法，其中，通过公式(1)计算大斜度井标准层的中子校正量：

ΔCNL＝CNL1–CNL2 (1)

其中，ΔCNL为中子校正量，CNL1为中子理论值，CNL2为中子实测值。

4.根据权利要求2所述的地层渗透率测井计算方法，其中，通过公式(2)计算大斜度井标准层的密度校正量：

ΔDEN＝DEN1–DEN2 (2)

其中，ΔDEN为密度校正量，DEN1为密度理论值，DEN2为密度实测值。

5.根据权利要求1所述的地层渗透率测井计算方法，其中，通过所述校正后曲线，采用中子-密度交会法进行储层孔隙度计算：

其中，φ_t为岩心地面孔隙度，CNL为环境校正后的中子值，DEN为环境校正后的密度值。

6.根据权利要求1所述的地层渗透率测井计算方法，其中，所述第一转换公式为：

其中，φ_t为岩心地面孔隙度，φ为覆压孔隙度，E和F为拟合系数。

7.根据权利要求1所述的地层渗透率测井计算方法，其中，所述第二转换公式为：

K＝G×e^H×φ (5)

其中，φ为覆压孔隙度，K为覆压渗透率，G和H为拟合系数。

8.一种地层渗透率测井计算装置，其特征在于，包括：

环境校正模块，针对大斜度井的中子和密度测井曲线进行环境校正，获得校正后曲线；

岩心地面孔隙度计算模块，通过所述校正后曲线，计算岩心地面孔隙度；

第一转换模块，建立所述岩心地面孔隙度和覆压孔隙度的第一转换公式，计算所述岩心地面孔隙度对应的覆压孔隙度；

第二转换模块，建立所述覆压孔隙度和覆压渗透率的第二转换公式；

预测模块，将所述岩心地面孔隙度对应的覆压孔隙度带入所述第二转换公式，获得所述岩心地面孔隙度对应的覆压渗透率，即为所述地层渗透率。

9.一种电子设备，其特征在于，所述电子设备包括：

存储器，存储有可执行指令；

处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-7中任一项所述的地层渗透率测井计算方法。

10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-7中任一项所述的地层渗透率测井计算方法。