CN102323198B - 岩心岩电实验参数m、n的校正方法及系统 - Google Patents
岩心岩电实验参数m、n的校正方法及系统 Download PDFInfo
- Publication number
- CN102323198B CN102323198B CN 201110147531 CN201110147531A CN102323198B CN 102323198 B CN102323198 B CN 102323198B CN 201110147531 CN201110147531 CN 201110147531 CN 201110147531 A CN201110147531 A CN 201110147531A CN 102323198 B CN102323198 B CN 102323198B
- Authority
- CN
- China
- Prior art keywords
- saturation
- rock core
- rock
- under
- core
- 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
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种岩心岩电实验参数m、n的校正方法及系统。其中所述校正方法包括在岩心在未达到完全饱和状态下，根据校正公式 对饱和度指数值n进行校正；在岩心未完全饱和状态下，根据校正公式对胶结指数值m进行校正。所述校正系统包括饱和度指数值校正模块及胶结指数值校正模块。根据本发明实施例提供的岩心岩电实验参数m、n的校正方法及系统，解决了低孔隙度低渗透率岩心未完全饱和状态下按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，致使饱和度指数n值降低、胶结指数m值增大的问题。
Description
技术领域
本发明涉及一种岩心岩电实验参数m、n的校正方法及系统。
背景技术
岩电实验首先要求用地层水完全饱和岩心，然而对于低孔隙度低渗透率岩心，由于受实验方法、实验设备以及实验周期的限制，岩心往往不能达到完全饱和，它会对地层因素F和电阻增大率I的实验测量结果有较明显的影响，进而影响到饱和度指数n值和胶结指数m值的求取，尤其是对于低孔隙度低渗透率地层的岩心，在不能达到完全饱和的情况下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，会导致饱和度指数n值降低、胶结指数m值增大。
岩心未达到完全饱和状态下求得的岩电实验参数n′值和m′值的绝对误差虽然较小，但相对误差却往往较大，应对未完全饱和状态下求得的岩电实验参数n′和m′加以校正。
发明内容
本发明的目的之一是为了解决低孔隙度低渗透率岩心未完全饱和状态下按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，致使饱和度指数n值降低、胶结指数m值增大的问题，提供一种岩心岩电实验参数m、n的校正方法及系统。
根据本发明的一个方面，提供一种岩心岩电实验参数m、n的校正方法包括：
在岩心在未达到完全饱和状态下，根据校正公式
对饱和度指数值n进行校正；其中，所述n′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的饱和度指数；所述n为n′经公式(1)校正后计算求取的饱和度指数；所述b′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的岩心系数，所述b′为大于1的小数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；
在岩心未完全饱和状态下，根据校正公式
对胶结指数值m进行校正；所述m′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的胶结指数；所述m为m′经公式(2)校正后计算求取的胶结指数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；n为n′经公式(1)校正后计算求取的饱和度指数；所述φ为岩心孔隙度，φ记为小数。
根据本发明的另一个方面，提供一种岩心岩电实验参数m、n的校正系统包括：
饱和度指数值校正模块，在岩心在未达到完全饱和状态下，根据校正公式
对饱和度指数值n进行校正；其中，所述n′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T5385-1991)》标准流程进行实验及计算求取的饱和度指数；所述n为n′经公式(1)校正后计算求取的饱和度指数；所述b′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的岩心系数，所述b′为大于1的小数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；
胶结指数值校正模块，在岩心未完全饱和状态下，根据校正公式
对胶结指数值m进行校正；所述m′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T5385-1991)》标准流程进行实验及计算求取的胶结指数；所述m为m′经公式(2)校正后计算求取的胶结指数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；n为n′经公式(1)校正后计算求取的饱和度指数；所述φ为岩心孔隙度，φ记为小数。
根据本发明提供的岩心岩电实验参数m、n的校正方法及系统，分别对未达到完全饱和状态下的岩心按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验求取的胶结指数m′值和饱和度指数n′值进行校正，校正后的胶结指数m值和饱和度指数n值更接近于岩心真实的胶结指数m值和饱和度指数n值，校正效果明显。
附图说明
图1是本发明实施例提供的岩心未达到完全饱和对胶结指数m的影响及校正示意图；
图2是本发明实施例提供的岩心未达到完全饱和对饱和度指数n的影响及校正的示意图；
图3是本发明实施例提供的岩心未达到完全饱和对计算含水饱和度的影响(校正前)的示意图；
图4岩心未达到完全饱和对计算含水饱和度的影响(校正后)的示意图。
具体实施方式
如图1所示，本发明实施例提供的一种岩心岩电实验参数m、n的校正方法包括：
步骤10、在岩心在未达到完全饱和状态下，根据校正公式 
对饱和度指数值n进行校正；其中，所述n′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的饱和度指数，对未达到完全饱和状态下的岩心，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，每块岩心测量5～7个不同的含水饱和度(含水饱和度介于0～100％之间)及其对应的电阻率，求得每块岩心在未达到完全饱和状态下的饱和度指数n′值。n为n′经公式(1)校正后计算求取的饱和度指数。b′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的岩心系数，b′为大于1的小数。Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数，例如，对每一块岩心，选取未达到完全饱和状态下按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程实验测量的5～7个不同的含水饱和度中最大的含水饱和度Sw。图1是岩心未达到完全饱和状态下实验的饱和度指数n校正前、校正后与完全饱和状态下实验结果对比图，从图1可以看出，校正后的饱和度指数n比校正前的饱和度指数n′明显更接近于岩心真实的饱和度指数值，校正效果明显。
步骤20、在岩心未完全饱和状态下，根据校正公式
对胶结指数值m进行校正；所述m′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的胶结指数。对未达到完全饱和状态下的岩心，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，测量每块岩心的地层因素F及其对应的孔隙度φ，求得未达到完全饱和状态下的胶结指数m′值。所述m为m′经公式(2)校正后计算求取的胶结指数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；n为n′经公式(1)校正后计算求取的饱和度指数；所述φ为岩心孔隙度，φ记为小数。图2是岩心未达到完全饱和状态下实验的胶结指数m校正前、校正后与完全饱和状态下实验结果对比图，从附图2可以看出，校正后的胶结指数m比校正前的胶结指数m′明显更接近于岩心真实的胶结指数值，校正效果明显。
例如，选择13块岩心样品(孔隙度在4％～18％之间)在岩心未达到完全饱和状态(测得的岩心含水饱和度Sw在85％～98％之间)条件下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，分别求得岩心未达到完全饱和状态下的胶结指数m′值和饱和度指数n′值。
通过本发明提供的校正公式(1)、(2)分别对岩心未达到完全饱和状态下的岩电实验参数m′值、n′值进行校正，求得岩石真实的胶结指数m值和饱和度指数n值。
上述方法还包括：
步骤30、利用阿尔奇饱和度模型公式，把岩心未达到完全饱和状态下计算求取的m′、n′代入阿尔奇饱和度模型公式，计算求得岩心未达到完全饱和状态下岩心含水饱和度S′w。图3提供了用m′、n′计算的岩心未达到完全饱和状态下的含水饱和度S′w。图3是13块岩心未达到完全饱和状态对计算含水饱和度的影响图(校正前)，从图3可以看出，岩心未达到完全饱和状态计算的每块岩心的含水饱和度(校正前)均高于实际含水饱和度，不对岩心未达到完全饱和状态下的胶结指数m′和饱和度指数n′进行校正，直接代入阿尔奇饱和度模型计算的含水饱和度与岩心实际含水饱和度会存在较大的误差。
步骤40、利用阿尔奇饱和度模型公式，把m′、n′经该发明提供的校正公式校正后的m、n值代入阿尔奇饱和度模型公式，计算求得饱和状态下岩心含水饱和度Sw″，更加接近于岩石的实际含水饱和度。图4提供了利用校正后的胶结指数m、饱和度指数n计算的饱和状态下的含水饱和度Sw″。图4是13块岩心未达到完全饱和状态对计算含水饱和度的影响图(校正后)，从图4可以看出，经对岩心未达到完全饱和状态下胶结指数m′和饱和度指数n′进行校正后的胶结指数m和饱和度指数n代入阿尔奇饱和度模型计算的含水饱和度在45度线附近，接近于岩心实际含水饱和度，校正效果明显。
本发明实施例还提供一种岩心岩电实验参数m、n的校正系统包括饱和度指数值校正模块、胶结指数值校正模块、第一计算模块及第二计算模块。其中，饱和度指数值校正模块，在岩心在未达到完全饱和状态下，根据校正公式
对饱和度指数值n进行校正；其中，所述n′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T5385-1991)》标准流程进行实验及计算求取的饱和度指数；所述n为n′经公式(1)校正后计算求取的饱和度指数；所述b′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的岩心系数，所述b′为大于1的小数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数。
胶结指数值校正模块，在岩心未完全饱和状态下，根据校正公式对胶结指数值m进行校正；所述m′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T5385-1991)》标准流程进行实验及计算求取的胶结指数；所述m为m′经公式(2)校正后计算求取的胶结指数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；n为n′经公式(1)校正后计算求取的饱和度指数；所述φ为岩心孔隙度，φ记为小数。
对每一块岩心，选取未达到完全饱和状态下按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程实验测量的5～7个不同的含水饱和度中最大的含水饱和度Sw。
第一计算模块，将所述m′、n′代入阿尔奇饱和度模型公式，计算求得岩心未达到完全饱和状态下岩心含水饱和度S′w。
第二计算模块，将所述胶结指数值m、饱和度指数值n代入阿尔奇饱和度模型公式，计算求得饱和状态下岩心含水饱和度Sw″。
根据本发明实施例提供的岩心岩电实验参数m、n的校正方法及系统，解决了低孔隙度低渗透率岩心未完全饱和状态下按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，致使饱和度指数n值降低、胶结指数m值增大的问题。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。
Claims (8)
1.一种岩心岩电实验参数m、n的校正方法，其特征在于，包括：
在岩心未达到完全饱和状态下，根据校正公式
对饱和度指数值n进行校正；其中，所述n′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的饱和度指数；所述n为n′经公式(1)校正后计算求取的饱和度指数；所述b′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的岩心系数，所述b′为大于1的小数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；
在岩心未完全饱和状态下，根据校正公式
对胶结指数值m进行校正；所述m′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的胶结指数；所述m为m′经公式(2)校正后计算求取的胶结指数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；n为n′经公式(1)校正后计算求取的饱和度指数；所述φ为岩心孔隙度，φ记为小数。
2.根据权利要求1所述的岩心岩电实验参数m、n的校正方法，其特征在于，还包括：
将所述m′、n′代入阿尔奇饱和度模型公式，计算求得岩心未达到完全饱和状态下岩心含水饱和度S′w。
3.根据权利要求1所述的岩心岩电实验参数m、n的校正方法，其特征在于，还包括：
将所述胶结指数值m、饱和度指数值n代入阿尔奇饱和度模型公式，计算求得饱和状态下岩心含水饱和度Sw″。
4.根据权利要求1-3任一项所述的方法，其特征在于：
对每一块岩心，选取未达到完全饱和状态下按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程实验测量的5～7个不同的含水饱和度中最大的含水饱和度Sw。
5.一种岩心岩电实验参数m、n的校正系统，其特征在于，包括：
饱和度指数值校正模块，在岩心在未达到完全饱和状态下，根据校正公式
对饱和度指数值n进行校正；其中，所述n′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T5385-1991)》标准流程进行实验及计算求取的饱和度指数；所述n为n′经公式(1)校正后计算求取的饱和度指数；所述b′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算求取的岩心系数，所述b′为大于1的小数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；
胶结指数值校正模块，在岩心未完全饱和状态下，根据校正公式对胶结指数值m进行校正；所述m′为岩心在未达到完全饱和状态下，按照《岩石电阻率参数实验室测量及计算方法(SY/T5385-1991)》标准流程进行实验及计算求取的胶结指数；所述m为m′经公式(2)校正后计算求取的胶结指数；所述Sw为岩心未达到完全饱和状态下实验测量所得含水饱和度中的最大含水饱和度，Sw记为小数；n为n′经公式(1)校正后计算求取的饱和度指数；所述φ为岩心孔隙度，φ记为小数。
6.根据权利要求5所述的校正系统，其特征在于，还包括：
第一计算模块，将所述m′、n′代入阿尔奇饱和度模型公式，计算求得岩心未达到完全饱和状态下岩心含水饱和度S′w。
7.根据权利要求5所述的校正系统，其特征在于，还包括：
第二计算模块，将所述胶结指数值m、饱和度指数值n代入阿尔奇饱和度模型公式，计算求得饱和状态下岩心含水饱和度Sw″。
8.根据权利要求5至7任一项所述的校正系统，其特征在于：
对每一块岩心，选取未达到完全饱和状态下按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程实验测量的5～7个不同的含水饱和度中最大的含水饱和度Sw。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110147531 CN102323198B (zh) | 2011-06-02 | 2011-06-02 | 岩心岩电实验参数m、n的校正方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110147531 CN102323198B (zh) | 2011-06-02 | 2011-06-02 | 岩心岩电实验参数m、n的校正方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102323198A CN102323198A (zh) | 2012-01-18 |
| CN102323198B true CN102323198B (zh) | 2013-07-31 |
Family
ID=45450984
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201110147531 Active CN102323198B (zh) | 2011-06-02 | 2011-06-02 | 岩心岩电实验参数m、n的校正方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102323198B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104101905B (zh) * | 2013-04-11 | 2017-02-08 | 中国石油天然气集团公司 | 一种基于岩电参数的储层分类方法 |
| CN104181090B (zh) * | 2013-05-27 | 2016-11-16 | 中国石油化工股份有限公司 | 一种基于岩电结构系数的含油气饱和度评价方法 |
| CN104675391B (zh) * | 2013-11-26 | 2018-04-10 | 中国石油化工股份有限公司 | 计算地层含油饱和度的方法 |
| CN109386281B (zh) * | 2017-08-02 | 2021-11-09 | 中国石油化工股份有限公司 | 一种获取裂缝性低孔低渗储层测井饱和度的方法 |
| CN109901238B (zh) * | 2019-02-28 | 2020-09-04 | 中国石油天然气集团有限公司 | 一种基于应力差电阻率实验的高应力地层电阻率校正方法 |
| CN112800646B (zh) * | 2021-01-08 | 2024-01-23 | 中海石油深海开发有限公司 | 一种修正阿尔奇模型中岩电参数的方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101358943A (zh) * | 2007-07-31 | 2009-02-04 | 中国石油天然气集团公司 | 一种欠饱和岩心岩电实验参数m、n的校正方法 |
-
2011
- 2011-06-02 CN CN 201110147531 patent/CN102323198B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101358943A (zh) * | 2007-07-31 | 2009-02-04 | 中国石油天然气集团公司 | 一种欠饱和岩心岩电实验参数m、n的校正方法 |
Non-Patent Citations (4)
| Title |
|---|
| 岩石孔隙度密度的自动测定;黄福堂;《大庆石油地质与开发》;19931231;第12卷(第4期);全文 * |
| 岩石气体孔隙度测量不确定度分析;柯式镇等;《计量学报》;20070430;第28卷(第2期);全文 * |
| 柯式镇等.岩石气体孔隙度测量不确定度分析.《计量学报》.2007,第28卷(第2期), |
| 黄福堂.岩石孔隙度密度的自动测定.《大庆石油地质与开发》.1993,第12卷(第4期), |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102323198A (zh) | 2012-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102323198B (zh) | 岩心岩电实验参数m、n的校正方法及系统 | |
| CN103617816B (zh) | 反应堆堆芯功率分布的测量方法 | |
| CN101358943A (zh) | 一种欠饱和岩心岩电实验参数m、n的校正方法 | |
| CN106053899B (zh) | 一种基准信号发生系统及方法 | |
| CN105258847A (zh) | 一种压力传感器校准的方法和装置 | |
| CN103983938B (zh) | 电能表小信号偏置的校正方法 | |
| CN107339938A (zh) | 一种用于单目立体视觉自标定的异型标定块及标定方法 | |
| CN105403800A (zh) | 一种电力系统参数的确定方法及装置 | |
| CN108594160A (zh) | 一种电能表的校表方法 | |
| TW201137373A (en) | Estimation method for residual discharging time of batteries | |
| CN108982949A (zh) | 一种基于线电压、相电压大小获取其相量的方法 | |
| CN103234647A (zh) | 一种嵌入式系统的温度校准方法及系统 | |
| CN101373211A (zh) | 三相电能表调试方法 | |
| CN105572191A (zh) | 一种电化学气体传感器的压力补偿方法 | |
| CN105372621A (zh) | 一种校准装置的方法、移动式校对装置以及现场测量设备 | |
| CN111198891B (zh) | 数据源融合方法、电子设备及非暂态计算机可读存储介质 | |
| CN109101466B (zh) | 基于分布函数取对数变换的威布尔分布参数估计方法 | |
| CN104373236A (zh) | 用于进气氧感测的静态压力和动态压力补偿 | |
| CN102288933B (zh) | 一种用于低于10us/cm的溶液电导率测量的校准方法 | |
| CN105259415B (zh) | 电池内阻计算方法 | |
| CN110657821B (zh) | 一种轨道检测惯性测量装置的标定方法 | |
| CN110133348A (zh) | 一种电能表自热误差的补偿方法、系统及存储介质 | |
| CN107196292B (zh) | 一种主动配电网静态等值方法及装置 | |
| CN113281692B (zh) | 一种直流互感器校验仪整检装置的误差自校准方法及系统 | |
| CN103926561B (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 | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |