CN111520134B - 原始含油饱和度的获取方法和装置、电子设备及存储介质 - Google Patents
原始含油饱和度的获取方法和装置、电子设备及存储介质 Download PDFInfo
- Publication number
- CN111520134B CN111520134B CN202010261122.1A CN202010261122A CN111520134B CN 111520134 B CN111520134 B CN 111520134B CN 202010261122 A CN202010261122 A CN 202010261122A CN 111520134 B CN111520134 B CN 111520134B
- Authority
- CN
- China
- Prior art keywords
- oil saturation
- reservoir
- resistivity
- original
- original oil
- 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
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000003860 storage Methods 0.000 title claims abstract description 16
- 238000013178 mathematical model Methods 0.000 claims abstract description 57
- 238000011084 recovery Methods 0.000 claims abstract description 43
- 230000008859 change Effects 0.000 claims abstract description 36
- 238000011161 development Methods 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims description 28
- 239000011148 porous material Substances 0.000 claims description 17
- 230000015572 biosynthetic process Effects 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 11
- 239000008398 formation water Substances 0.000 claims description 11
- 230000007704 transition Effects 0.000 claims description 10
- 239000011435 rock Substances 0.000 claims description 8
- 239000003921 oil Substances 0.000 description 265
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 20
- 230000000875 corresponding effect Effects 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 239000007924 injection Substances 0.000 description 7
- 238000002347 injection Methods 0.000 description 7
- 239000010779 crude oil Substances 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000004973 liquid crystal related substance Substances 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000706 filtrate Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本申请公开了一种原始含油饱和度的获取方法和装置、电子设备及存储介质，所述获取方法包括：根据储层的深电阻率与中电阻率的幅度差和储层含油性的变化规律，建立原始含油饱和度的数学模型，数学模型示出了原始含油饱和度和深电阻率与中电阻率的幅度差、以及储层的目前含油饱和度的数学关系；根据原始含油饱和度的数学模型计算得到原始含油饱和度。本申请的方法结合油藏储层的深‑中电阻率幅度差和含油性的变化，因此，能够提高计算注水开发油藏每口井的原始含油饱和度的准确性，并且通过建立的数学模型计算，因此操作快速、低成本，为油藏采收率的重新标定提供了准确的基础数据，保证高含水期油藏高效合理的开发。
Description
技术领域
本发明一般涉及油田开发检测技术领域，具体涉及一种原始含油饱和度的获取方法和装置、电子设备及存储介质。
背景技术
原始含油饱和度即油藏在原始未开发状态下孔隙体积内的含油百分比，它是储层含油性评价、油藏产能预测、地质储量计算的关键参数，在整个油藏生命周期的持续开发中起着重要的支撑作用。中国绝大多数油田为陆相沉积，主要采用注水方式开采，注水开发油田的储量占总储量的85％以上，经过几十年的注水开发，大多数油田已经进入高含水、高采出程度阶段；另外，随着常规油气藏开发潜力的不断下降，一些非均质性强、孔隙结构复杂的特殊油气藏逐渐成为增储挖潜的主要对象。因此，针对注水开发的强水淹油气藏以及具有复杂孔隙结构的特殊油气藏，原始含油饱和度的准确计算成为油气藏高效开发的前提。
目前常见的测量储层原始含油饱和度的方法主要有利用密闭取心资料确定原始含油饱和度、测定曲线和岩电参数利用Archie公式(阿尔奇公式)确定原始含油饱和度、利用毛细管压力曲线和J函数确定原始含油饱和度。但是，现有的三种方法存在成本比较大、操作不便、计算结果精度不高的等问题，同时对于非均质性强、孔隙结构复杂的特殊油气藏，无法得到准确的原始含油饱和度。
发明内容
鉴于现有技术中的上述缺陷或不足，期望提供一种原始含油饱和度的获取方法和装置、电子设备及存储介质，实现原始含油饱和度的准确获取。
第一方面，本发明提供的一种原始含油饱和度的获取方法，包括：
根据储层的深电阻率与中电阻率的幅度差和储层含油性的变化规律，建立原始含油饱和度的数学模型，数学模型示出了原始含油饱和度和深电阻率与中电阻率的幅度差、以及储层的目前含油饱和度的数学关系；
根据原始含油饱和度的数学模型计算得到原始含油饱和度。
进一步地，根据储层的深电阻率与中电阻率的幅度差和储层含油性的变化规律，建立原始含油饱和度的数学模型的步骤中，原始含油饱和度的数学模型表示为：
Sor＝[B*e-A*(Rt-Ri)+1]*So
式中，Sor表示储层的原始含油饱和度，％；
So表示储层的目前含油饱和度，％；
Rt表示储层的目前原状地层电阻率，即就是深电阻率，单位：Ω·m；
Ri表示储层的目前过渡带电阻率，即就是中电阻率，单位：Ω·m；
A和B分别表示含油饱和度的恢复系数，无量纲。
进一步地，根据原始含油饱和度的数学模型计算得到原始含油饱和度的步骤，还包括：确定恢复系数A和恢复系数B，具体包括：
根据需要获取原始含油饱和度的油藏储层已有的岩心实验数据和测井曲线，获取不同深电阻率与中电阻率的幅度差所对应的含油饱和度值；
根据不同深电阻率与中电阻率的幅度差与所对应的含油饱和度值，拟合含油饱和度比差值和深电阻率与中电阻率的幅度差的变化趋势图；其中，含油饱和度的比差值是储层开发初期储层的某个含油饱和度值与每一个深电阻率与中电阻率的所对应的含油饱和度值的比值再减去常数1；
根据变化趋势图通过最小二乘法确定储层的恢复系数A和恢复系数B。
进一步地，根据原始含油饱和度的数学模型计算得到原始含油饱和度的步骤，还包括：确定储层目前的含油饱和度，具体过程包括：
获取目前油藏状态下储层的工作参数，工作参数包括：岩性系数、与岩性有关的常数、孔隙指数、地层水电阻率、地层电阻率，储层孔隙度；
根据储层的工作参数，利用Archie公式计算得到储层目前的含油饱和度。
第二方面，本发明提供的一种原始含油饱和度的获取装置，包括：
建立模块，用于根据储层的深电阻率与中电阻率的幅度差和储层含油性的变化规律，建立原始含油饱和度的数学模型，数学模型示出了原始含油饱和度和深电阻率与中电阻率的幅度差、以及储层目前含油饱和度的数学关系；
计算模块，用于根据原始含油饱和度的数学模型计算得到原始含油饱和度。
进一步地，建立模块，执行用于根据储层的深电阻率与中电阻率的幅度差和储层含油性变化规律，建立原始含油饱和度的数学模型的过程，得到的原始含油饱和度的数学模型表示为：
Sor＝[B*e-A*(Rt-Ri)+1]*So
式中，Sor表示储层原始含油饱和度，％；
So表示储层目前含油饱和度，％；
Rt表示储层的目前原状地层电阻率，即就是深电阻率，单位：Ω·m；
Ri表示储层的目前过渡带电阻率，即就是中电阻率，单位：Ω·m；
A和B分别表示含油饱和度的恢复系数，无量纲。
进一步地，计算模块，执行用于根据原始含油饱和度的数学模型计算得到原始含油饱和度的过程，还包括：
获取单元，用于根据需要获取原始含油饱和度的油藏储层已有的岩心实验数据和测井曲线，获取不同深电阻率与中电阻率的幅度差所对应的含油饱和度值；
拟合单元，用于根据不同深电阻率与中电阻率的幅度差与所对应的含油饱和度值，拟合含油饱和度比差值和深电阻率与中电阻率的幅度差的变化趋势图；其中，含油饱和度的比差值是储层开发初期储层的某个含油饱和度值与每一个深电阻率与中电阻率的幅度差所对应的含油饱和度值的比值再减去常数1；
第一确定单元，用于根据变化趋势图通过最小二乘法确定储层的恢复系数A和恢复系数B。
进一步地，计算模块，执行用于根据原始含油饱和度的数学模型计算得到原始含油饱和度的过程，还包括第二确定单元，用于确定储层目前的含油饱和度。
第三方面，本发明提供一种电子设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，处理器执行程序时实现如第一方面所述的原始含有饱和度的获取方法。
第四方面，本发明提供的一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现如第一方面所述的原始含有饱和度的获取方法。
本申请的一种原始含油饱和度的获取方法和装置，电子设备及存储介质，通过电阻率径向幅度差与储层含油性变化的相关关系，建立原始含油饱和度的计算模型，通过现有的测井曲线和分析资料，确定不同储层类型的恢复系数，结合储层目前的含油饱和度，准确、快速地确定储层的原始含油饱和度，为注水开采油藏提供了可靠地技术支撑。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述，本申请的其它特征、目的和优点将会变得更明显：
图1为本发明的一实施例的原始含油饱和度的获取方法的流程示意图；
图2为深中电阻率幅度差与储层含油性变化特征对比图，其中，实线是深电阻率，虚线是中电阻率；
图3为本发明的一实施例中S20实施流程图；
图4为本发明的另一实施例的原始含油饱和度的获取装置的示意图；
图5本发明的另一实施例的计算模块200的示意图；
图6为七区八道组含油饱和度比差值与深中电阻率幅度差交会图；
图7为七区八道湾组XX井原始含油饱和度计算成果图；
图8为适用来实现本申请实施例的电子设备或伏缺期的计算机***600的结构示意图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释相关发明，而非对该发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与发明相关的部分。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
原始含油饱和度对整个油藏全生命周期的开发具有重要的意义。目前，储层原始含油饱和度的确定主要有三种方法：(1)利用密闭取心资料确定原始含油饱和度，该方法确定的原始含油饱和度数值最准确，可以作为储层的参考数据，但是其测量成本比较大，并且需要在油藏开发初期进行取心实验，当油藏注水开发以后，含油饱和度会发生变化，此时的测量值与原始值会存在差异。(2)利用测井曲线和岩电参数确定原始含油饱和度，该方法主要基于Archie公式，利用储层岩石物性、电阻率测井曲线及地层水资料对含水饱和度进行解释，其准确计算的前提是公式中的各个参数必须是油藏原始状态下获取的，包括孔隙度、岩电参数、地层水电阻率及原状地层电阻率；当储层注水开发以后，要利用Archie公式计算原始含油饱和度，虽然其他参数可以根据已有资料确定原始油藏状态下的数值，但是，由于电阻率的影响因素比较多且变化复杂，地层原始电阻率难以准确恢复，因此，计算的难度比较大。(3)利用毛管压力曲线和J函数确定原始含油饱和度，该方法认为油气运移的主要驱动力是浮力，成藏过程中油气通过克服毛细管压力而排驱孔隙体积内的水，最终达到平衡，此时含水饱和度随深度的变化具有毛管压力曲线的分布特征，对于同一油藏，该油藏条件下的J函数与毛细管压力具有相同的分布特征，通过实验数据建立含水饱和度与J函数的关系式，可以确定油藏原始含油饱和度；毛管压力曲线法适用于非均质性弱、储层连续性好的构造油气藏，对于孔隙结构复杂、非均质性强的特殊油气藏，由于毛管压力与不同油藏高度下的含水饱和度不具备对应关系，因而，不能确定该类油气藏的原始含油饱和度。
因此，可看出三种方法各有优缺点，密闭取心法测量数据最准确，但是成本比较大，而且需要在油藏原始状态下取心测试；测井曲线计算法可以准确计算每口井、每个层位的原始含油饱和度，但是对于注水开发油气藏，原始地层电阻率难以准确恢复，计算结果精度不高；毛管压力曲线法不用考虑储层电性对含油性的影响，适用于低阻、致密等构造油气藏，但是对于非均质性强、孔隙结构复杂的特殊油气藏，其计算公式不适用。
基于此，本申请提出一种原始含油饱和度的获取方法，根据原状地层电阻率(记为：深电阻率)与过渡带电阻率(记为：中电阻率)的径向幅度差(记为深-中电阻率的幅度差)可以表征含油性的变化，通常情况下，随着水淹程度增强，深-中电阻率的幅度差逐渐减小，储层孔隙原油减少，含油饱和度变小，超强水淹甚至会出现电阻率反转，幅度差呈现负值。因此，可以利用深-中电阻率的幅度差，结合原始油藏状态下与目前油藏状态下的含油饱和度值，建立数学模型，进而得到储层的原始含油饱和度，为油藏的合理高效开发提供技术支撑。
本申请的一个实施例的一种原始含油饱和度的获取方法，如图1所示，包括：
S10，根据深电阻率与中电阻率(深-中电阻率)的幅度差与储层含油性的变化规律，建立原始含油饱和度的数学模型，数学模型包括原始含油饱和度与深-中电阻率幅度差和储层的目前含油饱和度的数学关系；
在上述步骤中，深-中电阻率幅度差与储层含油性变化规律可以通过对油藏已有的测井曲线与岩心实验数据的综合分析得到，从现有测井曲线和岩心实验数据可得：在储层为油层且保持原始油藏状态的前提下，原状地层电阻率(深电阻率Rt)与过渡带电阻率(中电阻率Ri)的径向幅度差比较大，其主要原因是过渡带地层中部分原油被泥浆滤液所替代，导致电阻率降低，而原状地层中则保持原始的油水分布特征，因而深-中电阻率幅度差较大。图2示出了深中电阻率幅度差与储层含油性变化特征对比图，从图2中看出，当储层注水开发以后，原状地层中部分原油被注入水驱替，含油饱和度减小，电阻率降低，深-中电阻率幅度差减小(图中实线代表深电阻率Rt，虚线代表中电阻率Ri)；随着储层水淹程度的不断增强，孔隙原油持续减少，深-中电阻率幅度差也逐渐变小，当储层超强水淹时，甚至会出现电阻率反转，幅度差呈现负值。因此，深-中电阻率的幅度差与储层含油饱和度存在相关关系，能够通过数学模型进行定量表征。
因此，结合岩心实验数据和理论分析，得到数学模型，该数学模型很好的反映了原始含油饱和度与深-中电阻率幅度差和目前储层含油饱和度的数学关系。该数学模型的建立结合的油藏整个开发期的测井曲线以及分析实验数据，能够准确反映储层含油饱和度与幅度深-中电阻率幅度差的关系。方便计算原始含油饱和度，并且操作简单，不受油气藏性质大的影响，适用范围广，适用于非均质性弱、储层连续性好的构造油气藏，同样适用于孔隙结构复杂、非均质性强的特殊油气藏。
S20，根据原始含油饱和度的数学模型计算得到原始含油饱和度。
根据原始含油饱和度的数学模型可以得到，储层原始含油饱和度与深-中电阻率幅度差和目前含油饱和度有关。由于储层类型存在差异，该数学模型针对不同储层时需要确定恢复系数，以提高原始含油饱和度的计算精度。
本申请的原始含油饱和度获取方法，根据深-中电阻率幅度差与储层含油性变化规律，建立原始含油饱和度的数学模型，根据原始含油饱和度的数学模型，确定原始含油饱和度。本申请公开的方法，通过分析油藏开采以来的所有测井曲线和岩心实验数据，得到深-中电阻率幅度差与储层含油性的变化特征，从而建立了原始含油饱和度的计算模型，提高了计算原始含油饱和度的精确性，为油藏开采提供合理有效的技术支撑。
在本实施例的基础上，本实施例涉及的根据深-中电阻率幅度差与储层含油性变化规律，建议原始含油饱和度的数学模型的过程，其中所述数学模型可以表示如式(1)所示：
Sor＝[B*e-A*(Rt-Ri)+1]*So (1)
式中，Sor表示储层原始含油饱和度，％；
So表示储层目前含油饱和度，％；
Rt表示储层的目前原状地层电阻率，即就是深电阻率，单位：Ω·m；
Ri表示储层的目前过渡带电阻率，即就是中电阻率，单位：Ω·m；
A和B分别表示含油饱和度的恢复系数，无量纲。
由原始含油饱和度的计算公式(1)可以看出，储层原始含油饱和度包括深-中电阻率幅度差和目前含油饱和度两部分。其中，当幅度差较大时，含油饱和度的校正量比较小；当幅度差逐渐变小时，储层水淹强度增强，含油饱和度减少，校正量变大。其中油藏储层的深-中电阻率幅度差从测井曲线可以得到。另外，由于储层类型存在差异，恢复系数A和恢复系数B需要针对不同储层分别进行确定，以提高原始含油饱和度的计算精度。
在上述实施例的基础上，还需要进一步确定恢复系数A和恢复系数B，如图3所示，具体包括如下：
S201，根据需要获取原始含油饱和度的油藏储层已有的岩心实验数据和测井曲线，获取不同深-中电阻率幅度差所对应的含油饱和度值；
S202，根据不同深-中电阻幅度差与所对应的含油饱和度值，拟合含油饱和度比差值与深-中电阻率幅度差的变化趋势图，其中，含油饱和度的比差值是储层开发初期储层的某个含油饱和度值与每一个深-中电阻率幅度差所对应的含油饱和度值的比值再减去常数1；
S203，根据变化趋势图通过最小二乘法确定该类储层的恢复系数A和恢复系数B。
需要说明的是：
该实施方式中，恢复系数的计算方式对于其他类型的储层同样适用。
在S202中，绘制含油饱和度比差值与深-中电阻率幅度差变化的趋势图，其中含油饱和度比差值指的是，由于油藏自开发以来所有的含油饱和度分析资料、测井曲线以及岩心数据都是能够获取得到的，因此选择油藏在开发初期的某一个含油饱和度值记为Sor，将该数值与每个不同深-中电阻率幅度差对应的含油饱和度值So做比值后减去常数1，从而将该比差值(Sor/So-1)与不同深-中电阻率幅度差建立趋势图(也称为交会图)，即以Rt-Ri为横坐标，(Sor/So-1)为纵坐标，拟合出(Sor/So-1)与Rt-Ri曲线，通过对交互图进行处理得到恢复系数。此处，需要说明的是，含油饱和度的比差值还可以选择任意一段油藏开发以来的数据得到，但是为了更准确计算油藏原始含油饱和度，选择油藏开发初期的数据更合适。
在上述实施的基础上，根据建立的数学模型，可以看出，为了获取原始含油饱和度，还需确定储层目前的含油饱和度，储层目前的含油饱和度计算过程如下：
获取目前油藏状态下储层的工作参数，工作参数包括：岩性系数、与岩性有关的常数、孔隙指数、地层水电阻率、地层电阻率，储层孔隙度；
根据储层的工作参数，利用Archie公式计算得到储层目前的含油饱和度。
需要说明的是，储层的工作参数均可以通过岩心实验和测井曲线得到，其中储层孔隙度可以依据补偿密度测井曲线求取。
进一步地，Archie公式可以表示如式(2)所示：
其中，Sw为储层目前的含水饱和度；a为岩性系数，无量纲；b为与岩性有关的常数，记为第一拟合常数，无量纲；n为与岩性有关的常数，记为第二拟合常数，无量纲；m为孔隙指数，无量纲；Rw为地层水电阻率，单位为Ω.m；Rt为地层电阻率，单位为Ω.m；Φ为孔隙度，无量纲。
在本实施例中，上述公式中的第一拟合参数b和第二拟合参数n可以根据实验法获得。
本发明提供的原始含油饱和度的获取方法，通过分析储层深-中电阻率幅度差与储层含油性的变化，建立了原始含油饱和度的数学模型，综合考虑油藏开发以来的岩心实验数据和测井曲线数据，得到储层目前含油饱和度以及恢复系数A和恢复系数B，通过目前含油饱和度、恢复系数和深-中电阻率幅度差，得到原始含油饱和度值。本发明的方法结合油藏储层的深-中电阻率幅度差和含油性的变化，因此，能够提高计算注水开发油藏每口井的原始含油饱和度的准确性，并且通过建立的数学模型计算，因此操作快速、低成本。提供了油藏采收率的重新标定提供准确的基础数据，保证高含水期油藏高效合理的开发。
图1示出的本发明提供的原始含油饱和度的获取方法。上述实施例的执行主体可以是原始含油饱和度的获取装置，该原始含油饱和度的获取装置可以通过软件和/或硬件实现，该原始含油饱和度的获取装置可以设置在处理器中，或者为单独的处理器，为了便于阐述，以下将本实施例的执行主体简称为获取装置。上述实施例涉及的是获取装置根据深-中电阻率幅度差与储层含油性变化规律，建立原始含油饱和度的数学模型，根据原始含油饱和度的数学模型确定原始含油饱和度的具体过程。
第二方面，本发明的另一个实施例提供了一种原始含油饱和度的获取装置，如图4所示，可以包括：
建立模块100，根据深-中电阻率幅度差与储层含油性的变化规律，建立原始含油饱和度的数学模型，数学模型包括原始含有饱和度与深-中电阻率幅度差和目前储层含油饱和度的数学关系；
计算模块200，根据原始含油饱和度的数学模型计算得到原始含油饱和度。
本实施例的装置，可以用于执行上述所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
作为可选的实施方式，建立模块200具体建立的数学模型如下：
Sor＝[B*e-A*(Rt-Ri)+1]*So (1)
式中，Sor表示储层原始含油饱和度，％；
So表示储层目前含油饱和度，％；
Rt表示储层的目前原状地层电阻率，即就是深电阻率，单位：Ω·m；
Ri表示储层的目前过渡带电阻率，即就是中电阻率，单位：Ω·m；
A和B分别表示含油饱和度的恢复系数，无量纲。
进一步地，如图5所示，计算模块200还用于确定恢复系数A和恢复系数B，具体的执行功能单元可以包括：
获取单元201，用于根据需要获取原始含油饱和度的油藏储层已有的岩心实验数据和测井曲线，获取不同深-中电阻率幅度差所对应的含油饱和度值；
拟合单元202，用于根据不同深-中电阻幅度差与所对应的含油饱和度值，拟合含油饱和度比差值与深-中电阻率幅度差的变化趋势图，其中，含油饱和度的比差值是储层开发初期储层的某个含油饱和度值与每一个深-中电阻率幅度差所对应的含油饱和度值的比值再减去常数1；
第一确定单元203，用于根据趋势图通过最小二乘法确定该类储层的恢复系数A和恢复系数B。
此处，拟合单元202可以用于执行上述所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
进一步地，计算模块200还可以包括：
第二确定单元204，用于确定储层目前的含油饱和度，具体执行过程如下：
获取目前油藏状态下储层的工作参数，工作参数包括：岩性系数、与岩性有关的常数、孔隙指数、地层水电阻率、地层电阻率，储层孔隙度；
根据储层的工作参数，利用Archie公式计算得到储层目前的含油饱和度。
需要说明的是，储层的工作参数均可以通过岩心实验和测井曲线得到，其中储层孔隙度可以依据补偿密度测井曲线求取。
进一步地，Archie公式可以表示如式(2)所示：
其中，Sw为储层目前的含水饱和度；a为岩性系数，无量纲；b为与岩性有关的常数，无量纲；m为孔隙指数，无量纲；Rw为地层水电阻率，单位为Ω.m；Rt为地层电阻率，单位为Ω.m；Φ为孔隙度，无量纲。
本实施例的装置，可以用于执行上述所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
本发明的原始含油饱和度的获取方法，举例说明其具体过程和可行性。示例如下：
以克拉玛依油田七区八道湾组砾岩油藏为例，该油藏经过几十年的注水开发，目前已经进入高含水开发阶段，原始含油饱和度的准确恢复是油藏可持续开发的关键。八道湾组整体上属于冲积扇沉积体系中的辫流水道沉积，储层类型基本一致，通过对已有测井曲线与岩心实验数据的综合分析，发现储层开发初期深-中电阻率径向幅度差比较比较大；当储层注水开发以后，深-中电阻率幅度差减小，随着储层水淹程度的不断增强，孔隙原油持续减少，深-中电阻率幅度差也逐渐变小；在一些超强水淹段，甚至出现电阻率反转，幅度差呈现负值的情况。因此，对于八道湾组储层，深-中电阻率的幅度差与储层含油饱和度的变化具有很好的相关性，可以利用数学模型进行定量表征。
结合分析，根据七区八道湾组砾岩储层深-中电阻率的径向幅度差与含油性的变化规律，经建立建原始含油饱和度计算公式如下：
Sor＝[B*e-A*(Rt-Ri)+1]*So
式中，Sor表示储层原始含油饱和度，％；So表示储层目前含油饱和度，％；Rt表示目前储层原状地层电阻率，即就是深电阻率，单位：Ω·m；Ri表示目前储层过渡带电阻率，即就是中电阻率，单位：Ω·m；A和B分别表示含油饱和度的恢复系数，无量纲。
针对七区八道湾组砾岩储层，收集自储层开发以来所有的含油饱和度分析资料和测井曲线，确定不同深-中电阻率幅度差所对应的储层含油饱和度值，结合油藏初期开发的储层含油饱和度，绘制含油饱和度比差值与深中电阻率幅度差的变化趋势图，含油饱和度比差值(Sor/So-1)即对式(1)进行变形整理，把饱和度项放到等式一边；进而依据建立的交会图，如图6所示，确定七区八道湾组砾岩储层的恢复系数A和恢复系数B，其中A＝0.06，B＝0.2959。
目前含油饱和度利用Archie公式进行计算：
So＝100-Sw
式中，各个参数均取目前油藏状态下的数值，其中，a＝0.646；m＝1.919；b＝1.0148；n＝2.0819；地层水主要以NaHCO3型为主，电阻率为Rw＝0.42Ω·m；地层电阻率为Rt通过测量得到；储层总孔隙度Φ可以依据补偿密度测井曲线求取。
在确定七区八道湾组砾岩储层恢复系数A和B以及油藏目前含油饱和度的基础上，结合目前储层的测井曲线，利用上述数学模型可以准确计算每口井八道湾组储层的原始含油饱和度。
取七区八道湾组储层密闭取心分析资料的井，利用所建立的模型计算储层的原始含油饱和度。图7是七区八道湾组XX井的原始含油饱和度的计算成果图，在含油饱和度道中，深灰色表示储层目前含油饱和度，利用Archie公式及目前储层的原状地层电阻率测井曲线、岩电参数、地层水电阻率可以准确计算；浅灰色表示储层原始含油饱和度，依据深-中电阻率测井曲线的径向幅度差，利用建立的数学模型及系数A和B，可以有效地恢复储层的原始含油性。首先，从测井计算饱和度道中可以看出，杆状图表示岩心分析的目前含油性饱和度，曲线表示测井计算的目前含油性饱和度，测井计算结果与岩心分析数据吻合程度比较高，相关系数达到了95.6％。进而，把原始含油饱和度的计算结果与临近老井的原始含油饱和度岩心分析资料进行对比，解释精度也达到了93.2％，可以准确恢复八道湾组储层的原始含油饱和度，为油藏采收率的标定和地质储量的计算提供准确的基础数据。
下面参考图8，图8示出了适用来实现本申请实施例的电子设备或伏缺期的计算机***600的结构示意图。
如图8所示，计算机***600包括中央处理单元(CPU)601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM)603中的程序而执行各种适当的动作和处理。在RAM 603中，还存储有***600操作所需的各种程序和数据。CPU 601、ROM 602以及RAM 603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡、调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装入存储部分608。
特别地，根据本公开的实施例，上文参考图1或图3描述的过程可以被实现为计算机软件程序。例如，本公开的实施例包括一种计算机程序产品，其包括承载在机器上可读介质上的计算机程序，该计算机程序包含用于执行图1或图3的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。在该计算机程序中被中央处理单元(CPU)601执行时，执行本申请的***中限定的上述功能。
附图中的流程图和框图，图示了按照本发明各种实施例的***、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分，前述模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个接连地表示的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的***来实现，或者可以用专用硬件与计算机指令的组合来实现。
描述于本申请实施例中所涉及到的单元或模块可以通过软件的方式实现，也可以通过硬件的方式来实现。所描述的单元或模块也可以设置在处理器中，例如，可以描述为：一种处理器包括建立单元和计算单元。其中，这些单元或模块的名称在某种情况下并不构成对该单元或模块本身的限定，例如，计算单元还可以被描述为“用于根据建立单元得到数学模型计算得到原始含油饱和度的单元”。
作为另一方面，本申请还提供了一种计算机可读存储介质，该计算机可读存储介质可以是上述实施例中前述装置中所包含的计算机可读存储介质；也可以是单独存在，未装配入设备中的计算机可读存储介质。计算机可读存储介质存储有一个或者一个以上程序，前述程序被一个或者一个以上的处理器用来执行描述于本申请的原始含油饱和度的获取。
以上描述仅为本申请的较佳实施例以及对所运用技术原理的说明。本领域技术人员应当理解，本申请中所涉及的发明范围，并不限于上述技术特征的特定组合而成的技术方案，同时也应涵盖在不脱离所述发明构思的情况下，由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。
Claims (8)
1.一种原始含油饱和度的获取方法，其特征在于，包括：
根据储层的深电阻率与中电阻率的幅度差和储层含油性的变化规律，建立原始含油饱和度的数学模型，所述数学模型示出了原始含油饱和度和所述深电阻率与中电阻率的幅度差、以及储层的目前含油饱和度的数学关系；
所述原始含油饱和度的数学模型表示为：
Sor＝[B*e-A*(Rt-Ri)+1]*So
式中，Sor表示储层的原始含油饱和度；
So表示储层的目前含油饱和度；
Rt表示储层的目前原状地层电阻率，即就是深电阻率；
Ri表示储层的目前过渡带电阻率，即就是中电阻率；
A和B分别表示含油饱和度的恢复系数；
根据所述原始含油饱和度的数学模型计算得到原始含油饱和度。
2.根据权利要求1所述的方法，其特征在于，所述根据所述原始含油饱和度的数学模型计算得到原始含油饱和度的步骤，还包括：确定所述恢复系数A和所述恢复系数B，具体包括：
根据需要获取原始含油饱和度的油藏储层已有的岩心实验数据和测井曲线，获取不同深电阻率与中电阻率的幅度差所对应的含油饱和度值；
根据所述不同深电阻率与中电阻率的幅度差与所对应的含油饱和度值，拟合含油饱和度比差值和所述深电阻率与中电阻率的幅度差的变化趋势图；其中，所述含油饱和度的比差值是储层开发初期储层的某个含油饱和度值与每一个所述深电阻率与中电阻率的所对应的含油饱和度值的比值再减去常数1；
根据所述变化趋势图通过最小二乘法确定储层的所述恢复系数A和所述恢复系数B。
3.根据权利要求1所述的方法，其特征在于，所述根据所述原始含油饱和度的数学模型计算得到原始含油饱和度的步骤，还包括：确定储层目前的含油饱和度，具体过程包括：
获取目前油藏状态下储层的工作参数，所述工作参数包括：岩性系数、与岩性有关的常数、孔隙指数、地层水电阻率、地层电阻率，储层孔隙度；
根据储层的所述工作参数，利用Archie公式计算得到储层目前的含油饱和度。
4.一种原始含油饱和度的获取装置，其特征在于，包括：
建立模块，用于根据储层的深电阻率与中电阻率的幅度差和储层含油性的变化规律，建立原始含油饱和度的数学模型，所述数学模型示出了原始含油饱和度和所述深电阻率与中电阻率的幅度差、以及储层目前含油饱和度的数学关系；
所述原始含油饱和度的数学模型表示为：
Sor＝[B*e-A*(Rt-Ri)+1]*So
式中，Sor表示储层原始含油饱和度；
So表示储层目前含油饱和度；
Rt表示储层的目前原状地层电阻率，即就是深电阻率；
Ri表示储层的目前过渡带电阻率，即就是中电阻率；
A和B分别表示含油饱和度的恢复系数；
计算模块，用于根据所述原始含油饱和度的数学模型计算得到原始含油饱和度。
5.根据权利要求4所述的装置，其特征在于，所述计算模块，执行用于根据所述原始含油饱和度的数学模型计算得到原始含油饱和度的过程，还包括：
获取单元，用于根据需要获取原始含油饱和度的油藏储层已有的岩心实验数据和测井曲线，获取不同深电阻率与中电阻率的幅度差所对应的含油饱和度值；
拟合单元，用于根据所述不同深电阻率与中电阻率的幅度差与所对应的含油饱和度值，拟合含油饱和度比差值和所述深电阻率与中电阻率的幅度差的变化趋势图；其中，含油饱和度的比差值是储层开发初期储层的某个含油饱和度值与每一个所述深电阻率与中电阻率的幅度差所对应的含油饱和度值的比值再减去常数1；
第一确定单元，用于根据所述变化趋势图通过最小二乘法确定储层的所述恢复系数A和所述恢复系数B。
6.根据权利要求4所述的装置，其特征在于，所述计算模块，执行用于根据所述原始含油饱和度的数学模型计算得到原始含油饱和度的过程，还包括第二确定单元，用于确定储层目前的含油饱和度。
7.一种电子设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1-3中任一项所述的方法。
8.一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1-3中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010261122.1A CN111520134B (zh) | 2020-04-03 | 2020-04-03 | 原始含油饱和度的获取方法和装置、电子设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010261122.1A CN111520134B (zh) | 2020-04-03 | 2020-04-03 | 原始含油饱和度的获取方法和装置、电子设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111520134A CN111520134A (zh) | 2020-08-11 |
| CN111520134B true CN111520134B (zh) | 2023-06-06 |
Family
ID=71901832
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010261122.1A Active CN111520134B (zh) | 2020-04-03 | 2020-04-03 | 原始含油饱和度的获取方法和装置、电子设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111520134B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113327005A (zh) * | 2021-04-20 | 2021-08-31 | 中海油能源发展股份有限公司 | 基于试井资料计算储层含油饱和度的方法和电子设备 |
| CN115099991B (zh) * | 2022-07-20 | 2023-03-24 | 中国科学院大学 | 一种基于生产数据的非均质油藏动态含油饱和度计算方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1243958A (zh) * | 1998-08-03 | 2000-02-09 | 邓友明 | 储层流体性质测量方法及其设备 |
| CN101487390A (zh) * | 2009-02-23 | 2009-07-22 | 大庆油田有限责任公司 | 一种确定油层原始含油饱和度的阿尔奇模式方法 |
| CN103485770A (zh) * | 2013-06-19 | 2014-01-01 | 中国石油天然气集团公司 | 一种基于人工神经网络求取含油饱和度的方法及*** |
| CN104500055A (zh) * | 2014-12-16 | 2015-04-08 | 中国石油天然气股份有限公司 | 一种特低渗透油藏水淹层含水饱和度计算方法 |
| CN104847340A (zh) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | 一种水淹层测井定量评价方法 |
| CN106093083A (zh) * | 2016-07-01 | 2016-11-09 | 龙威 | 一种采用数字岩心模拟建立岩电关系的方法 |
| CN106869916A (zh) * | 2017-03-22 | 2017-06-20 | 中国石油天然气股份有限公司 | 一种碎屑岩稠油储层识别方法及装置 |
| CN107701180A (zh) * | 2017-09-19 | 2018-02-16 | 中国石油天然气集团公司 | 一种基于密闭取心的原始油藏含水饱和度计算方法 |
| CN107808032A (zh) * | 2017-09-28 | 2018-03-16 | 中国石油化工股份有限公司 | 页岩气储层含水饱和度计算方法及*** |
| CN110618170A (zh) * | 2019-09-26 | 2019-12-27 | 中国石油天然气股份有限公司 | 岩土的饱和度指数的确定方法及装置 |
-
2020
- 2020-04-03 CN CN202010261122.1A patent/CN111520134B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1243958A (zh) * | 1998-08-03 | 2000-02-09 | 邓友明 | 储层流体性质测量方法及其设备 |
| CN101487390A (zh) * | 2009-02-23 | 2009-07-22 | 大庆油田有限责任公司 | 一种确定油层原始含油饱和度的阿尔奇模式方法 |
| CN103485770A (zh) * | 2013-06-19 | 2014-01-01 | 中国石油天然气集团公司 | 一种基于人工神经网络求取含油饱和度的方法及*** |
| CN104500055A (zh) * | 2014-12-16 | 2015-04-08 | 中国石油天然气股份有限公司 | 一种特低渗透油藏水淹层含水饱和度计算方法 |
| CN104847340A (zh) * | 2015-03-26 | 2015-08-19 | 中国海洋石油总公司 | 一种水淹层测井定量评价方法 |
| CN106093083A (zh) * | 2016-07-01 | 2016-11-09 | 龙威 | 一种采用数字岩心模拟建立岩电关系的方法 |
| CN106869916A (zh) * | 2017-03-22 | 2017-06-20 | 中国石油天然气股份有限公司 | 一种碎屑岩稠油储层识别方法及装置 |
| CN107701180A (zh) * | 2017-09-19 | 2018-02-16 | 中国石油天然气集团公司 | 一种基于密闭取心的原始油藏含水饱和度计算方法 |
| CN107808032A (zh) * | 2017-09-28 | 2018-03-16 | 中国石油化工股份有限公司 | 页岩气储层含水饱和度计算方法及*** |
| CN110618170A (zh) * | 2019-09-26 | 2019-12-27 | 中国石油天然气股份有限公司 | 岩土的饱和度指数的确定方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 谭锋奇 ; 李洪奇 ; 武鑫 ; 范小秦 ; 于红岩 ; .砾岩油藏水淹层定量识别方法――以新疆克拉玛依油田六中区克下组为例.石油与天然气地质.2010,(第02期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111520134A (zh) | 2020-08-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111520134B (zh) | 原始含油饱和度的获取方法和装置、电子设备及存储介质 | |
| CN106869916B (zh) | 一种碎屑岩稠油储层识别方法及装置 | |
| US8626447B2 (en) | System and method for sweet zone identification in shale gas reservoirs | |
| CN110955982B (zh) | 一种变质岩储层渗透率计算方法、装置及计算机存储介质 | |
| RU2455483C2 (ru) | Способ оценки степени трещиноватости карбонатных пород через параметр диффузионно-адсорбционной активности | |
| CN111008482B (zh) | 一种变质岩储层含水饱和度计算方法及装置 | |
| Zhang et al. | Improved normalized full-gradient method and its application to the location of source body | |
| CN113775326B (zh) | 可动水饱和度的评价方法、装置、电子设备及介质 | |
| CN116504326B (zh) | 确定页岩有效含油率的方法和电子设备 | |
| CN110344825A (zh) | 一种低孔低渗低电阻率砂岩油层综合判识方法 | |
| CN110632654B (zh) | 断块圈闭含油边界确定方法及装置 | |
| CN115099014A (zh) | 一种基于随钻录井的天然气井地质探明储量计算方法 | |
| CN113945992B (zh) | 泥岩、油页岩识别方法、装置、电子设备及介质 | |
| CN113552653A (zh) | 一种确定碎屑岩古老地层剥蚀厚度的方法 | |
| CN114086938A (zh) | 一种非均质性砂岩储层的含气饱和度预测方法 | |
| CN111287739A (zh) | 一种基于地层原油粘度的剩余油分布预测方法 | |
| US11867862B2 (en) | Method for validating rock formations compaction parameters using geomechanical modeling | |
| RU2248444C2 (ru) | Способ определения остаточной нефтенасыщенности пластов | |
| WO2024099049A1 (zh) | 一种泥页岩储层的可动油含量评价方法、***、设备与存储介质 | |
| US20220275707A1 (en) | Volumetric treatment fluid distribution to remove formation damage | |
| CN112696197A (zh) | 一种油田储层指数曲线构建方法、***、设备及存储介质 | |
| Ye | Dominant Channels Identification Method for Multi-Fractured Horizontal Wells in Tight Reservoirs: Progress and Challenges | |
| Daae et al. | Reservoir Simulation Box Model to Provide Quick-Look Forecast of Remobilizing Trapped Gas in a Gas-Water Jurassic Reservoir with Large Aquifer | |
| Omrani et al. | Examining Key Factors Influencing the Accuracy of Water Saturation Calculations in Carbonate Reservoirs: A Case Study of Kangan and Dalan Formations in the Western Persian Gulf | |
| Brailovskaya et al. | A comprehensive approach to the re-evaluation of the geological structure of'low-pore'deposits in a limited set and low quality of initial data on the example of deposits in the Republic of Ingushetia (Russian) |
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 |