CN102220865B - 一种检测灰岩地层孔隙压力的方法 - Google Patents
一种检测灰岩地层孔隙压力的方法 Download PDFInfo
- Publication number
- CN102220865B CN102220865B CN201110124596.2A CN201110124596A CN102220865B CN 102220865 B CN102220865 B CN 102220865B CN 201110124596 A CN201110124596 A CN 201110124596A CN 102220865 B CN102220865 B CN 102220865B
- Authority
- CN
- China
- Prior art keywords
- formation pore
- pore pressure
- pressure
- model
- effective stress
- 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
Abstract
检测灰岩地层孔隙压力的方法,应用于油田钻井技术领域。包括:采集同一地质构造井岩样岩石力学参数数据、建立有效应力模型、建立上覆岩层压力模型、建立灰岩地层孔隙压力的检测模型、检测灰岩地层孔隙压力。效果是:与现场实测的地层孔隙压力数据进行对比,符合程度较高,趋势符合性较好,基本解决了传统方法不适用检测灰岩地层孔隙压力的难题。
Description
技术领域
本发明涉及油田钻井技术领域,特别涉及一种井下地层压力检测方法。
背景技术
在灰岩地层钻井过程中,如不能准确检测灰岩地层坍塌压力、破裂压力和地层孔隙压力,钻井过程中将会发生井漏、井塌等复杂事故。
灰岩地层的沉积压实机理不同于砂泥岩地层的压实理论;目前地层孔隙压力的检测模型基本都是采用砂泥岩地层孔隙压力检测模型,砂泥岩地层压力检测模型是基于压实理论建立的,利用砂泥岩地层压力检测模型就不能合理的解释灰岩地层孔隙压力,造成计算结果与实际检测的结果误差较大,有效应力定理不受沉积机理的限制,基于有效应力定理建立灰岩地层压力检测模型就是建立有效应力与某些岩石力学参数间的关系模型。
检索到相似的专利共4篇,其中与自己申报的技术相似的专利是:申请号是200610007768.7,专利名称是一种随钻预测钻头底下地层孔隙压力的方法。将自己申报的专利与该专利进行对比,主要有两个不同点:
(1)利用的测井数据不同:
(2)建立的地层孔隙压力模型和上覆岩层压力预测模型不同:
建立的上覆岩层压力模型和地层孔隙压力模型不同。该专利(申请号是200610007768.7)的方法包括以下主要内容:针对同处一区块的待钻井和已钻井,分别提取各自井旁若干道地震记录,作加权处理获得待钻井和已钻井的地震记录;利用已钻井不同层系地层的测井数据和地震记录,建立地震记录预测地层声波时差测井曲线的分层模型;预测待钻井钻头底下地层的声波时差测井曲线;利用已钻井的地层密度测井数据、声波时差测井数据和自然伽马测井数据建立该地质区块的上覆地层压力计算模型、孔隙度与泥质含量计算模型和分层段的地层孔隙压力计算模型;结合已建立的本区块的上覆地层压力计算模型、本区块的孔隙度与泥质含量计算模型和地层孔隙压力计算模型预测当前钻遇地层的孔隙压力。
经检索,灰岩地层压力检测模型主要两种:一是利用砂泥岩地层压力预测模型来预测;二是利用有效应力定理建立模型来检测。基于有效应力定理建立检测模型是检测灰岩地层孔隙压力是最有效的方法,但在基于有效应力定理建立的灰岩地层压力检测模型中没有考虑体积压缩系数对声波速度、地层孔隙压力和有效应力的影响。通过岩石力学参数实验研究发现,体积压缩系数对声波速度及地层孔隙压力的影响是不能忽略的。
发明内容
本发明的目的是:提供一种检测灰岩地层孔隙压力的方法,针对灰岩地层孔隙压力检测中遇到的一系列问题,解决砂泥岩地层孔隙压力解释模型难以检测灰岩地层孔隙压力这一难题,建立灰岩地层孔隙压力的解释模型;本发明采用有效应力定理建立灰岩地层孔隙压力检测模型,解决传统检测砂泥岩地层孔隙压力的解释模型不适用于检测灰岩地层孔隙压力的难题。
本发明中一种检测灰岩地层孔隙压力的方法包括下列步骤:包括:采集同一地质构造井岩样岩石力学参数数据、建立有效应力模型、建立上覆岩层压力模型、建立灰岩地层孔隙压力的检测模型、检测灰岩地层孔隙压力;
1)采集同一地质构造井岩样岩石力学参数数据:对处于同一地质构造的已经完成的某口取心井,根据钻井过程中的录井资料确定灰岩地层所处的深度;对取出的岩心在实验室加工成标准岩样,利用MTS实验装置,开展三轴实验,对岩石的力学参数进行测试,测试内容包括:骨架体积弹性模量Ema和纵波速度Vp,在三轴实验中,围压为20MPa,分别测定有效应力为10MPa,20MPa,30MPa,40MPa,50MPa,60MPa时的骨架体积弹性模量和纵波速度,由此得到两组岩石力学参数特征数据系列
2)对特征数据系列进行多元非线性回归,建立有效应力模型:将测得的两组岩石力学参数特征数据系列利用MATLAB软件对其分别进行多元非线性回归,得出有效应力与纵波速度,有效应力与骨架体积弹性模量间的函数关系,建立有效应力、纵波速度和骨架体积弹性模量之间的模型为:Vp=f(Pe,Ema),即
式中:a1,a2,a3,a4——回归系数,无量纲;
3)建立上覆岩层压力模型:利用同一地质构造已经完成的某口取心井的测井数据,通过地层密度与地层深度间的特征数据系列{ρ(i),h(i)},采用陈庭根主编的《钻井工程理论与技术》书中的积分公式
式中:Poi-一定深度的上覆岩层压力,MPa;
ρw,Hw-海水密度及水深,g/cm3,m;
ρ0,H0-上部无密度测井数据段平均密度及长度,g/cm3,m;
ρbi-一定深度的密度数据,g/cm3;
ΔH-积分步长,m;
由此得出上覆岩层压力Poi与地层深度间的特征数据系列{p0 (i),h(i)};
利用MATLAB软件对上覆岩层压力与地层深度间的特征数据系列{p0 (i),h(i)}进行多元非线性回归,建立上覆岩层压力模型为:P0=f(h),即Po=A+BH-Ce-DH,
式中:H——地层深度,m;
P0——上覆岩层压力,MPa;
A,B,C,D——回归系数,无量纲;
4)建立灰岩地层孔隙压力的检测模型:基于有效应力定理,建立灰岩地层孔隙压力的检测模型为:Pp=P0-Pe=f(h)-Pe,
式中:Pp——地层孔隙压力,MPa;
Pe——有效应力,MPa;
5)检测灰岩地层孔隙压力:利用建立的模型,应用计算机编程,计算已测井井段的地层孔隙压力值。
本发明专利的有益效果:本发明专利检测灰岩地层孔隙压力的方法,选取大港油田的某井,以此计算该井的灰岩地层孔隙压力剖面,并与现场实测的地层孔隙压力数据进行对比,符合程度较高,趋势符合性较好,基本解决了传统方法不适用检测灰岩地层孔隙压力的难题。
附图说明
图1是利用检测灰岩地层孔隙压力的方法,检测的井灰岩地层的密度与地层深度的测井数据图;
图2是利用检测灰岩地层孔隙压力的方法,检测的井灰岩地层的声波时差与地层深度的测井数据图;
图3是检测的地层孔隙压力剖面与实测数据对比图。
具体实施方式
实施例1:以一次检测灰岩地层孔隙压力的方法为例,对本发明专利作进一步详细说明。
本发明专利检测灰岩地层孔隙压力的方法
1.确定灰岩地层所处的深度,岩石力学参数实验
确定灰岩地层所处的深度,并在该深度范围内取心,将其加工成标准岩样,在实验室进行岩石力学参数实验,在三轴实验中,围压为20MPa,分别测定有效应力为10MPa,20MPa,30MPa,40MPa,50MPa,60MPa时的骨架体积弹性模量,和纵波速度,即得到两组岩石力学参数特征数据系列如下:
2.对特征数据系列进行多元非线性回归
根据测得的两组岩石力学参数特征数据系列,利用MATLAB软件对其分别进行多元非线性回归;
有效应力与骨架体积弹性模量间呈指数关系,基于有效应力与纵波速度建立的方程,将纵波速度与体积弹性模量间的关系方程表示为模型为:Vp=f(Ema);
3.建立上覆岩层压力模型
利用图1中的测井数据,地层密度与地层深度间的特征数据系列{ρ(i),h(i)},采用陈庭根主编的《钻井工程理论与技术》书中积分公式计算出上覆岩层压力,由此得到上覆岩层压力与地层深度间的特征数据系列{p0 (i),h(i)},如下
利用MATLAB软件对上覆岩层压力与地层深度间的特征数据系列{p0 (i),h(i)}进行多元非线性回归,建立上覆岩层压力模型为:Po=A+BH-Ce-DH,并回归出上覆岩层压力计算模型中的参数,即:A=2.1447,B=0.0434,C=0.3113,D=0.44927;
4.建立灰岩地层孔隙压力的检测模型
利用现场实际测量的某些点的地层孔隙压力值Pp,如下:
基于有效应力定理,P0=Pe+Pp,计算出对应地层深度的有效应力值,利用图2中的声波测井数据,通过利用MATLAB软件进行多元非线性回归,得到有效应力方程中的模型参数为:a1=-1.40811,a2=0.01766,a3=2.44991,a4=0.00014。
建立灰岩地层孔隙压力的检测模型为:
Po=2.1447+0.0434H-0.3113e-0.44927H
5.检测灰岩地层孔隙压力:
利用建立的模型,应用计算机编程,依据已测井段的测井数据,计算已测井井段的地层孔隙压力值,得到地层孔隙压力剖面,参阅附图3。
利用该方法,目前已完成了对该井的地层孔隙压力检测,并通过与现场实测的地层孔隙压力对比,该方法的检测精度都在12%以内,完全能够满足灰岩地层孔隙压力检测要求。参阅附图3。曲线代表该井灰岩地层孔隙压力的检测结果,两个圆点代表的是4678米和4735米的实测地层孔隙压力梯度。通过与检测结果对比发现,实测的地层压力梯度与检测的地层压力梯度的精度误差范围都在12%以内。由此证明该模型能很好的满足灰岩地层孔隙压力的检测要求。有了灰岩地层孔隙压力的检测模型后,本领域技术人员能熟练使用。
利用建立的模型,目前已完成了对该井等几口井的地层孔隙压力检测,并与现场实测的地层孔隙压力数据进行对比,较常规的地层孔隙压力检测方法相比有更高的准确度,基本解决了传统检测砂泥岩地层孔隙压力的方法不适用于检测灰岩地层孔隙压力的难题。
Claims (1)
1.一种检测灰岩地层孔隙压力的方法,其特征在于:包括:采集同一地质构造井岩样岩石力学参数数据、建立有效应力模型、建立上覆岩层压力模型、建立灰岩地层孔隙压力的检测模型、检测灰岩地层孔隙压力;
1)采集同一地质构造井岩样岩石力学参数数据:对处于同一地质构造的已经完成的某口取心井,根据钻井过程中的录井资料确定灰岩地层所处的深度;对取出的岩心在实验室加工成标准岩样,利用MTS实验装置,开展三轴实验,对岩石的力学参数进行测试,测试内容包括:骨架体积弹性模量Ema和纵波速度Vp,在三轴实验中,围压为20MPa,分别测定有效应力为10MPa,20MPa,30MPa,40MPa,50MPa,60MPa时的骨架体积弹性模量和纵波速度,由此得到两组岩石力学参数特征数据系列
2)对特征数据系列进行多元非线性回归,建立有效应力模型:将测得的两组岩石力学参数特征数据系列利用MATLAB软件对其分别进行多元非线性回归,得出有效应力与纵波速度,有效应力与骨架体积弹性模量间的函数关系,建立有效应力、纵波速度和骨架体积弹性模量之间的模型为:Vp=f(Pe,Ema),即
式中:a1,a2,a3,a4——回归系数,无量纲;
3)建立上覆岩层压力模型:利用同一地质构造已经完成的某口取心井的测井数据,通过地层密度与地层深度间的特征数据系列{ρ(i),h(i)},采用积分公式
式中:Poi-一定深度的上覆岩层压力,MPa;
ρw,Hw-海水密度及水深,g/cm3,m;
ρ0,H0-上部无密度测井数据段平均密度及长度,g/cm3,m;
ρbi-一定深度的密度数据,g/cm3;
ΔH-积分步长,m;
由此得出上覆岩层压力Poi与地层深度间的特征数据系列{p0 (i),h(i)};
利用MATLAB软件对上覆岩层压力与地层深度间的特征数据系列{p0 (i),h(i)}进行多元非线性回归,建立上覆岩层压力模型为:P0=f(h),即Po=A+BH-Ce-DH,
式中:H——地层深度,m;
P0——上覆岩层压力,MPa;
A,B,C,D——回归系数,无量纲;
4)建立灰岩地层孔隙压力的检测模型:基于有效应力定理,建立灰岩地层孔隙压力的检测模型为:Pp=P0-Pe=f(h)-Pe,
式中:Pp——地层孔隙压力,MPa;
Pe——有效应力,MPa;
5)检测灰岩地层孔隙压力:利用建立的模型,应用计算机编程,计算已测井井段的地层孔隙压力值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110124596.2A CN102220865B (zh) | 2011-05-13 | 2011-05-13 | 一种检测灰岩地层孔隙压力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110124596.2A CN102220865B (zh) | 2011-05-13 | 2011-05-13 | 一种检测灰岩地层孔隙压力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102220865A CN102220865A (zh) | 2011-10-19 |
CN102220865B true CN102220865B (zh) | 2014-04-02 |
Family
ID=44777547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110124596.2A Active CN102220865B (zh) | 2011-05-13 | 2011-05-13 | 一种检测灰岩地层孔隙压力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102220865B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10174613B2 (en) | 2014-01-14 | 2019-01-08 | Halliburton Energy Services, Inc. | Effective porosity determination for tight gas formations |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102608021B (zh) * | 2011-11-29 | 2014-04-02 | 中国石油大学(北京) | 一种砂泥岩薄互层储层压裂防砂堵工艺参数确定方法 |
CN103089253B (zh) * | 2013-01-22 | 2015-06-03 | 中国石油大学(北京) | 一种利用小波变换计算地层孔隙压力的方法 |
CN104863577B (zh) * | 2015-04-09 | 2017-07-21 | 中国石油大学(北京) | 利用地震纵波传播时间预测地层孔隙压力的方法 |
CN105484739B (zh) * | 2015-11-26 | 2018-07-06 | 中国科学院武汉岩土力学研究所 | 碳酸盐岩地层孔隙压力测试方法和装置 |
CN105569641B (zh) * | 2015-12-17 | 2021-09-28 | 西南石油大学 | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 |
CN106979006A (zh) * | 2017-05-17 | 2017-07-25 | 中国神华能源股份有限公司 | 地层压力的确定方法和装置 |
CN109509111B (zh) * | 2017-09-15 | 2021-09-28 | 中国石油化工股份有限公司 | 探井地层压力的预测方法及系统 |
CN109323954B (zh) * | 2018-11-15 | 2019-12-10 | 中国地质大学(武汉) | 一种碳酸盐岩地层孔隙压力预测方法 |
CN113128014A (zh) * | 2019-12-30 | 2021-07-16 | 中石化石油工程技术服务有限公司 | 基于灰色预测理论的钻头前方待钻地层孔隙压力预测方法 |
CN113536218B (zh) * | 2021-07-13 | 2023-06-06 | 中国石油大学(北京) | 一种计算岩石孔隙体积压缩系数的方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5937362A (en) * | 1998-02-04 | 1999-08-10 | Diamond Geoscience Research Corporation | Method for predicting pore pressure in a 3-D volume |
CN2527707Y (zh) * | 2002-02-07 | 2002-12-25 | 西南石油学院 | 高温高压孔隙介质模型 |
CN1828011A (zh) * | 2005-12-12 | 2006-09-06 | 中海油田服务股份有限公司 | 模拟油水两相电缆地层测试的方法 |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
CN101512100A (zh) * | 2006-08-07 | 2009-08-19 | 普拉德研究及开发股份有限公司 | 用于孔隙压力预测的方法及系统 |
CN201443393U (zh) * | 2009-05-15 | 2010-04-28 | 中国石油大学(华东) | 油藏流动模拟实验取压接头 |
-
2011
- 2011-05-13 CN CN201110124596.2A patent/CN102220865B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5937362A (en) * | 1998-02-04 | 1999-08-10 | Diamond Geoscience Research Corporation | Method for predicting pore pressure in a 3-D volume |
CN2527707Y (zh) * | 2002-02-07 | 2002-12-25 | 西南石油学院 | 高温高压孔隙介质模型 |
CN1828011A (zh) * | 2005-12-12 | 2006-09-06 | 中海油田服务股份有限公司 | 模拟油水两相电缆地层测试的方法 |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
CN101512100A (zh) * | 2006-08-07 | 2009-08-19 | 普拉德研究及开发股份有限公司 | 用于孔隙压力预测的方法及系统 |
CN201443393U (zh) * | 2009-05-15 | 2010-04-28 | 中国石油大学(华东) | 油藏流动模拟实验取压接头 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10174613B2 (en) | 2014-01-14 | 2019-01-08 | Halliburton Energy Services, Inc. | Effective porosity determination for tight gas formations |
US10301936B2 (en) | 2014-01-14 | 2019-05-28 | Halliburton Energy Service, Inc. | Tight gas formation pressure determination method |
Also Published As
Publication number | Publication date |
---|---|
CN102220865A (zh) | 2011-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102220865B (zh) | 一种检测灰岩地层孔隙压力的方法 | |
Lucier et al. | Geomechanical aspects of CO2 sequestration in a deep saline reservoir in the Ohio River Valley region | |
Warpinski et al. | In-situ stresses in low-permeability, nonmarine rocks | |
Clarkson et al. | Coalbed methane: current field-based evaluation methods | |
CN104453873A (zh) | 页岩油气经济有效层段的评价方法 | |
CN104278991A (zh) | 盐湖相烃源岩有机碳和生烃潜量的多元测井计算方法 | |
CN104500054B (zh) | 地层孔隙压力的确定方法及装置 | |
Djuraev et al. | A review on conceptual and practical oil and gas reservoir monitoring methods | |
Abbas et al. | Stability analysis of highly deviated boreholes to minimize drilling risks and nonproductive time | |
CN108301825A (zh) | 一种高温高压储层孔隙压力分析方法 | |
CN104374827A (zh) | 横观各向同性岩体原位动弹性模量的各向异性系数的测量方法 | |
CN104714252A (zh) | 分析流体因子敏感性的方法 | |
CN109555515A (zh) | 地层坍塌压力确定方法和装置 | |
Crawford et al. | Determining static elastic anisotropy in shales from sidewall cores: impact on stress prediction and hydraulic fracture modeling | |
CN111381292A (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
Dutta et al. | Integration of mini-frac stress results with core data to identify ductile behavior of overpressured sandstone reservoir in deep-water offshore Krishna-Godavari basin, India | |
CN110500089B (zh) | 基于阵列感应测井资料的页岩气水平井层理裂缝评价方法 | |
Mahmoud et al. | A Self-Adaptive Artificial Neural Network Technique to Estimate Static Young's Modulus Based on Well Logs | |
Roshan et al. | Borehole deformation based in situ stress estimation using televiewer data | |
CN104484573A (zh) | 一种地层的刚性系数的计算新方法 | |
Bakker et al. | Multiscale evaluation of potential damage in jetted lateral boreholes | |
Tan et al. | Leak-off mechanism and pressure prediction for shallow sediments in deepwater drilling | |
Elkhoury et al. | The First Pressuremeter Testing Campaign on Wireline Formation Testers in Deep Boreholes | |
Sharma et al. | Demarcating sweet spots in cambay shale by integrating rock eval pyrolysis, geomechanics and seismic data | |
Davies et al. | Determination of geomechanical properties of a typical Niger Delta reservoir rock using geophysical well logs |
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 |