CN101942992A - 一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法 - Google Patents
一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法 Download PDFInfo
- Publication number
- CN101942992A CN101942992A CN201010257163XA CN201010257163A CN101942992A CN 101942992 A CN101942992 A CN 101942992A CN 201010257163X A CN201010257163X A CN 201010257163XA CN 201010257163 A CN201010257163 A CN 201010257163A CN 101942992 A CN101942992 A CN 101942992A
- Authority
- CN
- China
- Prior art keywords
- pore pressure
- pressure
- curvature
- centerdot
- point
- 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.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法,该方法利用地质构造面曲率预测区域高压盐水层孔隙压力是通过地震资料和实钻录井资料确定盐膏层的区域分布规律和构造等值线图,利用调和趋势面法建立地质构造等值线图上的区域高程方程并确定构造区域任意点的主曲率,计算任意点的主应力,建立任意点的孔隙压力预测模型,确定高压盐水层的孔隙压力,以便在钻井设计确定现场施工时为确定安全钻井液密度提供科学依据,以有效阻止井壁坍塌、防止井下复杂情况的发生。
Description
技术领域
本发明涉及一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法。
背景技术
我国大部分油气资源集中在盐下构造,盐岩地层钻井是制约我国石油钻探的关键技术难题。由于地层受构造挤压变形的作用,形成异常高压盐水层,区域分布无规律,钻井过程中经常与高压盐水层发生遭遇战,当钻井液的密度不能平衡地层压力时,并会发生盐水将进入井筒,污染钻井液,产生钻井事故和复杂情况,给钻井作业的人力、物力带来不同程度的损失。为此如何在钻井前事先预测高压盐水层孔隙压力是非常必要的,若在钻井前能预测到高压盐水层孔隙压力,就可以为钻井设计确定现场施工时安全钻井液密度提供科学的依据,以有效阻止井壁失稳、防止井下复杂情况的发生。
为此本发明中的创作人凭借其多年从事相关行业的经验与实践,并经潜心研究与开发,终于创造出一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法
发明内容
本发明的目的在于提供一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法,利用该方法可以在钻井前预测到构造区域高压盐水层孔隙压力,以便在钻井设计确定现场施工时为确定安全钻井液密度提供科学依据,以有效阻止井壁失稳、防止井下复杂情况的发生。
本发明中利用地质构造面曲率预测区域高压盐水层孔隙压力的方法,包括有下列步骤:
1)根据地震反射剖面和实钻盐层层段,找出区域盐层的顶界和底界,获得区域分布;
2)利用趋势面法计算构造面底部的曲率,获得每一空间点的曲率,计算主曲率(ρ1,ρ2);
3)根据测井资料回归统计上覆岩层压力σv;
4)根据测井资料计算盐间砂岩夹层的弹性模量和泊松比;
5)利用已知盐层构造面底部的主曲率、弹性模量和泊松比,确定所在地层的主应力,建立孔隙压力预测模型;
6)计算孔隙压力,根据每一空间点的孔隙压力,作压力等值线,获得压力分布规律。
所述步骤2)中的方法是由盐层构造等值线图上的地理坐标(xi,yi)及相应点的高程数据wi利用调和趋势面法建立w=w(x,y)的构造图的高程方程,由方程w=w(x,y)计算构造图上任意一点(xi,yi)的主曲率(ρ1i,ρ2i)。
所述步骤5)中的方法是由计算的盐岩构造底部任一点的主曲率(ρ1i,ρ2i)结合相应地层的弹性模量、泊松比计算该点的主应力(σ1i,σ2i),即获得预测地层任意一点的三向主应力序列(σ1i,σ2i,σvi),利用岩石构造变形过程中任意点岩石体积变形Δvi近似为岩石孔隙变形Δvki,得出盐岩构造曲面上任意一点孔隙压力模型pi=f(ρ1i,ρ2i,h)。
本发明中利用地质构造面曲率预测区域高压盐水层孔隙压力是通过地震资料和实钻录井资料确定盐膏层的区域分布规律和构造等值线图,利用调和趋势面法建立地质构造等值线图上的区域高程方程并确定构造区域任意点的主曲率,计算任意点的主应力,建立任意点的孔隙压力预测模型,确定高压盐水层的孔隙压力,以便在钻井设计确定现场施工时为确定安全钻井液密度提供科学依据,以有效阻止井壁坍塌、防止井下复杂情况的发生。
附图说明
图1是预测区域地质构造图;
图2是预测区域高压盐水层孔隙压力分布图。
具体实施方式
下面将结合附图对本发明中的具体实施例作进一步详细说明。
在漫长的地质年代里,由于地壳的褶皱运动使盐层发生弯曲,或由于基底***使盐层上拱,盐间封闭的盐水层受构造挤压、抬升和盐岩塑形流动的作用形成异常高压,构造变形的程度与盐水层异常压力具有一定的联系。一般来说,构造变形越剧烈,异常高压盐水层孔隙压力越大,构造变形的程度可以用构造面上任意点曲率来描述,而曲率与岩石的力学本构变形有着内在联系。岩石本构关系是揭示岩石受力与变形的定量化描述,其中包含着如地层弹性参数、地层成分、密度、埋藏深度、地质年代、孔隙率、构造运动等因素。地质构造曲率能不同程度地反映地层所受外界作用的程度,因此利用地质构造面曲率可以进行钻前预测区域高压盐水层孔隙压力。
本发明中利用地质构造面曲率预测区域高压盐水层孔隙压力包括下列步骤:
1.区域盐层分布规律的确定
在预测区域根据地震反射剖面和实钻录井资料确定区域盐层分布规律和构造等值线分布图。
即首先确定区域地震反射剖面上盐顶和盐底的深度,同时根据实钻录井资料获得的盐顶和盐底的深度进行标定,建立区域的盐层分布规律。
2.确定盐层地质构造面的主曲率
根据划定的盐层区域的分布规律和获得的盐层构造等值线图1,获取地层构造等值线图上的任意点的地理坐标(xi,yi)及相应点的高程数据wi,利用调和趋势面法建立构造图的高程方程w=w(x,y),由方程w=w(x,y)计算构造图上任意点(xi,yi)的主曲率(ρ1i,ρ2i),具体方法如下:
为了计算上的方便,设构造图高程方程满足一阶傅立叶级数趋势面方程,具体形式如下:
w=a00+a10A1C0+a01A0C1+a11A1C1+b10B1C0 (2.1)
+b11B1C1+c01A0D1+c11A1D1+d11B1D1
一阶傅立叶级数趋势面方程中有9个特定系数;a00,a10,a01,a11,b10,b11,c01,c11,d11。式中:
为求出方程中的待定系数,可按最小二乘法原理,是每个待定系数对观测值与趋势值的离差平方和
的偏导数等于0,即
上述9个方程式经整理可得到一阶傅立叶级数趋势面的正规方程组,可写成如下矩阵形式
由式(2.1)可解出式(2.1)中的待定系数。
根据构造图上的各点坐标(xi,yi)及高程数据wi,可建立正归方程
式中N为总点数,w为海拔高程,解出c00,c10,c01,c20,c11…等系数,即得w(x,y)的近似表达式。
根据薄板小挠度弯曲理论,由于w是微小的,薄板中面在x和y方向的曲率及扭率可近似的表示为:
确定。其中,αi为主方向和x轴的夹角。构造曲面上任意点主曲率的计算公式为
即
3.根据测井资料回归统计上覆岩层压力σv
密度测井和声波测井可以直观地反映地层压实规律,并能获得岩石体积密度值。由于密度测井资料容易受到井径的大小及仪器检测程度影响,在利用密度测井数据之前,应结合井径测井资料过滤掉非真实数据,以获得比较可靠的密度,并利用这些密度散点数据,利用式(3.1)计算上覆岩层压力梯度:
式中为任意点i的上覆岩层压力梯度,g/cm3;ρw,Hw分别为地层水的密度及水深,g/cm3,m;ρ0,H0分别为上部无密度测井数据段的平均密度及井深,g/cm3,m;ρi,dhi为密度测井数据及与其对应的测井层段厚度,g/cm3,m。
由测井数据计算得到的密度散点求出上覆岩层压力梯度数据以后,通过数据回归,即可获得上覆岩层压力的统计规律。考虑到幂律和二项式回归模型精度较低的缺陷,采用式(3.2)形式的三次多项式进行回归:
Gz=a0+a1h+a2h2+a3h3 (3.2)
式中Gz为上覆岩层压力梯度,g/cm3;h为考察点的深度,m;a0,a1,a2,a3为待定回归系数。即构造曲面上任意点的上覆岩层压力为:
σvi=GZhi (3.3)
4.根据测井资料计算预测地层的弹性模量和泊松比
地层的弹性模量和泊松比由测井数据解释获得,具体步骤如下:
式中:
式中:a1、b1、a2和b2为系数,视具体区域取值。
5.建立高压盐水层孔隙压力预测模型
由计算的盐岩构造底部任一点的主曲率(ρ1i,ρ2i)结合相应盐岩地层的弹性模量、泊松比计算该点的主应力(σ1i,σ2i),即获得盐层任意一点的三向主应力序列(σ1i,σ2i,σvi),利用岩石构造变形过程中任意点岩石体积变形Δvi近似为岩石孔隙变形Δvki,得出盐岩构造曲面上任意一点孔隙压力模型pi=f(ρ1i,ρ2i,h)。
具体步骤如下:
(1)盐岩地质构造曲面上任意点三向主应力计算
根据薄板小挠度弯曲理论,板内一点的应变为:
式中:w=w(x,y)为薄板的挠度。
由于式(2.3)和(5.1)可得
则任意点的应力分量为
式中,E、μ分别是弹性模量和泊松比。
(2)高压盐水层孔隙压力预测模型的建立
岩石受应力作用后,一部分由岩石孔隙中的流体承受(即孔隙压力),一部分由岩石的骨架承受(即有效应力)。若岩石在封闭不排水情况下任意点受到平均应力为σi的三向应力作用,其平均有效应力σ0i与孔隙压力Ppi的关系为:
根据弹性理论,在应力作用下任意点岩石的体积变化为:
孔隙中流体发生的体积变化为:
式中,ni为岩石孔隙率,Kvi为孔隙的体积压缩系数。
由于成岩以后的岩石骨架体积压缩很小,岩石的体积变化近似等于孔隙的体积变化,则ΔVi=ΔVVi,于是:
式中,Bi为孔隙压力系数。
联立式(3.3)、(5.6)、(5.7)和(5.8)得高压盐水层孔隙压力预测模型:
6.计算高压盐水层孔隙压力
将待预测区域的地质构造面上的任意点的曲率带入步骤5)中的预测模型,如某并4464.72m构造曲率0.0208,水平最大、最小和上覆岩层地应力分别为116.06、83.48和112.94MPa,计算出该处孔隙压力1.86MPa/100m,实测该点压力数值压力数值1.92MPa/100m,满足工程需要。根据每一空间点的孔隙压力,作压力等值线,获得压力分布规律。
Claims (3)
1.一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法,包括下列步骤:
1)根据地震反射剖面和实钻盐层层段,找出区域盐层的顶界和底界,获得区域分布;
2)利用趋势面法计算构造面底部的曲率,获得每一空间点的曲率,计算主曲率(ρ1,ρ2);
3)根据测井资料回归统计上覆岩层压力σv;
4)根据测井资料计算盐间砂岩夹层的弹性模量和泊松比;
5)利用已知盐层构造面底部的主曲率、弹性模量和泊松比,确定所在地层的主应力,建立孔隙压力预测模型;
6)计算孔隙压力,根据每一空间点的孔隙压力,作压力等值线,获得压力分布规律。
2.根据权利要求1所述的利用地质构造面曲率预测区域高压盐水层孔隙压力的方法,其特征在于:所述步骤2)中的方法是由盐层构造等值线图上的地理坐标(xi,yi)及相应点的高程数据wi利用调和趋势面法建立w=w(x,y)的构造图的高程方程,由方程w=w(x,y)计算构造图上任意一点(xi,yi)的主曲率(ρ1i,ρ2i)。
3.根据权利要求1所述的利用地质构造面曲率预测区域高压盐水层孔隙压力的方法,其特征在于:所述步骤5)中的方法是由计算的盐岩构造底部任一点的主曲率(ρ1i,ρ2i)结合相应地层的弹性模量、泊松比计算该点的主应力(σ1i,σ2i),即获得预测地层任意一点的三向主应力序列(σ1i,σ2i,σvi),利用岩石构造变形过程中任意点岩石体积变形Δvi近似为岩石孔隙变形Δvki,得出盐岩构造曲面上任意一点孔隙压力模型pi=f(ρ1i,ρ2i,h)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010257163XA CN101942992B (zh) | 2010-08-19 | 2010-08-19 | 一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010257163XA CN101942992B (zh) | 2010-08-19 | 2010-08-19 | 一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101942992A true CN101942992A (zh) | 2011-01-12 |
CN101942992B CN101942992B (zh) | 2013-11-20 |
Family
ID=43435174
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010257163XA Expired - Fee Related CN101942992B (zh) | 2010-08-19 | 2010-08-19 | 一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101942992B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102128025A (zh) * | 2011-01-20 | 2011-07-20 | 中国石油天然气股份有限公司 | 一种盐膏层造斜钻井液密度确定方法 |
CN102182453A (zh) * | 2011-03-03 | 2011-09-14 | 中国石油集团钻井工程技术研究院 | 井壁坍塌分析方法 |
CN102253412A (zh) * | 2011-06-02 | 2011-11-23 | 中国石油天然气股份有限公司 | 一种利用趋势异常识别有利储层的方法及*** |
CN103390112A (zh) * | 2013-07-29 | 2013-11-13 | 中国石油大学(华东) | 考虑煤层地应力不确定性的钻井液密度设计方法 |
CN106368691A (zh) * | 2015-07-24 | 2017-02-01 | 中国石油化工股份有限公司 | 基于岩石物理地震信息三维异常孔隙压力预测方法 |
WO2017206158A1 (en) * | 2016-06-03 | 2017-12-07 | Schlumberger Technology Corporation | Pore pressure prediction |
CN108072901A (zh) * | 2016-11-10 | 2018-05-25 | 中国石油化工股份有限公司 | 一种获得准确的静水压力和上覆地层压力的方法及*** |
CN108645552A (zh) * | 2018-05-28 | 2018-10-12 | 中国地质大学(北京) | 一种基于双应力***的地应力测量方法 |
CN108801538A (zh) * | 2018-06-15 | 2018-11-13 | 中国石油大学(华东) | 一种计算压力衰减梯度的方法 |
CN110857626A (zh) * | 2018-08-14 | 2020-03-03 | 中国石油天然气股份有限公司 | 基于综合录井参数的随钻压力预测方法、装置及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3766994A (en) * | 1971-10-01 | 1973-10-23 | Continental Oil Co | Abnormal pressure detection during drilling of a well |
CN1966934A (zh) * | 2005-11-16 | 2007-05-23 | 中国石油大学(北京) | 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法 |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
WO2008052163A1 (en) * | 2006-10-27 | 2008-05-02 | Schlumberger Canada Limited | Sanding advisor |
-
2010
- 2010-08-19 CN CN201010257163XA patent/CN101942992B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3766994A (en) * | 1971-10-01 | 1973-10-23 | Continental Oil Co | Abnormal pressure detection during drilling of a well |
CN1966934A (zh) * | 2005-11-16 | 2007-05-23 | 中国石油大学(北京) | 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法 |
CN101025084A (zh) * | 2006-02-20 | 2007-08-29 | 中国石油大学(北京) | 一种随钻预测钻头底下地层孔隙压力的方法 |
WO2008052163A1 (en) * | 2006-10-27 | 2008-05-02 | Schlumberger Canada Limited | Sanding advisor |
Non-Patent Citations (3)
Title |
---|
曹润荣等: "面曲率法-趋势面拟合法在裂缝预测中的应用", 《石油工业计算机应用》, vol. 16, no. 3, 20 September 2008 (2008-09-20) * |
曾锦光等: "应用构造面主曲率研究油气藏裂缝问题", 《力学学报》, no. 2, 19 March 1982 (1982-03-19) * |
李志勇等: "裂缝预测主曲率法的新探索", 《石油勘探与开发》, vol. 30, no. 6, 30 December 2003 (2003-12-30) * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102128025A (zh) * | 2011-01-20 | 2011-07-20 | 中国石油天然气股份有限公司 | 一种盐膏层造斜钻井液密度确定方法 |
CN102128025B (zh) * | 2011-01-20 | 2013-03-13 | 中国石油天然气股份有限公司 | 一种盐膏层造斜钻井液密度确定方法 |
CN102182453A (zh) * | 2011-03-03 | 2011-09-14 | 中国石油集团钻井工程技术研究院 | 井壁坍塌分析方法 |
CN102182453B (zh) * | 2011-03-03 | 2014-01-29 | 中国石油集团钻井工程技术研究院 | 井壁坍塌分析方法 |
CN102253412A (zh) * | 2011-06-02 | 2011-11-23 | 中国石油天然气股份有限公司 | 一种利用趋势异常识别有利储层的方法及*** |
CN102253412B (zh) * | 2011-06-02 | 2013-07-10 | 中国石油天然气股份有限公司 | 一种利用趋势异常识别有利储层的方法及*** |
CN103390112A (zh) * | 2013-07-29 | 2013-11-13 | 中国石油大学(华东) | 考虑煤层地应力不确定性的钻井液密度设计方法 |
CN103390112B (zh) * | 2013-07-29 | 2016-08-10 | 中国石油大学(华东) | 考虑煤层地应力不确定性的钻井液密度设计方法 |
CN106368691A (zh) * | 2015-07-24 | 2017-02-01 | 中国石油化工股份有限公司 | 基于岩石物理地震信息三维异常孔隙压力预测方法 |
CN106368691B (zh) * | 2015-07-24 | 2021-05-28 | 中国石油化工股份有限公司 | 基于岩石物理地震信息三维异常孔隙压力预测方法 |
WO2017206158A1 (en) * | 2016-06-03 | 2017-12-07 | Schlumberger Technology Corporation | Pore pressure prediction |
CN108072901A (zh) * | 2016-11-10 | 2018-05-25 | 中国石油化工股份有限公司 | 一种获得准确的静水压力和上覆地层压力的方法及*** |
CN108645552A (zh) * | 2018-05-28 | 2018-10-12 | 中国地质大学(北京) | 一种基于双应力***的地应力测量方法 |
CN108801538A (zh) * | 2018-06-15 | 2018-11-13 | 中国石油大学(华东) | 一种计算压力衰减梯度的方法 |
CN108801538B (zh) * | 2018-06-15 | 2020-10-30 | 中国石油大学(华东) | 一种计算压力衰减梯度的方法 |
CN110857626A (zh) * | 2018-08-14 | 2020-03-03 | 中国石油天然气股份有限公司 | 基于综合录井参数的随钻压力预测方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN101942992B (zh) | 2013-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101942992B (zh) | 一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法 | |
CN105672971B (zh) | 一种储层裂缝开启压力、开启次序及注水压力预测方法 | |
CN107526892B (zh) | 一种海洋天然气水合物试采储层的稳定性评估方法 | |
CN104500050B (zh) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 | |
CN113821953B (zh) | 基于动态地应力和剩余油特征的暂堵重复压裂设计方法 | |
Carreón-Freyre et al. | Shearing along faults and stratigraphic joints controlled by land subsidence in the Valley of Queretaro, Mexico | |
CN104992468A (zh) | 缝洞型碳酸盐岩油气藏三维地质建模方法 | |
EP3084124A1 (en) | Method for determining hydraulic fracture orientation and dimension | |
Finol et al. | Numerical simulation of oil production with simultaneous ground subsidence | |
CN101806224A (zh) | 井下煤层气抽采产能预测方法 | |
Ringrose et al. | Building confidence in CO2 storage using reference datasets from demonstration projects | |
Harris | The role of geology in reservoir simulation studies | |
CN107742020A (zh) | 页岩气储层压裂改造体积预测方法 | |
CN103425877A (zh) | 基于多孔介质流固耦合模型的地质灾害预测***和方法 | |
CN105182412A (zh) | 煤矿地下水库煤柱坝体地震安全性的检测方法 | |
Ma et al. | An early warning method for water inrush in Dongjiahe coal mine based on microseismic moment tensor | |
CN106199712B (zh) | 一种确定压裂套管变形区域的方法及装置 | |
CN110705168A (zh) | 构造应力场的模拟方法 | |
Zhu et al. | Simulation analysis of influencing factors of subsidence based on mining under huge loose strata: a case study of Heze mining area, China | |
CN102562052A (zh) | 密井网浅层套损危害体识别方法 | |
Chengxuan et al. | An approach to the present-day three-dimensional (3D) stress field and its application in hydrocarbon migration and accumulation in the Zhangqiang depression, Liaohe field, China | |
CN114004157B (zh) | 一种基于盆地地质大数据的古构造特征动态恢复方法 | |
CN115935588A (zh) | 一种小井间距双井盐穴储气库的稳定性评价方法 | |
CN108459357A (zh) | 地层欠压实与流体膨胀超压的评价方法 | |
Zeng | Research on risk evaluation methods of groundwater bursting from aquifers underlying coal seams and applications to coalfields of North China |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131120 Termination date: 20160819 |
|
CF01 | Termination of patent right due to non-payment of annual fee |