CN103197348B - 利用各层内部样品进行加权编制测井交会图的方法 - Google Patents

利用各层内部样品进行加权编制测井交会图的方法 Download PDF

Info

Publication number
CN103197348B
CN103197348B CN201310098117.3A CN201310098117A CN103197348B CN 103197348 B CN103197348 B CN 103197348B CN 201310098117 A CN201310098117 A CN 201310098117A CN 103197348 B CN103197348 B CN 103197348B
Authority
CN
China
Prior art keywords
reservoir
sample
crossplot
rock core
chart
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.)
Expired - Fee Related
Application number
CN201310098117.3A
Other languages
English (en)
Other versions
CN103197348A (zh
Inventor
张小莉
冯乔
申怡博
程玉群
张涛
阎媛子
黄薇
巨银娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Northwest University Technology Engineering Research Institute Co Ltd
Northwest University
Shandong University of Science and Technology
Original Assignee
Xi'an Northwest University Technology Engineering Research Institute Co Ltd
Northwest University
Shandong University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xi'an Northwest University Technology Engineering Research Institute Co Ltd, Northwest University, Shandong University of Science and Technology filed Critical Xi'an Northwest University Technology Engineering Research Institute Co Ltd
Priority to CN201310098117.3A priority Critical patent/CN103197348B/zh
Publication of CN103197348A publication Critical patent/CN103197348A/zh
Application granted granted Critical
Publication of CN103197348B publication Critical patent/CN103197348B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

一种利用各层内部样品进行加权编制测井交会图的方法，充分考虑到油区目的层的每个储层的岩心分析样品数量的差异，以及测井资料分辨率与岩心分析采样分辨率的差异，对每个储层中的所有岩心分析样品的测试数值进行加权处理。通过上述处理，充分体现了各层厚度及其内部样品数的差异，符合统计学中对于样本空间的要求，明显提高测井交会图的绘制精度。

Description

利用各层内部样品进行加权编制测井交会图的方法
技术领域
本发明属于勘探地球物理领域，具体属于测井方法技术分支，尤其涉及到一种利用各层内部样品进行加权编制测井交会图的方法。
背景技术
随着油气田勘探开发程度的不断提高，所面临的地质目标复杂程度及测井解释难度也越来越大，因而有必要详细分析各类地质特征及其发生的频次。这样不仅需要详细分析具有各种属性的油藏地质特征，还需要有配套的技术方法对其属性进行提取和反映。测井交会图是一项常用来对油藏岩性、物性、含油性等特征表征的一种关键技术。如下所示通过编制测井交会图中最典型和最具代表性的声波时差—孔隙度测井交会图来反映出已有的制作测井交会图方法过程中存在的问题：
    设定根据测井资料、录井资料等综合分析，一个油区内井中目的层段具有岩心物性分析测试数据的储层共n个，在这n个储层中对应有m个岩心物性分析测试样品，m≥n，那么在n个储层中样品的数量分别为m1，m2，…，mn，则：m1+m2+…+mn=m，由于常规测井采样间隔是8点/m，即0.125m测量间隔，常规测井系列资料可以分辨岩层的最小厚度一般是0.5m，且存在界面效应，而实际岩心分析小柱样直径是2.5cm，长是3~5cm，可以看出两种资料的采样精度是不匹配的，界面附近测井资料无法逐一反映每个岩心样品分析值所对应的测井参数值，即无法严格按照岩心物性分析资料采用单点资料深度逐一对应的声波测井资料来编绘交会图，如何把岩心物性分析资料与测井资料有机对应起来，达到岩心刻度测井的目的是个问题。一般根据测井资料特征，目前采用层点分析方法，分别把n个储层内部的样品分析测试数值分别进行算术加权平均，如把第1、2、…，n层内的m1，m2，…，mn个样品的孔隙度值分别计算其算术平均值，把算术平均值分别记为p1，p2，…，pn，然后结合这n个平均值对应的各层声波时差值Δt1，Δt2，…，Δtn，编制成声波时差—孔隙度测井交会图，但是这样的声波时差—孔隙度测井交会图中所用样本数为n，在其编制交会图过程中，没有充分考虑到各层厚度差异以及各层样品数的差异所携带的样本空间信息，没有充分反映出样品的典型性和代表性。
发明内容
本发明提供一种利用各层内部样品进行加权编制测井交会图的方法，充分考虑到油区目的层的每个储层的岩心分析样品数量的差异，以及测井资料分辨率与岩心分析采样分辨率的差异，对每个储层中的所有岩心分析样品的测试数值进行加权处理。通过上述处理，充分体现了各层厚度及其内部样品数的差异，符合统计学中对于样本空间的要求，明显提高测井交会图的编制精度。
为实现上述目的，本发明的技术方案为：
一种利用各层内部样品进行加权编制测井交会图的方法，步骤如下：
步骤1：首先设定需要测试油区的目的层中，用取心器在需要测试油区的目的层预设的取样深度上取出m块岩心，将这m块岩心利用针对岩心的参数分析仪测试出其对应的参数，并结合对应岩心的测井资料和录井资料，在作为分析样本的岩心的取样深度及其以该取样深度为基准来设定的浮动范围条件下，按照岩层划分方法对井中目的层段进行岩层划分成一个以上的储层，并达到使被划分的每个储层内部的测井属性参数能够由一个测井参数值来反映，每个储层内部所述的一个测井参数值按照其对应的测井曲线的特征来确定，设定岩层划分成n个储层，n个储层的所述的一个测井参数值分别为Δt1，Δt2，…，Δtn，其中m和n均为大于等于1的整数，且m≥n；
步骤2：对于每个储层，分别统计出n个储层其各自对应的作为分析样本的岩心数量分别为m1，m2，…，mn，并对每个储层中的作为分析样本的岩心用对应的岩心分析测试工具进行分析测试，对每个储层中的作为分析样本的岩心分析测试后得到的样品分析测试数值按照预设的权值进行算术加权平均或几何加权平均，算术加权平均或几何加权平均后就得到了每个储层各自对应的样品分析测试数值的平均值，n个储层的样品分析测试数值的平均值分别设定为p1，p2，…，pn，其中m1，m2，…，mn均为大于等于1的整数，并且m1，m2，…，mn的总和为m；
步骤3：根据每个储层各自对应的作为分析样本的岩心数量、样品分析测试数值的平均值以及对应的测井参数值，确定每个储层参与编制交会图的样本，设定用各个储层的测井参数值做横坐标、各个储层的样品分析测试数值的平均值做纵坐标来用交会图绘制工具进行编绘测井交会图，由此得到第1储层参编资料为m1个坐标（Δt1，p1），第2储层参编资料为m2个坐标（Δt2，p2），…，第n储层参编资料为mn个坐标（Δtn，pn），n个储层中采用样本数为m1+m2+，…+mn=m个，进而利用交会图绘制工具在平面图纸上编制测井交会图。
所述的参数分析仪为物性分析仪。
所述的岩层划分方法为有序量最优分割的方法。
所述的测井参数值为声波时差测井参数值。
所述的样品分析测试数值为样品孔隙度分析测试数值。
本发明的技术方法优点为：
（1）充分体现了不同储层之间，由于各储层样品的差异所体现的样本代表性和权值差异；
（2）通过对每个储层中的所有岩心分析样品数值进行加权，这较常规层点编绘法拓展了样本空间，加大了具有典型性和代表性样品的权重，最终提高交会图精度。
总之，应用该加权法编绘交会图，发展了交会图编制方法，更有利于揭示岩—电之间关系，并有利于利用所建立的关系，实现岩—电参数转换，达到真实反映油藏参数及其特征的目的。
附图说明
图1为实施例1中利用各层内部样品进行加权编制测井交会图的方法编制的测井交会图及其回归相关参数R。
图2为实施例1中利用常规的层点法所编制的测井交会图及其回归相关参数R。
图3为实施例2中利用各层内部样品进行加权编制测井交会图的方法编制的测井交会图及其回归相关参数R。
图4为实施例2中利用常规的层点法所编制的测井交会图及其回归相关参数R。
图5为实施例3中利用各层内部样品进行加权编制测井交会图的方法编制的测井交会图及其回归相关参数R。
图6为实施例3中利用常规的层点法所编制的测井交会图及其回归相关参数R。
具体实施方式
下面通过具体实施实例对本发明做进一步说明：
实施例1：
利用各层内部样品进行加权编制测井交会图的方法，步骤如下：
步骤1：首先设定需要测试油区的目的层的5个井号分别为J33-22、J42-21、S231、W19-8和Z43的井中，用取心器在需要测试油区的目的层的5个井号分别为J33-22、J42-21、S231、W19-8和Z43的井中预设的取样深度上取出191块岩心，将这191块岩心利用针对岩心的物性分析仪测试出其对应的参数，并结合对应岩心的测井资料和录井资料，在作为分析样本的岩心的取样深度及其以该取样深度为基准来设定的浮动范围条件下，按照有序量最优分割法对井中目的层段进行岩层划分成层号分别为1、2、3…48的48个储层，并达到使被划分的每个储层内部的测井属性参数能够由一个声波时差测井参数值来反映，每个储层内部所述的一个声波时差测井参数值按照其对应的测井曲线的特征来确定并如表1统计所示；
步骤2：对于每个储层，分别统计出48个储层其各自对应的作为分析样本的岩心数量分别为如表1所示的6、4、…、1，并对每个储层中的作为分析样本的岩心用对应的岩心分析测试工具进行分析测试，对每个储层中的作为分析样本的岩心分析测试后得到的样品孔隙度分析测试数值按照预设的权值进行算术加权平均，算术加权平均后就得到了每个储层各自对应的样品孔隙度分析测试数值的平均值，48个储层的样品孔隙度分析测试数值的平均值分别为如表1所示的5.8、7.8、…、8.1；
步骤3：根据每个储层各自对应的作为分析样本的岩心数量、样品孔隙度分析测试数值的平均值以及对应的声波时差测井参数值，确定每个储层参与编制交会图的样本，设定用各个储层的测井参数值做横坐标、各个储层的样品孔隙度分析测试数值的平均值做纵坐标来用交会图绘制工具进行编绘测井交会图，由此得到第1储层参编资料为6个坐标（202.8，5.8），第2储层参编资料为4个坐标（221.8，7.8），…，第48储层参编资料为1个坐标（214.1，8.1），48个储层中采用样本数为191个，进而利用交会图绘制工具在平面图纸上编制测井交会图如图1所示。
为对比起见，图2为对应的按照常规的层点法所编制的交会图及其回归相关参数R。对比图1和图2可以看出，经过各层内部样品进行加权编制测井交会图，其回归关系的相关性有明显提高。
   实施例2：
利用各层内部样品进行加权编制测井交会图的方法，步骤如下：
步骤1：首先设定需要测试油区的目的层，用取心器在需要测试油区的目的层中的预设的取样深度上取出95块岩心，将这95块岩心利用针对岩心的物性分析仪测试出其对应的参数，并结合对应岩心的测井资料和录井资料，在作为分析样本的岩心的取样深度及其以该取样深度为基准来设定的浮动范围条件下，按照有序量最优分割的对井中目的层段进行岩层划分成层号分别为1、2、3…33的33个储层，并达到使被划分的每个储层内部的测井属性参数能够由一个声波时差测井参数值来反映，每个储层内部所述的一个声波时差测井参数值按照其对应的测井曲线的特征来确定并如表2统计所示；
步骤2：对于每个储层，分别统计出33个储层其各自对应的作为分析样本的岩心数量分别为如表2所示的2、4、…、3，并对每个储层中的作为分析样本的岩心用对应的岩心分析测试工具进行分析测试，对每个储层中的作为分析样本的岩心分析测试后得到的样品孔隙度分析测试数值按照预设的权值进行算术加权平均，算术加权平均后就得到了每个储层各自对应的样品孔隙度分析测试数值的平均值，33个储层的样品孔隙度分析测试数值的平均值分别为如表2所示的6.5、7.1、…、4.1；
步骤3：根据每个储层各自对应的作为分析样本的岩心数量、样品孔隙度分析测试数值的平均值以及对应的声波时差测井参数值，确定每个储层参与编制交会图的样本，设定用各个储层的测井参数值做横坐标、各个储层的样品孔隙度分析测试数值的平均值做纵坐标来用交会图绘制工具进行编绘测井交会图，由此得到第1储层参编资料为2个坐标（（228.6，6.5），第2储层参编资料为4个坐标（227.7，7.1），…，第33储层参编资料为3个坐标（206.0，4.1），33个储层中采用样本数为95个，进而利用交会图绘制工具在平面图纸上编制测井交会图如图3所示。
为对比起见，图4为对应的按照常规的层点法所编制的交会图及其回归相关参数R。对比图3和图4可以看出，经过各层内部样品进行加权编制测井交会图，其回归关系的相关性有明显提高。
实施例3：
利用各层内部样品进行加权编制测井交会图的方法，步骤如下：
步骤1：首先设定需要测试油区的目的层的18个井如表3所示的井号对应的井中，用取心器在需要测试油区的目的层的18个井如表3所示的井号对应的井中的预设的取样深度上取出251块岩心，将这251块岩心利用针对岩心的物性分析仪测试出其对应的参数，并结合对应岩心的测井资料和录井资料，在作为分析样本的岩心的取样深度及其以该取样深度为基准来设定的浮动范围条件下，按照有序量最优分割的对井中目的层段进行岩层划分成层号分别为1、2、3…152的152个储层，并达到使被划分的每个储层内部的测井属性参数能够由一个声波时差测井参数值来反映，每个储层内部所述的一个声波时差测井参数值按照其对应的测井曲线的特征来确定并如表3统计所示；
步骤2：对于每个储层，分别统计出152个储层其各自对应的作为分析样本的岩心数量分别为如表2所示的2、1、…、1，并对每个储层中的作为分析样本的岩心用对应的岩心分析测试工具进行分析测试，对每个储层中的作为分析样本的岩心分析测试后得到的样品孔隙度分析测试数值按照预设的权值进行算术加权平均，算术加权平均后就得到了每个储层各自对应的样品孔隙度分析测试数值的平均值，152个储层的样品孔隙度分析测试数值的平均值分别为如表3所示的11.0、10.4、…、9.6；
步骤3：根据每个储层各自对应的作为分析样本的岩心数量、样品孔隙度分析测试数值的平均值以及对应的声波时差测井参数值，确定每个储层参与编制交会图的样本，设定用各个储层的测井参数值做横坐标、各个储层的样品孔隙度分析测试数值的平均值做纵坐标来用交会图绘制工具进行编绘测井交会图，由此得到第1储层参编资料为2个坐为对比起见，图4为对应的按照常规的层点法所编制的交会图及其回归相关参数R。对比图3和图4可以看出，经过各层内部样品进行加权编制测井交会图，其回归关系的相关性有明显提高。
步骤3：根据每个储层各自对应的作为分析样本的岩心数量、样品孔隙度分析测试数值的平均值以及对应的声波时差测井参数值，确定每个储层参与编制交会图的样本，设定用各个储层的测井参数值做横坐标、各个储层的样品孔隙度分析测试数值的平均值做纵坐标来用交会图绘制工具进行编绘测井交会图，由此得到第1储层参编资料为2个坐标（229.2，11.0），第2储层参编资料为1个坐标（227.6，10.4），…，第152储层参编资料为1个坐标（223.9，9.6），152个储层中采用样本数为251个，进而利用交会图绘制工具在平面图纸上编制测井交会图如图5所示。
为对比起见，图6为对应的按照常规的层点法所编制的交会图及其回归相关参数R。对比图5和图6可以看出，经过各层内部样品进行加权编制测井交会图，其回归关系的相关性有明显提高。
从上述实施例可知，利用各层内部样品进行加权编制测井交会图的方法得到的回归关系的相关系数均在0.89以上。

Claims (2)

1.一种利用各层内部样品进行加权编制测井交会图的方法，其特征在于，步骤如下：
步骤1：首先设定需要测试油区的目的层中，用取心器在需要测试油区的目的层预设的取样深度上取出m块岩心，将这m块岩心利用针对岩心的参数分析仪测试出其对应的参数，并结合对应岩心的测井资料和录井资料，在作为分析样本的岩心的取样深度及其以该取样深度为基准来设定的浮动范围条件下，按照有序量最优分割这样一种岩层划分方法对井中目的层段进行岩层划分成一个以上的储层，并达到使被划分的每个储层内部的测井属性参数能够由一个声波时差测井参数值来反映，每个储层内部所述的一个声波时差测井参数值按照其对应的测井曲线的特征来确定，设定岩层划分成n个储层，n个储层的所述的一个声波时差测井参数值分别为Δt1，Δt2，…，Δtn，其中m和n均为大于等于1的整数，且m≥n；
步骤2：对于每个储层，分别统计出n个储层其各自对应的作为分析样本的岩心数量分别为m1，m2，…，mn，并对每个储层中的作为分析样本的岩心用对应的岩心分析测试工具进行分析测试，对每个储层中的作为分析样本的岩心分析测试后得到的样品的孔隙度分析测试数值按照预设的权值进行算术加权平均或几何加权平均，算术加权平均或几何加权平均后就得到了每个储层各自对应的样品的孔隙度分析测试数值的平均值，n个储层的样品的孔隙度分析测试数值的平均值分别设定为p1，p2，…，pn，其中m1，m2，…，mn均为大于等于1的整数，并且m1，m2，…，mn的总和为m；
步骤3：根据每个储层各自对应的作为分析样本的岩心数量、样品的孔隙度分析测试数值的平均值以及对应的声波时差测井参数值，确定每个储层参与编制交会图的样本，设定用各个储层的声波时差测井参数值做横坐标、各个储层的样品的孔隙度分析测试数值的平均值做纵坐标来用交会图绘制工具进行编绘测井交会图，由此得到第1储层参编资料为m1个坐标（Δt1，p1），第2储层参编资料为m2个坐标（Δt2，p2），…，第n储层参编资料为mn个坐标（Δtn，pn），n个储层中采用样本数为m1+m2+，…+mn=m个，进而利用交会图绘制工具在平面图纸上编制测井交会图。
2.根据权利要求1所述的利用各层内部样品进行加权编制测井交会图的方法，其特征在于，所述的参数分析仪为物性分析仪。
CN201310098117.3A 2013-03-26 2013-03-26 利用各层内部样品进行加权编制测井交会图的方法 Expired - Fee Related CN103197348B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310098117.3A CN103197348B (zh) 2013-03-26 2013-03-26 利用各层内部样品进行加权编制测井交会图的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310098117.3A CN103197348B (zh) 2013-03-26 2013-03-26 利用各层内部样品进行加权编制测井交会图的方法

Publications (2)

Publication Number Publication Date
CN103197348A CN103197348A (zh) 2013-07-10
CN103197348B true CN103197348B (zh) 2015-04-15

Family

ID=48720024

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310098117.3A Expired - Fee Related CN103197348B (zh) 2013-03-26 2013-03-26 利用各层内部样品进行加权编制测井交会图的方法

Country Status (1)

Country Link
CN (1) CN103197348B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104863572A (zh) * 2014-02-20 2015-08-26 中国石油化工股份有限公司 一种分析储层孔隙度的方法
CN105781536B (zh) * 2014-12-16 2019-01-18 中国石油天然气股份有限公司 一种石油储层孔隙度测算方法
CN105572747B (zh) * 2015-12-14 2018-08-14 中国石油大学(华东) 一种识别高盐富钾区注水开发油藏碎屑岩岩性的方法
CN107327294A (zh) * 2017-08-28 2017-11-07 西北大学 一种基于致密油储层变骨架参数条件下的孔隙度解释方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001033255A1 (en) * 1999-11-02 2001-05-10 Geco As Method and apparatus for generating a cross plot in attribute space from a plurality of attribute data sets and generating a class data set from the cross plot
CN102418518A (zh) * 2011-04-12 2012-04-18 北京师范大学 神经网络模拟交会图识别油层水淹级别的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130024122A1 (en) * 2011-07-18 2013-01-24 Schlumberger Technology Corporation Formation fluid detection
US9846256B2 (en) * 2011-08-09 2017-12-19 Schlumberger Technology Corporation Interactive display of results obtained from the inversion of logging data

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001033255A1 (en) * 1999-11-02 2001-05-10 Geco As Method and apparatus for generating a cross plot in attribute space from a plurality of attribute data sets and generating a class data set from the cross plot
CN102418518A (zh) * 2011-04-12 2012-04-18 北京师范大学 神经网络模拟交会图识别油层水淹级别的方法

Also Published As

Publication number Publication date
CN103197348A (zh) 2013-07-10

Similar Documents

Publication Publication Date Title
CN108931814B (zh) 一种基于多属性融合的基岩裂缝预测的方法
CN113759424B (zh) 基于频谱分解和机器学习的岩溶储层充填分析方法和系统
US8694261B1 (en) 3D-well log invention
CN102175832B (zh) 一种确定典型储层最佳饱和度计算模型的方法
CN113759425A (zh) 井震联合评价深层古岩溶储层充填特征的方法与系统
CN107748399B (zh) 利用重力界面反演识别山前带深部构造层方法
CN103777243A (zh) 砂泥岩薄互层储层厚度预测方法
CN103675907A (zh) 一种基于岩相约束的avo反演油气检测方法
CN105277982A (zh) 一种泥页岩总有机碳含量地震预测方法
Kianoush et al. Compressional and shear interval velocity modeling to determine formation pressures in an oilfield of SW Iran
CN103197348B (zh) 利用各层内部样品进行加权编制测井交会图的方法
CN103135131A (zh) 一种针对裂缝性储层预测的解释装置
CN106443770A (zh) 一种页岩气地质甜点的预测方法
CN105114067A (zh) 一种岩性测井相的方法
CN103628866A (zh) 获取地层的地层压力系数的方法
Cai et al. Intelligent calculation method of relative sonic attenuation and its application to fracture evaluation in tight sandstone reservoir
RU2011148308A (ru) Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления
CN112528106A (zh) 火山岩岩性识别方法
CN103529473A (zh) 一种地震反演过程中差异压实的补偿方法
CN103278852B (zh) 利用地震数据体波形结构特征模型预测油气的方法
CN104570104B (zh) 一种基于两步法avf的纵横波地震品质因子提取方法
CN103777244A (zh) 一种地震裂缝属性体的定量分析方法
Zhang et al. Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China
RU2255358C1 (ru) Способ геофизической разведки для выявления нефтегазопродуктивных типов геологического разреза в трехмерном межскважинном пространстве
Jiang et al. Quantitative identification of pseudofaults underlying igneous rocks through wave-equation-based velocity modeling: A case study of the Shuntuoguole uplift, Tarim Basin, China

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150415

Termination date: 20200326