一种基于测井资料进行岩屑回注选层的方法与流程

1.本发明涉及油气开发技术领域，具体涉及一种基于测井资料进行岩屑回注选层的方法。


背景技术：

2.近年来，油气开发越来越多地使用油基钻井液，这会产生大量的具有高污染性的油基岩屑。目前随着各个国家对环保的要求越来越严格，钻井废物必须达到环境要求才能排放，因此，如何高效低成本处理油基岩屑成为制约钻井进度的关键因素。
3.目前常用的处理方法有填埋法、机械分离法、生物法、热处理法及回注法。其中填埋法虽成本较低，但是由于填埋深度较浅，当遭遇雨季冲刷或地下水渗透可能导致油基废物泄漏，污染环境；机械分离法和热处理方法设备运行成本高；生物法处理周期长，且油基岩屑中的有毒物质会抑制生化降解过程；回注法是将油基岩屑研磨处理，通过高压泵将含有岩屑的浆液从环空或专门的回注井注入到目标地层，不仅可以从根本上解决油基岩屑污染环境的问题，同时处理成本低，时效性好。尤其是海上油气田，回注处理相比运回陆地处理，成本大幅降低，降低了海上平台空间的占用。但目前回注法在国外应用较多，国内尚未进行大规模应用，原因是回注目的层的筛选大都基于经验确定，缺少科学方法或指导方案，存在压力过高盖层破裂、孔隙不达标以及盖层稳定性差等风险，不仅影响回注岩屑的效率和回注量，更严重的是一旦发生泄漏，会对地层水或海洋环境造成不可逆的严重污染。


技术实现要素：

4.有鉴于此，本发明旨在提供一种基于测井资料进行岩屑回注选层的方法，以提供完整的、有科学依据的选层流程，以便有效解决经验选层的弊端。
5.本发明首先提出一种基于测井资料进行岩屑回注选层的方法，所述方法包括以下步骤：
6.s1、根据伽马测井资料计算待回注井对应地层泥质含量剖面；
7.s2、根据所述泥质含量和岩性判断标准，确定地层岩性纵向分布；
8.s3、选择所述地层岩性纵向分布中砂岩层为回注目标层，判断该目标层的深度是否符合要求；
9.s4、判断所述回注目标层的厚度是否符合要求；
10.s5、判断所述回注目标层盖层的厚度是否符合要求；
11.s6、根据测井资料，分别计算回注目标层孔隙度、渗透率、破裂压力和盖层破裂压力剖面；
12.s7、判断所述回注目标层和所述盖层的破裂压力差是否符合要求；
13.s8、根据所述目标层孔隙度、渗透率和破裂压力剖面分别计算目标层归一化孔隙度平均值，归一化渗透率平均值，归一化破裂压力平均值并确定所述回注目标层的回注潜力指数h；
14.s9：根据所述回注潜力指数h确定所述回注目标层是否合适。
15.根据本发明的一种实施方式，在步骤s8中，所述回注目标层的回注潜力指数h
16.按照下式确定：
[0017][0018]
其中a＝0.1365；b＝0.2385；c＝0.6250
[0019]
为归一化孔隙度平均值；
[0020]
为归一化渗透率平均值；
[0021]
为归一化破裂压力平均值；
[0022]
h指数是通过ahp层次分析法计算得到的一个无量纲数值，其代表的是目标地层回注潜能指数。
[0023]
根据本发明的一种实施方式，在步骤s1中，通过获取伽马测井曲线，根据以下公式计算地层泥质含量剖面，
[0024][0025]
其中：
[0026]
v
sh
‑
泥质含量，％
[0027]
gr
‑
测井伽马值
[0028]
gr
sand
‑
纯砂岩伽马值
[0029]
gr
mud
‑
纯泥岩伽马值。
[0030]
根据本发明的一种实施方式，在步骤s2中所述地层岩性纵向分布按如下标准判断，
[0031]
若泥质含量(v
sh
)≤40，且该层的厚度大于≥4m，则该层为砂岩层；
[0032]
若泥质含量(v
sh
)≤40，且该层的厚度大于<4m，则该层为砂岩夹层；
[0033]
若泥质含量(v
sh
)≥60，且该层的厚度大于≥4m，则该层为泥岩；
[0034]
若泥质含量(v
sh
)≥60，且该层的厚度大于<4m，则该层为泥岩夹层；
[0035]
若40<泥质含量(v
sh
)<60,且该层的厚度大于≥4m，则该层为泥质砂岩层。
[0036]
若40<泥质含量(v
sh
)<60,且该层的厚度大于<4m，则该层为泥质砂岩夹层。
[0037]
根据本发明的一种实施方式，在步骤s3中，结合回注作业区(陆地或海上)通过岩性纵向分布判断所选砂岩层顶部距离地表(陆地回注)或海底(海上平台回注)的距离是否大于600m，若大于600m则进行下一步选择，否则此目标层不符合回注深度要求。
[0038]
根据本发明的一种实施方式，在步骤s4中，通过所述地层岩性纵向分布判断所述回注目标层的厚度，若该回注目标层的厚度大于10m，则满足回注要求，否则此目标层不符合回注厚度要求。
[0039]
根据本发明的一种实施方式，在步骤s5中，通过所述地层岩性纵向分布判断所述回注目标层盖层的厚度，若盖层厚度大于20m，则满足回注盖层的要求，否则该盖层不满足要求。
[0040]
根据本发明的一种实施方式，在步骤s6中，通过伽马测井值、声波测井值和密度测井值分别计算所述回注目标层的孔隙度、渗透率、破裂压力以及所述盖层的破裂压力剖面。
[0041]
根据本发明的一种实施方式，在步骤s7中，若所述盖层底部破裂压力大于所述回注目标层顶部破裂压力，且两者的差值大于6.89mpa(即1000psi)，则所述盖层破裂压力满足要求，否则不符合回注要求。
[0042]
根据本发明的一种实施方式，在步骤s8中，分别将步骤s6获得的目标地层孔隙度、渗透率和破裂压力剖面做归一化处理，然后计算所述回注目标层的归一化孔隙度平均值、归一化渗透率平均值和归一化破裂压力平均值，再计算所述目标层回注潜力指数h；
[0043]
在步骤s9中，当h≥0.5，则所述回注目标层适合回注；当0.3<h<0.5
[0044]
，则所述回注目标层应谨慎回注；当h≤0.3，则所述回注目标层不适合回注；
[0045]
在步骤s3、s4、s5、s7中任一步不符合要求，则所述回注目标层不适合回注，需要重新返回第一步，选取另外目标层进行判断。
[0046]
本发明基于声波测井、伽马测井和密度测井资料，通过对回注目的层及盖层一系列参数的计算和判断，从而精确筛选回注目的层。该方法可提供完整的、有科学依据的岩屑回注选层流程，可有效解决经验选层的弊端，且计算成本低，筛选方便。
附图说明
[0047]
图1为本发明一实施例回注选层流程示意图。
具体实施方式
[0048]
以下将结合附图对本发明的较佳实施例进行详细说明，以便更清楚理解本发明的目的、特点和优点。应理解的是，附图所示的实施例并不是对本发明范围的限制，而只是为了说明本发明技术方案的实质精神。
[0049]
本发明提供一种基于常规测井资料的岩屑回注选层的方法，这种方法是基于声波测井、伽马测井和密度测井资料，通过对目的层及盖层一系列参数的计算和判断，从而精确筛选回注目的层。该方法可提供完整的、有科学依据的选层流程，可有效解决经验选层的弊端，且计算成本低，筛选方便。
[0050]
本实施方式一种基于测井资料进行岩屑回注选层的方法，通过声波测井、伽马测井和密度测井资料计算地层参数，根据筛选标准进行回注选层，如图1所示，主要包括如下步骤：
[0051]
步骤1：根据伽马测井资料计算地层泥质含量剖面；
[0052]
步骤2：根据岩性分类标准，确定地层岩性纵向分布；
[0053]
步骤3：选择砂岩层为回注目标层，判断目标层深度；
[0054]
步骤4：判断目标层厚度；
[0055]
步骤5：判断盖层厚度；
[0056]
步骤6：将测井资料代入计算公式，分别计算回注目标层孔隙度、渗透率、破裂压力和盖层破裂压力剖面；
[0057]
步骤7：判断目标层和盖层破裂压力差；
[0058]
步骤8：计算目标层归一化孔隙度平均值，归一化渗透率平均值，归一化破裂压力平均值，并代入公式计算回注潜力指数h；
[0059]
步骤9：根据选层判别标准选取回注层。
[0060]
下面通过实施方式将更具体地对每步骤进行详细介绍。
[0061]
根据本发明的一种实施方式，在具体执行上述步骤时，采取了以下实施方式。
[0062]
步骤1：获取伽马测井曲线，根据公式(1)计算地层泥质含量剖面。
[0063][0064]
其中：
[0065]
v
sh
‑
泥质含量，％
[0066]
gr
‑
测井伽马值
[0067]
gr
sand
‑
纯砂岩伽马值
[0068]
gr
mud
‑
纯泥岩伽马值
[0069]
步骤2：根据下表标准将地层进行分层，确定地层岩性纵向分布。
[0070]
地层类型泥质含量(vsh)厚度砂岩层≤40≥4m砂岩夹层≤40＜4m泥岩层≥60≥4m泥岩夹层≥60＜4m泥质砂岩层40＜vsh＜60≥4m泥质砂岩夹层40＜vsh＜60＜4m
[0071]
步骤3：根据岩性纵向分布选择厚度较大的砂岩层作为回注目标层，首先结合回注作业区(陆地或海上)通过岩性纵向分布判断所选砂岩层顶部距离地表(陆地回注)或海底(海上平台回注)的距离是否大于600m，若深度大于600m可进行下一步选择，否则此目标层不符合回注深度要求，需另选其它层。
[0072]
步骤4：当目标层深度满足要求后，通过岩性纵向分布判断目标层厚度，若目标地层厚度大于10m，满足回注层要求，否则此目标层不符合回注厚度要求，需另选其它层。
[0073]
步骤5：当目标层深度、厚度满足要求后，通过岩性纵向分布判断目标层的盖层厚度，若盖层厚度大于20m，满足回注盖层要求，否则盖层不满足要求，需另选其它层。
[0074]
步骤6：当目标层深度、厚度和盖层厚度都满足要求后，需将伽马测井、声波测井和密度测井带入公式(2
‑
9)，分别计算目标层孔隙度、渗透率、破裂压力和盖层破裂压力剖面。
[0075][0076]
σ
v
＝∫
0h
ρ
rock
*g*dh
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0077][0078][0079][0080]
[0081][0082][0083]
若取当计算结果s
wi
≤15％，取s
wi
＝15％
[0084]
其中：
[0085]
p
f
‑
破裂压力，mpa；
[0086]
σ
v
‑
上覆岩层压力，mpa；
[0087]
σ
h
‑
最小地应力，mpa；
[0088]
σ
h
‑
最大地应力，mpa；
[0089]
k
‑
渗透率，d；
[0090]
s
t
‑
岩石抗拉强度，实验获得，mpa；
[0091]
α
‑
biot系数，经验值φ＜α＜1；
[0092]
μ
‑
岩石静态泊松比，通过实验获得；
[0093]
φ
‑
孔隙度，％；
[0094]
ρ
rock
‑
岩石密度，g/cm3,由密度测井可得；
[0095]
p
p
‑
孔隙压力，mpa；
[0096]
δt
‑
测井声波时差，us/ft；
[0097]
δt
ma
‑
岩石骨架声波时差，砂岩182us/m,石灰岩156us/m，白云岩143us/m，此处取150us/m；
[0098]
δt
f
‑
泥浆声波时差，盐水608us/m，淡水620us/m；
[0099]
δt
sh
‑
泥质(泥岩)声波时差，283us/m；
[0100]
δt
n
‑
正常压实泥岩声波时差,将实际测井的声波取对数后，进行一元线性回归，求取声波的趋势线；
[0101]
h
‑
井深，m；
[0102]
s
wi
‑
束缚水饱和度，％；
[0103]
p
n
‑
静液柱压力，mpa；
[0104]
x
‑
eaton指数，由区域规律或实钻数据确定，取1.2；
[0105]
ε1、ε2‑
构造应力系数，根据实验数据反算；
[0106]
步骤:7：通过步骤6中计算出的破裂压力剖面判断盖层底部和目标层顶部破裂压力差值，若盖层底部破裂压力大于目标层顶部破裂压力，且差值大于6.89mpa(1000psi)，盖层破裂压力满足要求，否则需另选其他层。
[0107]
步骤8：当目标层及盖层满足步骤1到步骤7的要求后，分别将步骤6计算得到的目标层孔隙度、渗透率和破裂压力剖面按照公式(10
‑
11)进行归一化处理，然后计算目标层归一化孔隙度平均值、归一化渗透率平均值和归一化破裂压力平均值并代入公式(12)计算回注潜力指数h。
[0108]
对于孔隙度和渗透率越大越优型指标归一化处理公式：
[0109]
[0110]
对于破裂压力越小越优型指标归一化处理公式：
[0111][0112]
其中x
max
为样本数据的最大值，对应待求参数剖面中最大数据值，x
min
为样本数据的最小值，对应待求参数剖面中最小数据值，x
i
为待计算样本数据，对应剖面中每个数据点，y为归一值。
[0113][0114]
其中a＝0.1365；b＝0.2385；c＝0.6250
[0115]
为归一化孔隙度平均值；
[0116]
为归一化渗透率平均值；
[0117]
为归一化破裂压力平均值；
[0118]
其中，a、b、c是根据下述的ahp层次分析法得到，在公式(12)中a、b、c分别代表孔隙度、渗透率、破裂压力对回注层回注潜力指数的影响权重，无实际物理意义。
[0119]
需要说明的是，本发明根据ahp层次分析法得到孔隙度、渗透率和破裂压力对于回注潜力影响最大。定义回注潜力指数h等于回注目标层归一化孔隙度平均值、归一化渗透率平均值、归一化破裂压力平均值与其影响权重a、b、c的乘积之和，然后通过h大小判断该层是否适合回注。
[0120]
需要说明的是在本发明中“剖面”一词的定义是某参数与深度数据点的一一对应值组成的数据集合，表现为某参数值随深度的变化。
[0121]
需要说明的归一化的目的是将某参数实际数值转换为0
‑
1的标准取值范围内的数值，归一化平均值是指计算该参数归一化后的数据平均值。
[0122]
步骤9：根据步骤8计算所得h，参照表2确定目标层是否适合进行回注作业。
[0123]
回注建议适合回注谨慎回注不适合回注hh≥0.50.3<h<0.5h≤0.3
[0124]
需要说明的是h指数的取值范围并不是随意指定，而是根据多组室内实验效果确定。
[0125]
其中，公式(12)中系数a、b、c的确定方法可采用以下的ahp层次分析法：
[0126]
1.构建压裂效果判断矩阵。
[0127][0128]
通过比较渗透率k(c1)、孔隙度φ(c2)、破裂压力p
f
(c3)对回注效果的敏感性大小，确定影响回注效果的判断矩阵为：
[0129] c1c2c3c1121/3c21/211/4c3341
[0130]
2.各参数的权重计算与一致性检验步骤。
[0131]
(1)将矩阵各行向量进行几何平均：
[0132][0133]
(2)规范化处理
[0134][0135]
(3)各影响参数对上一层次的权重：
[0136]
w＝[w1,w2,
…
,w
n
]
t
[0137]
(4)判断矩阵a的最大特征根：
[0138][0139]
(5)判断a矩阵的一致性：
[0140][0141]
当cr<0.10时，即可认为判断矩阵具有符合的一致性，求得权重符合要求。若cr>0.10，则需要调整判断矩阵，直至cr<0.10。其中，一致性指标ri与n之间的关系见下表：
[0142]
n123456789ri0.000.000.580.901.121.241.321.411.45
[0143]
本发明实施方式中：
[0144]
λ
max
＝3.0183
[0145]
一致性符合要求，所以各参数对回注效果的影响权重如下：
[0146]
参数渗透率(k)孔隙度(φ)破裂压力(p
f
)权重0.23850.13650.6250
[0147]
3.根据多组实验结果对最终回注潜力指数h参数进行范围取值，最终结果如下：
[0148]
h≥0.5可作为回注层
[0149]
0.3<h<0.5谨慎回注
[0150]
h≤0.3不可回注。
[0151]
需要说明的是，在本文中，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的系统或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制；诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0152]
此外，在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
[0153]
上述各实施例仅用于说明本发明，其中实施例的各零部件、装置都是可以有所变化的，各实施方式都可根据需要进行组合或删减，附图中并非所有部件都是必要设置，本文中所定义的一般原理可以在不脱离本技术的精神或范围的情况下，在其它实施例中实现。因此，本技术将不会被限制于本文所述的这些实施例，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。