CN106570235B - 一种致密砂岩岩电参数获取方法及装置 - Google Patents

一种致密砂岩岩电参数获取方法及装置 Download PDF

Info

Publication number
CN106570235B
CN106570235B CN201610932948.XA CN201610932948A CN106570235B CN 106570235 B CN106570235 B CN 106570235B CN 201610932948 A CN201610932948 A CN 201610932948A CN 106570235 B CN106570235 B CN 106570235B
Authority
CN
China
Prior art keywords
spectrum
litho
parameter
magnetic resonance
nuclear magnetic
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
Application number
CN201610932948.XA
Other languages
English (en)
Other versions
CN106570235A (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.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
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 China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201610932948.XA priority Critical patent/CN106570235B/zh
Publication of CN106570235A publication Critical patent/CN106570235A/zh
Application granted granted Critical
Publication of CN106570235B publication Critical patent/CN106570235B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Theoretical Computer Science (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本申请实施例提供了一种致密砂岩岩电参数获取方法及装置。其中，所述方法包括：获取目的层的测井资料，以及预设数量岩心的岩电参数；所述岩心位于所述目的层内；确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系；根据所述关系和所述测井资料，确定所述目的层的岩电参数。本申请实施例可以获得高精度的岩电参数。

Description

一种致密砂岩岩电参数获取方法及装置
技术领域
本申请涉及石油勘探油气藏储层测井评价技术领域，尤其是涉及一种致密砂岩岩电参数获取方法及装置。
背景技术
近年来，常规油气资源产量逐年下降，非常规油气能源的勘探开发已引起国内外能源界的高度重视。非常规油气在全球油气产量所占的比重越来越大，并在全球能源结构中扮演着重要的角色。致密油一种重要的非常规石油资源，所谓致密油是指在低孔-低渗甚至超低孔渗砂岩储层中，没有自然产能或者自然产能很低，需要通过压裂或者特殊手段才能形成工业气流的油藏。但是由于致密油储层具有孔隙结构复杂、非均质性强等特点，使得致密油藏评价中存在很大的困难。
饱和度评价是油气储集层定量评价的核心，而岩电参数是计算饱和度过程中一个重要的参数。岩电参数的准确性直接影响油气饱和度的计算结果，而油气饱和度的计算精度一方面又影响着油气资源量的计算精度，从而进一步影响着决策者对于勘探目标的认识和决断；另一方面，在设计开发方案时，油藏中油水分布情况是开发人员必须考虑的因素。
目前测井评价领域中，使用最多的还是以电阻率为基础的饱和度评价，计算过程都是基于固定岩电参数的计算。但是致密砂岩储层等非均质储层逐渐成为油气勘探的重要领域，其储集空间具有如下两个典型特征：(1)油气储集空间致密，通常致密砂岩孔隙度在10％以下(常规砂岩为15～30％)；(2)孔隙类型多样，不仅包括粒间孔隙，还包括溶蚀孔隙等次生孔隙。在复杂孔隙结构的情况下，致密砂岩储层导电规律变得十分复杂，对于如致密砂岩储层，基于固定岩电参数的计算存在一定的不适用性。基于固定的胶结指数和饱和度指数计算的饱和度，计算结果准确率较低。
发明内容
本申请实施例的目的在于提供一种致密砂岩岩电参数获取方法及装置，可以获得高精度的岩电参数。
为达到上述目的，本申请实施例提供了一种致密砂岩岩电参数获取方法，该方法包括：
获取目的层的测井资料，以及预设数量岩心的岩电参数；所述岩心位于所述目的层内；
确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系；
根据所述关系和所述测井资料，确定所述目的层的岩电参数。
本申请实施例还提供了一种致密砂岩岩电参数获取装置，可以获得连续的高精度的岩电参数，该装置包括：
获取模块，用于获取目的层的测井资料，以及预设数量岩心的岩电参数；所述岩心位于所述目的层内；
关系确定模块，用于确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系；
岩电参数确定模块，用于根据所述关系和所述测井资料，确定所述目的层的岩电参数。
由上述本申请实施例所提供的技术方案可知，本申请实施例通过目的层内岩心样本相应数据，建立岩电参数与核磁共振T2谱和孔隙结构参数之间的关系，该关系反映了岩电参数随着核磁共振T2谱和孔隙结构参数变化而变化的关系。核磁共振T2谱和孔隙结构参数反映了储层孔隙度的变化，根据目的层的测井资料中相应数据并结合得到的所述关系，就可以获得随着孔隙度变化而变化的高精度的岩电参数。运用该岩电参数计算的饱和度相对于基于固定岩电参数得到的饱和度而言更加准确。
附图说明
此处所说明的附图用来提供对本申请实施例的进一步理解，构成本申请实施例的一部分，并不构成对本申请实施例的限定。在附图中：
图1为本申请实施例的一种致密砂岩岩电参数获取方法的流程示意图；
图2为本申请实施例的核磁共振T2谱划分的示意图；
图3为本申请实施例的另一种致密砂岩岩电参数获取方法的流程示意图；
图4为本申请实施例的胶结指数m与psm的拟合示意图；
图5为本申请实施例的饱和度指数n与psn的拟合示意图；
图6为本申请实施例的计算饱和度的成果示意图；
图7为本申请实施例的一种致密砂岩岩电参数获取装置示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本申请实施例做进一步详细说明。在此，本申请实施例的示意性实施例及其说明用于解释本申请实施例，但并不作为对本申请实施例的限定。
下面结合附图，对本申请实施例的具体实施方式作进一步的详细说明。
如图1所示为本申请实施例一种致密砂岩岩电参数获取方法的流程示意图，如图1所示，一种致密砂岩岩电参数获取方法可以包括以下步骤。
S101，获取目的层的测井资料，以及预设数量岩心的岩电参数；所述岩心位于所述目的层内。
测井资料是通过测井得到的地球物理相关资料。其中，测井是利用岩层的电化学特性、导电特性、声学特性、放射性等地球物理特性，测量地球物理参数的方法。所述测井资料中包括所述岩心的测井资料。岩心的岩电参数通常是通过岩电实验获得。岩电实验主要是通过钻取长度为2.5～6cm，直径为2.5cm的标准岩心，通过测量岩心的孔隙度、电阻率和饱和度等参数来求取岩电参数。
在本申请的一个具体实施方式中，获取预设数量的岩心的岩电参数通常指获取不少于10个岩心样品，对其进行岩电实验获得岩电参数。获取岩心的数量多少与实验结果精度有关，样本数量过少会导致结果准确性降低。
S102，确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系。
核磁共振测井是能测量到孔隙结构信息的测井方法，借助于核磁测井，测井人员能够获得地下储层孔隙结构参数。核磁共振T2谱以及孔隙结构参数都是能够从测井资料中获得的数据。其中，核磁共振T2谱可以包括：饱和状态的核磁共振T2谱和部分饱和状态的核磁共振T2谱，饱和状态的核磁共振T2谱里面的孔隙结构信息更加准确，部分饱和状态的T2谱进行油气校正，可以得到饱和状态的核磁共振T2谱。在本申请的各个实施例中，所用的核磁共振T2谱均为饱和状态的核磁共振T2谱。已知岩心样本的核磁共振T2谱、孔隙结构参数以及岩电参数，就可以确定岩电参数与核磁共振T2谱和孔隙结构参数之间的关系。
在本申请的各个实施例中，所述孔隙结构参数可以包括：T2谱几何均值和孔隙度。对应S102实施时可以为确定所述岩电参数与岩心测井资料中的核磁共振T2谱、T2谱几何均值和孔隙度之间的关系。
S103，根据所述关系和所述测井资料，确定所述目的层的岩电参数。
根据目的层测井资料中任意一处的核磁共振T2谱和孔隙结构参数，代入S102中得到的关系，就可以得到该出的岩电参数。
图1所示的流程图可知，本申请实施例通过目的层内岩心样本相应数据，确定了岩电参数与核磁共振T2谱和孔隙结构参数之间的关系，该关系反映了岩电参数随着核磁共振T2谱和孔隙结构参数变化而变化的关系，根据目的层的测井资料中任意一处的相应数据并结合得到的所述关系，就可以获得该处高精度的岩电参数。
在本申请的一个实施方式中，所述的核磁共振T2谱具体可以指按照一定预设规则划分的核磁共振T2谱中的三分量：第一预设分量，第二预设分量和第三预设分量。其划分预设规则为：
所述第一预设分量为所述核磁共振T2谱中弛豫时间位于第一预设区间内的孔隙度分量之和；
所述第二预设分量为所述核磁共振T2谱中弛豫时间位于第二预设区间内的孔隙度分量之和；
所述第三预设分量为所述核磁共振T2谱中弛豫时间位于第三预设区间内的孔隙度分量之和。
在本申请的一个具体实施例中，T2谱中的三分量划分依据如下所示。
第一预设分量为所述核磁共振T2谱中弛豫时间小于等于2的孔隙度分量之和；
第二预设分量为所述核磁共振T2谱中弛豫时间大于2且小于20的孔隙度分量之和；
第三预设分量为所述核磁共振T2谱中弛豫时间大于等于20的孔隙度分量之和。
在本实施例中，按照上述划分预设规则对核磁共振T2谱划分的示意图如图2所示。
上述实施例中，按照一定的预设划分规则将核磁共振T2谱划分成三个分量，从而可以将核磁共振T2谱的内容量化表示。核磁共振T2谱可以反映储层或岩心内部的孔隙结构，将该谱划分成三个分量可以准确的将储层或岩心内部的孔隙结构量化表示，从而为后续研究岩电参数随孔隙结构改变而改变的定量关系做准备。
在本申请各个实施方式中所述的岩电参数包括：胶结指数和饱和度指数。对应的，S102中所述确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系，包括：
确定所述胶结指数与所述核磁共振T2谱和孔隙结构参数之间的关系；
以及确定所述饱和度指数与所述核磁共振T2谱和孔隙结构参数之间的关系。
胶结指数和饱和度指数是两个重要的岩电参数，通常胶结指数用m表示，饱和度指数用n表示。在致密砂岩中，m、n的数值会随着孔隙结构的变化而变化。
在本申请的一个具体实施例中，S102具体实施时，可以包括以下步骤。
(1)按照预设规则将岩心的核磁共振T2谱划分成三分量，所述划分规则如下所示：
第一预设分量bin1，对应所述核磁共振T2谱中弛豫时间小于等于2的孔隙度分量之和；
第二预设分量bin2，对应所述核磁共振T2谱中弛豫时间大于2且小于20的孔隙度分量之和；
第三预设分量bin3，对应所述核磁共振T2谱中弛豫时间大于等于20的孔隙度分量之和。
(2)根据所述岩心的bin1、bin2、bin3、T2谱几何均值和孔隙度，得到第一参数psm(pore structure for m)和第二参数psn(pore structure for n)。其中，psm和psn为表征孔隙结构的参数，具体表达式如下所示。
psm＝0.05×bin1+0.65×bin2+0.2×bin3
式中，φ为孔隙度，T2lm为T2谱几何均值。
(3)将岩心的胶结指数与所述psm进行拟合，确定所述胶结指数与所述psm的关系式。
(4)将岩心的饱和度指数与所述psn进行拟合，确定所述饱和度指数与所述psn的关系式。
上述实施例中，首先构建了第一参数psm和第二参数psn，在确定psm和psn与岩电参数m和n的关系式，得到了m和n随psm和psn变化的连续关系式。
在本申请的一个实施例中，S102具体实施时，可以包括以下步骤。
(1)按照预设规则将岩心的核磁共振T2谱划分成三分量，所述划分规则如下所示：
bin1对应所述核磁共振T2谱中弛豫时间区间小于等于2的孔隙度分量；
bin2对应所述核磁共振T2谱中弛豫时间区间大于2且小于20的孔隙度分量；
bin3对应所述核磁共振T2谱中弛豫时间区间大于等于20的孔隙度分量。
(2)根据所述岩心的bin1、bin2、bin3、T2谱几何均值和孔隙度，得到表征孔隙结构的参数psm(pore structure for m)和psn(pore structure for n)。其中，psm和psn具体表达式如下所示。
psm＝0.05×bin1+0.65×bin2+0.2×bin3
式中，φ为孔隙度，T2lm为T2谱几何均值。
(3)将所述岩心的psm和胶结指数，按照第一关系式进行拟合，确定具体表达式；所述第一关系式为：m＝η×psm+λ；
将所述岩心的psn和饱和度指数，按照第二关系式进行拟合，确定具体表达式；所述第二关系式为：
上述实施例中，建立了目的层中岩心测井资料构建的参数psm和psn与岩电参数m和n的关系式，得到了m和n随psm和psn变化的连续关系式。
如图3所示，本申请的一个实施例中另一种致密砂岩岩电参数获取方法可以包括以下步骤。
S301，获取目的层的测井资料，以及预设数量岩心的岩电参数；所述岩心位于所述目的层内。
所述岩电参数包括胶结指数和饱和度指数。所述测井资料中包括所述岩心的测井资料。
S302，按照预设规则将岩心的核磁共振T2谱划分成三分量。
所述划分规则如下所示：
bin1对应所述核磁共振T2谱中弛豫时间区间小于等于2的孔隙度分量；
bin2对应所述核磁共振T2谱中弛豫时间区间大于2且小于20的孔隙度分量；
bin3对应所述核磁共振T2谱中弛豫时间区间大于等于20的孔隙度分量。
S303，根据所述岩心的bin1、bin2、bin3、T2谱几何均值和孔隙度，得到表征孔隙结构的参数psm(pore structure for m)和psn(pore structure for n)。
其中，psm和psn具体表达式如下所示。
psm＝0.05×bin1+0.65×bin2+0.2×bin3
式中，φ为孔隙度，T2lm为T2谱几何均值。
S304，将所述岩心的psm和胶结指数，按照第一关系式进行拟合，确定具体表达式。
所述第一关系式为：m＝η×psm+λ；
S305，将所述岩心的psn和饱和度指数，按照第二关系式进行拟合，确定具体表达式。
所述第二关系式为：
S306，根据所述目的层的测井资料中任意一处的核磁共振T2谱，得到该处的T2谱三分量bin1、bin2和bin3。
S307，根据该处的bin1、bin2、bin3、T2谱几何均值和孔隙度得到该处相应的psm和psn的数值。
S308，将该处的psm和psn的数值代入步骤S304拟合得到的关系式，计算相应的胶结指数和饱和度指数。
上述本申请实施例中，通过目的层内岩心样本相应数据，建立岩电参数中胶结指数和饱和度指数分别与核磁共振T2谱、T2谱几何均值和孔隙度关系式，该关系式反映了胶结指数和饱和度指数与随着核磁共振T2谱、T2谱几何均值和孔隙度变化而变化的定量关系，根据目的层的测井资料中相关数据并结合得到的所述定量关系，就可以获得连续的高精度的岩电参数。
运用上述各个实施例所提供的方法获得胶结指数m和饱和度指数n之后，可用于计算含油饱和度。含油饱和度能够反映岩石孔隙空间所含流体的性质和含量，因此是评价储层含油性的重要参数。目前在测井评价领域中，含油饱和度的计算通常是由电阻率测井资料通过饱和度模型计算得到，最常用的饱和度模型就是阿尔奇公式，为：
式中，Sw表示储层含水饱和度，Rw表示地层水电阻率，单位是Ω·m，φ表示储层孔隙度，Rt表示储层深电阻率，单位为Ω·m，a、b为岩性指数，m为胶结指数，n为饱和度指数，a、b、m、n统称为岩电参数。
阿尔奇公式中岩性指数a、b，通常认为是不变的，去定值带入计算即可，且大部分情况下，a和b都取1。
在本申请的一个实施例中，按照如图3所示的流程图获取目的层的岩电参数，步骤S304得到的m，n与psm和psn之间的关系进行拟合后得到的表达式如下所示：
m＝0.1065×psm+1.4125
m的拟合结果如图4所示，n的拟合结果如图5所示，图中R2表示拟合结果相关系数。由图4可知，拟合结果相关系数达到80.8％，说明两者拟合结果一致性较好。图5中拟合结果的平均相关系数达到89.1％，说明两者之间一致性较好。
在本实施例中，得到拟合表达式后，根据目的层的测井资料得到相应的三分量bin1、bin2和bin3；根据bin1、bin2、bin3、T2谱几何均值和孔隙度得到相应的psm和psn的数值；将所述psm和psn的数值代入拟合得到的关系式，计算相应的胶结指数和饱和度指数。
再根据得到的胶结指数和饱和度指数计算含油饱和度，按照岩性指数a和b都等于1的阿尔奇公式计算含油饱和度，计算结果如图6所示。图6中，第一道为深度道，第二道为自然伽马曲线和声波测井曲线，第三道为电阻率曲线(包含深、浅电阻率)，第四道为核磁共振T2谱以及计算得到的T2谱几何均值，第五道为计算得到的psm和psn，第六道为计算的胶结指数m和饱和度指数n，第七道为孔隙度，第八道为按照阿尔奇公式计算得到的含油饱和度以及取心分析测得的含油饱和度，第九道为解释结论，第十道为试油结论。
由图6中第八道所示可知，计算得到含油饱和度与取心实验测得含油饱和度一致性很高，说明根据本申请实施例所提供的方法得到的岩电参数计算得到的含油饱和度与实际情况吻合度较高，印证了本本申请实施例所获得岩电参数的准确性。
本申请实施例中还提供了一种致密砂岩岩电参数获取装置，如下面的实施例所述。由于该装置解决问题的原理与一种致密砂岩岩电参数获取方法相似，因此该装置的实施可以参见一种致密砂岩岩电参数获取方法的实施，重复之处不再赘述。
如图7所示，本申请实施例所提供的一种致密砂岩岩电参数获取装置可以包括：
获取模块701，用于获取目的层的测井资料，以及预设数量岩心的岩电参数；所述岩心位于所述目的层内。
关系确定模块702，用于确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系。
岩电参数确定模块703，岩电参数确定模块，用于根据所述关系和所述测井资料，确定所述目的层的岩电参数。
由上述装置实施例可知，本申请实施例通过目的层内岩心样本相应数据，建立岩电参数与核磁共振T2谱和孔隙结构参数之间的关系，该关系反映了岩电参数随着核磁共振T2谱和孔隙结构参数变化而变化的关系。核磁共振T2谱和孔隙结构参数反映了储层孔隙度的变化，根据目的层的测井资料中相应数据并结合得到的所述关系，就可以获得随着孔隙度变化而变化的高精度的岩电参数。运用该岩电参数计算的饱和度相对于基于固定岩电参数得到的饱和度而言更加准确。
本领域技术人员还可以了解到本申请实施例列出的各种说明性逻辑块、单元和步骤可以通过硬件、软件或两者的结合来实现。至于是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本申请实施例保护的范围。
本申请实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。
以上所述的具体实施例，对本申请的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本申请实施例的具体实施例而已，并不用于限定本申请的保护范围，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (11)

1.一种致密砂岩岩电参数获取方法，其特征在于，所述方法包括以下步骤：
获取目的层的测井资料，以及预设数量岩心的岩电参数；所述岩心位于所述目的层内；所述岩电参数包括：胶结指数和饱和度指数；
确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系；所述核磁共振T2谱包括：第一预设分量，第二预设分量和第三预设分量；所述孔隙结构参数包括：T2谱几何均值和孔隙度；所述确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系，包括：根据所述岩心的所述第一预设分量，第二预设分量、第三预设分量、T2谱几何均值和孔隙度，得到第一参数psm和第二参数psn：
psm＝0.05×bin1+0.65×bin2+0.2×bin3
式中，bin1表示第一预设分量，bin2表示第二预设分量，bin3表示第三预设分量，φ为孔隙度，T2lm为T2谱几何均值；将所述胶结指数与所述psm进行拟合，确定所述胶结指数与所述psm的关系式；将所述饱和度指数与所述psn进行拟合，确定所述饱和度指数与所述psn的关系式；
根据所述关系和所述测井资料，确定所述目的层的岩电参数。
2.如权利要求1所述的方法，其特征在于，所述第一预设分量为所述核磁共振T2谱中弛豫时间位于第一预设区间内的孔隙度分量之和；
所述第二预设分量为所述核磁共振T2谱中弛豫时间位于第二预设区间内的孔隙度分量之和；
所述第三预设分量为所述核磁共振T2谱中弛豫时间位于第三预设区间内的孔隙度分量之和。
3.如权利要求2所述的方法，其特征在于，所述确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系，包括：
确定所述胶结指数与所述核磁共振T2谱和孔隙结构参数之间的关系；
以及确定所述饱和度指数与所述核磁共振T2谱和孔隙结构参数之间的关系。
4.如权利要求1所述的方法，其特征在于，所述将所述胶结指数与所述psm进行拟合，确定所述胶结指数与所述psm的关系式，具体为：
将所述psm和胶结指数，按照第一关系式进行拟合，确定具体表达式；所述第一关系式为：m＝η×psm+λ，式中m表示胶结指数；
所述将所述饱和度指数与所述psn进行拟合，确定所述饱和度指数与所述psn的关系式，具体为：
将所述psn和饱和度指数，按照第二关系式进行拟合，确定具体表达式；所述第二关系式为：
式中，n表示饱和度指数。
5.如权利要求4所述的方法，其特征在于，所述预设数量至少为10。
6.一种致密砂岩岩电参数获取装置，其特征在于，所述装置包括：
获取模块，用于获取目的层的测井资料，以及预设数量岩心的岩电参数；所述岩心位于所述目的层内；所述岩电参数包括：胶结指数和饱和度指数；
关系确定模块，用于确定所述岩电参数与岩心测井资料中的核磁共振T2谱和孔隙结构参数之间的关系；所述核磁共振T2谱包括：第一预设分量，第二预设分量和第三预设分量；所述孔隙结构参数包括：T2谱几何均值和孔隙度；所述关系确定模块，包括：参数确定子模块，用于根据所述岩心的所述第一预设分量，第二预设分量、第三预设分量、T2谱几何均值和孔隙度，得到第一参数psm和第二参数psn：
psm＝0.05×bin1+0.65×bin2+0.2×bin3
式中，bin1表示第一预设分量，bin2表示第二预设分量，bin3表示第三预设分量，φ为孔隙度，T2lm为T2谱几何均值；第一拟合子模块，用于将所述胶结指数与所述psm进行拟合，确定所述胶结指数与所述psm的关系式；第二拟合子模块，用于将所述饱和度指数与所述psn进行拟合，确定所述饱和度指数与所述psn的关系式；
岩电参数确定模块，用于根据所述关系和所述测井资料，确定所述目的层的岩电参数。
7.如权利要求6所述的装置，其特征在于，所述第一预设分量为所述核磁共振T2谱中弛豫时间位于第一预设区间内的孔隙度分量之和；
所述第二预设分量为所述核磁共振T2谱中弛豫时间位于第二预设区间内的孔隙度分量之和；
所述第三预设分量为所述核磁共振T2谱中弛豫时间位于第三预设区间内的孔隙度分量之和。
8.如权利要求7所述的装置，其特征在于，所述关系确定模块用于确定所述胶结指数与所述核磁共振T2谱和孔隙结构参数之间的关系；以及确定所述饱和度指数与所述核磁共振T2谱和孔隙结构参数之间的关系。
9.如权利要求8所述的装置，其特征在于，所述孔隙结构参数包括：T2谱几何均值和孔隙度；
对应的，所述关系确定模块包括：
参数确定子模块，用于根据所述岩心的所述第一预设分量，第二预设分量、第三预设分量、T2谱几何均值和孔隙度，得到第一参数psm和第二参数psn：
psm＝0.05×bin1+0.65×bin2+0.2×bin3
式中，bin1表示第一预设分量，bin2表示第二预设分量，bin3表示第三预设分量，φ为孔隙度，T2lm为T2谱几何均值；
第一拟合子模块，用于将所述胶结指数与所述psm进行拟合，确定所述胶结指数与所述psm的关系式；
第二拟合子模块，用于将所述饱和度指数与所述psn进行拟合，确定所述饱和度指数与所述psn的关系式。
10.如权利要求9所述的装置，其特征在于，所述第一拟合子模块包括：
第一关系式拟合单元模块，用于将所述psm和胶结指数，按照第一关系式进行拟合，确定具体表达式；所述第一关系式为：m＝η×psm+λ，m表示胶结指数；
所述第二拟合子模块包括：
第二关系式拟合单元模块，用于将所述psn和饱和度指数，按照第二关系式进行拟合，确定具体表达式；所述第二关系式为：
式中，n表示饱和度指数。
11.如权利要求10所述的装置，其特征在于，所述预设数量至少为10。
CN201610932948.XA 2016-10-25 2016-10-25 一种致密砂岩岩电参数获取方法及装置 Active CN106570235B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610932948.XA CN106570235B (zh) 2016-10-25 2016-10-25 一种致密砂岩岩电参数获取方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610932948.XA CN106570235B (zh) 2016-10-25 2016-10-25 一种致密砂岩岩电参数获取方法及装置

Publications (2)

Publication Number Publication Date
CN106570235A CN106570235A (zh) 2017-04-19
CN106570235B true CN106570235B (zh) 2019-11-08

Family

ID=60414285

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610932948.XA Active CN106570235B (zh) 2016-10-25 2016-10-25 一种致密砂岩岩电参数获取方法及装置

Country Status (1)

Country Link
CN (1) CN106570235B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107300519B (zh) * 2017-05-11 2019-09-10 中国石油天然气股份有限公司 致密砂岩导电分析方法及系统
CN109386281B (zh) * 2017-08-02 2021-11-09 中国石油化工股份有限公司 一种获取裂缝性低孔低渗储层测井饱和度的方法
US11092714B2 (en) * 2018-11-21 2021-08-17 Baker Hughes, A Ge Company, Llc Fluid substitution method for T2 distributions of reservoir rocks

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101725344A (zh) * 2008-10-29 2010-06-09 中国石油天然气集团公司 一种确定岩电参数的方法
CN102434152A (zh) * 2011-12-05 2012-05-02 中国石油天然气股份有限公司 一种储层含油饱和度的计算方法
CN104453874A (zh) * 2014-10-23 2015-03-25 中国石油天然气集团公司 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101725344A (zh) * 2008-10-29 2010-06-09 中国石油天然气集团公司 一种确定岩电参数的方法
CN102434152A (zh) * 2011-12-05 2012-05-02 中国石油天然气股份有限公司 一种储层含油饱和度的计算方法
CN104453874A (zh) * 2014-10-23 2015-03-25 中国石油天然气集团公司 一种基于核磁共振的砂砾岩储层含油饱和度的计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
致密砂岩储层饱和度指数n计算方法研究;罗少成等;《西南石油大学学报（自然科学版）》;20140831;第36卷(第4期);第118页至119页 *
致密砂岩脆性指数测井评价方法--以鄂尔多斯盆地陇东地区长7段致密砂岩储集层为例;李华阳等;《新疆石油地质》;20141031;第35卷(第5期);第593页至597页 *

Also Published As

Publication number Publication date
CN106570235A (zh) 2017-04-19

Similar Documents

Publication Publication Date Title
Lai et al. Insight into the pore structure of tight sandstones using NMR and HPMI measurements
CA2903451C (en) Determining continuous capillary pressure curves for subsurface earth formations using saturation and nmr log data
CN104278991B (zh) 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法
CN110824556B (zh) 一种非常规致密砂岩储层的岩石物理模型建立方法及应用
Zielinski et al. Restricted diffusion effects in saturation estimates from 2D diffusion-relaxation NMR maps
CN105891089B (zh) 储层渗透率确定方法及装置
Cheng et al. Nuclear magnetic resonance features of low-permeability reservoirs with complex wettability
CN106570235B (zh) 一种致密砂岩岩电参数获取方法及装置
CN105525917B (zh) 识别储层的方法以及识别储层中流体类型的方法
Liang et al. Tight gas sandstone reservoirs evaluation from nuclear magnetic resonance (NMR) logs: case studies
CN105401937B (zh) 一种基于孔隙结构的饱和度指数预测方法
Masalmeh High oil recoveries from transition zones
Kuang et al. Predicting oil saturation of shale-oil reservoirs using nuclear magnetic resonance logs
CN109667576A (zh) 一种高矿化度成因低阻油层测井识别方法
Feng et al. A new empirical method for constructing capillary pressure curves from conventional logs in low-permeability sandstones
Qobi et al. Permeability determination from Stoneley waves in the Ara group carbonates, Oman
Gomes Characterization and modeling of a thick carbonate transition zone
CN106837318B (zh) 岩层稠油含量获取方法和装置
CN115901563A (zh) 用核磁共振和电阻率联合计算低孔低渗储层渗透率的方法
Lopez et al. Estimates of stress dependent properties in tight reservoirs: Their use with drill cuttings data
CN113294146B (zh) 含沥青质储层的有效孔隙度计算的方法和装置
Toumelin et al. Object-oriented approach for the pore-scale simulation of DC electrical conductivity of two-phase saturated porous media
Shi et al. Quantitative Evaluation of Water Flooding in a Low Resistivity Heavy Oil Reservoir with NMR and Conventional Logs
Li et al. Characterization of Tight Sandstone Reservoir Pore Structure and Validity from Geophysical Logging Data by Using Deep Learning Method
Harbi et al. A New NMR-Based Height Saturation Model of a Low Permeability Carbonate Reservoir

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant