CN107025358A - 一种基于数理统计的密闭取心井油水饱和度校正方法及系统 - Google Patents
一种基于数理统计的密闭取心井油水饱和度校正方法及系统 Download PDFInfo
- Publication number
- CN107025358A CN107025358A CN201710266068.8A CN201710266068A CN107025358A CN 107025358 A CN107025358 A CN 107025358A CN 201710266068 A CN201710266068 A CN 201710266068A CN 107025358 A CN107025358 A CN 107025358A
- Authority
- CN
- China
- Prior art keywords
- saturation
- correction
- oil
- water
- pore
- 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.)
- Pending
Links
- 238000012937 correction Methods 0.000 title claims abstract description 126
- 238000000034 method Methods 0.000 title claims abstract description 46
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 98
- 238000005259 measurement Methods 0.000 claims abstract description 54
- 238000002474 experimental method Methods 0.000 claims abstract description 21
- 238000011160 research Methods 0.000 claims abstract description 12
- 230000002159 abnormal effect Effects 0.000 claims abstract description 7
- 239000012530 fluid Substances 0.000 claims description 51
- 230000008569 process Effects 0.000 claims description 20
- 238000013329 compounding Methods 0.000 claims description 17
- 239000002002 slurry Substances 0.000 claims description 17
- 230000000704 physical effect Effects 0.000 claims description 16
- 238000009966 trimming Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 abstract description 12
- 239000011435 rock Substances 0.000 abstract description 10
- 238000009826 distribution Methods 0.000 abstract description 5
- 238000011161 development Methods 0.000 abstract description 4
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000012217 deletion Methods 0.000 abstract 1
- 230000037430 deletion Effects 0.000 abstract 1
- 238000005553 drilling Methods 0.000 description 10
- 238000007789 sealing Methods 0.000 description 7
- 238000011545 laboratory measurement Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 230000009897 systematic effect Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007654 immersion Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 229920002472 Starch Polymers 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000007872 degassing Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002332 oil field water Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 235000019698 starch Nutrition 0.000 description 1
- 239000008107 starch Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于数理统计的密闭取心井油水饱和度校正方法及系统,涉及油气田开发技术领域。所述方法包括：S1、选取研究区密闭取心井不同层位、不同岩性的岩心样品，实验测量每个样品的油、水饱和度数值，删除测量数据中的异常数据；S2、根据测量获得的每个样品的油、水饱和度数值确定造成油水饱和度损失的影响因素；S3、对确定的影响因素进行校正。本发明基于密闭取心井油水饱和度损失的影响因素分析，确定导致实测饱和度与原始地层饱和度存在偏差的主要原因，从数理统计的角度探讨了岩心油、水挥发的变化规律，可以有效地恢复原始地层条件下油水饱和度的真实值，为储层含油性的评价及剩余油分布规律的研究提供重要的基础数据。
Description
技术领域
本发明涉及油气田开发技术领域，尤其涉及一种基于数理统计的密闭取心井油水饱和度校正方法及系统。
背景技术
随着国内各油田勘探开发的不断深入以及老油田综合含水率的不断上升，水淹层的准确评价成为当前各个油田控水稳油和增储挖潜需要解决的关键技术难题，而原始地层条件下储层油水饱和度资料则是判断油层水淹级别、分析剩余油分布规律以及制定合理开发方案的重要依据。为了获取真实地层条件下的油水饱和度资料，通常采用密闭取心方式,使之尽量保持地层原始状态，然后在实验室测量油水饱和度数据。但是，由于在钻井、取心、制样、样品保管以及分析化验过程中，泥浆滤液侵入、脱气、轻质组分挥发及系统误差等工程和地质因素的影响，最终测量的密闭取心井油水饱和度资料与地层真实的饱和度有明显的偏差，一般约为5％～35％。因此，实验室测量的饱和度数据不能直接应用于油藏地质研究。为了消除实测油、水饱和度受降压脱气油水溢出的影响，需要对岩心分析测量的油水饱和度进行校正，使其恢复到地层的真实情况。
目前，常用的油水饱和度校正方法主要有物理模拟实验和经验公式校正两种。前者实验操作复杂，对测量样品的质量要求比较高，一般只考虑单因素的影响，无法模拟饱和度损失的全过程，并且针对大量岩心样品的饱和度校正，其高昂的测试费用与较低的校正效率是该方法最大的不足；后者针对特定的油气藏类型，确定影响油水饱和度损失的地质及油藏因素，优选表征参数建立区域经验公式进行饱和度校正，该方法一方面缺少理论依据，针对特定油气藏适用性较好，但是可移植性与推广性不强；另一方面校正结果的误差较大，不能满足高含水期或复杂油气藏饱和度校正的精度要求。
发明内容
本发明为了解决上述技术问题，提出一种基于数理统计的密闭取心井油水饱和度校正方法及系统。
本发明解决上述技术问题的技术方案如下：一种基于数理统计的密闭取心井油水饱和度校正方法，所述方法包括：
S1、选取研究区密闭取心井不同层位、不同岩性的岩心样品，实验测量每个样品的油、水饱和度数值，删除测量数据中的异常数据；
S2、根据测量获得的每个样品的油、水饱和度数值确定造成油水饱和度损失的影响因素；
S3、对确定的影响因素进行校正。
本发明的有益效果是：本发明基于密闭取心井油水饱和度损失的影响因素分析，确定导致实测饱和度与原始地层饱和度存在偏差的主要原因，从数理统计的角度探讨了岩心油、水挥发的变化规律，可以有效地恢复原始地层条件下油水饱和度的真实值，为储层含油性的评价及剩余油分布规律的研究提供重要的基础数据。
在上述技术方案的基础上，本发明还可以做如下改进。
进一步，S2中所述的影响因素包括泥浆滤液侵入量、孔隙流体体积变化、轻质组分挥发以及饱和度数值测量过程中的误差。
采用上述进一步方案的有益效果是：通过对密闭取心井油水饱和度损失的影响因素进行综合分析，确定导致实测饱和度与原始地层饱和度存在偏差的主要原因，然后分别从泥浆滤液侵入、孔隙流体体积变化、轻质组分挥发、测试系统误差等几个方面分析造成油水饱和度存在偏差的成因机理，并且建立相应的校正模型，准确恢复原始地层条件下油水饱和度数值，为储层含油性的综合评价提供准确的地质参数。
进一步，对所述泥浆滤液侵入量的校正过程具体为：模拟真实钻井过程中各种地质与工程条件，确定S1中所述的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量，并用测量获得的不同层位、不同物性的密闭取心井的岩心样品的含水饱和度减去所对应的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量，从而实现泥浆滤液侵入量的校正。
进一步，对所述孔隙流体体积变化的校正过程具体为：
根据公式：
Sol＝Sos*Bo
获得孔隙流体体积校正后含油饱和度，其中，Sol为孔隙流体体积校正后含油饱和度；Sos为实验测量含油饱和度；Bo为油的体积系数；
根据公式：
Swl＝Sws*Bw
获得孔隙流体体积校正后含水饱和度，其中，Swl为孔隙流体体积校正后含水饱和度；Sws为实验测量含水饱和度；Bw为水的体积系数。
进一步，对所述轻质组分挥发的校正过程具体为：
根据公式
实现对所述轻质组分挥发的校正，其中，η1为轻质组分挥发后的油的剩余率，η2为轻质组分挥发后的水的剩余率。
进一步，对所述饱和度数值测量过程中的误差进行校正的具体过程为：
设饱和度总损失量中油的损失百分比为Y，水的损失百分比为1-Y，则有：
Y＝(1-η1)/[(1-η1)+(1-η2)]；
因此，校正后的油饱和度为：
So＝(1-Sol-Swl)*Y+Sol；
校正后的水饱和度为：
Sw＝(1-Sol-Swl)*(1-Y)+Swl。
为了解决上述技术问题，本发明还提出一种基于数理统计的密闭取心井油水饱和度校正系统，所述系统包括：
数据选取模块，用于选取研究区密闭取心井不同层位、不同岩性的岩心样品，实验测量每个样品的油、水饱和度数值，删除测量数据中的异常数据；
因素确定模块，用于根据测量获得的每个样品的油、水饱和度数值确定造成油水饱和度损失的影响因素；
校正模块，用于对确定的影响因素进行校正。
本发明的有益效果是：本发明基于密闭取心井油水饱和度损失的影响因素分析，确定导致实测饱和度与原始地层饱和度存在偏差的主要原因，从数理统计的角度探讨了岩心油、水挥发的变化规律，可以有效地恢复原始地层条件下油水饱和度的真实值，为储层含油性的评价及剩余油分布规律的研究提供重要的基础数据。
进一步，所述因素确定模块中影响因素包括泥浆滤液侵入量、孔隙流体体积变化、轻质组分挥发以及饱和度数值测量过程中的误差。
采用上述进一步方案的有益效果是：通过对密闭取心井油水饱和度损失的影响因素进行综合分析，确定导致实测饱和度与原始地层饱和度存在偏差的主要原因，然后分别从泥浆滤液侵入、孔隙流体体积变化、轻质组分挥发、测试系统误差等几个方面分析造成油水饱和度存在偏差的成因机理，并且建立相应的校正模型，准确恢复原始地层条件下油水饱和度数值，为储层含油性的综合评价提供准确的地质参数。。
进一步，所述校正模块包括泥浆滤液侵入量校正模块，所述泥浆滤液侵入量校正模块用于模拟真实钻井过程中各种地质与工程条件，确定数据选取模块中所述的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量，并用测量获得的不同层位、不同物性的密闭取心井的岩心样品的含水饱和度减去所对应的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量。
进一步，所述校正模块还包括孔隙流体体积变化校正模块，所述孔隙流体体积变化校正模块用于根据公式
Sol＝Sos*Bo
获得孔隙流体体积校正后含油饱和度，其中，Sol为孔隙流体体积校正后含油饱和度；Sos为实验测量含油饱和度；Bo为油的体积系数；
根据公式：
Swl＝Sws*Bw
获得孔隙流体体积校正后含水饱和度，其中，Swl为孔隙流体体积校正后含水饱和度；Sws为实验测量含水饱和度；Bw为水的体积系数；
所述校正模块还包括轻质组分挥发校正模块，所述轻质组分挥发校正模块用于根据公式
实现对所述轻质组分挥发的校正，其中，η1为轻质组分挥发后的油的剩余率，η2为轻质组分挥发后的水的剩余率；
所述校正模块还包括误差校正模块，所述误差校正模块用于设饱和度总损失量中油的损失百分比为Y，水的损失百分比为1-Y，则有：
Y＝(1-η1)/[(1-η1)+(1-η2)]；
因此，校正后的油饱和度为：
So＝(1-Sol-Swl)*Y+Sol；
校正后的水饱和度为：
Sw＝(1-Sol-Swl)*(1-Y)+Swl。
附图说明
图1为本发明实施例所述的基于数理统计的密闭取心井油水饱和度校正方法的流程图；
图2为本发明实施例所述的孔隙流体体积校正后的油、水饱和度交会图；
图3为本发明实施例所述的实验测量油水饱和度与理论值的对比示意图；
图4为本发明实施例所述的基于数理统计的密闭取心井油水饱和度校正系统原理图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
如图1所示，本实施例提出一种基于数理统计的密闭取心井油水饱和度校正方法，所述方法包括：
S101、选取研究区密闭取心井不同层位、不同岩性的岩心样品，实验测量每个样品的油、水饱和度数值，删除测量数据中的异常数据。
具体的，以克拉玛依油田八区530井区八道湾组砾岩油藏为例，该油藏目前正在进行老油田“二三结合”的试点工程开发，饱和度资料是整个开发方案设计的关键地质参数。因此，研究区设计钻取3口密闭取心井以获得各种储层参数，以T88724井为例说明油水饱和度校正的技术细节和整个过程。该井共计钻取饱和度样品129个，在实验室中分别测量了油、水饱和度数值，129个样品油水饱和度之和的平均值为87.2％，整体损失量达到12.8％，需要进行饱和度校正。
S102、根据测量获得的每个样品的油、水饱和度数值确定造成油水饱和度损失的影响因素，并对确定的影响因素进行校正。
具体的，八道湾组砾岩储层不含游离气，当储层不存在游离气时，岩石孔隙内油水饱和度之和等于100％，但实验室测量的油水饱和度之和通常小于100％，大约存在5％～35％的偏差，造成偏差的主要影响因素包括以下几方面：①钻井过程中泥浆滤液侵入造成油水饱和度异常；②孔隙流体体积变化引起油水饱和度变化；③压力和温度变化引起孔隙中轻质组分挥发，进而导致油水饱和度变化；④测量方法的系统误差也会导致油水饱和度出现偏差。在以上四个影响因素中，压力和温度变化引起孔隙中轻质组分挥发是导致岩心饱和度存在偏差的最主要原因，孔隙流体体积变化的影响次之，泥浆滤液侵入和测量系统误差造成的饱和度变化量虽然比较小，但是也不可忽略，同样需要进行校正。以下对每种影响因素分别进行校正，以获取原始地层条件下的油、水饱和度值。
泥浆滤液侵入量校正：在钻井过程中为了保证施工安全，会利用人工配制的泥浆来平衡地层压力，密闭取心样品也会侵入少量泥浆，导致孔隙中含水饱和度增加，因此，首先需要对泥浆侵入量进行校正。泥浆侵入量与储层岩性、渗透性、泥浆滤液矿化度等内部因素和泥浆柱压力、地层孔隙压力等外部因素有关，选择不同岩性、不同物性的密闭取心样品，在实验室中模拟真实钻井过程中各种地质与工程条件，确定不同岩性和物性样品的泥浆侵入量。用实验室测量的岩心含水饱和度减去泥浆侵入量，消除泥浆侵入对饱和度校正的影响。
选择密闭取心井的岩心样品，在实验室中模拟真实钻井过程中各种地质与工程条件，确定不同样品的泥浆浸入量，如表1所示，用实验测量的含水饱和度减去泥浆滤液侵入量，消除泥浆侵入对饱和度的影响。对于泥浆侵入量为0的岩心样品则不需要进行校正。
表1
孔隙流体体积变化的校正：密闭取心样品从地下取出至地面，由于地层条件发生了很大变化，流体体积也会发生较大变化，油水体积均会变小。因此，实验室测量的油水饱和度，由于体积变化，转化成原始地层条件下的饱和度时将有所增大，两者可以通过油、水体积系数进行换算，校正公式如下：
根据公式：
Sol＝Sos*Bo
获得孔隙流体体积校正后含油饱和度，其中，Sol为孔隙流体体积校正后含油饱和度；Sos为实验测量含油饱和度；Bo为油的体积系数。
根据公式：
Swl＝Sws*Bw
获得孔隙流体体积校正后含水饱和度，其中，Swl为孔隙流体体积校正后含水饱和度；Sws为实验测量含水饱和度；Bw为水的体积系数。
利用实验模拟技术确定八道湾组砾岩油藏的油、水体积系数，其数值分别为1.078和1.012。在通过上述公式进行校正后的油、水饱和度数值如表1所示。
轻质组分挥发校正：在泥浆滤液浸入量和孔隙流体体积变化的校正后，基于油、水剩余率的概念，利用数理统计的方式制作孔隙流体体积校正后的油、水饱和度交会图，如图2所示，从图中可以看出，油水饱和度呈现出很好的线性关系，相关系数的平方达到0.8857，并且其斜率为-0.8404，说明该井取心段的含油饱和度损失大于含水饱和度，即油的挥发大于水的挥发。依据图2建立的油水饱和度相关关系。
当储层孔隙中不含游离气而只有油水两相流体时，油水饱和度之间的关系满足下式：
So+Sw＝1
式中，So为原始地层条件下储层的含油饱和度，Sw为原始地层条件下储层的含水饱和度。
由上式可知，原始地层条件下孔隙中的含油饱和度和含水饱和度之和等于1，两者呈线性关系，其斜率为-1，截距为1。当地层中的岩心被钻取之后，由于温度和压力的变化，轻烃组分和空隙中的水就会挥发，假设挥发以后油的剩余率为η1，水的剩余率为η2，则有：
Sol＝So*η1
Swl＝Sw*η2
综合上式可以得到经过泥浆侵入与孔隙流体体积校正后的油水饱和度满足以下关系：
Sol/η1+Swl/η2＝1
对上式进行变形后可以得到一个线性关系式：
由上式可以看出，经过挥发后岩心测量的油水饱和度也呈线性关系，其斜率和截距与油、水的剩余率有关，即油、水的挥发率决定着岩心测量油水饱和度的线性关系。对于岩性与物性相同的密闭取心样品，油水损失率基本一致，那么就可以通过数理统计回归的方法来确定油水的剩余率η1、η2，进而进行轻质组分挥发校正。
饱和度测试方法的系统误差校正：密闭取心过程中，由于油、水饱和度是分两次独立测量，校正后得到的油、水饱和度之和不一定为100％，通常会有少量浮动。因此，还需要对数据进行系统误差校正，校正方法采用的是将饱和度数据的校正系数转换成总损失量的百分含量，则校正后油水饱和度数据之和为100％。
设饱和度总损失量中油的损失百分比为Y，水的损失百分比为1-Y，则有：
Y＝(1-η1)/[(1-η1)+(1-η2)]；
因此，校正后的油饱和度为：
So＝(1-Sol-Swl)*Y+Sol；
校正后的水饱和度为：
Sw＝(1-Sol-Swl)*(1-Y)+Swl。
基于上述计算过程可以确定饱和度总损失量中油的损失百分比为0.8653，水的损失百分比为0.1347。然后进行饱和度校正，恢复原始地层条件下含油饱和度和含水饱和度值，校正后的油水饱和度之和都等于100％。
本实施例中所述的饱和度校正方法在实际应用时需要具备两个条件：①储层孔隙流体中不含游离气体；②岩心经过挥发后实验室测量的油水饱和度呈线性关系。依据大量的实验统计数据发现，对于同一口密闭取心井，其油、水饱和度之间都存在较好的线性关系，但这种线性关系的截距和斜率不再是1和-1，说明基于密闭取心井测量的油、水饱和度资料存在误差。从图3中实验测量油水饱和度与理论值的对比示意图中可以看出，如果密闭取心资料不受任何外界因素的影响，油水饱和度分析数据会落在MN线上，如果分析数据存在系统误差，则分析数据会落在MN线的两侧，当油水饱和度之和大于1时，落在MN线的右上方；当油水饱和度之和小于1时，落在MN线的左下方。一般情况下，由于油水挥发的影响，分析数据会落在MN线的左下方。另外，如果线性关系的斜率小于-1，即线性关系如CD，则反映了含水饱和度的损失大于含油饱和度；如果斜率大于-1，即线性关系如EF，则反映了含油饱和度的损失大于含水饱和度。因此，依据饱和度的数理统计关系可以分析油水损失的系统误差。
本实施例所述的校正方法是对密闭取心井岩心饱和度损失的总校正，方法分别从泥浆滤液侵入、孔隙流体体积变化、轻质组分挥发、测试系统误差几方面明确了造成油水饱和度存在偏差的成因机理，并且建立了相应的校正模型。在实际资料处理中，该方法具有快速、准确及易操作的特点，能够准确恢复原始地层条件下油、水饱和度数值，为储层含油性的综合评价及剩余油分布规律的研究提供可靠的地质参数。特别是对于高含水期老油田的二次调整开发具有重要的技术支撑作用，可以快速、准确地校正大量密闭取心井的油水饱和度资料，具有重要的油田参考价值与广泛的市场应用前景。
对应的，如图4所示，本实施例还提出了一种基于数理统计的密闭取心井油水饱和度校正系统，所述系统包括：
数据选取模块，用于选取研究区密闭取心井不同层位、不同岩性的岩心样品，实验测量每个样品的油、水饱和度数值，删除测量数据中的异常数据；
因素确定模块，用于根据测量获得的每个样品的油、水饱和度数值确定造成油水饱和度损失的影响因素，所述影响因素包括泥浆滤液侵入量、孔隙流体体积变化、轻质组分挥发以及饱和度数值测量过程中的误差；
校正模块，用于对确定的影响因素进行校正。
优选的，所述校正模块包括：泥浆滤液侵入量校正模块、孔隙流体体积变化校正模块、轻质组分挥发校正模块和误差校正模块。
所述泥浆滤液侵入量校正模块用于模拟真实钻井过程中各种地质与工程条件，确定数据选取模块中所述的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量，并用测量获得的不同层位、不同物性的密闭取心井的岩心样品的含水饱和度减去所对应的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量。
所述校正模块还包括孔隙流体体积变化校正模块，所述孔隙流体体积变化校正模块用于根据公式
Sol＝Sos*Bo
获得孔隙流体体积校正后含油饱和度，其中，Sol为孔隙流体体积校正后含油饱和度；Sos为实验测量含油饱和度；Bo为油的体积系数；
根据公式：
Swl＝Sws*Bw
获得孔隙流体体积校正后含水饱和度，其中，Swl为孔隙流体体积校正后含水饱和度；Sws为实验测量含水饱和度；Bw为水的体积系数；
所述校正模块还包括轻质组分挥发校正模块，所述轻质组分挥发校正模块用于根据公式
实现对所述轻质组分挥发的校正，其中，η1为轻质组分挥发后的油的剩余率，η2为轻质组分挥发后的水的剩余率；
所述校正模块还包括误差校正模块，所述误差校正模块用于设饱和度总损失量中油的损失百分比为Y，水的损失百分比为1-Y，则有：
Y＝(1-η1)/[(1-η1)+(1-η2)]；
因此，校正后的油饱和度为：
So＝(1-Sol-Swl)*Y+Sol；
校正后的水饱和度为：
Sw＝(1-Sol-Swl)*(1-Y)+Swl。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于数理统计的密闭取心井油水饱和度校正方法，其特征在于，所述方法包括：
S1、选取研究区密闭取心井不同层位、不同岩性的岩心样品，实验测量每个样品的油、水饱和度数值，删除测量数据中的异常数据；
S2、根据测量获得的每个样品的油、水饱和度数值确定造成油水饱和度损失的影响因素；
S3、对确定的影响因素进行校正。
2.根据权利要求1所述的一种基于数理统计的密闭取心井油水饱和度校正方法，其特征在于，S2中所述的影响因素包括泥浆滤液侵入量、孔隙流体体积变化、轻质组分挥发以及饱和度数值测量过程中的误差。
3.根据权利要求2所述的一种基于数理统计的密闭取心井油水饱和度校正方法，其特征在于，对所述泥浆滤液侵入量的校正过程具体为：模拟真实钻井过程中各种地质与工程条件，确定S1中所述的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量，并用测量获得的不同层位、不同物性的密闭取心井的岩心样品的含水饱和度减去所对应的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量，从而实现泥浆滤液侵入量的校正。
4.根据权利要求3所述的一种基于数理统计的密闭取心井油水饱和度校正方法，其特征在于，对所述孔隙流体体积变化的校正过程具体为：
根据公式：
Sol＝Sos*Bo
获得孔隙流体体积校正后含油饱和度，其中，Sol为孔隙流体体积校正后含油饱和度；Sos为实验测量含油饱和度；Bo为油的体积系数；
根据公式：
Swl＝Sws*Bw
获得孔隙流体体积校正后含水饱和度，其中，Swl为孔隙流体体积校正后含水饱和度；Sws为实验测量含水饱和度；Bw为水的体积系数。
5.根据权利要求4所述的一种基于数理统计的密闭取心井油水饱和度校正方法，其特征在于，对所述轻质组分挥发的校正过程具体为：
根据公式
实现对所述轻质组分挥发的校正，其中，η1为轻质组分挥发后的油的剩余率，η2为轻质组分挥发后的水的剩余率。
6.根据权利要求5所述的一种基于数理统计的密闭取心井油水饱和度校正方法，其特征在于，对所述饱和度数值测量过程中的误差进行校正的具体过程为：
设饱和度总损失量中油的损失百分比为Y，水的损失百分比为1-Y，则有：
Y＝(1-η1)/[(1-η1)+(1-η2)]；
因此，校正后的油饱和度为：
So＝(1-Sol-Swl)*Y+Sol；
校正后的水饱和度为：
Sw＝(1-Sol-Swl)*(1-Y)+Swl。
7.一种基于数理统计的密闭取心井油水饱和度校正系统，其特征在于，所述系统包括：
数据选取模块，用于选取研究区密闭取心井不同层位、不同岩性的岩心样品，实验测量每个样品的油、水饱和度数值，删除测量数据中的异常数据；
因素确定模块，用于根据测量获得的每个样品的油、水饱和度数值确定造成油水饱和度损失的影响因素；
校正模块，用于对确定的影响因素进行校正。
8.根据权利要求7所述的一种基于数理统计的密闭取心井油水饱和度校正系统，其特征在于，所述因素确定模块中影响因素包括泥浆滤液侵入量、孔隙流体体积变化、轻质组分挥发以及饱和度数值测量过程中的误差。
9.根据权利要求8所述的一种基于数理统计的密闭取心井油水饱和度校正系统，其特征在于，所述校正模块包括泥浆滤液侵入量校正模块，所述泥浆滤液侵入量校正模块用于模拟真实钻井过程中各种地质与工程条件，确定数据选取模块中所述的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量，并用测量获得的不同层位、不同物性的密闭取心井的岩心样品的含水饱和度减去所对应的不同层位、不同物性的密闭取心井的岩心样品的泥浆侵入量。
10.根据权利要求9所述的一种基于数理统计的密闭取心井油水饱和度校正系统，其特征在于，所述校正模块还包括孔隙流体体积变化校正模块，所述孔隙流体体积变化校正模块用于根据公式
Sol＝Sos*Bo
获得孔隙流体体积校正后含油饱和度，其中，Sol为孔隙流体体积校正后含油饱和度；Sos为实验测量含油饱和度；Bo为油的体积系数；
根据公式：
Swl＝Sws*Bw
获得孔隙流体体积校正后含水饱和度，其中，Swl为孔隙流体体积校正后含水饱和度；Sws为实验测量含水饱和度；Bw为水的体积系数；
所述校正模块还包括轻质组分挥发校正模块，所述轻质组分挥发校正模块用于根据公式
实现对所述轻质组分挥发的校正，其中，η1为轻质组分挥发后的油的剩余率，η2为轻质组分挥发后的水的剩余率；
所述校正模块还包括误差校正模块，所述误差校正模块用于设饱和度总损失量中油的损失百分比为Y，水的损失百分比为1-Y，则有：
Y＝(1-η1)/[(1-η1)+(1-η2)]；
因此，校正后的油饱和度为：
So＝(1-Sol-Swl)*Y+Sol；
校正后的水饱和度为：
Sw＝(1-Sol-Swl)*(1-Y)+Swl。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710266068.8A CN107025358A (zh) | 2017-04-21 | 2017-04-21 | 一种基于数理统计的密闭取心井油水饱和度校正方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710266068.8A CN107025358A (zh) | 2017-04-21 | 2017-04-21 | 一种基于数理统计的密闭取心井油水饱和度校正方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107025358A true CN107025358A (zh) | 2017-08-08 |
Family
ID=59527688
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710266068.8A Pending CN107025358A (zh) | 2017-04-21 | 2017-04-21 | 一种基于数理统计的密闭取心井油水饱和度校正方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107025358A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112710806A (zh) * | 2019-10-25 | 2021-04-27 | 中国石油化工股份有限公司 | 一种砂岩油藏取心饱和度的校正方法 |
| WO2021108877A1 (pt) * | 2019-12-02 | 2021-06-10 | Petróleo Brasileiro S.A. - Petrobras | Metodologia do mapa de qualidade de hphiso filtrado por transmisisbilidade e qualidade de varrido para pósprocessamento de simulações de fluxo de reservatórios de petróleo |
| CN114965140A (zh) * | 2021-02-24 | 2022-08-30 | 中国石油天然气股份有限公司 | 活油储层密闭取心油水饱和度校正方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105301220A (zh) * | 2014-06-06 | 2016-02-03 | 中国石油化工股份有限公司 | 密闭取心分析饱和度的校正方法 |
| CN106153854A (zh) * | 2015-03-30 | 2016-11-23 | 中国石油化工股份有限公司 | 密闭取心井饱和度校正方法 |
-
2017
- 2017-04-21 CN CN201710266068.8A patent/CN107025358A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105301220A (zh) * | 2014-06-06 | 2016-02-03 | 中国石油化工股份有限公司 | 密闭取心分析饱和度的校正方法 |
| CN106153854A (zh) * | 2015-03-30 | 2016-11-23 | 中国石油化工股份有限公司 | 密闭取心井饱和度校正方法 |
Non-Patent Citations (2)
| Title |
|---|
| 王谦 等: "油基泥浆密闭取心饱和度校正方法", 《测井技术》 * |
| 谭锋奇 等: "密闭取心井油水饱和度校正方法", 《地球科学(中国地质大学学报)》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112710806A (zh) * | 2019-10-25 | 2021-04-27 | 中国石油化工股份有限公司 | 一种砂岩油藏取心饱和度的校正方法 |
| WO2021108877A1 (pt) * | 2019-12-02 | 2021-06-10 | Petróleo Brasileiro S.A. - Petrobras | Metodologia do mapa de qualidade de hphiso filtrado por transmisisbilidade e qualidade de varrido para pósprocessamento de simulações de fluxo de reservatórios de petróleo |
| CN114965140A (zh) * | 2021-02-24 | 2022-08-30 | 中国石油天然气股份有限公司 | 活油储层密闭取心油水饱和度校正方法 |
| CN114965140B (zh) * | 2021-02-24 | 2024-05-28 | 中国石油天然气股份有限公司 | 活油储层密闭取心油水饱和度校正方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bourbiaux | Fractured reservoir simulation: a challenging and rewarding issue | |
| US20070112518A1 (en) | Method and Apparatus for Measuring the Wettability of Geological Formations | |
| CA2692425A1 (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
| US20200378249A1 (en) | Systems and Methods for Identifying Two or More Charges into Reservoir Using Downhole Fluid Analysis | |
| CN108194077B (zh) | 一种气测全烃校正方法 | |
| CN105257285A (zh) | 油气由烃源岩层向岩溶储层充注能力的评价方法 | |
| CN107025358A (zh) | 一种基于数理统计的密闭取心井油水饱和度校正方法及系统 | |
| CN110593855B (zh) | 一种测井电阻率的校正方法、油层的识别方法 | |
| Woodhouse | Accurate reservoir water saturations from oil-mud cores: Questions and answers from Prudhoe Bay and beyond | |
| CN116047602B (zh) | 基于生烃数值模拟的ii型水合物饱和度预测方法 | |
| CN112560246A (zh) | 一种目标井散点地层压力系数的预测方法 | |
| CN102134994B (zh) | 一种基于基岩油藏油水层电阻率的地层数据处理方法 | |
| Brummett Jr et al. | Reservoir description by simulation at SACROC-a case history | |
| Carnegie | Understanding the pressure gradients improves production from oil/water transition carbonate zones | |
| Oedietram et al. | Integrating petrophysics and production performance to determine variable initial water salinity conditions impacting water saturation calculations | |
| Seth et al. | Saturation height function in a field under imbibition: a case study | |
| CN113777655A (zh) | 一种盐湖盆地小层烃源岩平面分布预测方法 | |
| Al-Mozan et al. | Reservoir Modeling for Mishrif Formation in Nasiriyah Oilfield | |
| Underschultz et al. | Estimating formation water salinity from wireline pressure data: Case study in the Vulcan sub-basin | |
| Baker et al. | History matching standards; quality control and risk analysis for simulation | |
| CN113720745A (zh) | 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法 | |
| Wach et al. | Well placement, cost reduction, and increased production using reservoir models based on outcrop, core, well-log, seismic data, and modern analogs: Onshore and offshore Western Trinidad | |
| Kefeng et al. | Reserves evaluation based on reservoir geological knowledge database | |
| CN111827989B (zh) | 一种沥青削减指数识别含沥青砾岩油层的方法 | |
| Junjie et al. | COMPREHENSIVE CORRECTION METHOD OF AIRTIGHT CORING SATURATION BASED ON CORE CLASSIFICATION. |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170808 |