一种基于邻井测井数据钻前预测当前井目标地层孔隙压力的方法技术

本发明专利技术属于石油与天然气工程领域，提供一种基于邻井测井数据钻前预测当前井目标地层孔隙压力的方法，首先依据邻井测井数据计算邻井地层孔隙压力，然后划分当前井目标地层的深度范围，最后结合邻井与当前井的几何方位图，计算当前井目标地层各深度范围内的地层孔隙压力，并通过与现场试验对比可知该方法预测精度超过95％，较常规方法预测精度提高约20％。本发明专利技术还提供一种测试当前井目标地层孔隙压力的现场试验方法，通过与开采阶段的试油测试方法相比，吻合率为95.9％，解决了钻井阶段难以测试目标地层孔隙压力的问题。以测试目标地层孔隙压力的问题。以测试目标地层孔隙压力的问题。

【技术实现步骤摘要】
一种基于邻井测井数据钻前预测当前井目标地层孔隙压力的方法


[0001]本专利技术属于石油与天然气工程领域，涉及一种基于邻井测井数据实现钻前预测当前井目标地层孔隙压力的方法。

技术介绍

[0002]地层孔隙压力是油气钻井中的关键参数，通常简称为地层压力。准确预测地层孔隙压力可以极大降低井喷、井漏等钻井风险，有效缩短工程周期节省钻井成本。然而，现阶段地层孔隙压力的钻前预测体系尚不成熟，严重依赖于地震层速度数据，考虑到地震层速度分辨精度不够(远低于测井的0.1m分辨精度)，使得地层孔隙压力的钻前预测始终难以达到现场要求。在同一井场内，已钻井测井数据分辨率足够高(一般为0.1m的分辨精度)且对当前井目标地层孔隙压力的钻前预测具有重要价值，鉴于此迫切需要提出一种基于邻井测井数据钻前预测当前井目标地层孔隙压力的方法。
[0003]现阶段已提出不少方法预测当前井目标地层孔隙压力，如一种适用于碳酸盐岩地层的地层孔隙流体压力预测方法(CN202210040987.4)、基于地震数据的三维地层压力预测方法及装置(CN202111674358.9)、一种地层压力预测方法及岩性油藏评价方法(CN201910996100.7)、一种使用地震数据预测地层压力系数的方法、装置及设备(CN202110292026.8)、地层压力预测分析方法、装置、介质及设备(CN202110763391.2)、一种基于地震资料的地层孔隙压力预测方法(CN201510246679.7)。然而，上述方法可分成两类：一类是利用测井资料在钻后分析地层孔隙压力，显然钻后分析地层孔隙压力难以指导现场实际钻井；另一类是利用地震资料对地层孔隙压力进行钻前分析，然而这类方法未考虑地震层速度的分辨精度，因此获取的地层孔隙压力准确性存在问题。为此，提出一种基于邻井测井数据钻前预测当前井目标地层孔隙压力的方法已势在必行。

技术实现思路

[0004]本专利技术提供一种基于邻井测井数据钻前预测当前井目标地层孔隙压力的方法，首先依据邻井测井数据计算邻井地层孔隙压力，然后划分当前井目标地层的深度范围，最后结合邻井与当前井的几何方位图，计算当前井目标地层各深度范围内的地层孔隙压力，并通过现场试验验证所提方法的有效性。
[0005]为实现上述目的，本专利技术第一方面提供一种基于邻井测井数据钻前预测当前井目标地层孔隙压力的方法，其操作步骤是：
[0006](1)从同一井场内筛选当前井的所有有效邻井，有效邻井数目应大于等于5；
[0007]有效邻井的筛选应同时满足下述5个条件：
[0008]1)每口有效邻井应与当前井在深度方向上地层层序一致；
[0009]2)每口有效邻井应包含目标地层且在目标地层上具有密度测井数据和声波时差数据；
[0010]3)每口有效邻井的目标地层深度范围应与当前井的目标地层深度范围有交集；由于地质构造原因，有效邻井的目标地层深度范围与当前井的目标地层深度范围一般不同，每口有效邻井的目标地层深度范围应与当前井目标地层的深度范围相交，优选的，实际预测时，以相交5m及以上来预测当前井的地层孔隙压力；
[0011]4)所有有效邻井在目标地层上的深度范围构成一个集合，该集合的上下限应包含当前井目标地层的深度范围；
[0012]5)任意两口有效邻井的连线不能通过当前井，如若通过，应舍弃掉离当前井水平距离较远的有效邻井；
[0013](2)有效邻井的地层孔隙压力预测
[0014]1)建立有效邻井正常压实趋势线
[0015]根据每口有效邻井的岩性确认各有效邻井的纯泥岩地层，筛选出纯泥岩地层在正常压实状态下的声波时差，即从纯泥岩地层中找出“随深度增加而声波时差降低”的各深度位置处的声波时差，然后利用Matlab将筛选出的声波时差散点数据拟合成下式(1)所述的对数曲线，其中，因变量为Δt
n
，自变量为H，建立每口有效邻井的正常压实趋势线，见下式：
[0016]lnΔt
n
＝lnΔt0‑
CH
ꢀꢀ
(1)
[0017]式中，Δt
n
为某深度处的声波时差，μs/m；Δt0为深度为零时的声波时差，μs/m；C为压实系数，无量纲；H为深度，m；
[0018]2)获取每口有效邻井的地层孔隙压力
[0019]将密度测井数据和声波时差数据代入Eaton公式，即可获取每口有效邻井的地层孔隙压力，见式(2)：
[0020][0021]式中：P
P
为地层孔隙压力，MPa；P0为上覆地层垂向压力，对密度沿深度方向积分获取，MPa；ρ为密度，由密度测井获取，g/cm3；P
h
为静水压力，对水密度沿深度方向积分获取，MPa；Δt
n
为正常压实趋势线上的声波时差，通过公式(1)的正常压实趋势线方程获取，μs/m；Δt为实际测井资料中的声波时差，μs/m；N为Eaton指数，依托同一井场内有实测地层孔隙压力的指定井获取，将指定井的实测地层孔隙压力及相对应深度处的上覆地层垂向压力、静水压力、正常趋势线上的声波时差和实测声波时差代入下式(3)计算Eaton指数：
[0022][0023]式中：P
p
'为指定井的地层孔隙压力实测值，MPa；P0'为指定井地层孔隙压力实测深度处的上覆地层垂向压力，由指定井的密度沿深度方向积分获取，MPa；P
h
'为静水压力，由水密度沿深度方向积分获取，MPa；Δt'
n
为指定井正常压实趋势线上的声波时差，通过公式(1)的正常压实趋势线方程获取，μs/m；Δt'为指定井地层孔隙压力实测深度处的声波时差实测值，μs/m；
[0024](3)划分当前井目标地层的深度范围
[0025]假设当前井目标地层的深度范围为[h,h']，划分方法如下：
[0026]1)定义有效邻井与当前井在目标地层上重合的深度范围为有效深度范围；
[0027]2)假设共有i口有效邻井，标记各有效邻井的有效深度范围为[h1,h2]、[h3,h4]、
……
、[h
2i
‑1,h
2i
]，将深度h1、h2、
……
、h
2i
按大小排序(其中相等的变量仅保留一个)并重新标记为a1(a1＝h)、a2、
……
、a
u
(a
u
＝h')；
[0028]3)将深度最小值a1和第二小值a2构成一个新的深度范围[a1,a2]，观察[a1,a2]与各有效邻井的有效深度范围[h1,h2]、[h3,h4]、
……
、[h
2i
‑1,h
2i
]的包含关系，将有效深度范围包含[a1,a2]的各有效邻井筛选出，构成第1井组；
[0029]4)将深度第二小值a2和第三小值a3构成一个新的深度范围[a2,a3]，观察[a2,a3]与各有效邻本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于邻井测井数据钻前预测当前井目标地层孔隙压力的方法，其特征在于，具体操作步骤是：(1)从同一井场内筛选当前井的所有有效邻井，共有i口有效邻井，i≥5；(2)有效邻井的地层孔隙压力预测；(3)划分当前井目标地层的深度范围假设当前井目标地层的深度范围为[h,h
′
]，划分方法如下：1)定义有效邻井与当前井在目标地层上重合的深度范围为有效深度范围；2)标记各有效邻井的有效深度范围为[h1,h2]、[h3,h4]、
……
、[h
2i
‑1,h
2i
]，将深度h1、h2、
……
、h
2i
按大小排序并重新标记为a1、a2、
……
、a
u
，其中深度数值相等的变量仅保留一个，a1＝h，a
u
＝h
′
；3)将深度最小值a1和第二小值a2构成一个新的深度范围[a1,a2]，观察[a1,a2]与各有效邻井的有效深度范围[h1,h2]、[h3,h4]、
……
、[h
2i
‑1,h
2i
]的包含关系，将有效深度范围包含[a1,a2]的所有有效邻井筛选出，构成第1井组；4)将深度第二小值a2和第三小值a3构成一个新的深度范围[a2,a3]，观察[a2,a3]与各有效邻井的有效深度范围[h1,h2]、[h3,h4]、
……
、[h
2i
‑1,h
2i
]的包含关系，将有效深度范围包含[a2,a3]的所有有效邻井筛选出，构成第2井组；5)以此类推，直到将深度第二大值a
u
‑1和最大值a
u
构成一个新的深度范围[a
u
‑1,a
u
]，观察[a
u
‑1,a
u
]与各有效邻井的有效深度范围[h1,h2]、[h3,h4]、
……
、[h
2i
‑1,h
2i
]的包含关系，将有效深度范围包含[a
u
‑1,a
u
]的所有有效邻井筛选出，构成第u
‑
1井组；(4)构建各井组内有效邻井与当前井的几何方位图以第j井组为例，假设第j井组由k口有效邻井构成，这k口有效邻井与当前井的几何方位图的构建方法如下：1)从有效邻井中选取关键邻井关键邻井选取的条件为：关键邻井与当前井的连线应能将井场划分为两个区域，每个区域内有效邻井的数量相等，将这两个区域命名为区域1、区域2；2)命名各区域内的有效邻井区域1中按照与关键邻井水平距离的远近，由近及远将各有效邻井依次命名为1号邻井、3号邻井
…
k
‑
2号邻井；区域2中按照与关键邻井水平距离的远近，由近及远将各有效邻井依次命名为2号邻井、4号邻井、
…
、k
‑
1号邻井；区域1和区域2内的k
‑
1口有效邻井与关键邻井共同构成第j井组的k口有效邻井；也就是说，区域1中与关键邻井水平距离最近的为1号邻井，与关键邻井水平距离第2近的为3号邻井，依此类推；区域2中与关键邻井水平距离最近的为2号邻井，与关键邻井水平距离第2近的为4号邻井，依此类推；3)将区域1与区域2内的相邻有效邻井连线并测量连线的距离将1号邻井与2号邻井、3号邻井与4号邻井、
…
、k
‑
2号邻井与k
‑
1号邻井分别连线，共有条连线，这条连线将分别与关键邻井和当前井的连线相交，交点定义为1号虚拟井、2号虚拟井、
…
、号虚拟井，测量各虚拟井与相应两口邻井间的距离，以1号虚拟井为例，需要测量1号虚拟井与1号邻井的距离L1、1号虚拟井与2号邻井的距离L2，依次类推，分别
记为L1、L2…
L
k
‑1；4)测量关键邻井与当前井连线上的各相邻井之间的距离测量关键邻井与距离关键邻井最近的虚拟井之间的距离为d1、各相邻虚拟井之间的距离为d2、d3…
当前井与距离当前井最近的虚拟井之间的距离为(5)计算当前井目标地层各深度范围的地层孔隙压力以第j井组为例，假设第j井组由k口有效邻井构成，对应的当前井目标地层的深度范围为[a
j
,a
j+1
]，则当前井目标地层[a
j
,a
j+1
]深度范围内地层孔隙压力的计算方法如下：
①
由步骤(2)获取第j井组内有效邻井的地层孔隙压力，然后结合虚拟井与相应两口邻井之间的距离，计算各虚拟井在[a
j
,a
j+1
]深度范围的地层孔隙压力，见下式：式中：D
虚拟1
为1号虚拟井的地层孔隙压力，MPa；D
虚拟2
为2号虚拟井的地层孔隙压力，MPa；
……
；为号虚拟井的地层孔隙压力，MPa；D1为1号邻井的地层孔隙压力，MPa；D2为2号邻井的地层孔隙压力，MPa；
……
；D
k
‑1为k
‑
1号邻井的地层孔隙压力，MPa；L1为1号虚拟井与1号邻井的距离，m；L2为1号虚拟井与2号邻井的距离，m；
……
；L
k
‑1为号虚拟井与k
‑
1号邻井的距离，m；
②
考虑虚拟井与当前井之间的距离关系，当前井目标地层在[a
j
,a
j+1
]深度范围的地层孔隙压力表示为下式：式中：P
j
为当前井目标地层在[a
j
,a
j+1
]深度范围内的地层孔隙压力，MPa；D
a
为关键邻井的地层孔隙压力，MPa；D
虚拟1
为1号虚拟井的地层孔隙压力，MPa；D
虚拟2
为2号虚拟井的地层孔隙压力，MPa；
……
；为号虚拟井的地层孔隙压力，MPa；d1为关键邻井与1号虚拟井之间的距离，m；d2为1号虚拟井与2号虚拟井之间的距离，m；
……
；为号虚拟井与当前井之间的距离，m；重复上述方法，即可获得当前井目标地层在各深度范围[a1,a2]、[a2,a3]、
……
[a
j
,a
j+1
]...

【专利技术属性】
技术研发人员：张立松，杨清纯，王伟，山丽洁，孙致学，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：