CN105240001B - 核磁共振测井孔隙度校正方法及装置 - Google Patents
核磁共振测井孔隙度校正方法及装置 Download PDFInfo
- Publication number
- CN105240001B CN105240001B CN201510614235.4A CN201510614235A CN105240001B CN 105240001 B CN105240001 B CN 105240001B CN 201510614235 A CN201510614235 A CN 201510614235A CN 105240001 B CN105240001 B CN 105240001B
- Authority
- CN
- China
- Prior art keywords
- mrow
- porosity
- pebble
- magnetic resonance
- msub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例提供一种核磁共振测井孔隙度校正方法及装置。该方法包括：获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度；依据目标岩石大小、砾石组分大小、非砾石组分大小以及孔隙度的定义建立核磁共振测井孔隙度的校正模型；依据井壁电成像测井图像获取砾石组分在目标岩石中的比例；依据砾石组分在目标岩石中的比例和校正模型计算核磁共振校正孔隙度，核磁共振校正孔隙度是核磁共振测井孔隙度的校正值。本发明实施例利用电成像测井技术获得的井壁电成像测井图像对核磁共振测井孔隙度进行校正，获得核磁共振测井孔隙度的校正值，避免了核磁共振测井确定的储层孔隙度比地层实际的储层孔隙度偏小，提高了对油气藏的评价精度。
Description
技术领域
本发明实施例涉及石油探测领域，尤其涉及一种核磁共振测井孔隙度校正方法及装置。
背景技术
随着石油勘探技术的不断发展，以及勘探对象的日益复杂，核磁共振测井作为一种新的测井技术，成为确定储层孔隙度的重要的测井方法之一。储层孔隙度是指岩样中所有孔隙空间体积之和与该岩样体积的比值，是衡量储层岩石中所含孔隙体积多少的一种参数，反映着储层岩石储存流体的能力。因此，储层孔隙度的准确计算对油气藏的精细评价具有十分重要的意义。
核磁共振测井的测量对象为储层孔隙中流体如氢核，其测量结果基本不受岩性、岩石骨架等因素的影响，因此，在确定储层孔隙度方面具有其他测井方法无法比拟的优势。
但是，当储层岩性为砂砾岩时，利用核磁共振测井确定的储层孔隙度比地层实际的储层孔隙度偏小，导致无法对油气藏进行精细评价。
发明内容
本发明实施例提供一种核磁共振测井孔隙度校正方法及装置，以校正核磁共振测井孔隙度，提高对油气藏的评价精度。
本发明实施例的一个方面是提供一种核磁共振测井孔隙度校正方法，包括：
获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，所述核磁共振测井孔隙度表示孔隙空间在所述目标岩石中的比例，所述目标岩石包括砾石组分和所述非砾石组分，所述非砾石组分包括所述孔隙空间；
依据所述目标岩石的大小、所述砾石组分的大小和所述非砾石组分的大小以及孔隙度的定义建立所述核磁共振测井孔隙度的校正模型；
依据所述井壁电成像测井图像获取所述砾石组分在所述目标岩石中的比例；
依据所述砾石组分在所述目标岩石中的比例和所述校正模型计算核磁共振校正孔隙度，所述核磁共振校正孔隙度是所述核磁共振测井孔隙度的校正值。
本发明实施例的另一个方面是提供一种核磁共振测井孔隙度校正装置，包括：
获取模块，用于获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，所述核磁共振测井孔隙度表示孔隙空间在所述目标岩石中的比例，所述目标岩石包括砾石组分和非砾石组分，所述非砾石组分包括所述孔隙空间；依据所述井壁电成像测井图像获取所述砾石组分在所述目标岩石中的比例；
建模模块，用于依据所述目标岩石的大小、所述砾石组分的大小和所述非砾石组分的大小以及孔隙度的定义建立所述核磁共振测井孔隙度的校正模型；
计算模块，用于依据所述砾石组分在所述目标岩石中的比例和所述校正模型计算核磁共振校正孔隙度，所述核磁共振校正孔隙度是所述核磁共振测井孔隙度的校正值。
本发明实施例提供的核磁共振测井孔隙度校正方法及装置，通过获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，利用电成像测井技术获得的井壁电成像测井图像对核磁共振测井孔隙度进行校正，获得核磁共振测井孔隙度的校正值，避免了核磁共振测井确定的储层孔隙度比地层实际的储层孔隙度偏小，提高了对油气藏的评价精度。
附图说明
图1为本发明实施例提供的核磁共振测井孔隙度校正方法流程图；
图2为本发明实施例提供的井壁微电阻率图像示意图；
图3为本发明实施例提供的核磁共振测井的原始数据的示意图；
图4为本发明实施例提供的核磁共振测井孔隙度校正方法适用的检测区域；
图5为本发明另一实施例提供的核磁共振测井孔隙度校正方法适用的检测区域；
图6为本发明实施例提供的核磁共振测井孔隙度与岩心分析孔隙度校正前的交会图；
图7为本发明实施例提供的核磁共振测井孔隙度与岩心分析孔隙度校正后的交会图；
图8为本发明实施例提供的核磁共振测井孔隙度校正装置的结构；
图9为本发明另一实施例提供的核磁共振测井孔隙度校正装置的结构图。
具体实施方式
图1为本发明实施例提供的核磁共振测井孔隙度校正方法流程图；图2为本发明实施例提供的井壁微电阻率图像示意图；图3为本发明实施例提供的核磁共振测井的原始数据的示意图。本发明实施例针对储层岩性为砂砾岩时，利用核磁共振测井确定的储层孔隙度比地层实际的储层孔隙度偏小，提供了核磁共振测井孔隙度校正方法，该方法的具体步骤如下：
步骤S101、获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，所述核磁共振测井孔隙度表示孔隙空间在所述目标岩石中的比例，所述目标岩石包括砾石组分和非砾石组分，所述非砾石组分包括所述孔隙空间；
所述获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，包括：利用电成像测井仪测量所述目标岩石对应的地层电阻率，依据所述地层电阻率获得所述井壁电成像测井图像；利用所述核磁共振测井仪测量所述目标岩石对应的核磁共振信号，依据所述核磁共振信号获得所述核磁共振测井孔隙度。
电成像测井仪利用井下仪器极板上密集排列的钮扣状电极测量井壁圆周上几十条甚至上百条微电阻率信息，对该微电阻率信息经过高分辨率数学成像处理产生一幅高清晰度的井壁微电阻率图像，该井壁微电阻率图像如图2所示，对该井壁微电阻率图像经过现有的测井解释软件进行分析处理获得井壁电成像测井图像。
在核磁共振测井之前，地层中的质子是随机取向排列的，当核磁共振测井仪经过地层时，核磁共振测井仪的磁场使质子磁化。首先核磁共振测井仪的永久磁场使质子的旋转轴朝向一致。然后，核磁共振测井仪发射交变磁场使这些被极化的质子从新的平衡位置翻转。当核磁共振测井仪停止发射交变磁场后，质子就开始进动回到静磁场使之极化的位置，这一过程称为弛豫。核磁共振测井仪利用特定的脉冲序列采集弛豫过程中一系列自旋回波，该自旋回波即为核磁共振测井的原始数据，核磁共振测井的原始数据如图3所示，对该核磁共振测井的原始数据经过现有的测井解释软件进行分析处理获得核磁共振测井孔隙度。
步骤S102、依据所述目标岩石的大小、所述砾石组分的大小和所述非砾石组分的大小以及孔隙度的定义建立所述核磁共振测井孔隙度的校正模型；
依据以下三个特性建立核磁共振测井孔隙度的校正模型：1)目标岩石的大小等于砾石组分的大小和非砾石组分的大小的和(孔隙存在于非砾石组分中)；2)核磁共振测井孔隙度表示孔隙空间在所述目标岩石中的比例；3)砾石组分的孔隙度为0，即砾石组分中没有孔隙。
步骤S103、依据所述井壁电成像测井图像获取所述砾石组分在所述目标岩石中的比例；
所述井壁电成像测井图像包括多个图像点；所述依据所述井壁电成像测井图像获取所述砾石组分在所述目标岩石中的比例，包括：依据各个图像点的亮暗程度分别对每个图像点进行标识，且所述图像点的亮暗程度与标识值成正比；统计所述标识值大于阈值的图像点的个数，所述个数与所述图像点的总个数的比值是所述砾石组分在所述目标岩石中的比例。
在本发明实施例中井壁电成像测井图像包括多个图像点，每个图像点亮暗程度不同，依据各个图像点的亮暗程度分别对每个图像点进行标识，标识值在0-256之间，图像点的亮度越大，标识值越大，图像点的亮度越小，标识值越小，将砾石的电阻率临界值作为井壁电成像测井图像对应的图像亮度临界值，并将该图像亮度临界值作为标识值的阈值，则标识值大于阈值的图像点对应砾石组分，标识值小于阈值的图像点对应非砾石组分即孔隙空间。标识值大于阈值的图像点的个数与图像点的总个数的比值是所述砾石组分在所述目标岩石中的比例。
步骤S104、依据所述砾石组分在所述目标岩石中的比例和所述校正模型计算核磁共振校正孔隙度，所述核磁共振校正孔隙度是所述核磁共振测井孔隙度的校正值。
将S103获得的砾石组分在所述目标岩石中的比例和步骤S101中获取到的核磁共振测井孔隙度带入到校正模型中即可计算出核磁共振校正孔隙度，核磁共振校正孔隙度即是所述核磁共振测井孔隙度的校正值。
本发明实施例通过获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，利用电成像测井技术获得的井壁电成像测井图像对核磁共振测井孔隙度进行校正，获得核磁共振测井孔隙度的校正值，避免了核磁共振测井确定的储层孔隙度比地层实际的储层孔隙度偏小，提高了对油气藏的评价精度。
图4为本发明实施例提供的核磁共振测井孔隙度校正方法适用的检测区域。在上述实施例的基础上，所述校正模型的检测对象是长度为2πR、宽度为H、厚度为D的目标岩石；
所述目标岩石的大小为所述目标岩石的体积，所述砾石组分的大小为所述砾石组分的体积，所述孔隙空间的大小为所述孔隙空间的体积，且所述目标岩石的体积是所述砾石组分的体积与所述非砾石组分的体积的和；
所述校正模型对应公式(1)、(2)、(3)：
Vgravel+Vnon＝2πR·H·D (1)
φgravel＝0 (3)
其中，Vgravel表示所述砾石组分的体积，Vnon表示所述非砾石组分的体积，2πR·H·D表示所述目标岩石的体积，φgravel表示所述砾石组分的孔隙度，φreal表示所述核磁共振校正孔隙度，φNMR表示所述核磁共振测井孔隙度。
在本发明实施例中R具体为目标岩石对应的目标井的井眼半径，取值范围为20cm-25cm；H为核磁共振测井仪的纵向分辨率，取值范围为1m-1.2m；D为核磁共振测井仪的敏感区域的厚度，取值范围为1mm-2mm。
将公式(3)带到公式(2)中获得进一步获得由公式(1)可知Vnon＝2πR·H·D-Vgravel，将Vnon＝2πR·H·D-Vgravel和公式(1)带入公式中，可获得Vgravel/2πR·H·D表示所述砾石组分的体积占所述目标岩石的体积的比例(其值是小于1)，上述步骤S103中标识值大于阈值的图像点的个数与图像点的总个数的比值等于所述砾石组分的体积占所述目标岩石的体积的比例，即将上述步骤S103中获得的标识值大于阈值的图像点的个数与图像点的总个数的比值作为Vgravel/2πR·H·D便可计算获得核磁共振校正孔隙度φreal。
图5为本发明另一实施例提供的核磁共振测井孔隙度校正方法适用的检测区域。由于核磁探测区域的圆柱壳的厚度很小(D很小)，认为径向均质(D近似等于0)，故将所述的长度为2πR、宽度为H、厚度为D的目标岩石近似看做长度为2πR、宽度为H的平面，体积计算等效为面积计算，在上述实施例的基础上，所述校正模型的检测对象是长度为2πR、宽度为H的横截面；
所述目标岩石的大小为所述横截面的面积，所述砾石组分的大小为所述横截面中所述砾石组分的面积，所述非砾石组分的大小为所述横截面中所述非砾石组分的面积，且所述横截面的面积等于所述砾石组分的面积和所述非砾石组分的面积的和；
所述校正模型对应公式(4)、(5)、(6)：
Sgravel+Snon＝2πR·H (4)
φgravel＝0 (6)
其中，Sgravel表示所述砾石组分的面积，Snon表示所述非砾石组分的面积，2πR·H表示所述横截面的面积，φgravel表示所述砾石组分的孔隙度，φreal表示所述核磁共振校正孔隙度，φNMR表示所述核磁共振测井孔隙度。
在本发明实施例中R具体为目标岩石对应的目标井的井眼半径，取值范围为20cm-25cm；H为核磁共振测井仪的纵向分辨率，取值范围为1m-1.2m。
将公式(6)带到公式(5)中获得进一步获得由公式(4)可知Snon＝2πR·H-Sgravel，将Snon＝2πR·H-Sgravel和公式(4)带入公式中，可获得sgravel/2πR·H表示所述砾石组分的面积占所述横截面的面积的比例，上述步骤S103中标识值大于阈值的图像点的个数与图像点的总个数的比值等于所述砾石组分的面积占所述横截面的面积的比例，即将上述步骤S103中获得的标识值大于阈值的图像点的个数与图像点的总个数的比值作为sgravel/2πR·H便可计算获得核磁共振校正孔隙度φreal。
本发明实施例通过检测对象是长度为2πR、宽度为H、厚度为D的目标岩石，以及检测对象是长度为2πR、宽度为H、厚度为D的目标岩石沿厚度方向的横截面，分别提供了两种校正模型，通过两种校正模型获得两种核磁共振测井孔隙度的校正方法，实现了对核磁共振测井孔隙度的校正。
图6为本发明实施例提供的核磁共振测井孔隙度与岩心分析孔隙度校正前的交会图；图7为本发明实施例提供的核磁共振测井孔隙度与岩心分析孔隙度校正后的交会图。在上述实施例的基础上，所述依据所述目标岩石的大小、所述砾石组分的大小和所述非砾石组分的大小以及孔隙度的定义建立所述核磁共振测井孔隙度的校正模型之前，还包括：获取所述目标岩石对应的岩心分析孔隙度；依据交会图分析方法对所述岩心分析孔隙度和所述核磁共振测井孔隙度进行误差分析，判断获知所述核磁共振测井孔隙度小于所述岩心分析孔隙度。
本发明实施例对目标岩石中的同一个检测点获取其核磁共振测井孔隙度与岩心分析孔隙度，将每个检测点的核磁共振测井孔隙度对应到X轴，将每个检测点的岩心分析孔隙度对应到Y轴，获得如图6所示的坐标图，并在所示坐标图中画一条45度角的线，多数检测点落在45度-90度范围内，表示核磁共振测井孔隙度小于岩心分析孔隙度，以岩心分析孔隙度为地层实际的储层孔隙度，则核磁共振测井孔隙度比地层实际的储层孔隙度偏小，需要对核磁共振测井孔隙度进行校正。图7所示为利用本发明实施例提供的核磁共振测井孔隙度校正方法校正后的核磁共振测井孔隙度即核磁共振校正孔隙度与岩心分析孔隙度的交会图，如图7所示，检测点均匀分布在45度线的周围，表示通过本发明实施例提供的核磁共振测井孔隙度校正方法精确校正了核磁共振测井孔隙度。
本发明实施例通过交会图分析方法对目标岩石的岩心分析孔隙度和核磁共振测井孔隙度进行误差分析，通过预先判断核磁共振测井孔隙度是否小于岩心分析孔隙度，进一步确定是否需要对核磁共振测井孔隙度进行校正，避免多次校正操作，提高了校正效率。
图8为本发明实施例提供的核磁共振测井孔隙度校正装置的结构图。本发明实施例提供的核磁共振测井孔隙度校正装置可以执行核磁共振测井孔隙度校正方法实施例提供的处理流程，如图8所示，核磁共振测井孔隙度校正装置80包括获取模块81、建模模块82和计算模块83，其中，获取模块81用于获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，所述核磁共振测井孔隙度表示孔隙空间在所述目标岩石中的比例，所述目标岩石包括砾石组分和非砾石组分，所述非砾石组分包括所述孔隙空间；依据所述井壁电成像测井图像获取所述砾石组分在所述目标岩石中的比例；建模模块82用于依据所述目标岩石的大小、所述砾石组分的大小和所述非砾石组分的大小以及孔隙度的定义建立所述核磁共振测井孔隙度的校正模型；计算模块83用于依据所述砾石组分在所述目标岩石中的比例和所述校正模型计算核磁共振校正孔隙度，所述核磁共振校正孔隙度是所述核磁共振测井孔隙度的校正值。
本发明实施例通过获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，利用电成像测井技术获得的井壁电成像测井图像对核磁共振测井孔隙度进行校正，获得核磁共振测井孔隙度的校正值，避免了核磁共振测井确定的储层孔隙度比地层实际的储层孔隙度偏小，提高了对油气藏的评价精度。
图9为本发明另一实施例提供的核磁共振测井孔隙度校正装置的结构图。在上述实施例的基础上，所述校正模型的检测对象是长度为2πR、宽度为H、厚度为D的目标岩石；所述目标岩石的大小为所述目标岩石的体积，所述砾石组分的大小为所述砾石组分的体积，所述孔隙空间的大小为所述孔隙空间的体积，且所述目标岩石的体积是所述砾石组分的体积与所述非砾石组分的体积的和；所述校正模型对应公式(1)、(2)、(3)：
Vgravel+Vnon＝2πR·H·D (1)
φgravel＝0 (3)
其中，Vgravel表示所述砾石组分的体积，Vnon表示所述非砾石组分的体积，2πR·H·D表示所述目标岩石的体积，φgravel表示所述砾石组分的孔隙度，φreal表示所述核磁共振校正孔隙度，φNMR表示所述核磁共振测井孔隙度。
所述校正模型的检测对象是长度为2πR、宽度为H的横截面；所述目标岩石的大小为所述横截面的面积，所述砾石组分的大小为所述横截面中所述砾石组分的面积，所述非砾石组分的大小为所述横截面中所述非砾石组分的面积，且所述横截面的面积等于所述砾石组分的面积和所述孔隙空间的面积的和；所述校正模型对应公式(4)、(5)、(6)：
Sgravel+Snon＝2πR·H (4)
φgravel＝0 (6)
其中，Sgravel表示所述砾石组分的面积，Snon表示所述非砾石组分的面积，2πR·H表示所述横截面的面积，φgravel表示所述砾石组分的孔隙度，φreal表示所述核磁共振校正孔隙度，φNMR表示所述核磁共振测井孔隙度。
获取模块81具体用于利用电成像测井仪测量所述目标岩石对应的地层电阻率，依据所述地层电阻率获得所述井壁电成像测井图像；利用所述核磁共振测井仪测量所述目标岩石对应的核磁共振信号，依据所述核磁共振信号获得所述核磁共振测井孔隙度。
获取模块81还用于获取所述目标岩石对应的岩心分析孔隙度；核磁共振测井孔隙度校正装置80还包括误差分析模块84，误差分析模块84用于依据交会图分析方法对所述岩心分析孔隙度和所述核磁共振测井孔隙度进行误差分析，判断获知所述核磁共振测井孔隙度小于所述岩心分析孔隙度。
所述井壁电成像测井图像包括多个图像点；获取模块81包括标识单元811和统计单元812，其中，标识单元811用于依据各个图像点的亮暗程度分别对每个图像点进行标识，且所述图像点的亮暗程度与标识值成正比；统计单元812用于统计所述标识值大于阈值的图像点的个数，所述个数与所述图像点的总个数的比值是所述砾石组分在所述目标岩石中的比例。
本发明实施例提供的核磁共振测井孔隙度校正装置可以具体用于执行上述图1所提供的方法实施例，具体功能此处不再赘述。
本发明实施例将检测对象是长度为2πR、宽度为H、厚度为D的目标岩石近演化为长度为2πR、宽度为H的横截面，依据该横截面提供了校正模型，通过该校正模型获得了核磁共振测井孔隙度的校正方法，实现了对核磁共振测井孔隙度的校正；通过交会图分析方法对目标岩石的岩心分析孔隙度和核磁共振测井孔隙度进行误差分析，通过预先判断核磁共振测井孔隙度是否小于岩心分析孔隙度，进一步确定是否需要对核磁共振测井孔隙度进行校正，避免多次校正操作，提高了校正效率。
本发明实施例通过获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，利用电成像测井技术获得的井壁电成像测井图像对核磁共振测井孔隙度进行校正，获得核磁共振测井孔隙度的校正值，避免了核磁共振测井确定的储层孔隙度比地层实际的储层孔隙度偏小，提高了对油气藏的评价精度；通过检测对象是长度为2πR、宽度为H、厚度为D的目标岩石，以及检测对象是长度为2πR、宽度为H、厚度为D的目标岩石沿厚度方向的横截面，分别提供了两种校正模型，通过两种校正模型获得两种核磁共振测井孔隙度的校正方法，实现了对核磁共振测井孔隙度的校正；通过交会图分析方法对目标岩石的岩心分析孔隙度和核磁共振测井孔隙度进行误差分析，通过预先判断核磁共振测井孔隙度是否小于岩心分析孔隙度，进一步确定是否需要对核磁共振测井孔隙度进行校正，避免多次校正操作，提高了校正效率。
综上所述，本发明实施例通过发送设备向接收设备发送NDM文件对应的完整性信息，使得接收设备能够依据所述完整性信息验证所述NDM文件是否被篡改，提高了NDM文件的安全性；通过向接收设备发送NDM文件的有效性信息，使得接收设备能够确定所述NDM文件的有效期限，提高了NDM文件的有效性；通过发送设备向接收设备发送预设算法信息和证书，使得接收设备准确验证所述NDM文件是否被篡改。
在本发明所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元，可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或处理器(processor)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(Random Access Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
本领域技术人员可以清楚地了解到，为描述的方便和简洁，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。上述描述的装置的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (12)
1.一种核磁共振测井孔隙度校正方法，其特征在于，包括：
获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，所述核磁共振测井孔隙度表示孔隙空间在所述目标岩石中的比例，所述目标岩石包括砾石组分和非砾石组分，所述非砾石组分包括所述孔隙空间；
依据所述目标岩石的大小、所述砾石组分的大小和所述非砾石组分的大小以及孔隙度的定义建立所述核磁共振测井孔隙度的校正模型；
依据所述井壁电成像测井图像获取所述砾石组分在所述目标岩石中的比例；
依据所述砾石组分在所述目标岩石中的比例和所述校正模型计算核磁共振校正孔隙度，所述核磁共振校正孔隙度是所述核磁共振测井孔隙度的校正值。
2.根据权利要求1所述的方法，其特征在于，所述校正模型的检测对象是长度为2πR、宽度为H、厚度为D的目标岩石；
所述目标岩石的大小为所述目标岩石的体积，所述砾石组分的大小为所述砾石组分的体积，所述孔隙空间的大小为所述孔隙空间的体积，且所述目标岩石的体积是所述砾石组分的体积与所述非砾石组分的体积的和；
所述校正模型对应公式(1)、(2)、(3)：
Vgravel+Vnon＝2πR·H·D (1)
<mrow>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>N</mi>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
φgravel＝0 (3)
其中，Vgravel表示所述砾石组分的体积，Vnon表示所述非砾石组分的体积，2πR·H·D表示所述目标岩石的体积，φgravel表示所述砾石组分的孔隙度，φreal表示所述核磁共振校正孔隙度，φNMR表示所述核磁共振测井孔隙度。
3.根据权利要求1所述的方法，其特征在于，所述校正模型的检测对象是长度为2πR、宽度为H的横截面；
所述目标岩石的大小为所述横截面的面积，所述砾石组分的大小为所述横截面中所述砾石组分的面积，所述非砾石组分的大小为所述横截面中所述非砾石组分的面积，且所述横截面的面积等于所述砾石组分的面积和所述非砾石组分的面积的和；
所述校正模型对应公式(4)、(5)、(6)：
Sgravel+Snon＝2πR·H (4)
<mrow>
<mfrac>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>N</mi>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
φgravel＝0 (6)
其中，Sgravel表示所述砾石组分的面积，Snon表示所述非砾石组分的面积，2πR·H表示所述横截面的面积，φgravel表示所述砾石组分的孔隙度，φreal表示所述核磁共振校正孔隙度，φNMR表示所述核磁共振测井孔隙度。
4.根据权利要求1-3任一项所述的方法，其特征在于，所述获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，包括：
利用电成像测井仪测量所述目标岩石对应的地层电阻率，依据所述地层电阻率获得所述井壁电成像测井图像；
利用核磁共振测井仪测量所述目标岩石对应的核磁共振信号，依据所述核磁共振信号获得所述核磁共振测井孔隙度。
5.根据权利要求4所述的方法，其特征在于，所述依据所述目标岩石的大小、所述砾石组分的大小和所述非砾石组分的大小以及孔隙度的定义建立所述核磁共振测井孔隙度的校正模型之前，还包括：
获取所述目标岩石对应的岩心分析孔隙度；
依据交会图分析方法对所述岩心分析孔隙度和所述核磁共振测井孔隙度进行误差分析，判断获知所述核磁共振测井孔隙度小于所述岩心分析孔隙度。
6.根据权利要求5所述的方法，其特征在于，所述井壁电成像测井图像包括多个图像点；
所述依据所述井壁电成像测井图像获取所述砾石组分在所述目标岩石中的比例，包括：
依据各个图像点的亮暗程度分别对每个图像点进行标识，且所述图像点的亮暗程度与标识值成正比；
统计所述标识值大于阈值的图像点的个数，所述个数与所述图像点的总个数的比值是所述砾石组分在所述目标岩石中的比例。
7.一种核磁共振测井孔隙度校正装置，其特征在于，包括：
获取模块，用于获取目标岩石对应的井壁电成像测井图像和核磁共振测井孔隙度，所述核磁共振测井孔隙度表示孔隙空间在所述目标岩石中的比例，所述目标岩石包括砾石组分和非砾石组分，所述非砾石组分包括所述孔隙空间；依据所述井壁电成像测井图像获取所述砾石组分在所述目标岩石中的比例；
建模模块，用于依据所述目标岩石的大小、所述砾石组分的大小和所述非砾石组分的大小以及孔隙度的定义建立所述核磁共振测井孔隙度的校正模型；
计算模块，用于依据所述砾石组分在所述目标岩石中的比例和所述校正模型计算核磁共振校正孔隙度，所述核磁共振校正孔隙度是所述核磁共振测井孔隙度的校正值。
8.根据权利要求7所述的核磁共振测井孔隙度校正装置，其特征在于，所述校正模型的检测对象是长度为2πR、宽度为H、厚度为D的目标岩石；
所述目标岩石的大小为所述目标岩石的体积，所述砾石组分的大小为所述砾石组分的体积，所述孔隙空间的大小为所述孔隙空间的体积，且所述目标岩石的体积是所述砾石组分的体积与所述非砾石组分的体积的和；
所述校正模型对应公式(1)、(2)、(3)：
Vgravel+Vnon＝2πR·H·D (1)
<mrow>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>N</mi>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
φgravel＝0 (3)
其中，Vgravel表示所述砾石组分的体积，Vnon表示所述非砾石组分的体积，2πR·H·D表示所述目标岩石的体积，φgravel表示所述砾石组分的孔隙度，φreal表示所述核磁共振校正孔隙度，φNMR表示所述核磁共振测井孔隙度。
9.根据权利要求7所述的核磁共振测井孔隙度校正装置，其特征在于，所述校正模型的检测对象是长度为2πR、宽度为H的横截面；
所述目标岩石的大小为所述横截面的面积，所述砾石组分的大小为所述横截面中所述砾石组分的面积，所述非砾石组分的大小为所述横截面中所述非砾石组分的面积，且所述横截面的面积等于所述砾石组分的面积和所述非砾石组分的面积的和；
所述校正模型对应公式(4)、(5)、(6)：
Sgravel+Snon＝2πR·H (4)
<mrow>
<mfrac>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>g</mi>
<mi>r</mi>
<mi>a</mi>
<mi>v</mi>
<mi>e</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>n</mi>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>N</mi>
<mi>M</mi>
<mi>R</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
φgravel＝0 (6)
其中，Sgravel表示所述砾石组分的面积，Snon表示所述非砾石组分的面积，2πR·H表示所述横截面的面积，φgravel表示所述砾石组分的孔隙度，φreal表示所述核磁共振校正孔隙度，φNMR表示所述核磁共振测井孔隙度。
10.根据权利要求7-9任一项所述的核磁共振测井孔隙度校正装置，其特征在于，所述获取模块具体用于利用电成像测井仪测量所述目标岩石对应的地层电阻率，依据所述地层电阻率获得所述井壁电成像测井图像；
利用核磁共振测井仪测量所述目标岩石对应的核磁共振信号，依据所述核磁共振信号获得所述核磁共振测井孔隙度。
11.根据权利要求10所述的核磁共振测井孔隙度校正装置，其特征在于，所述获取模块还用于获取所述目标岩石对应的岩心分析孔隙度；
所述核磁共振测井孔隙度校正装置还包括：
误差分析模块，用于依据交会图分析方法对所述岩心分析孔隙度和所述核磁共振测井孔隙度进行误差分析，判断获知所述核磁共振测井孔隙度小于所述岩心分析孔隙度。
12.根据权利要求11所述的核磁共振测井孔隙度校正装置，其特征在于，所述井壁电成像测井图像包括多个图像点；
所述获取模块包括：
标识单元，用于依据各个图像点的亮暗程度分别对每个图像点进行标识，且所述图像点的亮暗程度与标识值成正比；
统计单元，用于统计所述标识值大于阈值的图像点的个数，所述个数与所述图像点的总个数的比值是所述砾石组分在所述目标岩石中的比例。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510614235.4A CN105240001B (zh) | 2015-09-23 | 2015-09-23 | 核磁共振测井孔隙度校正方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510614235.4A CN105240001B (zh) | 2015-09-23 | 2015-09-23 | 核磁共振测井孔隙度校正方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105240001A CN105240001A (zh) | 2016-01-13 |
| CN105240001B true CN105240001B (zh) | 2018-03-09 |
Family
ID=55037794
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510614235.4A Active CN105240001B (zh) | 2015-09-23 | 2015-09-23 | 核磁共振测井孔隙度校正方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105240001B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105651805B (zh) * | 2016-03-29 | 2018-01-02 | 西南石油大学 | 利用核磁共振测量岩石孔隙度的方法 |
| CN105954170B (zh) * | 2016-05-04 | 2017-02-22 | 中国石油大学(华东) | 一种考虑背景噪声信号量的核磁孔隙度计算方法 |
| CN106066494B (zh) * | 2016-05-24 | 2018-05-18 | 中国地质大学（北京） | 一种火成岩核磁共振孔隙度校正方法和t2分布校正方法 |
| CN106066492B (zh) * | 2016-06-06 | 2018-05-22 | 中国石油大学(华东) | 一种核磁共振测井孔隙度回波间隔影响的校正方法 |
| CN107829731B (zh) * | 2017-11-06 | 2020-10-09 | 陈国军 | 一种黏土蚀变的火山岩孔隙度校正方法 |
| CN108873082B (zh) * | 2018-05-08 | 2019-11-05 | 中国石油大学(华东) | 一种考虑弛豫组分区间的致密岩心核磁共振孔隙度校正方法 |
| CN111075439B (zh) * | 2018-10-19 | 2023-02-10 | 中国石油天然气股份有限公司 | 确定砾岩储层有效厚度的方法、装置和存储介质 |
| CN110057727B (zh) * | 2019-04-17 | 2021-08-31 | 中国石油大学(华东) | 一种针对砂砾岩储层中岩屑含量统计的方法 |
| CN111220522A (zh) * | 2019-10-28 | 2020-06-02 | 中国石油大学(华东) | 一种高泥质细粒沉积物水合物饱和度的岩心刻度测井计算方法 |
| CN111337408B (zh) * | 2020-03-27 | 2022-02-08 | 西南石油大学 | 一种利用低场核磁共振设备测试岩石裂缝孔隙度的方法 |
| CN114509466B (zh) * | 2022-02-17 | 2024-08-16 | 中国海洋石油集团有限公司 | 稠油含氢指数校正方法及装置 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3780302A (en) * | 1971-09-20 | 1973-12-18 | Texaco Inc | Pulsed neutron logging system |
| US5086643A (en) * | 1990-09-18 | 1992-02-11 | Mobil Oil Corporation | System and method for determining multi-phase relative permeability of a subterranean reservoir |
| MXPA01006730A (es) * | 1998-12-30 | 2003-07-14 | Baker Hughes Inc | Determinacion de la fraccion de arena y saturacion de agua por medio de una herramienta de mapeo de la resistividad de pozos, el registro de induccion transversal y un modelo tensorial de saturacion de agua. |
| CN101929973B (zh) * | 2009-06-22 | 2012-10-17 | 中国石油天然气股份有限公司 | 裂缝储层含油气饱和度定量计算方法 |
| CN102262041B (zh) * | 2011-04-20 | 2013-03-13 | 中国石油天然气股份有限公司 | 一种基于多谱孔隙结构分析的饱和度确定方法 |
| CN103291287A (zh) * | 2013-05-09 | 2013-09-11 | 中国石油天然气股份有限公司 | 一种孔洞型储层有效性级别测定方法 |
-
2015
- 2015-09-23 CN CN201510614235.4A patent/CN105240001B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105240001A (zh) | 2016-01-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105240001B (zh) | 核磁共振测井孔隙度校正方法及装置 | |
| US12025576B2 (en) | Method for assessment of pore-throat size distribution and permeability in porous media | |
| EP2965126B1 (en) | Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data | |
| Schöniger et al. | Finding the right balance between groundwater model complexity and experimental effort via Bayesian model selection | |
| CN107250828A (zh) | 用于改进的井底nmr特性描述的切屑分析 | |
| CN104775811B (zh) | 一种地层各向异性信息的提取及校正方法和系统 | |
| CN105593703A (zh) | 根据基于核心的岩石定型构建并通过岩石置换建模增强的模拟到地震工作流程 | |
| CN105074456A (zh) | 用于改进由岩石样品的材料特性的直接数值模拟以及确定材料特性的不确定性的系统和方法 | |
| Stuart et al. | A two-stage Markov chain Monte Carlo method for seismic inversion and uncertainty quantification | |
| WO2006048596A2 (en) | Method and system for predicting production of a well | |
| CN106370817A (zh) | 一种基于岩心分析和电成像测井的缝洞定量表征方法 | |
| CN101978290A (zh) | 慢度-频率投影显示和动画 | |
| WO2019204004A1 (en) | Determining permeability of porous media based on nuclear magnetic resonance measurement | |
| CN108897066A (zh) | 碳酸盐岩裂缝密度定量预测方法及装置 | |
| US20140270393A1 (en) | Systems and methods for improving direct numerical simulation of material properties from rock samples and determining uncertainty in the material properties | |
| CN106154322A (zh) | 测井曲线校正方法和装置 | |
| NO344398B1 (no) | Fremgangsmåte og anordning for inkorporering av intern gradient og begrenset diffusjon i NMR-inversjon ved brønnlogging | |
| SA520412674B1 (ar) | تحديد تأثيرات أساسها أداة في بيانات تسجيل أداء بالرنين المغناطيسي النووي | |
| Kumar et al. | Measurement of aperture and multiphase flow in fractures with NMR imaging | |
| WO2014142976A1 (en) | Systems and methods for improving direct numerical simulation of material properties from rock samples and determining uncertainty in the material properties | |
| Garcia et al. | Integrated workflow to estimate permeability through quantification of rock fabric using joint interpretation of nuclear magnetic resonance and electric measurements | |
| Chugunova et al. | Explicit fracture network modelling: from multiple point statistics to dynamic simulation | |
| CN105241911B (zh) | 基于lbm模拟低场核磁共振分析流体的方法及装置 | |
| CN109118019A (zh) | 一种煤层气含量预测方法及装置 | |
| CN109031437A (zh) | 随钻核磁共振自旋回波校正方法与装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |