CN105223116B - 一种基于核磁共振谱系数法计算束缚水饱和度的方法 - Google Patents
一种基于核磁共振谱系数法计算束缚水饱和度的方法 Download PDFInfo
- Publication number
- CN105223116B CN105223116B CN201510541975.XA CN201510541975A CN105223116B CN 105223116 B CN105223116 B CN 105223116B CN 201510541975 A CN201510541975 A CN 201510541975A CN 105223116 B CN105223116 B CN 105223116B
- Authority
- CN
- China
- Prior art keywords
- water saturation
- irreducible water
- weight
- rock sample
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于核磁共振谱系数法计算束缚水饱和度的方法，属于储层评价领域。所述方法包括：首先对岩心进行核磁共振标准T2测试实验，获取岩心核磁T2分布谱，利用称重法测量计算岩样的实际束缚水饱和度STURE；利用核磁T2谱和岩样的实际束缚水饱和度STURE，采用最小二乘法确定不同孔隙分量下的待定系数，建立计算束缚水饱和度的新模型，最后得到岩样计算的束缚水饱和度。本发明实现了利用核磁共振谱系数法确定储层束缚水饱和度，且通过对比采用新模型计算的结果，能够更准确、更可靠的获取储层束缚水饱和度，为石油勘探和油藏开发评价提供可靠的储层参数。
Description
技术领域
本发明涉及储层评价领域，储层束缚水饱和度评价领域，特别涉及一种基于核磁共振谱系数法计算束缚水饱和度的方法。
背景技术
利用核磁共振测井计算储层束缚水饱和度尤为重要，它直接或间接影响到有效孔隙度、渗透率、自由流体饱和度等储层参数的计算精度，目前利用核磁共振测井计算储层束缚水饱和度的方法主要有两种：T2截止值法与谱系数法。
T2截止值法是目前最常用的方法，计算模型如下：
式中：SBVI为束缚水饱和度，T2min为T2分布最小值，T2max为T2分布最大值，φ为T2对应的幅度值，T2cutoff为束缚流体的T2截止值，具有地区经验性，需要由大量岩心实验确定，也是计算精度的关键所在，在砂岩储层通常取33ms，在碳酸岩储层通常取92ms。
谱系数法的计算模型如下：
式中：SSSBVI为计算的束缚水饱和度，Wi为计算的第i个权系数分量，为T2分布第i个孔隙分量，m、b为待定系数，Coates等人通过大量实验在砂岩储层m为0.0618，b为1。
在实现本发明的过程中，发明人发现现有技术至少存在以下问题：
T2截止值法模型简单，截止值地区性依赖强取值差异大，在特定区域计算结果精度难以满足要求，尤其在稠油储层计算结果偏大；谱系数法模型合理但是在小孔占优的储层计算结果往往偏小，权系数模型有待改进。
发明内容
本发明的目的是在谱系数法基础上通过对权系数的改进，从而得到新的计算模型，使其计算束缚水饱和度更加接近实际情况，满足对储层评价的精度要求。
为达到上述目的，本发明采用的技术方案为：
一种基于核磁共振谱系数法计算束缚水饱和度的方法，包括以下步骤：
步骤(1)，对岩样岩心进行核磁共振标准T2测试实验，获取核磁T2谱分布；
步骤(2)，利用称重法测量计算岩样的实际束缚水饱和度STURE；
步骤(3)，利用核磁T2谱和岩样的实际束缚水饱和度STURE，采用最小二乘法确定待定系数a、b；
步骤(4)，根据步骤(3)确定的待定系数a、b，建立新模型：
其中，SSSBVI为计算的束缚水饱和度；
Wi为权系数分量；
φi为T2分布孔隙分量；
T2,i为核磁测量的横向弛豫时间；
a、b为待定系数；
利用该模型得到计算束缚水饱和度SSSBVI。
在进行步骤(1)之时，按照《岩样核磁共振参数实验室测量规范SY/T 6490-2007》标准规定的流程对岩样进行核磁共振测量获取岩样T2分布谱。
在进行步骤(2)之前，需要测量岩样的干重、湿重、离心后的重量。
具体的，利用测量岩样干重、湿重、离心重的数据采用称重法计算岩样束缚水饱和度作为岩样实际的束缚水饱和度STURE值，通过以下公式计算得到，
其中，GL为岩样离心后重量，即离心重；GG为岩样烘干后的重量，即干重；GS为岩样饱和水状态下的重量，即湿重。
具体的，利用步骤(1)得到的岩样T2谱与步骤(2)称重法得到的岩样的实际束缚水饱和度STURE，采用最小二乘法利用excel中规划求解来确定新模型中的待定系数a、b值。
具体的，所述步骤(3)中Wi为权系数分量的求取，所述步骤(4)中权系数分量Wi的求取，每一个权系数分量与其所对应的T2时间值为指数关系，而且随着T2时间的减小权系数Wi逐渐趋于1，具体的是将步骤(3)中确定的系数a、b带入到公式求取。
具体的，所述步骤(4)中计算的束缚水饱和度SSSBVI是根据确定的权系数得到计算的束缚水饱和度SSSBVI。
与实际束缚水饱和度相比，所述束缚水饱和度SSSBVI绝对误差在4个单位以内，该单位以百分比计；相对误差在10％以内。
相对于现有技术，本发明的有益效果为：
本发明提供了一种确定储层束缚水饱和度的新方法，通过对现有谱系数法计算模型的改进，使得权系数更加容易求取；通过实际研究发现a、b值与T2截止值具有相关性，认为改进模型是原谱系数法和T2截止值法两者的综合考虑。使得利用核磁共振计算储层束缚水饱和度更加准确、可加可靠，为利用核磁共振测井计算其他相关的储层参数提供有力保证。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的确定储层束缚水饱和度的方法流程图；
图2是本发明实施例提供的新模型在不同a值下权系数Wi变化曲线；
图3是本发明实施例提供的原模型计算的束缚水饱和度与称重法测量束缚水饱和度对比图。
图4是本发明实施例提供的新模型计算的束缚水饱和度与称重法测量束缚水饱和度对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
参见图1，本发明实施例提供了一种确定储层束缚水饱和度的方法，按照如下步骤进行操作：
步骤一，本实施例中选取华北油田某区域的30块岩样，按照《岩样核磁共振参数实验室测量规范SY/T 6490-2007》标准规定的流程对这些岩样进行核磁共振测量获取岩样T2分布谱。
步骤二，通过测量每一块岩样干重、湿重、离心重的数据采用称重法计算岩样束缚水饱和度作为岩样真实的束缚水饱和度值，见表1，计算公式如下：
其中，GL为岩样离心后重量，即离心重；GG为岩样烘干后的重量，即干重；GS为岩样饱和水状态下的重量，即湿重。
步骤三，结合步骤一中获得的核磁共振T2谱与步骤二获得的束缚水饱和度，采用最小二乘法求取a值为0.032，b值为1.05。
步骤四，利用步骤三得到的待定系数a、b，最终得到计算的束缚水饱和度模型：
其中，SSSBVI为计算的束缚水饱和度；
Wi为权系数分量；
为T2分布孔隙分量；
T2,i为核磁测量的横向弛豫时间。
图2是本发明实施例提供的新模型在不同a值下权系数Wi变化曲线，随着T2值的减小，权系数逐渐趋近1，其物理意义为随着孔隙半径的减小，其束缚水所占空间比例逐渐增大；而且随着a值的增大，束缚水所占空间比例增大的速率更快。
由上述新模型计算每一块岩样的计算的束缚水饱和度SSSBVI，计算结果见表1。
表1
从该实施例中可以看出该批岩样束缚水饱和度的计算采用新模型与原模型相比绝对误差平均降低1.56，相对误差评价降低4.92％，与实际束缚水饱和度相比绝对误差是3.04(％)，相对误差为8.51％。其对比见图3、图4所示。图中充分说明了模型的改进有助于束缚水饱和度计算精度的提高。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
以上对本发明实施例所提供的技术方案进行了详细介绍，本文中应用了具体个例对本发明实施例的原理以及实施方式进行了阐述，以上实施例的说明只适用于帮助理解本发明实施例的原理；同时，对于本领域的一般技术人员，依据本发明实施例，在具体实施方式以及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种基于核磁共振谱系数法计算束缚水饱和度的方法，其特征在于，所述方法按照以下步骤进行操作：
步骤(1)，对岩样岩心进行核磁共振标准T2测试实验，获取核磁T2谱分布；
步骤(2)，利用称重法测量计算岩样的实际束缚水饱和度STURE；
步骤(3)，利用核磁T2谱和岩样的实际束缚水饱和度STURE，采用最小二乘法确定待定系数a、b；
步骤(4)，根据步骤(3)确定的待定系数a、b，建立束缚水饱和度计算模型：
<mrow>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>a</mi>
<mo>&CenterDot;</mo>
<msup>
<msub>
<mi>T</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mi>b</mi>
</msup>
</mrow>
</msup>
</mrow>
其中，SSSBVI为计算的束缚水饱和度；
Wi为权系数分量；
为T2分布孔隙分量；
T2,i为核磁测量的横向弛豫时间；
a、b为待定系数；
利用该模型得到计算束缚水饱和度SSSBVI。
2.如权利要求1所述的一种基于核磁共振谱系数法计算束缚水饱和度的方法，其特征在于，所述步骤(1)中岩心核磁共振T2测试实验按照《岩样核磁共振参数实验室测量规范SY/T 6490-2007》标准规定的流程进行。
3.如权利要求1所述的一种基于核磁共振谱系数法计算束缚水饱和度的方法，其特征在于，在进行步骤(2)之前，需要测量岩样的干重、湿重、离心后的重量。
4.如权利要求3所述的一种基于核磁共振谱系数法计算束缚水饱和度的方法，其特征在于，利用测量岩样干重、湿重、离心重的数据，利用称重法测量计算岩样的实际束缚水饱和度，通过以下公式计算：
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>T</mi>
<mi>U</mi>
<mi>R</mi>
<mi>E</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>G</mi>
<mi>L</mi>
</msub>
<mo>-</mo>
<msub>
<mi>G</mi>
<mi>G</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>G</mi>
<mi>S</mi>
</msub>
<mo>-</mo>
<msub>
<mi>G</mi>
<mi>G</mi>
</msub>
</mrow>
</mfrac>
<mo>*</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
其中，STURE为岩样的实际束缚水饱和度；
GL为岩样离心后重，即离心重；
GG为岩样烘干后的重量，即干重；
GS为岩样饱和水状态下的重量，即湿重。
5.如权利要求1所述的一种基于核磁共振谱系数法计算束缚水饱和度的方法，其特征在于，所述步骤(3)中待定系数a、b是利用步骤(1)得到的岩样T2谱与步骤(2)称重法得到的岩样的实际束缚水饱和度STURE，采用最小二乘法原理利用excel中规划求解来确定模型中的待定系数a、b值。
6.如权利要求1所述的一种基于核磁共振谱系数法计算束缚水饱和度的方法，其特征在于，所述步骤(4)中权系数分量Wi的求取，每一个权系数分量与其所对应的T2时间值为指数关系，而且随着T2时间的减小权系数分量Wi逐渐趋于1。
7.如权利要求1所述的一种基于核磁共振谱系数法计算束缚水饱和度的方法，其特征在于，所述步骤(4)中计算的束缚水饱和度SSSBVI是根据确定的权系数分量得到计算的束缚水饱和度SSSBVI。
8.如权利要求1所述的一种基于核磁共振谱系数法计算束缚水饱和度的方法，其特征在于，所述束缚水饱和度SSSBVI与实际束缚水饱和度相比，绝对误差在4个单位以内，该单位以百分比计；相对误差在10％以内。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510541975.XA CN105223116B (zh) | 2015-08-28 | 2015-08-28 | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510541975.XA CN105223116B (zh) | 2015-08-28 | 2015-08-28 | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105223116A CN105223116A (zh) | 2016-01-06 |
| CN105223116B true CN105223116B (zh) | 2017-09-01 |
Family
ID=54992191
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510541975.XA Active CN105223116B (zh) | 2015-08-28 | 2015-08-28 | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105223116B (zh) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105804734B (zh) * | 2016-04-12 | 2018-10-16 | 中国石油天然气集团公司 | 一种利用核磁共振测井识别稠油储层的方法 |
| CN106706686B (zh) * | 2016-12-06 | 2019-04-26 | 东北大学 | 一种确定富有机质页岩水饱和度与水化损伤关系的方法 |
| US10338014B2 (en) * | 2017-07-27 | 2019-07-02 | Saudi Arabian Oil Company | Estimating formation properties using saturation profiles |
| CN108007954A (zh) * | 2017-11-29 | 2018-05-08 | 西南石油大学 | 一种致密岩心束缚水环境建立的实验装置及测试方法 |
| CN108956417B (zh) * | 2018-06-09 | 2019-09-24 | 中国石油大学(华东) | 分析岩石孔隙无效吸水量的同位素核磁方法 |
| CN109781765A (zh) * | 2019-01-18 | 2019-05-21 | 西南石油大学 | 一种计算致密储层束缚水液膜厚度的新方法 |
| CN110146940B (zh) * | 2019-04-04 | 2021-03-26 | 中国石油大学(华东) | 一种气藏束缚水饱和度、气藏可动水范围确定方法和装置 |
| CN112147172B (zh) * | 2019-06-27 | 2024-07-02 | 中国石油化工股份有限公司 | 基于核磁共振t2谱评估含水饱和度的方法和装置 |
| CN112782477B (zh) * | 2019-11-11 | 2024-05-14 | 中国石油化工股份有限公司 | 一种测定岩心不同润湿状态下电性响应特征的方法及系统 |
| CN111220540B (zh) * | 2020-01-22 | 2020-12-18 | 长安大学 | 一种土壤毛细孔隙中土-水粘附力的测定方法 |
| CN113294146B (zh) * | 2020-02-05 | 2023-09-26 | 中国石油天然气股份有限公司 | 含沥青质储层的有效孔隙度计算的方法和装置 |
| CN112540096B (zh) * | 2020-11-27 | 2022-03-04 | 武汉大学 | 饱和冻结岩石未冻束缚水和未冻自由水含量的获取方法 |
| CN113153284B (zh) * | 2021-04-30 | 2023-06-30 | 中国石油天然气股份有限公司 | 束缚水饱和参数的确定方法、装置、设备和存储介质 |
| CN113447987B (zh) * | 2021-06-24 | 2022-03-29 | 中国石油大学(北京) | 一种确定地层流体饱和度的方法、装置及设备 |
| CN113586043B (zh) * | 2021-07-19 | 2023-10-27 | 中国石油天然气股份有限公司 | 束缚水饱和参数的确定方法、装置和计算机设备 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103225506B (zh) * | 2013-04-19 | 2015-06-03 | 中国石油大学(华东) | 三组分自动混联导电饱和度模型的建立方法 |
| CN103279651B (zh) * | 2013-05-10 | 2016-12-28 | 中国石油天然气股份有限公司 | 一种储层参数预测方法及装置 |
| CN103884633B (zh) * | 2014-03-05 | 2016-05-11 | 中国石油天然气股份有限公司 | 一种确定岩石渗透率的方法及装置 |
| RU2558838C1 (ru) * | 2014-07-02 | 2015-08-10 | Общество с ограниченной ответственностью "Научно-исследовательский институт природных газов и газовых технологий-Газпром ВНИИГАЗ" | Способ моделирования и оценки активного объема подземного хранилища газа в водоносных трещиновато-поровых структурах |
| CN104636819A (zh) * | 2014-12-31 | 2015-05-20 | 中国石油天然气集团公司 | 一种利用储层有效厚度加权系数进行储层定量产能预测的方法 |
-
2015
- 2015-08-28 CN CN201510541975.XA patent/CN105223116B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105223116A (zh) | 2016-01-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105223116B (zh) | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 | |
| CN109838230B (zh) | 油藏水淹层的定量评价方法 | |
| CN105866009B (zh) | 一种计算致密油储层有效孔隙度的方法及装置 | |
| CN102253069B (zh) | 一种根据核磁共振t2谱确定渗透率的方法和装置 | |
| CN101892837B (zh) | 地层因数确定方法及含油饱和度确定方法 | |
| CN105891089B (zh) | 储层渗透率确定方法及装置 | |
| CN103422853B (zh) | 一种确定地层含水饱和度的方法及装置 | |
| CN104634804B (zh) | 一种利用核磁共振t2谱确定储层相对渗透率的方法 | |
| CN103884633B (zh) | 一种确定岩石渗透率的方法及装置 | |
| CN106050225B (zh) | 一种核磁共振测井t2谱100％纯水谱的确定方法 | |
| CN105114064B (zh) | 确定致密砂岩储层饱和度的方法 | |
| CN109443867B (zh) | 一种对致密岩石的物性参数进行连续检测的方法 | |
| Bagarello et al. | Estimating field‐saturated soil hydraulic conductivity by a simplified Beerkan infiltration experiment | |
| CN109725016B (zh) | 一种用于含重油、沥青质岩心的核磁共振实验室测量方法 | |
| CN105804734A (zh) | 一种利用核磁共振测井识别稠油储层的方法 | |
| CN109932297B (zh) | 一种致密砂岩储层渗透率的计算方法 | |
| CN109580689B (zh) | 一种核磁共振测井t2截止值的逐点计算方法 | |
| CN104374679B (zh) | 一种确定储层渗透率贡献值的方法 | |
| Fityus et al. | Water content measurement in expansive soils using the neutron probe | |
| CN105547957B (zh) | 土壤下渗湿润峰吸力和降雨初损量计算方法 | |
| CN102052074A (zh) | 一种岩心数字化确定岩电参数的方法 | |
| CN105547959B (zh) | 致密基岩中基质孔隙储集天然气能力的分析方法 | |
| CN112049623B (zh) | 一种中基性火山岩束缚水出水的判定方法和装置 | |
| CN103046925B (zh) | 基于二项式的凝析气藏绝对无阻流量的获取方法及系统 | |
| CN113552036A (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 |