CN114165226B - 岩电参数的确定方法、确定装置与计算机可读存储介质 - Google Patents
岩电参数的确定方法、确定装置与计算机可读存储介质 Download PDFInfo
- Publication number
- CN114165226B CN114165226B CN202010845683.6A CN202010845683A CN114165226B CN 114165226 B CN114165226 B CN 114165226B CN 202010845683 A CN202010845683 A CN 202010845683A CN 114165226 B CN114165226 B CN 114165226B
- Authority
- CN
- China
- Prior art keywords
- experimental data
- equivalent
- determining
- rock
- corrected
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 207
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000003860 storage Methods 0.000 title claims abstract description 24
- 230000033558 biomineral tissue development Effects 0.000 claims abstract description 181
- 239000008398 formation water Substances 0.000 claims abstract description 124
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 120
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 79
- 238000012937 correction Methods 0.000 claims abstract description 46
- 230000005611 electricity Effects 0.000 claims abstract description 14
- 238000002474 experimental method Methods 0.000 claims description 39
- 230000000007 visual effect Effects 0.000 claims description 38
- 229920006395 saturated elastomer Polymers 0.000 claims description 33
- 238000011426 transformation method Methods 0.000 claims description 12
- 230000001502 supplementing effect Effects 0.000 abstract description 5
- 238000005755 formation reaction Methods 0.000 description 66
- 239000000243 solution Substances 0.000 description 51
- 238000010586 diagram Methods 0.000 description 14
- 238000004458 analytical method Methods 0.000 description 13
- 238000009826 distribution Methods 0.000 description 13
- 239000011148 porous material Substances 0.000 description 13
- 238000004590 computer program Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000011156 evaluation Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000000295 complement effect Effects 0.000 description 4
- 238000012216 screening Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 239000007864 aqueous solution Substances 0.000 description 2
- 238000003556 assay Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000012797 qualification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000013589 supplement Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 229930195734 saturated hydrocarbon Natural products 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
-
- 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
- E21B47/00—Survey of boreholes or wells
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)
- Geophysics (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
本申请提供了一种岩电参数的确定方法、确定装置与计算机可读存储介质。该确定方法包括：获取第一实验数据，第一实验数据为非真实地层水矿化度下的实验数据，第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；对部分第一实验数据进行校正，得到校正实验数据，校正实验数据包括校正地层因素和校正电阻增大率；至少根据第一实验数据和校正实验数据确定真实地层水矿化度下的岩电参数。该方案在无真实地层水矿化度下实验数据的情况下，仍然实现了对真实地层水矿化度下的岩电参数的确定，无需补选真实矿化度下的实验数据，即无需增加额外的工作量，简化了岩电参数的确定流程，缩短了岩电参数的确定周期。
Description
技术领域
本申请涉及石油测井技术领域，具体而言，涉及一种岩电参数的确定方法、确定装置、计算机可读存储介质与处理器。
背景技术
储层含油饱和度是储层流体性质判别的重要指标，也是油气藏研究中评价储层的关键特征参数，储层含油饱和度的测井评价依据是阿尔奇公式，计算含油饱和度要用到岩电参数(a，m，b，n；其中m为孔隙度结构指数，与储层的孔喉结构及孔隙迂曲程度有关，孔隙结构越复杂，m值越大。理想条件下，当储层的骨架与孔隙空间完全独立时，m＝1；孔隙迂曲程度高时，则m>1，反之孔喉连通程度较好，则m<1，一般m值变化范围为1.5～3；n为饱和度指数，其值为孔隙结构以及油水在孔隙空间中微观分布特征共同作用的结果，一般情况下n值接近2；a，b为岩性相关的系数，a值变化范围为0.6～1.5，b值接近1，分别与m、n值为互补关系)，岩电参数确定的准确度就影响着含油饱和度和储量计算的精度。因此，研究清岩电参数的影响因素及变化规律对准确计算含油饱和度至关重要。实验及分析表明地层水矿化度对于岩电参数存在以下影响规律：①岩石导电能力主要与岩石中孔隙流体带电离子的浓度有关，即与孔隙流体的矿化度有关。随着矿化度增加，地层可移动的导电离子数量增加，使得R0、Rt、Rw都减小。②对于纯砂岩储层，由于只有孔隙内地层水导电，当地层水矿化度增加时，岩石100％含水电阻率R0与地层水电阻率Rw同步降低，地层因素F基本不变；而对于含泥质的储层，当地层水矿化度较低(一般低于上万mg/L)时，岩石中泥质附加导电性将与孔隙水共同作用，增强岩石导电性，地层因素F降低；当地层水矿化度较高时，岩石附加导电作用基本可忽略，地层因素F与孔隙度有较好的相关性。③由于烃类的存在，电流在烃类与孔隙水共存的岩石中流动比只有岩石中只有孔隙水时更加曲折，当矿化度增加时，R0降低速度必然比部分饱和烃类、部分饱和地层水的岩石电阻率Rt减小要快，则电阻增大率I增加。因此，随着矿化度增大而增大，地层因素F、电阻增大率I曲线越来越陡，斜率m、n值增大。
实践中由于勘探初期资料少，岩电实验矿化度参考邻区或邻井试油证实的地层水矿化度值，一般与真实矿化度有明显差异，相应的实验数据与真实矿化度下的实验点之间存在规律性的偏差。与本发明相关的现有技术已知真实的地层水矿化度后，采用筛选法及补选实验法来甄选或者补充实验数据：即在处理实验数据时一般筛选出与真实矿化度接近的实验点确定岩电参数，或者在确定真实矿化度后补选实验样品来补充数据点。以上方法在数据样本充足的情况下能够达到测井解释和储层评价的精度，但均舍去了其他矿化度下实验数据点，当数据样本较少或岩心客观条件不满足补选实验时，则不能更加准确地评价储层，并且补选实验法涉及额外工作量，会产生额外费用、岩电参数确定周期延长。
在背景技术部分中公开的以上信息只是用来加强对本文所描述技术的背景技术的理解，因此，背景技术中可能包含某些信息，这些信息对于本领域技术人员来说并未形成在本国已知的现有技术。
发明内容
本申请的主要目的在于提供一种岩电参数的确定方法、确定装置、计算机可读存储介质与处理器，以解决现有技术中受地层水矿化度影响的岩电参数的确定方法的确定效率较低的问题。
为了实现上述目的，根据本申请的一个方面，提供了一种岩电参数的确定方法，包括：获取第一实验数据，所述第一实验数据为非真实地层水矿化度下的实验数据，所述第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；对部分所述第一实验数据进行校正，得到校正实验数据，所述校正实验数据包括校正地层因素和校正电阻增大率；至少根据所述第一实验数据和所述校正实验数据确定所述真实地层水矿化度下的岩电参数。
可选地，至少根据所述第一实验数据和所述校正实验数据确定所述真实地层水矿化度下的岩电参数，包括：获取第二实验数据，所述第二实验数据为真实地层水矿化度下的实验数据，所述第二实验数据包括第二地层因素、第二有效孔隙度、第二电阻增大率和第二含水饱和度；根据所述第二实验数据、所述第一实验数据和所述校正实验数据确定所述真实地层水矿化度下的所述岩电参数。
可选地，获取第一实验数据，包括：选取多个岩心样品；根据所述非真实地层水矿化度配置第一溶液；将多个所述岩心样品中的部分所述岩心样品放置在所述第一溶液中进行饱和处理，得到多个第一饱和岩心样品；对所述第一饱和岩心样品进行实验，得到所述第一实验数据。
可选地，对部分所述第一实验数据进行校正，得到校正实验数据，包括：根据所述真实地层水矿化度配置第二溶液，所述第二溶液为等效真实地层水溶液；将多个所述岩心样品中的另一部分所述岩心样品放置在所述第二溶液中进行饱和处理，得到多个第二饱和岩心样品；对所述第二饱和岩心样品进行实验，得到等效实验数据，所述等效实验数据包括等效地层因素、等效有效孔隙度、等效电阻增大率和等效含水饱和度；根据所述第一实验数据确定非真实地层水矿化度下的视岩电参数，所述视岩电参数包括第一视岩性系数、第二视岩性系数、视胶结指数和视饱和度指数；根据所述等效实验数据确定等效岩电参数，所述等效岩电参数包括第一等效岩性系数、第二等效岩性系数、等效胶结指数和等效饱和度指数；根据所述视岩电参数、所述等效岩电参数、所述第一地层因素和所述第一电阻增大率确定所述校正地层因素和所述校正电阻增大率。
可选地，根据所述视岩电参数、所述等效岩电参数、所述第一地层因素和所述第一电阻增大率确定所述校正地层因素和所述校正电阻增大率，包括：采用双对数线性变换法构建第一公式，所述第一公式表示为F′＝a'·(F/a)(m′/m)，其中，F′为所述校正地层因素，a'为所述第一等效岩性系数，F为所述第一地层因素，a为所述第一视岩性系数，m′为所述等效胶结指数，m为视胶结指数；采用所述双对数线性变换法构建第二公式，所述第二公式表示为I′＝b'·(I/b)(n′/n)，其中，I′为所述校正电阻增大率，b'为所述第二等效岩性系数，I为所述第一电阻增大率，b为所述第二视岩性系数，n′为所述等效饱和度指数，n为所述视饱和度指数；根据所述第一公式和所述第二公式确定所述校正地层因素和所述校正电阻增大率。
可选地，至少根据所述第一实验数据和所述校正实验数据确定所述真实地层水矿化度下的岩电参数，包括：根据所述校正地层因素、所述校正电阻增大率、所述第一有效孔隙度和所述第一含水饱和度，确定所述真实地层水矿化度下的岩电参数，所述岩电参数包括第一岩性系数、第二岩性系数、胶结指数和饱和度指数。
可选地，根据所述校正地层因素、所述校正电阻增大率、所述第一有效孔隙度和所述第一含水饱和度，确定所述真实地层水矿化度下的岩电参数，包括：构建第三公式，所述第三公式表示为F′＝a″·φ-m″，其中，φ为所述第一有效孔隙度，a″为所述第一岩性系数，m″为所述胶结指数；构建第四公式，所述第四公式表示为I′＝b″·Sw-n″，其中，Sw为所述第一含水饱和度，b″为所述第二岩性系数，n″为所述饱和度指数；根据所述第三公式和所述第四公式，确定所述真实地层水矿化度下的岩电参数。
根据本申请的一个方面，提供了一种岩电参数的确定装置，包括：获取单元，用于获取第一实验数据，所述第一实验数据为非真实地层水矿化度下的实验数据，所述第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；校正单元，用于对部分所述第一实验数据进行校正，得到校正实验数据，所述校正实验数据包括校正地层因素和校正电阻增大率；确定单元，用于至少根据所述第一实验数据和所述校正实验数据确定所述真实地层水矿化度下的岩电参数。
根据本申请的又一个方面，提供了一种计算机可读存储介质，所述计算机可读存储介质包括存储的程序，其中，在所述程序运行时控制所述计算机可读存储介质所在设备执行任意一种所述的岩电参数的确定方法。
根据本申请的再一个方面，提供了一种处理器，所述处理器用于运行程序，其中，所述程序运行时执行任意一种所述的岩电参数的确定方法。
应用本申请的技术方案，通过获取非真实地层水矿化度下的实验数据，再对非真实地层水矿化度下的实验数据进行校正得到校正实验数据，然后利用校正实验数据和第一实验数据确定上述真实地层水矿化度下的岩电参数，在无真实地层水矿化度下实验数据的情况下，仍然实现了对真实地层水矿化度下的岩电参数的确定，无需补选真实矿化度下的实验数据，即无需增加额外的工作量，简化了岩电参数的确定流程，缩短了岩电参数的确定周期，实现了利用非真实地层水矿化度下的实验数据对真实地层水矿化度下的岩电参数的确定。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1是根据本申请实施例的岩电参数的确定方法的流程图；
图2是根据本申请实施例的岩电参数的确定装置示意图；
图3是根据本申请实施例的双对数线性变换方法推导原理示意图；
图4(a)是根据本申请实施例的第一地层因素与第一有效孔隙度的分布图；
图4(b)是根据本申请实施例的第一电阻增大率和第一含水饱和度的分布图；
图5(a)是根据本申请实施例的第一地层因素与第一有效孔隙度、第二地层因素与第二有效孔隙度的分布图；
图5(b)是根据本申请实施例的第一电阻增大率和第一含水饱和度、第二电阻增大率和第二含水饱和度的分布图；
图6(a)是根据本申请实施例的校正地层因素与第一有效孔隙度的分布图；
图6(b)是根据本申请实施例的校正电阻增大率与第一含水饱和度的分布图；
图7(a)是根据本申请实施例的一种取心分析饱和度、真实岩电参数解释饱和度与岩电参数校正后解释饱和度的分布关系图；
图7(b)是根据本申请实施例的另一种取心分析饱和度、真实岩电参数解释饱和度与岩电参数校正后解释饱和度的分布关系图。
具体实施方式
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本申请的实施例。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
应该理解的是，当元件(诸如层、膜、区域、或衬底)描述为在另一元件“上”时，该元件可直接在该另一元件上，或者也可存在中间元件。而且，在说明书以及权利要求书中，当描述有元件“连接”至另一元件时，该元件可“直接连接”至该另一元件，或者通过第三元件“连接”至该另一元件。
正如背景技术中所介绍的，现有技术中的受地层水矿化度影响的岩电参数的确定方法的确定效率较低，为解决如上受地层水矿化度影响的岩电参数的确定方法的确定效率较低的问题，本申请的实施例提供了一种岩电参数的确定方法、确定装置、计算机可读存储介质与处理器。
根据本申请的实施例，提供了一种岩电参数的确定方法。
图1是根据本申请实施例的岩电参数的确定方法的流程图。如图1所示，该方法包括以下步骤：
步骤S101，获取第一实验数据，上述第一实验数据为非真实地层水矿化度下的实验数据，上述第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；
步骤S102，对部分上述第一实验数据进行校正，得到校正实验数据，上述校正实验数据包括校正地层因素和校正电阻增大率；
步骤S103，至少根据上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的岩电参数。
上述方案中，通过获取非真实地层水矿化度下的实验数据，再对非真实地层水矿化度下的实验数据进行校正得到校正实验数据，然后利用校正实验数据和第一实验数据确定上述真实地层水矿化度下的岩电参数，在无真实地层水矿化度下实验数据的情况下，仍然实现了对真实地层水矿化度下的岩电参数的确定，无需补选真实矿化度下的实验数据，即无需增加额外的工作量，简化了岩电参数的确定流程，缩短了岩电参数的确定周期，实现了利用非真实地层水矿化度下的实验数据对真实地层水矿化度下的岩电参数的确定。
具体地，根据上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的岩电参数，可以按照国标SY/T 5385—2007《岩石电阻率参数实验室测量及计算方法》确定岩电参数，以进一步地保证岩电参数的精确确定。
具体地，上述真实地层水矿化度可以为20000mg/L，非真实地层水矿化度可以为10000mg/L、15000mg/L以及25000mg/L等，可以获取多种非真实地层水矿化度下的第一实验数据，例如，获取地层水矿化度为10000mg/L下的多个第一实验数据，获取地层水矿化度为15000mg/L下的多个第一实验数据，获取地层水矿化度为25000mg/L下的多个第一实验数据。然后，对地层水矿化度为10000mg/L下的多个第一实验数据进行校正得到校正实验数据，对地层水矿化度为15000mg/L下的多个第一实验数据进行校正得到校正实验数据，对地层水矿化度为25000mg/L下的多个第一实验数据进行校正得到校正实验数据，然后，根据所有的第一实验数据和所有的校正实验数据确定真实矿化度下的岩电参数。
需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
本申请的一种实施例中，至少根据上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的岩电参数，包括：获取第二实验数据，上述第二实验数据为真实地层水矿化度下的实验数据，上述第二实验数据包括第二地层因素、第二有效孔隙度、第二电阻增大率和第二含水饱和度；根据上述第二实验数据、上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的上述岩电参数。即在有真实地层水矿化度下的实验数据的情况下，联合第一实验数据、校正实验数据和第二实验数据确定真实地层水矿化度下的上述岩电参数，由于综合了真实地层水矿化度下的实验数据，实现了岩电参数的更精确的确定。克服了筛选法和补选实验法对实验数据点利用不充分、测井解释精度不足以及可能产生额外化验成本、确定周期长的问题。
本申请的又一种实施例，获取第一实验数据，包括：选取多个岩心样品；根据上述非真实地层水矿化度配置第一溶液；将多个上述岩心样品中的部分上述岩心样品放置在上述第一溶液中进行饱和处理，得到多个第一饱和岩心样品；对上述第一饱和岩心样品进行实验，得到上述第一实验数据。根据非真实地层水矿化度配置第一溶液，具体地，可以根据矿化度10000mg/L配置第一溶液，然后将选取的多个岩心样品放置在第一溶液中进行饱和处理，再对得到的多个第一饱和岩心样品进行实验，得到多组第一实验数据。第一溶液即为具有非真实地层水矿化度的溶液，实现了第一实验数据的准确获取。且本方案在增加少量额外选样分析工作、节约时间和费用的前提下，达到将所有的第一实验数据校正到等效真实地层水矿化度下的状态，得出的岩电参数包含更加丰富的储层信息，更加准确合理，显著提高了油水饱和度解释和储层评价的精度。
本申请的又一种实施例，对部分上述第一实验数据进行校正，得到校正实验数据，包括：根据上述真实地层水矿化度配置第二溶液，上述第二溶液为等效真实地层水溶液；将多个上述岩心样品中的另一部分上述岩心样品放置在上述第二溶液中进行饱和处理，得到多个第二饱和岩心样品；对上述第二饱和岩心样品进行实验，得到等效实验数据，上述等效实验数据包括等效地层因素、等效有效孔隙度、等效电阻增大率和等效含水饱和度；根据上述第一实验数据确定非真实地层水矿化度下的视岩电参数，上述视岩电参数包括第一视岩性系数、第二视岩性系数、视胶结指数和视饱和度指数；根据上述等效实验数据确定等效岩电参数，上述等效岩电参数包括第一等效岩性系数、第二等效岩性系数、等效胶结指数和等效饱和度指数；根据上述视岩电参数、上述等效岩电参数、上述第一地层因素和上述第一电阻增大率确定上述校正地层因素和上述校正电阻增大率。即选取一批岩心样品，选取的岩心样品具有相同或者近乎相同的参数，将选取的岩心样品中的一部分放置在第一溶液中进行饱和处理，得到第一饱和岩心样品，然后对第一饱和岩心样品进行实验得到第一实验数据，第一实验数据即为非真实地层水矿化度下的实验数据；将选取的岩心样品中的另一部分放置在第二溶液中进行饱和处理，得到第二饱和岩心样品，然后对第二饱和岩心样品进行实验得到等效实验数据，由于选取的一批岩心样品具有相同的参数，再根据等效实验数据和第一实验数据就确定了校正实验数据，实现了对非真实地层水矿化度下的实验数据的校正。
优选地，根据真实地层水矿化度配置等效真实地层水溶液，然后从已经在第一溶液中进行过饱和处理的岩心样品中选取一些代表性的岩心样品，进行清洗处理后，再将选取的代表性的岩心样品放置在等效真实地层水溶液中再次进行饱和处理，得到第二饱和岩心样品；再对第二饱和岩心样品进行实验得到等效实验数据；然后根据第一实验数据可以确定非真实地层水矿化度下的视岩电参数，根据等效实验数据可以确定等效岩电参数，再根据视岩电参数、等效岩电参数、第一地层因素和上述第一电阻增大率就可以确定校正地层因素和上述校正电阻增大率，即对非真实地层水矿化度下的第一实验数据的校正，具体是通过配置等效真实地层水溶液，然后将在第一溶液中做过实验的岩心样品放置在真实地层水溶液重新做实验，获取等效实验数据，再根据等效实验数据确定等效岩电参数，再利用等效岩电参数、视岩电参数、第一地层因素和上述第一电阻增大率确定上述校正地层因素和上述校正电阻增大率，实现对非真实地层水矿化度下的实验数据的校正。
本申请的再一种实施例，根据上述视岩电参数、上述等效岩电参数和上述第一地层因素和上述第一电阻增大率确定上述校正地层因素和上述校正电阻增大率，包括：采用双对数线性变换法构建第一公式，上述第一公式表示为F′＝a'·(F/a)(m′/m)，其中，F′为上述校正地层因素，a'为上述第一等效岩性系数，F为上述第一地层因素，a为上述第一视岩性系数，m′为上述等效胶结指数，m为视胶结指数；采用上述双对数线性变换法构建第二公式，上述第二公式表示为I′＝b'·(I/b)(n′/n)，其中，I′为上述校正电阻增大率，b'为上述第二等效岩性系数，I为上述第一电阻增大率，b为上述第二视岩性系数，n′为上述等效饱和度指数，n为上述视饱和度指数；根据上述第一公式和上述第二公式确定上述校正地层因素和上述校正电阻增大率。第一公式中的a'、F、a、m′和m，通过对岩心样品在第一溶液和等效真实地层水溶液做实验可以得到，在a'、F、a、m′和m已经求得的情况下，通过第一公式可以得到F′；同理，第二公式中的b'、I、b、n′和n，通过对岩心样品在第一溶液和等效真实地层水溶液做实验可以得到，在b'、I、b、n′和n已经求得的情况下，通过第二公式可以得到I′，实现了对校正地层因素和校正电阻增大率的精确确定。
本申请的一种具体的实施例中，至少根据上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的岩电参数，包括：根据上述校正地层因素、上述校正电阻增大率、上述第一有效孔隙度和上述第一含水饱和度，确定上述真实地层水矿化度下的岩电参数，上述岩电参数包括第一岩性系数、第二岩性系数、胶结指数和饱和度指数。在确定校正地层因素和上述校正电阻增大率之后，结合根据第一实验数据确定的第一有效孔隙度和上述第一含水饱和度，实现了对真实地层水矿化度下的岩电参数的确定。
本申请的另一种实施例，根据上述校正地层因素、上述校正电阻增大率、上述第一有效孔隙度和上述第一含水饱和度，确定上述真实地层水矿化度下的岩电参数，包括：构建第三公式，上述第三公式表示为F′＝a″·φ-m″，其中，φ为上述第一有效孔隙度，a″为上述第一岩性系数，m″为上述胶结指数；构建第四公式，上述第四公式表示为I′＝b″·Sw-n″，其中，Sw为上述第一含水饱和度，b″为上述第二岩性系数，n″为上述饱和度指数；根据上述第三公式和上述第四公式，确定上述真实地层水矿化度下的岩电参数。在得到F′和φ的前提下，通过对多组F′和φ进行拟合得到第三公式，然后根据第三公式确定a″和m″，同理，在得到I′和Sw的前提下，通过对多组I′和Sw进行拟合得到第四公式，然后根据第四公式确定b″和n″，最终实现了对a″、m″、b″和n″的确定，即实现了对真实地层水矿化度下的岩电参数的精确确定。
本申请实施例还提供了一种岩电参数的确定装置，需要说明的是，本申请实施例的岩电参数的确定装置可以用于执行本申请实施例所提供的用于岩电参数的确定方法。以下对本申请实施例提供的岩电参数的确定装置进行介绍。
图2是根据本申请实施例的岩电参数的确定装置的示意图。如图2所示，该装置包括：
获取单元10，用于获取第一实验数据，上述第一实验数据为非真实地层水矿化度下的实验数据，上述第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；
校正单元20，用于对部分上述第一实验数据进行校正，得到校正实验数据，上述校正实验数据包括校正地层因素和校正电阻增大率；
确定单元30，用于至少根据上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的岩电参数。
上述方案中，获取单元获取非真实地层水矿化度下的实验数据，校正单元对非真实地层水矿化度下的实验数据进行校正得到校正实验数据，确定单元利用校正实验数据和第一实验数据确定上述真实地层水矿化度下的岩电参数，在无真实地层水矿化度下实验数据的情况下，仍然实现了对真实地层水矿化度下的岩电参数的确定，无需补选真实矿化度下的实验数据，即无需增加额外的工作量，简化了岩电参数的确定流程，缩短了岩电参数的确定周期，实现了利用非真实地层水矿化度下的实验数据对真实地层水矿化度下的岩电参数的确定。
具体地，上述真实地层水矿化度可以为20000mg/L，非真实地层水矿化度可以为10000mg/L、15000mg/L以及25000mg/L等，可以获取多种非真实地层水矿化度下的第一实验数据，例如，获取地层水矿化度为10000mg/L下的多个第一实验数据，获取地层水矿化度为15000mg/L下的多个第一实验数据，获取地层水矿化度为25000mg/L下的多个第一实验数据。然后，对地层水矿化度为10000mg/L下的多个第一实验数据进行校正得到校正实验数据，对地层水矿化度为15000mg/L下的多个第一实验数据进行校正得到校正实验数据，对地层水矿化度为25000mg/L下的多个第一实验数据进行校正得到校正实验数据，然后，根据所有的第一实验数据和所有的校正实验数据确定真实矿化度下的岩电参数。
本申请的一种实施例中，确定单元包括获取模块和第一确定模块，获取模块用于获取第二实验数据，上述第二实验数据为真实地层水矿化度下的实验数据，上述第二实验数据包括第二地层因素、第二有效孔隙度、第二电阻增大率和第二含水饱和度；第一确定模块用于根据上述第二实验数据、上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的上述岩电参数。即在有真实地层水矿化度下的实验数据的情况下，联合第一实验数据、校正实验数据和第二实验数据确定真实地层水矿化度下的上述岩电参数，由于综合了真实地层水矿化度下的实验数据，实现了岩电参数的更精确的确定。克服了筛选法和补选实验法对实验数据点利用不充分、测井解释精度不足以及可能产生额外化验成本、确定周期长的问题。
本申请的又一种实施例，获取单元包括选取模块、第一配置模块和第一实验模块，选取模块用于选取多个岩心样品；第一配置模块用于根据上述非真实地层水矿化度配置第一溶液；第一实验模块用于将多个上述岩心样品中的部分上述岩心样品放置在上述第一溶液中进行饱和处理，得到多个第一饱和岩心样品；对上述第一饱和岩心样品进行实验，得到上述第一实验数据。根据非真实地层水矿化度配置第一溶液，具体地，可以根据矿化度10000mg/L配置第一溶液，然后将选取的多个岩心样品放置在第一溶液中饱和处理，再对得到的多个第一饱和岩心样品进行实验，得到多组第一实验数据。第一溶液即为具有非真实地层水矿化度的溶液，实现了第一实验数据的准确获取。且本方案在增加少量额外选样分析工作、节约时间和费用的前提下，达到将所有的第一实验数据校正到等效真实地层水矿化度下的状态，得出的岩电参数包含更加丰富的储层信息，更加准确合理，显著提高了油水饱和度解释和储层评价的精度。
本申请的又一种实施例，校正单元包括第二配置模块、第二实验模块、第二确定模块、第三确定模块和第四确定模块，第二配置模块用于根据上述真实地层水矿化度配置第二溶液，上述第二溶液为等效真实地层水溶液；第二实验模块由于将多个上述岩心样品中的另一部分上述岩心样品放置在上述第二溶液中进行饱和处理，得到多个第二饱和岩心样品；对上述第二饱和岩心样品进行实验，得到等效实验数据，上述等效实验数据包括等效地层因素、等效有效孔隙度、等效电阻增大率和等效含水饱和度；第二确定模块用于根据上述第一实验数据确定非真实地层水矿化度下的视岩电参数，上述视岩电参数包括第一视岩性系数、第二视岩性系数、视胶结指数和视饱和度指数；第三确定模块用于根据上述等效实验数据确定等效岩电参数，上述等效岩电参数包括第一等效岩性系数、第二等效岩性系数、等效胶结指数和等效饱和度指数；第四确定模块用于根据上述视岩电参数、上述等效岩电参数、上述第一地层因素和上述第一电阻增大率确定上述校正地层因素和上述校正电阻增大率。即选取一批岩心样品，选取的岩心样品具有相同或者近乎相同的参数，将选取的岩心样品中的一部分放置在第一溶液中进行饱和处理，得到第一饱和岩心样品，然后对第一饱和岩心样品进行实验得到第一实验数据，第一实验数据即为非真实地层水矿化度下的实验数据；将选取的岩心样品中的另一部分放置在第二溶液中进行饱和处理，得到第二饱和岩心样品，然后对第二饱和岩心样品进行实验得到等效实验数据，由于选取的一批岩心样品具有相同的参数，再根据等效实验数据和第一实验数据就确定了校正实验数据，实现了对非真实地层水矿化度下的实验数据的校正。
本申请的再一种实施例，第四确定模块还用于采用双对数线性变换法构建第一公式，上述第一公式表示为F′＝a'·(F/a)(m′/m)，其中，F′为上述校正地层因素，a'为上述第一等效岩性系数，F为上述第一地层因素，a为上述第一视岩性系数，m′为上述等效胶结指数，m为视胶结指数；采用上述双对数线性变换法构建第二公式，上述第二公式表示为I′＝b'·(I/b)(n′/n)，其中，I′为上述校正电阻增大率，b'为上述第二等效岩性系数，I为上述第一电阻增大率，b为上述第二视岩性系数，n′为上述等效饱和度指数，n为上述视饱和度指数；根据上述第一公式和上述第二公式确定上述校正地层因素和上述校正电阻增大率。第一公式中的a'、F、a、m′和m，通过对岩心样品在第一溶液和等效真实地层水溶液做实验可以得到，在a'、F、a、m′和m已经求得的情况下，通过第一公式可以得到F′；同理，第二公式中的b'、I、b、n′和n，通过对岩心样品在第一溶液和等效真实地层水溶液做实验可以得到，在b'、I、b、n′和n已经求得的情况下，通过第二公式可以得到I′，实现了对校正地层因素和校正电阻增大率的精确确定。
本申请的一种具体的实施例中，确定单元还用于根据上述校正地层因素、上述校正电阻增大率、上述第一有效孔隙度和上述第一含水饱和度，确定上述真实地层水矿化度下的岩电参数，上述岩电参数包括第一岩性系数、第二岩性系数、胶结指数和饱和度指数。在确定校正地层因素和上述校正电阻增大率之后，结合根据第一实验数据确定的第一有效孔隙度和上述第一含水饱和度，实现了对真实地层水矿化度下的岩电参数的确定。
本申请的另一种实施例，确定单元还用于构建第三公式，上述第三公式表示为F′＝a″·φ-m″，其中，φ为上述第一有效孔隙度，a″为上述第一岩性系数，m″为上述胶结指数；构建第四公式，上述第四公式表示为I′＝b″·Sw-n″，其中，Sw为上述第一含水饱和度，b″为上述第二岩性系数，n″为上述第一含水饱和度；根据上述第三公式和上述第四公式，确定上述真实地层水矿化度下的岩电参数。在得到F′和φ的前提下，通过对多组F′和φ进行拟合得到第三公式，然后根据第三公式确定a″和m″，同理，在得到I′和Sw的前提下，通过对多组I′和Sw进行拟合得到第四公式，然后根据第四公式确定b″和n″，最终实现了对a″、m″、b″和n″的确定，即实现了对真实地层水矿化度下的岩电参数的精确确定。
上述岩电参数的确定装置包括处理器和存储器，上述获取单元、校正单元和确定单元等均作为程序单元存储在存储器中，由处理器执行存储在存储器中的上述程序单元来实现相应的功能。
处理器中包含内核，由内核去存储器中调取相应的程序单元。内核可以设置一个或以上，通过调整内核参数来提高岩电参数的确定效率。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)，存储器包括至少一个存储芯片。
本发明实施例提供了一种计算机可读存储介质，上述计算机可读存储介质包括存储的程序，其中，在上述程序运行时控制上述计算机可读存储介质所在设备执行上述岩电参数的确定方法。
本发明实施例提供了一种处理器，上述处理器用于运行程序，其中，上述程序运行时执行上述岩电参数的确定方法。
本发明实施例提供了一种设备，设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序，处理器执行程序时实现至少以下步骤：
步骤S101，获取第一实验数据，上述第一实验数据为非真实地层水矿化度下的实验数据，上述第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；
步骤S102，对部分上述第一实验数据进行校正，得到校正实验数据，上述校正实验数据包括校正地层因素和校正电阻增大率；
步骤S103，至少根据上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的岩电参数。
本文中的设备可以是服务器、PC、PAD、手机等。
本申请还提供了一种计算机程序产品，当在数据处理设备上执行时，适于执行初始化有至少如下方法步骤的程序：步骤S101，获取第一实验数据，上述第一实验数据为非真实地层水矿化度下的实验数据，上述第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；
步骤S102，对部分上述第一实验数据进行校正，得到校正实验数据，上述校正实验数据包括校正地层因素和校正电阻增大率；
步骤S103，至少根据上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的岩电参数。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
实施例1
本实施例涉及一种具体的岩电参数的确定方法，具体包括如下步骤：
步骤1：选取多个岩心样品，根据上述非真实地层水矿化度配置第一溶液，将多个上述岩心样品放置在上述第一溶液中进行实验，得到上述第一实验数据，上述第一实验数据为非真实地层水矿化度下的实验数据，上述第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；
步骤2：根据上述第一实验数据确定非真实地层水矿化度下的视岩电参数，上述视岩电参数包括第一视岩性系数、第二视岩性系数、视胶结指数和视饱和度指数，具体的实现方式是获取多组第一地层因素和第一有效孔隙度，然后对多组第一地层因素和第一有效孔隙度进行拟合得到第一地层因素和第一有效孔隙度的关系式，然后根据第一地层因素和第一有效孔隙度的关系式确定第一视岩性系数和视胶结指数；获取多组第一电阻增大率和第一含水饱和度，然后对多组第一电阻增大率和第一含水饱和度进行拟合得到第一电阻增大率和第一含水饱和度的关系式，然后根据第一电阻增大率和第一含水饱和度的关系式确定第一电阻增大率和第一含水饱和度；
步骤3：根据上述真实地层水矿化度配置第二溶液，上述第二溶液为等效真实地层水溶液；从多个上述岩心样品中选取部分上述岩心样品；将选取的部分上述岩心样品放置在上述第二溶液中进行实验，得到等效实验数据，上述等效实验数据包括等效地层因素、等效有效孔隙度、等效电阻增大率和等效含水饱和度；根据上述等效实验数据确定等效岩电参数，上述等效岩电参数包括第一等效岩性系数、第二等效岩性系数、等效胶结指数和等效饱和度指数；确定第一等效岩性系数、第二等效岩性系数、等效胶结指数和等效饱和度指数的方法与确定第一视岩性系数、第二视岩性系数、视胶结指数和视饱和度指数的方法相同；
步骤4：根据上述视岩电参数、上述等效岩电参数、上述第一地层因素和上述第一电阻增大率确定上述校正地层因素和上述校正电阻增大率；
步骤4的具体的实现方式为：采用双对数线性变换法构建第一公式，上述第一公式表示为F′＝a'·(F/a)(m′/m)，其中，F′为上述校正地层因素，a'为上述第一等效岩性系数，F为上述第一地层因素，a为上述第一视岩性系数，m′为上述等效胶结指数，m为视胶结指数，采用上述双对数线性变换法构建第二公式，上述第二公式表示为I′＝b'·(I/b)(n′/n)，其中，I′为上述校正电阻增大率，b'为上述第二等效岩性系数，I为上述第一电阻增大率，b为上述第二视岩性系数，n′为上述等效饱和度指数，n为上述视饱和度指数；根据上述第一公式和上述第二公式确定上述校正地层因素和上述校正电阻增大率。
采用双对数线性变换法构建第一公式和第二公式的具体的方法是：
根据阿尔奇公式，第一地层因素F与第一有效孔隙度在双对数坐标中的关系如公式1所示，第一电阻增大率I与第一含水饱和度Sw的关系如公式2所示：
lgI＝lgb-n·lgsw (公式2)
如图3所示，假设双对数坐标中已知一条直线L1，斜率为b1，截距为a1；按顺时针旋转一个角度θ变为L2，斜率为b2，截距为a2，由此可知θ＝arctan(b2)-arctan(b1)。L1上的一点坐标为(X1，Y1)，旋转后在L2坐标为(X1，Y2)，两条直线的交点为(X0，Y0)，则可以求出(X0，Y0)，Y2的表达式分别为：
并且已知X1＝X2＝lgx1，Y1＝lgy1，Y2＝lgy’(x1，y1，y’分别为X1或X2，Y1，Y2算术坐标系中的对应值)则将双对数关系变换为算术坐标关系，可推导出：
即：
旋转后的假设X2＝X1是基于岩电实验的地质意义考虑，如图3所示，X轴代表岩样的孔隙度或含水饱和度，为固定值，Y轴代表实验测得的地层因素或电阻增大率值。校正过程相当于对样品的一次虚拟实验，实验结果受矿化度的影响发生变化，但样品应始终为同一块，因此假设X2＝X1。若直线L1逆时针旋转成L2，变换公式不变。
应用到阿尔奇公式中则可得到公式7和公式8，公式7即为第一公式，公式8即为第二公式：
F′＝a'·(F/a)(m′/m) (公式7)
I′＝b'·(I/b)(n′/n) (公式8)
其中：
a、m、b、n、F、I为各组非真实地层水矿化度下的实验点确定的视岩电参数及测得的第一地层因素及第一电阻增大率；
a'、m′、b'、n′为各组非真实地层水矿化度下选取有代表性的岩心样品后在等效真实地层水矿化度下补做岩电实验后确定的等效岩电参数；
F′为上述校正地层因素，I′为上述校正电阻增大率。
步骤5：根据上述校正地层因素、上述校正电阻增大率、上述第一有效孔隙度和上述第一含水饱和度，确定上述真实地层水矿化度下的岩电参数，上述岩电参数包括第一岩性系数、第二岩性系数、胶结指数和饱和度指数，具体的实现方式是：构建第三公式，上述第三公式表示为F′＝a″·φ-m″，其中，φ为上述第一有效孔隙度，a″为上述第一岩性系数，m″为上述胶结指数；构建第四公式，上述第四公式表示为I′＝b″·Sw-n″，其中，Sw为上述第一含水饱和度，b″为上述第二岩性系数，n″为上述饱和度指数；根据上述第三公式和上述第四公式，确定上述真实地层水矿化度下的岩电参数。
步骤6：根据真实地层水矿化度下的岩电参数进行油水饱和度解释(即解释饱和度)和储层评价，由于通过本方案的方法获得了较准确的岩电参数，进一步地保证了油水饱和度解释和储层评价的准确性。
实施例2
本实施例涉及一种具体的岩电参数的确定方法，具体包括如下步骤：
该实施例在实施例1的基础上联合了真实地层水矿化度下的实验数据，即第二实验数据，上述第二实验数据包括第二地层因素、第二有效孔隙度、第二电阻增大率和第二含水饱和度；
根据上述第二实验数据、上述第一实验数据和上述校正实验数据确定上述真实地层水矿化度下的上述岩电参数。
采用本方案确定的岩电参数的实验结果分析如下：
如图4(a)所示，表示第一地层因素F与第一有效孔隙度φ的分布图；
如图4(b)所示，表示第一电阻增大率I和第一含水饱和度Sw的分布图；
由图4(a)和图4(b)可知，第一地层因素和第一有效孔隙度数据点较分散，且还有个别的第一电阻增大率与其他的数据点有显著的偏差，由多个第一地层因素与多个第一有效孔隙度拟合得到的关系式的相关系数较低；
如图5(a)所示，表示第一地层因素F与第一有效孔隙度φ、第二地层因素F1与第二有效孔隙度φ1的分布图；
如图5(b)所示，表示第一电阻增大率I和第一含水饱和度Sw、第二电阻增大率I1和第二含水饱和度Sw1的分布图；
图5(a)中的黑色圆点表示第一地层因素和第一有效孔隙度分布的点，三角形表示第二地层因素和第二有效孔隙度分布的点；
图5(b)中黑色圆点表示第一电阻增大率和第一含水饱和度分布的点，三角形表示第二电阻增大率和第二含水饱和度分布的点；
由图4(a)和图4(b)可知，真实地层水矿化度下的实验数据只占总数据点很小的比例，由第二实验数据得到的岩电参数较符合阿尔奇公式的要求，相关系数提高；
如图6(a)所示，为根据第二实验数据、上述第一实验数据和上述校正实验数据确定的校正地层因素F′与第一有效孔隙度φ的分布图；
如图6(b)所示，为根据第二实验数据、上述第一实验数据和上述校正实验数据确定的校正电阻增大率I′与第一含水饱和度Sw的分布图；
由图6(a)和图6(b)可知，所有的数据点明显收敛，相关系数更高，即联合第二实验数据、上述第一实验数据和上述校正实验数据确定的岩电参数符合阿尔奇公式的要求，得到了较准确的岩电参数。
如图7(a)所示和图7(b)所示，为取心分析饱和度、根据第二实验数据确定的岩电参数的解释饱和度(即真实岩电参数解释饱和度)与利用本方案得到的岩电参数的解释饱和度(即岩电参数校正后解释饱和度)的分布关系图，其中，横坐标表示样品深度，纵坐标表示含油饱和度。
由图7(a)和图7(b)可知，根据本方案得到的岩电参数的解释饱和度与取心分析饱和度几乎重合，然而，仅仅利用第二实验数据确定的岩电参数的解释饱和度与取心分析饱和度有明显的偏差，进行点对点误差分析得到：真实岩电参数点对点误差分析合格率为8.33％，平均绝对误差为9.13％～10.55％，超出石油天然气储量估算规范(DZ/T0217—2020)中关于“大型以上油(气)田(藏)用测井解释资料确定探明储量含油(气)饱和度(％)时，应有油基泥浆取心或密闭取心分析验证，绝对误差不超过±5个百分点，特殊情况除外”的要求；校正后的点对点误差分析合格率91.67％，平均绝对误差为0.26％～2.73％，体现出本发明方法的优越性，校正后测井解释得到的含油饱和度更加准确。
从以上的描述中，可以看出，本申请上述的实施例实现了如下技术效果：
1)、本申请的岩电参数的确定方法，通过获取非真实地层水矿化度下的实验数据，再对非真实地层水矿化度下的实验数据进行校正得到校正实验数据，然后利用校正实验数据和第一实验数据确定上述真实地层水矿化度下的岩电参数，在无真实地层水矿化度下实验数据的情况下，仍然实现了对真实地层水矿化度下的岩电参数的确定，无需补选真实矿化度下的实验数据，即无需增加额外的工作量，简化了岩电参数的确定流程，缩短了岩电参数的确定周期，实现了利用非真实地层水矿化度下的实验数据对真实地层水矿化度下的岩电参数的确定。
2)、本申请的岩电参数的确定装置，获取单元获取非真实地层水矿化度下的实验数据，校正单元对非真实地层水矿化度下的实验数据进行校正得到校正实验数据，确定单元利用校正实验数据和第一实验数据确定上述真实地层水矿化度下的岩电参数，在无真实地层水矿化度下实验数据的情况下，仍然实现了对真实地层水矿化度下的岩电参数的确定，无需补选真实矿化度下的实验数据，即无需增加额外的工作量，简化了岩电参数的确定流程，缩短了岩电参数的确定周期，实现了利用非真实地层水矿化度下的实验数据对真实地层水矿化度下的岩电参数的确定。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (7)
1.一种岩电参数的确定方法，其特征在于，包括：
获取第一实验数据，所述第一实验数据为非真实地层水矿化度下的实验数据，所述第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；
对部分所述第一实验数据进行校正，得到校正实验数据，所述校正实验数据包括校正地层因素和校正电阻增大率；
至少根据所述第一实验数据和所述校正实验数据确定真实地层水矿化度下的岩电参数；
其中，获取第一实验数据，包括：
选取多个岩心样品；
根据非真实地层水矿化度配置第一溶液；
将多个所述岩心样品中的部分所述岩心样品放置在所述第一溶液中进行饱和处理，得到多个第一饱和岩心样品；
对所述第一饱和岩心样品进行实验，得到所述第一实验数据；
其中，对部分所述第一实验数据进行校正，得到校正实验数据，包括：
根据真实地层水矿化度配置第二溶液，所述第二溶液为等效真实地层水溶液；
将多个所述岩心样品中的另一部分所述岩心样品放置在所述第二溶液中进行饱和处理，得到多个第二饱和岩心样品；
对所述第二饱和岩心样品进行实验，得到等效实验数据，所述等效实验数据包括等效地层因素、等效有效孔隙度、等效电阻增大率和等效含水饱和度；
根据所述第一实验数据确定非真实地层水矿化度下的视岩电参数，所述视岩电参数包括第一视岩性系数、第二视岩性系数、视胶结指数和视饱和度指数；
根据所述等效实验数据确定等效岩电参数，所述等效岩电参数包括第一等效岩性系数、第二等效岩性系数、等效胶结指数和等效饱和度指数；
根据所述视岩电参数、所述等效岩电参数、所述第一地层因素和所述第一电阻增大率确定所述校正地层因素和所述校正电阻增大率；
至少根据所述第一实验数据和所述校正实验数据确定真实地层水矿化度下的岩电参数，包括：
获取第二实验数据，所述第二实验数据为真实地层水矿化度下的实验数据，所述第二实验数据包括第二地层因素、第二有效孔隙度、第二电阻增大率和第二含水饱和度；
根据所述第二实验数据、所述第一实验数据和所述校正实验数据确定真实地层水矿化度下的所述岩电参数。
2.根据权利要求1所述的确定方法，其特征在于，根据所述视岩电参数、所述等效岩电参数、所述第一地层因素和所述第一电阻增大率确定所述校正地层因素和所述校正电阻增大率，包括：
采用双对数线性变换法构建第一公式，所述第一公式表示为F′＝a'·(F/a)(m′/m)，其中，F′为所述校正地层因素，a'为所述第一等效岩性系数，F为所述第一地层因素，a为所述第一视岩性系数，m′为所述等效胶结指数，m为视胶结指数；
采用所述双对数线性变换法构建第二公式，所述第二公式表示为I′＝b'·(I/b)(n′/n)，其中，I′为所述校正电阻增大率，b'为所述第二等效岩性系数，I为所述第一电阻增大率，b为所述第二视岩性系数，n′为所述等效饱和度指数，n为所述视饱和度指数；
根据所述第一公式和所述第二公式确定所述校正地层因素和所述校正电阻增大率。
3.根据权利要求1所述的确定方法，其特征在于，至少根据所述第一实验数据和所述校正实验数据确定真实地层水矿化度下的岩电参数，包括：
根据所述校正地层因素、所述校正电阻增大率、所述第一有效孔隙度和所述第一含水饱和度，确定真实地层水矿化度下的岩电参数，所述岩电参数包括第一岩性系数、第二岩性系数、胶结指数和饱和度指数。
4.根据权利要求3所述的确定方法，其特征在于，根据所述校正地层因素、所述校正电阻增大率、所述第一有效孔隙度和所述第一含水饱和度，确定真实地层水矿化度下的岩电参数，包括：
构建第三公式，所述第三公式表示为F′＝a″·φ-m″，其中，φ为所述第一有效孔隙度，a″为所述第一岩性系数，m″为所述胶结指数；
构建第四公式，所述第四公式表示为I′＝b″·Sw-n″，其中，Sw为所述第一含水饱和度，b″为所述第二岩性系数，n″为所述饱和度指数；
根据所述第三公式和所述第四公式，确定真实地层水矿化度下的岩电参数。
5.一种岩电参数的确定装置，其特征在于，包括：
获取单元，用于获取第一实验数据，所述第一实验数据为非真实地层水矿化度下的实验数据，所述第一实验数据包括第一地层因素、第一有效孔隙度、第一电阻增大率和第一含水饱和度；
校正单元，用于对部分所述第一实验数据进行校正，得到校正实验数据，所述校正实验数据包括校正地层因素和校正电阻增大率；
确定单元，用于至少根据所述第一实验数据和所述校正实验数据确定真实地层水矿化度下的岩电参数；
其中，所述获取单元包括：
选取模块用于选取多个岩心样品；
第一配置模块用于根据非真实地层水矿化度配置第一溶液；
第一实验模块用于将多个所述岩心样品中的部分所述岩心样品放置在所述第一溶液中进行饱和处理，得到多个第一饱和岩心样品；对所述第一饱和岩心样品进行实验，得到所述第一实验数据；
其中，所述校正单元包括：
第二配置模块用于根据真实地层水矿化度配置第二溶液，所述第二溶液为等效真实地层水溶液；
第二实验模块用于将多个所述岩心样品中的另一部分所述岩心样品放置在所述第二溶液中进行饱和处理，得到多个第二饱和岩心样品；对所述第二饱和岩心样品进行实验，得到等效实验数据，所述等效实验数据包括等效地层因素、等效有效孔隙度、等效电阻增大率和等效含水饱和度；
第二确定模块用于根据所述第一实验数据确定非真实地层水矿化度下的视岩电参数，所述视岩电参数包括第一视岩性系数、第二视岩性系数、视胶结指数和视饱和度指数；
第三确定模块用于根据所述等效实验数据确定等效岩电参数，所述等效岩电参数包括第一等效岩性系数、第二等效岩性系数、等效胶结指数和等效饱和度指数；
第四确定模块用于根据所述视岩电参数、所述等效岩电参数、所述第一地层因素和所述第一电阻增大率确定所述校正地层因素和所述校正电阻增大率；
其中，所述确定单元包括：
获取模块，用于获取第二实验数据，所述第二实验数据为真实地层水矿化度下的实验数据，所述第二实验数据包括第二地层因素、第二有效孔隙度、第二电阻增大率和第二含水饱和度；
第一确定模块，用于根据所述第二实验数据、所述第一实验数据和所述校正实验数据确定真实地层水矿化度下的所述岩电参数。
6.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质包括存储的程序，其中，在所述程序运行时控制所述计算机可读存储介质所在设备执行权利要求1至4中任意一项所述的岩电参数的确定方法。
7.一种处理器，其特征在于，所述处理器用于运行程序，其中，所述程序运行时执行权利要求1至4中任意一项所述的岩电参数的确定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010845683.6A CN114165226B (zh) | 2020-08-20 | 2020-08-20 | 岩电参数的确定方法、确定装置与计算机可读存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010845683.6A CN114165226B (zh) | 2020-08-20 | 2020-08-20 | 岩电参数的确定方法、确定装置与计算机可读存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114165226A CN114165226A (zh) | 2022-03-11 |
| CN114165226B true CN114165226B (zh) | 2024-03-01 |
Family
ID=80475351
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010845683.6A Active CN114165226B (zh) | 2020-08-20 | 2020-08-20 | 岩电参数的确定方法、确定装置与计算机可读存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114165226B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4752882A (en) * | 1986-05-05 | 1988-06-21 | Mobil Oil Corporation | Method for determining the effective water saturation in a low-resistivity hydrocarbon-bearing rock formation based upon rock matrix conductance |
| CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
| CN108802073A (zh) * | 2017-05-05 | 2018-11-13 | 中国石油化工股份有限公司 | 一种基于数字岩心的岩石电性参数获取方法及装置 |
| CN109117505A (zh) * | 2018-07-10 | 2019-01-01 | 中国海洋石油集团有限公司 | 一种基于介电实验的孔隙结构储层冲洗带含水饱和度计算方法 |
-
2020
- 2020-08-20 CN CN202010845683.6A patent/CN114165226B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4752882A (en) * | 1986-05-05 | 1988-06-21 | Mobil Oil Corporation | Method for determining the effective water saturation in a low-resistivity hydrocarbon-bearing rock formation based upon rock matrix conductance |
| CN102434152A (zh) * | 2011-12-05 | 2012-05-02 | 中国石油天然气股份有限公司 | 一种储层含油饱和度的计算方法 |
| CN108802073A (zh) * | 2017-05-05 | 2018-11-13 | 中国石油化工股份有限公司 | 一种基于数字岩心的岩石电性参数获取方法及装置 |
| CN109117505A (zh) * | 2018-07-10 | 2019-01-01 | 中国海洋石油集团有限公司 | 一种基于介电实验的孔隙结构储层冲洗带含水饱和度计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114165226A (zh) | 2022-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Soulaine et al. | The impact of sub-resolution porosity of X-ray microtomography images on the permeability | |
| Bauer et al. | Improving the estimations of petrophysical transport behavior of carbonate rocks using a dual pore network approach combined with computed microtomography | |
| Glover et al. | Grain-size to effective pore-size transformation derived from electrokinetic theory | |
| Glover et al. | Streaming-potential coefficient of reservoir rock: A theoretical model | |
| Dakhelpour-Ghoveifel et al. | Capillary-based method for rock typing in transition zone of carbonate reservoirs | |
| Almroth-Rosell et al. | A new approach to model oxygen dependent benthic phosphate fluxes in the Baltic Sea | |
| CN111475920B (zh) | 一种深水盆地古水深的获取方法、系统、电子设备及存储介质 | |
| Tagliabue et al. | A ventilation‐based framework to explain the regeneration‐scavenging balance of iron in the ocean | |
| Iwasaki et al. | Micro‐CT Scanning of Tests of Three Planktic Foraminiferal Species to Clarify DissolutionProcess and Progress | |
| Xue et al. | Temporal variation and stoichiometric ratios of organic matter remineralization in bottom waters of the northern G ulf of M exico during late spring and summer | |
| Wang et al. | Anchoring multi‐scale models to micron‐scale imaging of multiphase flow in rocks | |
| CN109630090B (zh) | 储盖组合测井评价方法和装置 | |
| CN108681793B (zh) | 深层油藏采油指数预测方法及装置 | |
| CN114165226B (zh) | 岩电参数的确定方法、确定装置与计算机可读存储介质 | |
| CN110160927A (zh) | 一种确定束缚水饱和度的方法及装置 | |
| CN116930023A (zh) | 致密砂岩相控分类孔隙度测井精细解释方法和装置 | |
| WO2021208578A1 (zh) | 测井参数的获取方法、装置和介质 | |
| CN112698401A (zh) | 一种岩石物理模板构建方法、装置、设备及存储介质 | |
| CN112231881A (zh) | 基质渗透率计算方法及系统 | |
| CN111485874A (zh) | 确定水平井储层含水饱和度的方法及装置 | |
| CN112782778A (zh) | 碳酸盐岩台地咸化环境烃源岩发育确定方法及装置 | |
| CN113818870B (zh) | 储层流体性质识别方法及系统 | |
| Zhang et al. | Hydrodynamics, Diagenesis and Hypoxia Variably Drive Benthic Oxygen Flux in a River‐Reservoir System | |
| CN113914853B (zh) | 深水沉积体系砂泥岩薄互层砂岩储层含水饱和度确定方法 | |
| Liu et al. | Correction of Inexactitude in Surficial Sediment Calcium Carbonate Content Data of the Global Ocean |
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 |