CN101025084B - 一种随钻预测钻头底下地层孔隙压力的方法 - Google Patents
一种随钻预测钻头底下地层孔隙压力的方法 Download PDFInfo
- Publication number
- CN101025084B CN101025084B CN2006100077687A CN200610007768A CN101025084B CN 101025084 B CN101025084 B CN 101025084B CN 2006100077687 A CN2006100077687 A CN 2006100077687A CN 200610007768 A CN200610007768 A CN 200610007768A CN 101025084 B CN101025084 B CN 101025084B
- Authority
- CN
- China
- Prior art keywords
- degree
- drilling
- model
- porosity
- pore pressure
- 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.)
- Expired - Fee Related
Links
Images
Abstract
一种利用地震资料随钻预测钻头底下地层孔隙压力的方法,包括下列步骤:针对同处一区块的待钻井和已钻井,分别提取各自井旁若干道地震记录,作加权处理获得待钻井和已钻井的地震记录;利用已钻井不同层系地层的测井数据和地震记录,建立地震记录预测地层声波时差测井曲线的分层模型;预测待钻井钻头底下地层的声波时差测井曲线;利用已钻井的地层密度测井数据、声波时差测井数据和自然伽马测井数据建立该地质区块的上覆地层压力计算模型、孔隙度与泥质含量计算模型和分层段的地层孔隙压力计算模型;结合已建立的本区块的上覆地层压力计算模型、本区块的孔隙度与泥质含量计算模型和地层孔隙压力计算模型预测当前钻遇地层的孔隙压力。
Description
技术领域
本发明涉及一种随钻预测地层孔隙压力的方法,更具体地说,涉及一种利用地震资料随钻预测地层孔隙压力的方法。
背景技术
在石油钻井过程中,不确定性的异常地层孔隙压力经常造成井下作业复杂,且容易发生事故,给钻井作业的人力、物力带来不同程度的损失。为此如何在钻井前事先预测地层孔隙压力是非常必要的,因若在钻井前能预测到地层孔隙压力就可以科学设计套管程序,合理选择钻井液密度,以实现快速、安全、高效、低成本钻井和完井。
为此本发明中的创作人凭借其多年从事相关行业的经验与实践,并经潜心研究与开发,终于创造出一种利用地震资料来随钻预测地层孔隙压力的方法。
发明内容
本发明的主要目的在于提供一种利用地震资料随钻预测钻头底下地层孔隙压力的方法。
本发明中利用地震资料随钻预测邻井钻头底下地层孔隙压力的方法,包括下列步骤:
步骤1,针对同处一区块的待钻井和已钻井,分别提取各自井旁若干道地震记录,作加权处理获得待钻井和已钻井的地震记录;
步骤2,利用已钻井不同层系地层的测井数据和地震记录,建立地震记录预测地层声波时差测井曲线的分层模型;
步骤3,预测待钻井钻头底下地层的声波时差测井曲线;
步骤4,利用已钻井的地层密度测井数据、声波时差测井数据和自然伽马测井数据建立该地质区块的上覆地层压力计算模型、孔隙度与泥质含量计算模型和分层段的地层孔隙压力计算模型;
步骤5,结合已建立的本区块的上覆地层压力计算模型、该层系的孔隙度与泥质含量计算模型和地层孔隙压力计算模型预测当前钻遇地层的孔隙压力。
所述步骤2包括:
步骤21,提取已钻井的声波时差测井数据、地层密度测井数据和自然伽马测井数据;
步骤22,提取已钻井相邻若干道地震记录,作加权处理获得已钻井地震记录;
步骤23,查阅已钻井的地质录井资料,获取准确的地质分层数据(地层地质分层的顶部深度和底部深度);
步骤24,取处理完毕的已钻井每一地质分层层段的地震记录;
步骤25,分层段求取已钻井的关联维数、R/S分数维、振幅谱维、最大李雅普洛夫指数、突变参数和综合反射地震参数序列;
步骤26,结合已钻井地质资料和测井资料,得到已钻井分层段的测井曲线;
步骤27,以关联维数、R/S分数维、振幅谱维、最大李雅普洛夫指数、突变参数和综合反射地震参数为输入层的神经元,以声波时差值为输出层的神经元,分层段创建地震参数与声波时差关系的反向传播神经网络模型;
步骤28,分层段训练神经网络模型,获得分层段预测测井曲线的反向传播神经网络模型。
所述步骤3包括:
步骤31,根据地质录井数据,确定当前钻井钻头所处的地质层段;
步骤32,选取该层段的声波时差测井曲线预测的神经网络模型;
步骤33,将该层段的地震参数代入已训练的神经网络模型,预测出该层系井段的声波时差测井曲线。
所述步骤4包括:
步骤41,建立区域上覆地层压力梯度计算模型;
步骤42,建立区域孔隙度和泥质含量计算模型;
步骤43,建立区域分层孔隙压力计算模型。
所述步骤41包括:
1)利用密度测井资料计算不同深度的上覆地层压力梯度
其中G0i为某深度的上覆地层压力梯度(单位:g/ml),ρi,Δhi分别为该深度上部各层的密度平均值(单位:g/ml)和层间深度间隔(单位:km),ρ0,h0分别为上部无密度测井数据层段的平均密度(单位:g/ml)和厚度(单位:km),
2)利用下式建立上覆地层压力梯度和深度的非线性关系模型
其中,G0为上覆地层压力梯度值(单位:g/ml),h为测井深度值(单位: km),a1,a2,a3,a4为模型参数。
所述步骤42包括:
利用下式建立孔隙度、泥质含量与纵横波速度的非线性关系模型
φ=a1+a2vp+a3vs
Vsh=b1+b2vp+b3vs
其中φ为孔隙度,由声波测井资料计算得到;Vsh为泥质含量,由自然伽马测井资料计算得到;a1,a2,a3,b1,b2,b3为模型参数,vs为横波速度(单位km/s),可以利用纵波速度vp由下面的经验关系确定
纵波速度可以直接由声波时差测井数据计算出
式中,vp为纵波速度(单位:km/s),Δt为声波时差测井值(单位:μs/ft)。
所述步骤43包括:
1)建立纵波速度、泥质含量、孔隙度以及垂直有效应力之间非线性关系的分层模型
上式中vp为纵波速度(单位:km/s),Vsh为泥质含量,φ为孔隙度,Pe为垂直有效应力(单位MPa),a1,a2,a3,a4,a5为模型参数。
2)利用有效应力定律预测地层孔隙压力
Pp=Po-Pe
其中,Pe为垂直有效应力(单位MPa),Po为上覆地层压力(单位MPa),Pp为地层孔隙压力(单位MPa)。
本发明中利用地震资料随钻预测邻井地层孔隙压力的方法主要是找出同 一井眼中不同层系地层的已钻井段测井结果与地震特征参数的关联,而根据地质分层随钻预测地层孔隙压力,以有效防止井下复杂情况的发生。
附图说明
图1是已钻井通过相邻若干道地震记录加权处理获得的地震记录;
图2是利用已钻井资料建立的该地质区块地震双程旅行时和测井深度之间的对应关系;
图3是已钻井第三系安居安组井段关联维的计算结果图;
图4为已钻井第三系安居安组井段R/S分数维的计算结果图;
图5为已钻井第三系安居安组井段振幅谱维的计算结果图;
图6为已钻井第三系安居安组井段最大李雅普洛夫指数的计算结果图;
图7为已钻井第三系安居安组井段突变参数的计算结果图;
图8为已钻井第三系安居安组井段综合反射地震参数的计算结果图;
图9为根据已钻井测井资料和地质资料提取的已钻井第三系安居安组井段的声波时差测井曲线;
图10是待钻井通过相邻若干道地震记录加权处理获得的地震记录;
图11是待钻井第三系安居安组井段关联维的计算结果图;
图12为待钻井第三系安居安组井段R/S分数维的计算结果图;
图13为待钻井第三系安居安组井段振幅谱维的计算结果图;
图14为待钻井第三系安居安组井段最大李雅普洛夫指数的计算结果图;
图15为待钻井第三系安居安组井段突变参数的计算结果图;
图16为已钻井第三系安居安组井段综合反射地震参数的计算结果图;
图17为待钻井所预测的第三系安居安组井段的声波时差测井曲线;
图18是根据已钻井密度测井资料建立的该地质区块的上覆地层压力计算模型;
图19为待钻井第三系安居安组井段的地层孔隙压力;
图20为本发明预测钻头底下地层孔隙压力的方法流程图。
具体实施方式
以下参照附图20对本发明中利用地震记录随钻预测井眼待钻井段孔隙压力的方法进行详细说明,主要包括下列步骤:
步骤1,对已钻井和待钻井分别提取各自井旁若干道地震记录,作加权处理获得已钻井和待钻井地震记录(图1)。
步骤2,对已钻井进行声波时差、中子密度测井和自然伽马测井,获得声波时差测井曲线、地层密度测井曲线和自然伽马测井曲线。
步骤3,建立地震特征参数与测井数据间的关系,即建立利用地震记录构建测井曲线的模型:
地震波的反射特征与地下介质的波阻抗有关,波阻抗的自然对数对时间增量的导数定义为反射率函数,地震记录的反射率函数深刻反映了地震的线性与非线性特征;而地层的声波速度直接决定着反射率函数。因此利用地震的线性与非线性特征可以映射地震的反射率函数,也就是说地震的线性与非线性特征通过适当的映射模型可以反映地层的声波特征。
其中包括:
1)提取已钻井的测井数据:声波时差测井数据和地层密度测井数据;
2)查阅已钻井的地质录井资料,获取准确的地质分层数据:地层地质分层的顶部深度和底部深度;
3)取已处理完毕的已钻井的每一地质分层层段的地震记录;
4)利用所述每一地质分层层段的地震记录,分层段求取地震特征参数;其中,所述地震特征参数包括:关联维数、R/S分数维、振幅谱维、最大李雅普洛夫指数、突变参数和综合反射地震参数;其中,
a.用声波测井时差进行深时转换,其转换公式为
式中为H0起始深度,H为终止深度,T(H)为深度H所对应的时间,Δt(h)为深度h处的声波测井的时差值。将测井的深度采样点变成了时间序列,然后根据地震记录时间采样间隔对声波、密度测井值进行重采样,将深度采样值转换为时间采样值,以此时深关系(图2)为标准来实现地震和测井数据在时间和深度上的对应,将地震记录的时间刻度转换为深度刻度,得到分层段的地震记录,用于提取地震特征参数。
b.关联维数
从时间序列直接计算关联维数的算法,即G-P算法。
第一步是将地震记录的时间序列xi=x(ti)i=1,2,......,N重新排列创立一个m维的向量相空间
Xn(m,τ)=(xn,xn+τ,......,xn+(m-1)τ) n=1,2,......,Nm (2)
式中τ=kΔt是固定时间间隔,即时间延迟,Δt是两次相邻采样的间隔,k是整数
Nm=N-(m-1)τ (3)
从这Nm个点中任意选定一个参考点Xi,计算其余Nm-1个点到Xi的距离
对所有点重复这一过程,得到关联积分函数
式中H是Heaviside函数
在lgCm(r)-lgr的坐标图上,如点分布在一条直线附近,则分形存在,采用直线拟合的方法,拟合直线的斜率就是关联维D2。图3为关联维的计算结果图。
c.R/S分数维
重标度极差分析方法(R/S分析法)是Hurst于1965年提出的一种新的统计方法,它在分形理论中有着广泛的应用。
设已知地震记录的时间序列为xi=x(ti),i=1,2,......,N,则τ个时间数据的均值为
由此可以求得累积离差
极差
标准差
根据Hurst分析的统计规律的关系式为
R/S∝(τ/2)H (10)
具体计算时,计算出R/S的一组值。在lg(R/S)-lg(τ/2)的坐标图上拟合 出一条直线,其斜率就是R/S分数维。
图4为R/S分数维的计算结果图。
d.振幅谱维
对地震记录的时间序列xi=x(ti),i=1,2,......,N作功率谱分析,功率谱S(f)随着f的变化可以表示为S(f)∝fβ。因为振幅谱 所以有A(f)∝fβ/2。令Df=β/2,则有A(f)∝fDf。绘制一张lgA(f)-lgf的图,然后用直线拟合,其斜率即为振幅谱维Df。
图5为振幅谱维的计算结果图。
e.李雅普洛夫指数
李雅普洛夫指数是指相空间中邻近轨道发散或收敛的平均指数率,反映了系统的性态对于初值的敏感程度。
一般计算最大李雅普洛夫指数采用A.Wolf重构法,其步骤如下:
A.重建相空间
设已知地震记录的时间序列为xi=x(ti),i=1,2,......,N,i为采样点的序号,xi为在该时间下的振幅值,输入嵌入空间维数m及延迟时间τ得相空间:
Xn(m,τ)=(xn,xn+τ,......,xn+(m-1)τ) n=1,2,......,Nm (11)
Nm=N-(m-1)τ (12)
B.在相空间中,以初始点A(t0)为参考点,选取A(t0)的最近邻点B(t0),设在t1=t0+kΔt时,A(t0)和B(t0)分别演化到A(t1)和B(t1),计算从t0到t1时的指数增长率
C.在A(t1)的若干邻近点中,找出一个与A(t1)的夹角θ1很小的邻近点C(t1),如果找不到,仍然选取B(t1),设在t2=t1+kΔt时,A(t1)和C(t1)分别演化到A(t2)和C(t2),则
将这一过程一直进行到点集的终点,然后取的平均值作为最大李雅普洛夫指数的估计值LE(m)。
D.增加嵌入空间维数m,重复A-C步,直到LE(m)保持平稳为止,此时的LE(m)即为所求的最大李雅普洛夫指数。
图6为最大李雅普洛夫指数的计算结果图。
f.突变参数
将地震信号看成对时间变量的连续函数x(t),x(t)可展成级数形式如下
y=x(t)=a0+a1 t+a2t2+…+antn+… (14)
式中,t为时间,y为对应t的位移,a0,a1,a2,…an为待定的系数。实际分析发现,对具有一定趋势规律的时间序列,截取到4次项时,精度已足够高。则对上式可近似表示为
y=x(t)=a0+a1 t+a2t2+a3t3+a4t4 (15)
对上式作变量代换,化为尖点突变的标准形式,先令
t=zt-q (16)
将其代入上式,可得
b0=a4q4-a3q3+a2q2-a1q+a0
式中, b1=-4a4q3+3a3q2-2a2q+a1
b2=6a4q2-3a3q+a2
b4=a4
上面仍不是尖点突变的标准形式,作进一步变量代换,令
则
式中:c=b0
z为状态变量,a,b为控制变量,由突变理论可知,平衡曲面方程为
z3+az+b=0 (20)
分叉集方程为
4a3+27b2=0 (21)
只有在控制变量满足分叉集方程时,系统才是不稳定的,才有可能从一个平衡态突变到另一个平衡态。
图7为突变参数的计算结果图。
g.综合反射地震参数
将常规的反射地震参数,即瞬时振幅、瞬时频率、瞬时相位、自相关函数极小值、自相关函数极大值以及五阶自回归模型系数共十种参数进行综合评判,优选出一种综合反射参数参与建模。
图8为综合反射地震参数的计算结果图。
5)根据已钻井的地质资料和测井资料,提取已钻井分层段的声波时差测井曲线(图9)。
6)将以上提取的各种地震特征参数以及对应的声波时差测井数据输入神经网络,分层创建反向传播神经网络模型。
步骤4,利用待钻井段的地震特征参数预测待钻井段的声波时差测并曲线:
1)根据地质录井数据,确定待钻井钻头所处的地质层段;
2)选取该层段的声波时差测井曲线神经网络模型;
3)将根据待钻井的地震记录(图10)计算出的待钻井段的关联维数(图11)、R/S分数维(图12)、振幅谱维(图13)、最大李雅普洛夫指数(图14)、突变参数(图15)和综合反射地震参数(图16)代入已训练好的对应层段神经网络模型,预测出待钻井在不同深度地层的声波时差测井曲线(图17)。
步骤5,利用已钻井的地层密度测井数据建立该地质区块的上覆地层压力计算模型:
上覆地层压力计算模型首先要求确定上覆地层压力梯度,上覆地层压力梯度可以由密度测井的散点数据计算得到
上式中,G0i为某深度的上覆地层压力梯度(单位:g/ml),ρi,Δhi分别为该深度上部各层的密度平均值(单位:g/ml)和层间深度间隔(单位:km),ρ0,h0分别为上部无密度测井数据层段的平均密度(单位:g/ml)和厚度(单位:km)。
由上覆地层压力梯度可计算各深度的上覆地层压力
P0i=0.01G0ihi (23)
上式中,P0i为地下某点上覆地层压力(单位:MPa),G0i为该点上覆地层压力梯度(单位:g/ml),hi为该点深度(单位:m)。
通常利用已钻井密度测井的散点数据将上覆地层压力梯度回归为深度的函数,研究认为将上覆地层压力梯度回归为下式效果最好
为了克服通常密度测井数据不足的问题,在只有n个密度测井散点数据的情况下,可以建立变量W和G0的关系
Wi=0 (i=1)
根据以上关系,将W回归为深度h的函数
利用上式进行非线性回归,得到a1,a2,a3,a4代入(24)式中进行计算即可以得到随深度变化的上覆地层压力梯度剖面(图18),利用该剖面可以方便地计算该地质区块任一深度的上覆地层压力。
步骤6,利用已钻井的声波时差测井曲线和自然伽马测井曲线建立该地质区块的孔隙度和泥质含量计算模型:
大量实验证明,声波速度、孔隙度和泥质含量之间存在良好的线性关系,这种关系随地质区块的差异而有所不同,所以可以将声波速度、孔隙度和泥质含量代入下式,建立适合于本区块的孔隙度和泥质含量计算模型
φ=a1+a2vp+a3vs (27)
Vsh=b1+b2vp+b3vs
其中φ为孔隙度,Vsh为泥质含量,vs为横波速度(单位km/s),可以利用纵波速度vp由下面的经验关系确定
纵波速度可以直接由声波时差测井数据计算出
式中,vp为纵波速度(单位:km/s),Δt为声波时差测井值(单位:μs/ft),
在获取已钻井的声波时差测井数据和自然伽马测井数据的基础上,就可以利用测井解释的结论来求取孔隙度和泥质含量。通常选取声波时差测井数据由下式可以计算孔隙度
式中,Δt,Δtma,Δtf分别为目的层、岩石骨架和地层流体的声波时差测井值(单位:μs/ft)。
在获取自然伽马测井数据的基础上利用相对值法计算泥质含量
式中,Vsh为泥质含量,GCUR为希尔奇指数,与地质时代有关,可根据取芯分析资料和自然伽马测井值统计确定,一般对于第三系地层取3.7,老地层取2。IGR为泥质含量指数, GR,GRmin,GRmax分别表示目的层、纯泥岩层和纯砂岩层的自然伽马测井值(API)。
将计算出的声波速度、孔隙度和泥质含量按照(27)式通过多元线性回归的方法确定a1,a2,a3,b1,b2,b3,从而得到该区块的孔隙度和泥质含量计算模型。
步骤7,利用已钻井资料建立该区块分层段的地层孔隙压力计算模型:
地层孔隙压力预测的理论基础是饱和多孔介质的有效应力定律
Pe=Po-Pp (32)
其中,Pe为垂直有效应力(单位MPa),Po为上覆地层压力(单位MPa),Pp为地层孔隙压力(单位MPa)。从上式可以知道,在已知上覆地层压力的情况下,只要求得垂直有效应力就可以计算出地层孔隙压力,所以建立地层孔隙压力计算模型实际上就是建立垂直有效应力计算模型。
大量的岩心测试实验已经证明,影响声波在岩石中传播速度的主要因素
是岩性、孔隙度和垂直有效应力,所以可以将纵波速度表示为岩性、孔隙度和垂直有效应力的函数
vp=f(Vsh,φ,Pe) (33)
上式中vp为纵波速度(单位:km/s),Vsh为泥质含量,φ为孔隙度,Pe为垂直有效应力(单位MPa)。
实验工作表明建立以下模型来表征上面的函数关系最为理想
在建立分层垂直有效应力计算模型时,vp的值由声波时差测井资料直接计算,孔隙度和泥质含量的计算利用已建立的孔隙度和泥质含量分层计算模型。
利用已建立的该区块的上覆地层压力模型和静液压力或实测压力就确定用于建模的分层垂直有效应力
Pe=Po-gw h (35)
式中,gw表示静液压力梯度,一般取0.0105MPa/m。
将vp,φ,Vsh,Pe分层数据代入式(34)中,利用多元非线性回归的方法计算出a1,a2,a3,a4,a5,从而建立起该地质区块分层的地层孔隙压力计算模型。
在实际应用中还经常使用以下简化的模型,即忽略孔隙度和泥质含量对垂直有效应力的影响,而只考虑声波速度对垂直有效应力的影响,在这种情况下,就无需建立孔隙度和泥质含量计算模型,这就大大简化了计算程序。这种简化计算模型如下式所示
步骤8,利用预测得到待钻井段的声波时差测井数据结合已建立的区域上覆地层压力计算模型、孔隙度与泥质含量计算模型和分层段的地层孔隙压力 计算模型计算当前钻遇层段的地层孔隙压力:
实际钻进中,通过岩屑录井资料和分层神经网络预测模型随钻预测出钻头下地层的声波时差测井数据,将其转换为声波速度,结合分层孔隙度与泥质含量计算模型计算待钻地层的孔隙度与泥质含量,再利用分层地层孔隙压力计算模型反算待钻地层的垂直有效应力,最后利用区域上覆地层压力计算模型计算待钻地层的孔隙压力(图19)。
上述实施例仅用于说明本发明,而并非用于限定本发明。
Claims (8)
1.一种利用地震资料随钻预测钻头底下地层孔隙压力的方法,其特征在于,包括下列步骤:
步骤1,针对同处一区块的待钻井和已钻井,分别提取各自井旁若干道地震记录,作加权处理获得待钻井和已钻井的地震记录,提取已钻井的声波时差测井数据、地层密度测井数据和自然伽马测井数据;
步骤2,利用已钻井不同层系地层的测井数据和地震记录,建立地震记录预测地层声波时差测井曲线的分层模型;
步骤3,预测待钻井钻头底下地层的声波时差测井曲线;
步骤4,利用已钻井的地层密度测井数据、声波时差测井数据和自然伽马测井数据建立该地质区块的上覆地层压力计算模型、孔隙度与泥质含量计算模型和分层段的地层孔隙压力计算模型;
步骤5,结合上述本区块的上覆地层压力计算模型、本区块的孔隙度与泥质含量计算模型和地层孔隙压力计算模型预测当前钻遇地层的孔隙压力。
2.根据权利要求1所述的随钻预测钻头底下地层孔隙压力的方法,其特征在于,所述步骤2包括:
所述步骤2包括:
步骤21,查阅已钻井的地质录井资料,获取准确的地质分层数据,包括地层地质分层的顶部深度和底部深度;
步骤22,将地震记录的时间刻度转换为深度刻度,求得已钻井每一地质分层段的地震记录;
步骤23,根据分层段的地震记录求取已钻井的地震特征参数;
步骤24,结合已钻井地质资料和测井资料,得到已钻井分层段的声波时差测井曲线;
步骤25,以地震特征参数为输入层的神经元,以声波时差值为输出层的神经元,分层段创建地震特征参数与声波时差关系的反向传播神经网络模型。
3.根据权利要求2所述的随钻预测钻头底下地层孔隙压力的方法,其特征在于,所述步骤3包括:
步骤31,根据地质录井数据,确定当前钻井钻头所处的地质层段;
步骤32,选取该层段的声波时差测井曲线预测的神经网络模型;
步骤33,将该层段的地震特征参数代入上述的神经网络模型,预测出该层系井段的声波时差测井曲线。
4.根据权利要求2或3所述的随钻预测钻头底下地层孔隙压力的方法,其特征在于,所述地震特征参数包括:关联维数、R/S分数维、振幅谱维、最大李雅普洛夫指数、突变参数和综合反射地震参数。
5.根据权利要求3所述的随钻预测钻头底下地层孔隙压力的方法,其特征在于,所述步骤4包括:
步骤41,利用已钻井的地质密度测井数据将上覆地层压力梯度回归为深度的函数,建立地质区域上覆地层压力梯度计算模型;
步骤42,利用已钻井的声波时差测井曲线和自然伽马测井曲线确定声波速度、孔隙度和泥质含量的双线形关系模型,建立区域孔隙度和泥质含量计算模型;
步骤43,利用已钻井的测井资料、上述区域孔隙度和泥质含量计算模型以及纵波速度与泥质含量、孔隙度和垂直有效应力之间的函数关系,建立垂直有效应力计算模型,再利用饱和多孔介质的有效应力定律、所述垂直有效应力计算模型以及上述上覆地层压力梯度计算模型,建立区域分层孔隙压力计算模型。
7.根据权利要求5所述的随钻预测钻头底下地层孔隙压力的方法,其特征在于,所述步骤42包括:
1)利用下式建立孔隙度、泥质含量与纵横波速度之间的双线性关系模型
φ=a1+a2vp+a3vs
Vsh=b1+b2vp+b3vs
式中,φ为孔隙度,由声波测井资料计算得到;Vsh为泥质含量,由自然伽马测井资料计算得到;a1,a2,a3,b1,b2,b3为模型参数,vs为横波速度,单位:km/svp为纵波速度单位:km/s;
2)vs根据经验关系确定
vp由声波时差测井数据计算出
式中,vp为纵波速度,单位:km/s,Δt为声波时差测井值,单位:μs/ft;由声波时差测井数据计算孔隙度
式中,Δt,Δtma,Δtf分别为目的层、岩石骨架和地层流体的声波时差测井值,单位:μs/ft;
由自然伽马测井数据计算泥质含量
3)利用计算出的声波速度、孔隙度和泥质含量通过多元线性回归的方法确定a1,a2,a3,b1,b2,b3,得到该地质区域的孔隙度和泥质含量计算模型。
8.根据权利要求5所述的随钻预测钻头底下地层孔隙压力的方法,其特征在于,所述步骤43包括:
1)建立纵波速度、泥质含量、孔隙度以及垂直有效应力之间非线性关系的分层模型
式中,vp为纵波速度,单位:km/s,Vsh为泥质含量,φ为孔隙度,Pe为垂直有效应力,单位:MPa,a1,a2,a3,a4,a5为模型参数;
其中,利用上述的上覆地层压力模型和静液压力求得分层垂直有效应力值
Pe=Po-gw h
式中,Po为上覆地层压力,h为深度,gw表示静液压力梯度,取0.0105MPa/m;
由声波时差测井资料计算得出纵波速度,利用上述孔隙度和泥质含量计算模型计算得出孔隙度和泥质含量;
将计算得出的分层垂直有效应力、纵波速度、孔隙度和泥质含量代入上式中利用多元非线性回归的方法计算得出a1,a2,a3,a4,a5;
2)利用上述垂直有效应力计算模型和上述上覆地层压力模型,根据有效应力定律
Pp=Po-Pe
其中,Pe为垂直有效应力,单位:MPa,Po为上覆地层压力,单位:MPa,Pp为地层孔隙压力,单位:MPa,得出该地质区域的地层孔隙压力计算模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2006100077687A CN101025084B (zh) | 2006-02-20 | 2006-02-20 | 一种随钻预测钻头底下地层孔隙压力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2006100077687A CN101025084B (zh) | 2006-02-20 | 2006-02-20 | 一种随钻预测钻头底下地层孔隙压力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101025084A CN101025084A (zh) | 2007-08-29 |
CN101025084B true CN101025084B (zh) | 2011-04-13 |
Family
ID=38743667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2006100077687A Expired - Fee Related CN101025084B (zh) | 2006-02-20 | 2006-02-20 | 一种随钻预测钻头底下地层孔隙压力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101025084B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2739875C1 (ru) * | 2020-05-12 | 2020-12-29 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Уфимский государственный нефтяной технический университет" | Способ определения работоспособности породоразрушающего инструмента |
Families Citing this family (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101942992B (zh) * | 2010-08-19 | 2013-11-20 | 中国石油大学(北京) | 一种利用地质构造面曲率预测区域高压盐水层孔隙压力的方法 |
CN101963056B (zh) * | 2010-08-19 | 2014-04-09 | 中国石油大学(北京) | 一种利用测井资料预测碳酸盐岩地层孔隙压力的方法 |
CN101936157B (zh) * | 2010-08-19 | 2013-06-05 | 中国石油大学(北京) | 一种利用测井资料检测高压盐水层孔隙压力的方法 |
CN102220865B (zh) * | 2011-05-13 | 2014-04-02 | 中国石油天然气股份有限公司 | 一种检测灰岩地层孔隙压力的方法 |
CN102252948B (zh) * | 2011-06-10 | 2014-05-14 | 中国石油天然气股份有限公司 | 泥页岩孔隙度测定方法 |
CN102426390B (zh) * | 2011-10-21 | 2013-07-03 | 中国石油大学(北京) | 一种非均质泥砂岩储层储量确定方法 |
CN103018778B (zh) * | 2012-12-12 | 2015-10-14 | 中国石油天然气股份有限公司 | 校正声波测井曲线进行储层预测的方法及设备 |
CN103257363B (zh) * | 2013-05-21 | 2015-09-30 | 西南石油大学 | 一种探测地下裂缝性储层中裂缝倾角的方法 |
CN103306672B (zh) * | 2013-05-24 | 2016-04-06 | 中国石油大学(北京) | 一种预测页岩地层不同钻进方向的研磨性的方法 |
CN103573250B (zh) * | 2013-07-22 | 2016-08-10 | 中国石油天然气股份有限公司 | 一种计算水平井井眼到地层上下界面距离的方法 |
CN103485758A (zh) * | 2013-08-12 | 2014-01-01 | 中国石油天然气股份有限公司 | 一种利用多种资料拟合反演复杂岩性解释方法及装置 |
CN103615236A (zh) * | 2013-11-26 | 2014-03-05 | 中国石油集团川庆钻探工程有限公司 | 一种远程录井信息实时监测地层压力的方法 |
CN104036119B (zh) * | 2014-05-16 | 2017-05-03 | 陕西延长石油(集团)有限责任公司研究院 | 划分沉积地层的方法 |
CN104863577B (zh) * | 2015-04-09 | 2017-07-21 | 中国石油大学(北京) | 利用地震纵波传播时间预测地层孔隙压力的方法 |
CN106285642B (zh) * | 2015-05-14 | 2020-01-10 | 中国石油化工股份有限公司 | 一种基于地震资料的地层孔隙压力预测方法 |
CN106368691B (zh) * | 2015-07-24 | 2021-05-28 | 中国石油化工股份有限公司 | 基于岩石物理地震信息三维异常孔隙压力预测方法 |
CN105116466B (zh) * | 2015-07-30 | 2017-12-19 | 中国石油天然气股份有限公司 | 一种确定油田储层物性特征的方法和装置 |
CN105298478A (zh) * | 2015-09-08 | 2016-02-03 | 中国石油大学(北京) | 一种用于断裂构造的地层孔隙压力确定方法 |
CN106934075B (zh) * | 2015-12-29 | 2020-07-10 | 中国石油化工股份有限公司 | 钻井液密度确定方法及静态当量密度确定方法 |
CN106321090B (zh) * | 2016-08-25 | 2019-10-29 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | 一种盐间地层孔隙压力的预测方法 |
CN107025360A (zh) * | 2017-04-21 | 2017-08-08 | 中国石油大港油田勘探开发研究院 | 一种利用测井资料恢复古地层压力方法 |
CN107575219B (zh) * | 2017-09-15 | 2020-08-07 | 中石化石油工程技术服务有限公司 | 一种页岩气储层地层破裂压力梯度计算方法 |
CN108089227B (zh) * | 2017-12-12 | 2020-05-29 | 中国石油天然气集团有限公司 | 新的基于三维地震资料的地层孔隙压力预测方法 |
CN108304959B (zh) * | 2017-12-21 | 2022-09-16 | 中国石油化工股份有限公司 | 提高地层流体压力预测精度的方法 |
CN109736784B (zh) * | 2018-04-27 | 2019-11-05 | 长江大学 | 沉积岩地层孔隙压力预测计算方法 |
CN110857626B (zh) * | 2018-08-14 | 2022-11-04 | 中国石油天然气股份有限公司 | 基于综合录井参数的随钻压力预测方法、装置及存储介质 |
CN109061729B (zh) * | 2018-08-22 | 2019-11-12 | 西安石油大学 | 一种高温高压气藏含气敏感曲线重构方法 |
CN109165737B (zh) * | 2018-08-29 | 2020-09-01 | 电子科技大学 | 基于条件随机场和bp神经网络的孔隙度预测方法 |
US11852010B2 (en) | 2018-09-14 | 2023-12-26 | Eni S.P.A. | Method for estimating a pore pressure value in geological formations to be drilled by a drilling apparatus |
CN109184644B (zh) * | 2018-09-28 | 2021-06-15 | 中国海洋石油集团有限公司 | 一种考虑聚合物非牛顿性和渗流附加阻力的早期注聚效果评价方法 |
CN111257945B (zh) * | 2018-11-30 | 2021-03-16 | 中国石油化工股份有限公司 | 一种更新已钻开井段井周地层地震速度的方法 |
CN109521221B (zh) * | 2018-12-10 | 2020-09-29 | 东北大学 | 一种钻爆法施工硬岩隧道微震波波速实时获取方法 |
CN109356573A (zh) * | 2018-12-12 | 2019-02-19 | 中法渤海地质服务有限公司 | 一种地层声波时差的测试提取方法 |
CN109498370B (zh) * | 2018-12-26 | 2021-02-19 | 杭州电子科技大学 | 基于肌电小波关联维的下肢关节角度预测方法 |
CN109458176A (zh) * | 2018-12-28 | 2019-03-12 | 西南石油大学 | 碳酸盐岩储层压力的预测方法及其应用 |
CN110376645A (zh) * | 2019-08-02 | 2019-10-25 | 中国石油大学(北京) | 一种薄层砂体顶界面位置的确定方法、装置及设备 |
CN113027427B (zh) * | 2019-12-25 | 2024-04-05 | 中石化石油工程技术服务有限公司 | 一种基于贝叶斯理论的含可信度地层压力随钻修正方法 |
CN113128014A (zh) * | 2019-12-30 | 2021-07-16 | 中石化石油工程技术服务有限公司 | 基于灰色预测理论的钻头前方待钻地层孔隙压力预测方法 |
CN111502647B (zh) * | 2020-03-27 | 2021-02-02 | 中国石油化工股份有限公司石油工程技术研究院 | 一种钻井地质环境因素确定方法、装置及存储介质 |
CN114075972A (zh) * | 2020-08-05 | 2022-02-22 | 中石化石油工程技术服务有限公司 | 基于随钻资料的钻井过程中动态井身结构优化设计方法 |
CN112230283B (zh) * | 2020-10-12 | 2021-08-10 | 北京中恒利华石油技术研究所 | 一种基于测井曲线支持向量机建模的地震孔隙度预测方法 |
CN112560246B (zh) * | 2020-12-07 | 2024-04-26 | 中海石油(中国)有限公司 | 一种目标井散点地层压力系数的预测方法 |
CN112664188B (zh) * | 2020-12-31 | 2023-09-22 | 中国海洋石油集团有限公司 | 一种基于砂体传递模式的地层孔隙压力预测方法 |
CN113553780B (zh) * | 2021-09-22 | 2021-11-26 | 西南石油大学 | 一种基于机器学习的地层孔隙压力预测方法 |
CN114492757B (zh) * | 2022-01-29 | 2023-11-10 | 中国石油化工股份有限公司 | 模型训练方法、测井信息预测方法、装置和设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5128866A (en) * | 1989-09-20 | 1992-07-07 | Chevron Corporation | Pore pressure prediction method |
US5282384A (en) * | 1992-10-05 | 1994-02-01 | Baroid Technology, Inc. | Method for calculating sedimentary rock pore pressure |
CN1588127A (zh) * | 2004-10-22 | 2005-03-02 | 石油大学(北京) | 利用地震层速度钻前预测坍塌压力与破裂压力的方法 |
CN1588128A (zh) * | 2004-10-22 | 2005-03-02 | 石油大学(北京) | 利用地震记录预测井眼待钻井段坍塌压力和破裂压力的方法 |
-
2006
- 2006-02-20 CN CN2006100077687A patent/CN101025084B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5128866A (en) * | 1989-09-20 | 1992-07-07 | Chevron Corporation | Pore pressure prediction method |
US5282384A (en) * | 1992-10-05 | 1994-02-01 | Baroid Technology, Inc. | Method for calculating sedimentary rock pore pressure |
CN1588127A (zh) * | 2004-10-22 | 2005-03-02 | 石油大学(北京) | 利用地震层速度钻前预测坍塌压力与破裂压力的方法 |
CN1588128A (zh) * | 2004-10-22 | 2005-03-02 | 石油大学(北京) | 利用地震记录预测井眼待钻井段坍塌压力和破裂压力的方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2739875C1 (ru) * | 2020-05-12 | 2020-12-29 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Уфимский государственный нефтяной технический университет" | Способ определения работоспособности породоразрушающего инструмента |
Also Published As
Publication number | Publication date |
---|---|
CN101025084A (zh) | 2007-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101025084B (zh) | 一种随钻预测钻头底下地层孔隙压力的方法 | |
CN1966934B (zh) | 一种随钻预测钻头底下地层坍塌压力和破裂压力的方法 | |
Islam et al. | Advanced petroleum reservoir simulation: Towards developing reservoir emulators | |
Abbas et al. | Drilling rate of penetration prediction of high-angled wells using artificial neural networks | |
Magara | Compaction and fluid migration | |
Zerrouki et al. | Prediction of natural fracture porosity from well log data by means of fuzzy ranking and an artificial neural network in Hassi Messaoud oil field, Algeria | |
RU2669948C2 (ru) | Оптимизация многоступенчатого проекта нефтяного месторождения в условиях неопределенности | |
US8515720B2 (en) | Determine field fractures using geomechanical forward modeling | |
Alpay | A practical approach to defining reservoir heterogeneity | |
US8768672B2 (en) | Method for predicting time-lapse seismic timeshifts by computer simulation | |
US7933750B2 (en) | Method for defining regions in reservoir simulation | |
CN106368691A (zh) | 基于岩石物理地震信息三维异常孔隙压力预测方法 | |
El Ouahed et al. | Application of artificial intelligence to characterize naturally fractured zones in Hassi Messaoud Oil Field, Algeria | |
CN112746837A (zh) | 基于分布式光纤传感的页岩油藏勘探数据采集系统及方法 | |
CN108138555A (zh) | 预测储层性质的方法、系统及设备 | |
CN102220865A (zh) | 一种检测灰岩地层孔隙压力的方法 | |
Burbey et al. | Modeling three-dimensional deformation in response to pumping of unconsolidated aquifers | |
Harris | The role of geology in reservoir simulation studies | |
CN104714252A (zh) | 分析流体因子敏感性的方法 | |
Abbas et al. | Stability analysis of highly deviated boreholes to minimize drilling risks and nonproductive time | |
Hui et al. | Investigation on two Mw 3.6 and Mw 4.1 earthquakes triggered by poroelastic effects of hydraulic fracturing operations near crooked lake, Alberta | |
CN110426751A (zh) | 一种利用测井资料预测横波时差的方法 | |
CN106979006A (zh) | 地层压力的确定方法和装置 | |
Chakravarty et al. | Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures | |
Delorme et al. | An efficient model to simulate reservoir stimulation and induced microseismic events on 3D discrete fracture network for unconventional reservoirs |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110413 Termination date: 20140220 |