一种含油饱和度的确定方法及装置

1.本说明书实施例涉及地质勘探技术领域，特别涉及一种含油饱和度的确定方法及装置。


背景技术：

2.含油饱和度是致密砂岩储层表征和储量计算的重要参数之一。在地层水矿化度变化较大、孔隙结构复杂、泥质含量较高等特征的储层，仅利用电阻率测井评价含油饱和度较为困难。除电阻率测井之外，介电测井资料也用于计算储层含油饱和度。介电测井测量地层兆赫兹频段复介电常数，特别是斯伦贝谢公司的阵列介电测井(adt)可提供兆赫兹频段介电频散曲线，因而蕴含丰富的地层信息。
3.现有技术中常基于介电测井资料利用正演模型评价含油饱和度，正演模型可以包括crim模型、泥质砂岩模型等，但这些正演模型有固定的假设条件，使得正演结果与实际情况会存在差异。并且利用介电测井评价储层的含油饱和度时，储层的泥质含量、地层水矿化度、矿物类型、孔隙结构等因素会影响模型性能，进而会影响确定的含油饱和度的精确度。因此，采用现有技术中的技术方案无法准确的确定储层含油饱和度。
4.针对上述问题，目前尚未提出有效的解决方案。


技术实现要素：

5.本说明书实施例提供了一种含油饱和度的确定方法及装置，以解决现有技术中无法准确的确定储层含油饱和度的问题。
6.本说明书实施例提供了一种含油饱和度的确定方法，包括：获取目标井的测井资料集；其中，所述测井资料集中包括：阵列感应测井电阻率曲线、介电测井资料中四种频率的介电常数曲线和电导率曲线；利用多种响应方程构建最小平方目标函数；基于所述目标井的测井资料集，利用粒子群优化算法对所述目标函数进行寻优反演，得到所述目标井的含水饱和度曲线；根据所述目标井的含水饱和度曲线，确定所述目标井的含油饱和度曲线。
7.本说明书实施例还提供了一种含油饱和度的确定装置，包括：获取模块，用于获取目标井的测井资料集；其中，所述测井资料集中包括：阵列感应测井电阻率曲线、介电测井资料中四种频率的介电常数曲线和电导率曲线；构建模块，用于利用多种响应方程构建最小平方目标函数；反演模块，用于基于所述目标井的测井资料集，利用粒子群优化算法对所述目标函数进行寻优反演，得到所述目标井的含水饱和度曲线；确定模块，用于根据所述目标井的含水饱和度曲线，确定所述目标井的含油饱和度曲线。
8.本说明书实施例提供了一种含油饱和度的确定方法，可以获取目标井的阵列感应测井电阻率曲线、介电测井资料中四种频率的介电常数曲线和电导率曲线。由于测井资料集中包含多种测井曲线，而不同的测井曲线适用的响应方程会存在差异，因此，可以利用多种响应方程构建最小平方目标函数，并基于所述目标井的测井资料集，利用粒子群优化算法对所述目标函数进行寻优反演，得到所述目标井的含水饱和度曲线。进一步的，可以根据
所述目标井的含水饱和度曲线，确定所述目标井的含油饱和度曲线。利用多种响应模型可以综合考虑储层的泥质含量、地层水矿化度、岩石孔隙结构等因素的影响，进而可以利用构建得到的最小平方目标函数实现准确反演得到含油饱和度的目的。
附图说明
9.此处所说明的附图用来提供对本说明书实施例的进一步理解，构成本说明书实施例的一部分，并不构成对本说明书实施例的限定。在附图中：
10.图1是根据本说明书实施例提供的含油饱和度的确定方法的步骤示意图；
11.图2是根据本说明书实施例提供的l1井中目标层段理论测井响应和实际测井响应的对比示意图；
12.图3是根据本说明书实施例提供的l1井的常规测井曲线、介电测井曲线和储层参数计算结果的示意图；
13.图4是根据本说明书实施例提供的l2井中目标层段常规测井、介电测井及其反演结果的示意图；
14.图5是根据本说明书实施例提供的含油饱和度的确定装置的结构示意图；
15.图6是根据本说明书实施例提供的含油饱和度的确定设备的结构示意图。
具体实施方式
16.下面将参考若干示例性实施方式来描述本说明书实施例的原理和精神。应当理解，给出这些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本说明书实施例，而并非以任何方式限制本说明书实施例的范围。相反，提供这些实施方式是为了使本说明书实施例公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
17.本领域的技术人员知道，本说明书实施例的实施方式可以实现为一种系统、装置设备、方法或计算机程序产品。因此，本说明书实施例公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
18.虽然下文描述流程包括以特定顺序出现的多个操作，但是应该清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
19.请参阅图1，本实施方式可以提供一种含油饱和度的确定方法，可以应用于致密砂岩等非常规油气储层。该含油饱和度的确定方法可以用于利用构建得到的最小平方目标函数准确的反演得到含油饱和度。上述含油饱和度的确定方法可以包括以下步骤。
20.s101：获取目标井的测井资料集；其中，测井资料集中包括：阵列感应测井电阻率曲线、介电测井资料中四种频率的介电常数曲线和电导率曲线。
21.在本实施方式中，为了定量评估目标井的地层情况，可以获取目标井的测井资料集，上述测井资料集可以为目标井预先测井得到的资料。由于致密砂岩等非常规油气储层岩石致密，大多数情况下泥浆侵入较浅，因此，相对于常规油气储层而言，介电测井在非常规油气储层饱和度评价中具有较大优势。基于此，上述测井资料集中可以包含常规测井资料和介电测井资料。
22.在本实施方式中，上述阵列感应测井电阻率曲线可以是基于常规测井资料获取的，上述阵列感应测井电阻率曲线可以为阵列感应测井10in(0.25m)的探测深度电阻率测井响应。当然可以理解的是，在一些实施例中还可以为其它探测深度的电阻率测井响应，例如：0.5米等，具体的可以根据实际情况确定，本说明书实施例对此不作限定。
23.在本实施方式中，上述四种频率所属的频段范围可以为20mhz(兆赫兹)～1ghz(千兆赫兹)之间，或者1ghz左右以及20mhz～500mhz。当然，四种频率所属的频段范围不限于上述举例，所属领域技术人员在本说明书实施例技术精髓的启示下，还可能做出其它变更，但只要其实现的功能和效果与本说明书实施例相同或相似，均应涵盖于本说明书实施例保护范围内。
24.s102：利用多种响应方程构建最小平方目标函数。
25.在本实施方式中，由于测井资料集中包含多种测井曲线，而不同的测井曲线适用的响应方程会存在差异，因此，可以利用多种响应方程构建最小平方目标函数。
26.在本实施方式中，上述响应方程可以包含复折射率(crim)模型、泥质砂岩(shsd)模型、西门杜方程(一种泥质砂岩电导模型)、stroud
–
milton
–
de(smd，一种考虑介电频散的纯岩石饱和流体响应方程)模型、waxman
‑
smith方程(w
‑
s方程，一种泥质砂岩电导率方程)等。当然可以理解的是，上述响应方程还可以包含更多或者更少的模型方程，具体的可以根据实际需求进行选择，本说明书实施例对此不作限定。
27.在本实施方式中，利用多种响应模型可以综合考虑储层的泥质含量、地层水矿化度、岩石孔隙结构等因素的影响，进而可以利用构建得到的最小平方目标函数实现准确评价致密砂岩储层含油饱和度的目的。
28.在本实施方式中，最小平方目标函数可以为利用最小二乘法构建的函数，最小二乘法(又称最小平方方法)是一种数学优化技术，通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。其中，上述实际数据在本实施例中可以是指测井资料集中的实测数据。
29.s103：基于目标井的测井资料集，利用粒子群优化算法对目标函数进行寻优反演，得到目标井的含水饱和度曲线。
30.在本实施方式中，可以基于目标井的测井资料集，利用粒子群优化算法对目标函数进行寻优反演，得到目标井的含水饱和度曲线。其中，上述目标井的测井资料集可以作为粒子群优化算法的输入数据。
31.在本实施方式中，粒子群算法是一种经典的仿生全局随机搜索算法，它的基本思路是将优化问题的每一个解视为一个代表人工鸟的粒子，每个粒子都有自己的位置和速度并根据一定的法则调整位置和速度，整体飞向最优目标。粒子的速度变化主要由粒子前一次迭代的速度、个体最佳适应值位置和种群最佳适应值位置三方面所决定。
32.粒子的速度和位置迭代公式分别如下所示：
[0033][0034][0035]
其中，是第i个粒子在第r次迭代的速度；ω是指数衰减的惯性权重系数；是
第i个粒子在第r次迭代的位置；为寻优过程中迭代至第r次为止第i个粒子最佳适应值的位置；gbest
r
为粒子群迭代至第r次最优适应值的位置；r1、r2是两个随机函数，取0～1间的随机数来增加寻优过程的随机性；c1、c2是学习因子。
[0036]
在本实施方式中，粒子群优化算法的实现步骤如下所示：
[0037]
①
将含水饱和度s
w
、水相曲折度指数mn、冲洗带矿化度sal以及泥质频散项体积含量p
d
四个反演参数设置为粒子，随机初始化位置和速度，设置迭代次数、粒子个数、迭代终止条件等参数；
[0038]
②
根据粒子的初始位置，利用目标函数、四组不同频率的相对介电常数和电导率、阵列感应测井电阻率曲线计算每个粒子的适应度值；
[0039]
③
更新每个粒子个体在寻优过程中的pbest
ir
和gbest
r
；
[0040]
④
根据速度迭代公式(1)赋予粒子新的速度；
[0041]
⑤
根据位置迭代公式(2)更新粒子的位置；
[0042]
⑥
重复步骤
②
至步骤
⑤
，直至满足迭代终止条件，输出四个反演参数的最优值。由于测井曲线是连续的点构成的曲线，在对测井曲线进行反演计算时，得到一系列的最优点构成的最优曲线。
[0043]
在本实施方式中，可以假设mn＝m＝n，m和n分别为胶结指数和饱和度指数。由于在含水纯砂岩或碳酸盐岩中，mn相当于archie公式(阿尔奇公式)中的m，因此，为了方便计算，可以假设mn＝m＝n。
[0044]
在本实施方式中，输入参数除了包含四组不同频率的相对介电常数和电导率、阵列感应测井电阻率曲线外，还可以包含孔隙度曲线、油气的相对介电常数、地层温度和岩石骨架相对介电常数等参数。具体的可以根据实际情况确定，本说明书实施例对此不作限定。
[0045]
s104：根据目标井的含水饱和度曲线，确定目标井的含油饱和度曲线。
[0046]
在本实施方式中，由于同一深度地层的含水饱和度与含油饱和度之和为1，因此，可以根据目标井的含水饱和度曲线，确定目标井的含油饱和度曲线。
[0047]
从以上的描述中，可以看出，本说明书实施例实现了如下技术效果：可以获取目标井的阵列感应测井电阻率曲线、介电测井资料中四种频率的介电常数曲线和电导率曲线。由于测井资料集中包含多种测井曲线，而不同的测井曲线适用的响应方程会存在差异，因此，可以利用多种响应方程构建最小平方目标函数，并基于所述目标井的测井资料集，利用粒子群优化算法对所述目标函数进行寻优反演，得到所述目标井的含水饱和度曲线。进一步的，可以根据所述目标井的含水饱和度曲线，确定所述目标井的含油饱和度曲线。利用多种响应模型可以综合考虑储层的泥质含量、地层水矿化度、岩石孔隙结构等因素的影响，进而可以利用构建得到的最小平方目标函数实现准确反演得到含油饱和度的目的。
[0048]
在一个实施方式中，所述多种响应方程包括：复折射率模型、泥质砂岩响应模型、西门杜方程。
[0049]
在本实施方式中，上述响应方程可以包含复折射率(crim)模型、泥质砂岩(shsd)模型、西门杜方程(一种泥质砂岩电导模型)等。当然可以理解的是，上述响应方程还可以包含更多或者更少的模型方程，例如：stroud
–
milton
–
de(smd，一种考虑介电频散的纯岩石饱和流体响应方程)模型、waxman
‑
smith方程(w
‑
s方程，一种泥质砂岩电导率方程)等，具体的可以根据实际需求进行选择，本说明书实施例对此不作限定。
[0050]
在本实施方式中，crim模型是解释1ghz左右介电常数的常用模型。由于在1ghz附近的界面极化效应可以忽略不计，而界面极化与岩石结构相关，因此，crim模型不包括岩石结构的变化。对于一个包含三种介质的crim模型，其表达式如下所示：
[0051][0052]
其中，ε
*
是岩石复介电常数；φ是岩石孔隙度；s
w
是含水饱和度；ε
m
和ε
h
分别是岩石骨架和油气的相对介电常数；ε
w*
是地层水的复介电常数。
[0053]
在本实施方式中，地层水的复介电常数ε
w*
由实部和虚部组成，其表达式为：
[0054][0055]
其中，ε
rw
为水的相对介电常数；ε0是自由空间的介电常数；σ
w
为水的直流电导率，s/m；ω为角频率。
[0056]
在本实施方式中，相对介电常数ε
rw
和电导率σ
w
均为储层温度和地层水矿化度的函数，采用如下公式计算：
[0057]
ε
rw
(0,t)＝94.88
‑
0.2317t+0.000217t2ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0058][0059][0060]
其中，ε
rw
(s,t)为在矿化度为s，地层温度为t时的相对介电常数；σ
w
(s,t)为在矿化度为s，地层温度为t时的电导率；t为地层温度，
°
f(华氏度)；s为矿化度，kppm(1000体积浓度)。上述公式(7)是根据斯伦贝谢公司的地层水电阻率与矿化度关系图版获得的。
[0061]
在一个实施方式中，泥质砂岩响应模型的表达式为：
[0062][0063]
其中，ε
*
为岩石复介电常数；φ为岩石孔隙度；s
w
为含水饱和度；p
d
为泥质频散项体积含量；mn为水相曲折度指数；为复合导电相；ε
rm
为岩石骨架和油气的混合介电常数；γ为伽马函数。
[0064]
在本实施方式中，泥质砂岩响应(shsd)模型是在考虑介电频散的纯岩石饱和流体smd模型的基础上，将湿粘土对介电频散增加的影响考虑进来。smd模型利用岩石孔隙度、流体饱和度、岩石骨架相对介电常数及岩石结构指数等计算饱和流体岩石兆赫兹频段复介电常数及其频散曲线，反过来可根据岩石复介电常数及频散曲线反演岩石流体饱和度。smd模型的表达式如下：
[0065][0066]
其中，ε
*
是岩石复介电常数；是地层水的复介电常数；φ为岩石孔隙度；s
w
为含水饱和度；mn是水相曲折度指数；γ是伽马函数；ε
rm
是岩石骨架和油气的混合介电常数。
[0067]
在本实施方式中，ε
rm
按照以下公式计算：
[0068][0069]
其中，ε
rm
是岩石骨架和油气的混合介电常数；ε
m
和ε
h
分别是岩石骨架和油气的相对介电常数；φ为岩石孔隙度；s
w
为含水饱和度。
[0070]
在本实施方式中，shsd模型在smd模型的充水孔隙中引入一种代表粘土诱发分散的频散相，粘土频散相按照以下公式计算：
[0071][0072]
其中，ε
s
是频散相的低频介电常数，它应该足够高，以表征粘土表面效应的频散性质；ε
∞
是频散项的高频介电常数，定义为干粘土介电常数；σ
d
是粘土束缚水的电导率；f
c
是临界频率，大约200mhz；ω为角频率。
[0073]
在本实施方式中，粘土相与游离水相混合，复合导电相被定义为：
[0074][0075]
其中，φ
w
是含水孔隙度；是地层水的复介电常数；是粘土频散相；p
d
是量化频散相的分数，即泥质频散项体积含量，它表征额外的粘土诱导效应对总频散的贡献。
[0076]
最后，将复合导电相代入smd模型中，可以得到一种新的shsd模型：
[0077][0078]
其中，ε
*
为岩石复介电常数；φ为岩石孔隙度；s
w
为含水饱和度；p
d
为泥质频散项体积含量；mn为水相曲折度指数；为复合导电相；ε
rm
为岩石骨架和油气的混合介电常数；γ为伽马函数。
[0079]
在一个实施方式中，西门杜方程的表达式为：
[0080][0081]
其中，r
xo
为冲洗带电阻率；s
xo
为冲洗带含水饱和度；v
sh
为泥质含量；r
sh
为泥岩电阻率；r
mf
为冲洗带流体电阻率；φ为岩石孔隙度；m为胶结指数；n为饱和度指数。
[0082]
在本实施方式中，经典西门杜方程是一种泥质砂岩电导模型，采用的是深电阻率r
t
，其不考虑泥质的分布形式，具有较好的应用效果，表达式如下：
[0083][0084]
其中，r
t
为深电阻率，ω
·
m(电阻率)；v
sh
为泥质含量，v/v(体积比)；r
sh
为泥岩电阻率，ω
·
m；r
w
为地层水电阻率，ω
·
m；a为阿尔奇公式指数，通常为1；m为胶结指数；n为饱和度指数。
[0085]
在本实施方式中，由于介电测井探测深度相对较浅，为了配合crim模型和shsd模型，需要利用冲洗带电阻率r
xo
代替深电阻率r
t
，这样反演计算得到的是冲洗带含水饱和度，与crim模型和shsd模型反演计算得到的冲洗带含水饱和度一致。因此，采用探测深度较浅
的电阻率曲线与其匹配，得到的表达式为：
[0086][0087]
其中，r
xo
为冲洗带电阻率，ω
·
m；s
xo
为冲洗带含水饱和度；v
sh
为泥质含量；r
sh
为泥岩电阻率；r
mf
为冲洗带流体电阻率，ω
·
m，与公式(7)中电导率互为倒数；φ为岩石孔隙度；m为胶结指数；n为饱和度指数。
[0088]
在一个实施方式中，按照以下公式，利用多种响应方程构建最小平方目标函数：
[0089][0090]
其中，minf为最小平方目标函数值；j为变量，j＝1,2,3,4分别表示四种不同的频率；ε
logj
为介电测井实际测量的相对介电常数；ε
calj
为利用响应方程计算得到的相对介电常数的理论值；σ
logj
为介电测井实际测量的电导率值；σ
calj
为利用响应方程计算得到的电导率的理论值；r
log
为电阻率测井实测值；r
cal
为利用响应方程计算得到的电阻率的理论值。
[0091]
在本实施方式中，ε
calj
、σ
calj
和r
cal
表示在对含水饱和度、水相曲折度指数、冲洗带矿化度以及泥质频散项体积含量这四个参数进行寻优过程中利用响应方程(3)、(12)、(14)计算的理论测井值。具体的，ε
calj
为反演参数带入公式(3)、(12)，基于四种频率计算得到的相对介电常数；σ
calj
为反演参数带入公式(3)、(12)，基于四种频率计算得到的电导率值；r
cal
为反演参数带入公式(14)，计算得到的冲洗带电阻率值。
[0092]
在一个实施例中，在粒子群算法的每次迭代中，均会得到四个反演参数。将冲洗带含水饱和度和冲洗带矿化度代入公式(3)可计算出960mhz测量的相对介电常数和电导率，即ε
calj
和σ
calj
；将冲洗带含水饱和度、冲洗带矿化度、水相曲折度指数、和泥质频散项体积含量代入公式(12)可计算出其它三个频段下的相对介电常数和电导率，即ε
calj
和σ
calj
；将冲洗带含水饱和度和水相曲折度指数(mn＝m＝n)代入公式(14)可以计算出冲洗带电阻率r
xo
，即r
cal
。
[0093]
在本实施例中，根据计算出的ε
calj
、σ
calj
和r
cal
分别与其测井实测值进行最小平方计算(公式15)，可以得到最小平方目标函数的值。其中，上述ε
logj
、σ
logj
和r
log
可以是根据测井资料集中的各个测井曲线确定的实测值。
[0094]
在本实施方式中，在利用多种响应方程构建最小平方目标函数之前，还可以包括：基于所述目标井的测井资料集选择多个响应方程；其中，阵列感应测井电阻率曲线对应西门杜方程，1ghz测量的介电常数曲线和电导率曲线对应复折射率模型，20mhz至500mhz频段的介电常数曲线和电导率曲线对应泥质砂岩响应模型。
[0095]
在本实施方式中，由于不同的测井曲线适用的响应方程不同，因此，可以基于所述目标井的测井资料集选择多个响应方程。响应方程和测井曲线之间的对应关系可以包括：阵列感应测井电阻率曲线对应西门杜方程，1ghz测量的介电常数曲线和电导率曲线对应复折射率模型，20mhz至500mhz频段的介电常数曲线和电导率曲线对应泥质砂岩响应模型等。当然，响应方程和测井曲线之间的对应关系不限于上述举例，所属领域技术人员在本说明书实施例技术精髓的启示下，还可能做出其它变更，但只要其实现的功能和效果与本说明书实施例相同或相似，均应涵盖于本说明书实施例保护范围内。
[0096]
在本实施方式中，20mhz至500mhz之间的频段，优选的可以为：22mhz、100mhz和350mhz频段等，具体的可以根据实际情况确定，本说明书实施例对此不作限定。
[0097]
在一个实施方式中，基于所述目标井的测井资料集，利用粒子群优化算法对所述目标函数进行寻优反演，得到所述目标井的含水饱和度曲线，可以包括：将含水饱和度、水相曲折度指数、冲洗带矿化度、泥质频散项体积含量作为反演参数，并获取所述目标井的特征信息集；其中，所述特征信息集中包含所述目标井的孔隙度曲线、油气的相对介电常数、地层温度和岩石骨架相对介电常数。进一步的，可以根据所述目标井的测井资料集和特征信息集，利用粒子群优化算法对所述目标函数进行寻优反演，得到各个反演参数的最优值，从而可以基于所述各个反演参数的最优值确定所述目标井的含水饱和度曲线。
[0098]
在本实施方式中，在进行寻优反演之前可以获取目标井的孔隙度曲线、油气的相对介电常数、地层温度和岩石骨架相对介电常数，这些特征参数可以用于公式(3)、(12)、(14)的计算。其中，孔隙度曲线利用常规测井测量得到的曲线基于经验公式计算得到的；油气的相对介电常数、地层温度和岩石骨架相对介电常数均为经验值，在一些实施例中可以分别为：2.2、160华氏度、5.0。
[0099]
在本实施方式中，利用粒子群优化算法对所述目标函数进行寻优反演的步骤可以包括：
[0100]
①
将含水饱和度s
w
、水相曲折度指数mn、冲洗带矿化度sal以及泥质频散项体积含量p
d
四个反演参数设置为粒子，随机初始化位置和速度，设置迭代次数、粒子个数、迭代终止条件等参数；
[0101]
②
根据粒子的初始位置，根据四组不同频率的相对介电常数和电导率、阵列感应测井电阻率曲线和特征信息集，利用公式(15)计算每个粒子的适应度值；
[0102]
③
更新每个粒子个体在寻优过程中的和gbest
r
；
[0103]
④
根据速度迭代公式(1)赋予粒子新的速度；
[0104]
⑤
根据位置迭代公式(2)更新粒子的位置；
[0105]
⑥
重复步骤
②
至步骤
⑤
，直至满足迭代终止条件(适应度小于预设阈值)，输出四个反演参数的最优值。由于测井曲线是连续的点构成的曲线，在对测井曲线进行反演计算时，得到一系列的最优点构成的最优曲线，从而可以得到含水饱和度曲线。
[0106]
在本实施方式中，可以假设mn＝m＝n，m和n分别为胶结指数和饱和度指数。由于在含水纯砂岩或碳酸盐岩中，mn相当于archie公式中的m，因此，为了方便计算，可以假设mn＝m＝n。
[0107]
在本实施方式中，输入参数除了包含四组不同频率的相对介电常数和电导率、阵列感应测井电阻率曲线外，还可以包含孔隙度曲线、油气的相对介电常数、地层温度和岩石骨架相对介电常数等特征参数。
[0108]
在本实施方式中，寻优反演得到的水相曲折度指数mn、冲洗带矿化度sal以及泥质频散项体积含量p
d
这三个参数，均具有一定的岩石物理含义，可用于辅助评估储层，水相曲折度指数表征岩石孔隙结构，与archie公式中的胶结指数m有关；冲洗带矿化度表征地层水矿化度的分布；泥质频散项体积含量与储层粘土含量有关。这三个参数是模型中的参数，与含水饱和度是同时输出的。
[0109]
在一个实施方式中，根据所述目标井的含水饱和度曲线，确定所述目标井的含油
饱和度曲线，可以包括：将1与所述目标井的含水饱和度曲线的差值作为所述目标井的含油饱和度曲线。
[0110]
在本实施方式中，由于同一深度地层的含水饱和度与含油饱和度之和为1，因此，可以将1与所述目标井的含水饱和度曲线的差值作为所述目标井的含油饱和度曲线。
[0111]
在一个实施例中，可以将反演得到的最优参数的理论介电常数和电导率测井曲线响应值与实际测井作误差分析。电阻率测井曲线的相对误差定义为：
[0112][0113]
其中，err
res
为电阻率测井曲线的相对误差；r
xo
为冲洗带电阻率；meas和sys分别表示实际测量的测井响应和理论计算的测井响应。
[0114]
介电测井每种频率的平均相对误差定义为：
[0115][0116]
其中，err
dielec
为介电测井每种频率的平均相对误差；perm为相对介电常数；cond为电导率；meas和sys分别表示实际测量的测井响应和理论计算的测井响应。
[0117]
在本实施例中，图2为l1井中目标层段理论测井响应和实际测井响应的对比示意图。图2中第一道为自然伽马测井曲线(gr)和井径曲线(cal)，第2道至第9道显示理论和现场测量的介电常数和电导率测井响应的对比。其中，epsi_xo_f0、epsi_xo_f1、epsi_xo_f2和epsi_xo_f3分别代表介电测井仪器在22mhz、100mhz、350mhz和960mhz频率测量的相对介电常数；epsi0、epsi1、epsi2和epsi3是合成的对应频率的理论介电常数。类似地，cond_xo_f0、cond_xo_f1、cond_xo_f2和cond_xo_f3分别代表在22mhz、100mhz、350mhz和960mhz频率测量的电导率；cond0、cond1、cond2和cond3是合成的相应的电导率测井响应。第10道为阵列感应测井af10(表示探测深度10in)和合成的理论电阻率曲线rxos的对比，第11道和第12道中err_f0、err_f1、err_f2和err_f3分别是在22mhz、100mhz、350mhz和960mhz频率下正演合成测井值和实际测井值的平均相对误差，err_rxo表示计算的和实际测量的冲洗带电阻率相对误差。从第2
‑
10道可以看出，合成的测井曲线和实际测量的曲线吻合良好，二者数值接近，变化趋势也相同。最后两道中，定量计算的相对误差仅xx55m处较大，该深度段可能的原因是泥质含量较重，孔隙度较小(如图3所示)。其它深度段相对误差均较小，大多数数值小于10％，而少部分高于10％的数值也基本小于20％，这验证了本说明书实施方式中提出的致密砂岩含油饱和度介电测井和电阻率曲线处理方法的准确性和可靠性。
[0118]
在本实施例中，图3为l1井的常规测井曲线、介电测井曲线和储层参数计算结果的示意图。图3中第1道包括为自然伽马测井曲线(gr)、自然电位(sp)和井径曲线(cal)。第2道为阵列感应电阻率测井曲线af90、af60、af30、af10，其中，af90是探测深度最深的电阻率，af10是探测深度最浅的电阻率。第3道为孔隙度测井曲线，包括密度(den)、声波(ac)和中子(cnl)测井曲线。第4道中cond_xo_f0、cond_xo_f1、cond_xo_f2和cond_xo_f3分别代表介电测井仪器在22mhz、100mhz、350mhz和960mhz频率测量的电导率，第5道中epsi_xo_f0、epsi_xo_f1、epsi_xo_f2和epsi_xo_f3分别代表介电测井仪器在22mhz、100mhz、350mhz和960mhz频率测量的相对介电常数。第6道为利用常规经验公式计算出的泥质含量(sh)和孔隙度(por)，sand为砂岩含量。计算出的孔隙度(por)和岩心分析的孔隙度(porc)如第7道所示，二者吻合较好。第8道为反演计算的含水饱和度(sxo)，以及岩心分析含水饱和度(swc)和岩
心分析含油饱和度(shc)，第9道包含反演计算得到的参数泥质频散项体积含量p
d
、水相曲折度指数mn和冲洗带矿化度sal。
[0119]
在本实施例中，图3中的目标层段试油测试为水层，岩心分析含油饱和度也表明该层含油性较差。从图3中第8道中可以看到，计算的整个目标层段冲洗带含水饱和度大多为1，只有两个小层有一定的含油性显示。第9道计算的mn基本上等于2，频散项体积含量p
d
值介于0.06～0.09之间；冲洗带矿化度分布于15～40kppm，该值并不代表地层水矿化度，而是介于地层水矿化度和泥浆滤液矿化度间，进而验证了本说明书实施例中技术方案的有效性和可靠性。
[0120]
在本实施例中，图4为l2井中目标层段常规测井、介电测井及其反演结果的示意图，图4中所有曲线道与图3中对应，重复之处不再赘述。其中，l1井中目标层段试油测试为水层，l2井中目标层段试油测试为油层，图4中第8道岩心分析的含油饱和度也表明其含油性较好。图4中第2道显示原状地层电阻率(af90)高于冲洗带电阻率，因此，实际原状地层含油饱和度应高于该值，说明目标层段为油层。而反演计算结果显示，冲洗带含水饱和度除xx90m外，介于60％～70％之间，显示较好的含油性；计算的p
d
和冲洗带矿化度与l1井水层的计算值相近；mn值介于2～2.4之间，略大于l1井水层mn值，说明油层中水的导电路径比水层的复杂，这也与实际情况相符合。
[0121]
基于同一发明构思，本说明书实施例中还提供了一种含油饱和度的确定装置，如下面的实施例所述。由于含油饱和度的确定装置解决问题的原理与含油饱和度的确定方法相似，因此含油饱和度的确定装置的实施可以参见含油饱和度的确定方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图5是本说明书实施例的含油饱和度的确定装置的一种结构框图，如图5所示，可以包括：获取模块501、构建模块502、反演模块503、确定模块504，下面对该结构进行说明。
[0122]
获取模块501，可以用于获取目标井的测井资料集；其中，所述测井资料集中包括：阵列感应测井电阻率曲线、介电测井资料中四种频率的介电常数曲线和电导率曲线；
[0123]
构建模块502，可以用于利用多种响应方程构建最小平方目标函数；
[0124]
反演模块503，可以用于基于所述目标井的测井资料集，利用粒子群优化算法对所述目标函数进行寻优反演，得到所述目标井的含水饱和度曲线；
[0125]
确定模块504，可以用于根据所述目标井的含水饱和度曲线，确定所述目标井的含油饱和度曲线。
[0126]
本说明书实施例实施方式还提供了一种电子设备，具体可以参阅图6所示的基于本说明书实施例提供的含油饱和度的确定方法的电子设备组成结构示意图，所述电子设备具体可以包括输入设备61、处理器62、存储器63。其中，所述输入设备61具体可以用于输入目标井的测井资料集。所述处理器62具体可以用于获取目标井的测井资料集；其中，所述测井资料集中包括：阵列感应测井电阻率曲线、介电测井资料中四种频率的介电常数曲线和电导率曲线；利用多种响应方程构建最小平方目标函数；基于所述目标井的测井资料集，利用粒子群优化算法对所述目标函数进行寻优反演，得到所述目标井的含水饱和度曲线；根据所述目标井的含水饱和度曲线，确定所述目标井的含油饱和度曲线。所述存储器63具体
可以用于存储含油饱和度曲线等数据。
[0127]
在本实施方式中，所述输入设备具体可以是用户和计算机系统之间进行信息交换的主要装置之一。所述输入设备可以包括键盘、鼠标、摄像头、扫描仪、光笔、手写输入板、语音输入装置等；输入设备用于把原始数据和处理这些数的程序输入到计算机中。所述输入设备还可以获取接收其他模块、单元、设备传输过来的数据。所述处理器可以按任何适当的方式实现。例如，处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(application specific integrated circuit，asic)、可编程逻辑控制器和嵌入微控制器的形式等等。所述存储器具体可以是现代信息技术中用于保存信息的记忆设备。所述存储器可以包括多个层次，在数字系统中，只要能保存二进制数据的都可以是存储器；在集成电路中，一个没有实物形式的具有存储功能的电路也叫存储器，如ram、fifo等；在系统中，具有实物形式的存储设备也叫存储器，如内存条、tf卡等。
[0128]
在本实施方式中，该电子设备具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
[0129]
本说明书实施例实施方式中还提供了一种基于含油饱和度的确定方法的计算机存储介质，所述计算机存储介质存储有计算机程序指令，在所述计算机程序指令被执行时可以实现：获取目标井的测井资料集；其中，所述测井资料集中包括：阵列感应测井电阻率曲线、介电测井资料中四种频率的介电常数曲线和电导率曲线；利用多种响应方程构建最小平方目标函数；基于所述目标井的测井资料集，利用粒子群优化算法对所述目标函数进行寻优反演，得到所述目标井的含水饱和度曲线；根据所述目标井的含水饱和度曲线，确定所述目标井的含油饱和度曲线。
[0130]
在本实施方式中，上述存储介质包括但不限于随机存取存储器(random access memory，ram)、只读存储器(read
‑
only memory，rom)、缓存(cache)、硬盘(hard disk drive，hdd)或者存储卡(memory card)。所述存储器可以用于存储计算机程序指令。网络通信单元可以是依照通信协议规定的标准设置的，用于进行网络连接通信的接口。
[0131]
在本实施方式中，该计算机存储介质存储的程序指令具体实现的功能和效果，可以与其它实施方式对照解释，在此不再赘述。
[0132]
显然，本领域的技术人员应该明白，上述的本说明书实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本说明书实施例不限制于任何特定的硬件和软件结合。
[0133]
虽然本说明书实施例提供了如上述实施例或流程图所述的方法操作步骤，但基于常规或者无需创造性的劳动在所述方法中可以包括更多或者更少的操作步骤。在逻辑性上不存在必要因果关系的步骤中，这些步骤的执行顺序不限于本说明书实施例提供的执行顺序。所述的方法的在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
[0134]
应该理解，以上描述是为了进行图示说明而不是为了进行限制。通过阅读上述描述，在所提供的示例之外的许多实施方式和许多应用对本领域技术人员来说都将是显而易见的。因此，本说明书实施例的范围不应该参照上述描述来确定，而是应该参照前述权利要求以及这些权利要求所拥有的等价物的全部范围来确定。
[0135]
以上所述仅为本说明书实施例的优选实施例而已，并不用于限制本说明书实施例，对于本领域的技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的保护范围之内。