CN113404451A - 一种基于测井资料进行岩屑回注选层的方法 - Google Patents
一种基于测井资料进行岩屑回注选层的方法 Download PDFInfo
- Publication number
- CN113404451A CN113404451A CN202110154307.7A CN202110154307A CN113404451A CN 113404451 A CN113404451 A CN 113404451A CN 202110154307 A CN202110154307 A CN 202110154307A CN 113404451 A CN113404451 A CN 113404451A
- Authority
- CN
- China
- Prior art keywords
- reinjection
- layer
- target layer
- target
- thickness
- 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
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000011435 rock Substances 0.000 title claims abstract description 42
- 230000035699 permeability Effects 0.000 claims abstract description 29
- 239000010410 layer Substances 0.000 claims description 163
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 239000011229 interlayer Substances 0.000 claims description 9
- 238000005553 drilling Methods 0.000 claims description 5
- 238000005520 cutting process Methods 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- 239000011148 porous material Substances 0.000 claims description 3
- 235000019738 Limestone Nutrition 0.000 claims description 2
- 239000010459 dolomite Substances 0.000 claims description 2
- 229910000514 dolomite Inorganic materials 0.000 claims description 2
- 239000013505 freshwater Substances 0.000 claims description 2
- 239000006028 limestone Substances 0.000 claims description 2
- 238000012417 linear regression Methods 0.000 claims description 2
- 239000002002 slurry Substances 0.000 claims description 2
- 230000003068 static effect Effects 0.000 claims description 2
- 238000002347 injection Methods 0.000 claims 2
- 239000007924 injection Substances 0.000 claims 2
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims 1
- 239000011780 sodium chloride Substances 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 13
- 238000010606 normalization Methods 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 239000002699 waste material Substances 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000010170 biological method Methods 0.000 description 1
- 239000012267 brine Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 1
- 239000003440 toxic substance Substances 0.000 description 1
- 239000010913 used oil Substances 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B21/00—Methods or apparatus for flushing boreholes, e.g. by use of exhaust air from motor
- E21B21/01—Arrangements for handling drilling fluids or cuttings outside the borehole, e.g. mud boxes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mining & Mineral Resources (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Geology (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Databases & Information Systems (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种基于测井资料进行岩屑回注选层的方法，包括：S1、根据伽马测井资料计算待回注井的泥质含量随深度的变化；S2、确定地层岩性纵向分布；S3、选择砂岩所在层为回注目标层，判断该目标层的深度，若符合要求，执行下一步，S4、判断回注目标层的厚度是否符合要求，S5、判断回注目标层的盖层的厚度是否符合要求，S6、根据测井资料，分别求得目标层孔隙度、渗透率、破裂压力和盖层破裂压力剖面；S7、判断回注目标层和盖层的破裂压力差是否符合要求，S8、根据目标层孔隙度、渗透率和破裂压力值确定回注目标层的回注潜力指数H；S9：根据回注潜力指数确定回注目标层是否合适。本发明提供有科学依据的选层流程，可有效解决经验选层的不足问题。
Description
技术领域
本发明涉及油气开发技术领域，具体涉及一种基于测井资料进行岩屑回注选层的方法。
背景技术
近年来，油气开发越来越多地使用油基钻井液，这会产生大量的具有高污染性的油基岩屑。目前随着各个国家对环保的要求越来越严格，钻井废物必须达到环境要求才能排放，因此，如何高效低成本处理油基岩屑成为制约钻井进度的关键因素。
目前常用的处理方法有填埋法、机械分离法、生物法、热处理法及回注法。其中填埋法虽成本较低，但是由于填埋深度较浅，当遭遇雨季冲刷或地下水渗透可能导致油基废物泄漏，污染环境；机械分离法和热处理方法设备运行成本高；生物法处理周期长，且油基岩屑中的有毒物质会抑制生化降解过程；回注法是将油基岩屑研磨处理，通过高压泵将含有岩屑的浆液从环空或专门的回注井注入到目标地层，不仅可以从根本上解决油基岩屑污染环境的问题，同时处理成本低，时效性好。尤其是海上油气田，回注处理相比运回陆地处理，成本大幅降低，降低了海上平台空间的占用。但目前回注法在国外应用较多，国内尚未进行大规模应用，原因是回注目的层的筛选大都基于经验确定，缺少科学方法或指导方案，存在压力过高盖层破裂、孔隙不达标以及盖层稳定性差等风险，不仅影响回注岩屑的效率和回注量，更严重的是一旦发生泄漏，会对地层水或海洋环境造成不可逆的严重污染。
发明内容
有鉴于此，本发明旨在提供一种基于测井资料进行岩屑回注选层的方法，以提供完整的、有科学依据的选层流程，以便有效解决经验选层的弊端。
本发明首先提出一种基于测井资料进行岩屑回注选层的方法，所述方法包括以下步骤：
S1、根据伽马测井资料计算待回注井对应地层泥质含量剖面；
S2、根据所述泥质含量和岩性判断标准，确定地层岩性纵向分布；
S3、选择所述地层岩性纵向分布中砂岩层为回注目标层，判断该目标层的深度是否符合要求；
S4、判断所述回注目标层的厚度是否符合要求；
S5、判断所述回注目标层盖层的厚度是否符合要求；
S6、根据测井资料，分别计算回注目标层孔隙度、渗透率、破裂压力和盖层破裂压力剖面；
S7、判断所述回注目标层和所述盖层的破裂压力差是否符合要求；
S8、根据所述目标层孔隙度、渗透率和破裂压力剖面分别计算目标层归一化孔隙度平均值，归一化渗透率平均值，归一化破裂压力平均值并确定所述回注目标层的回注潜力指数H；
S9：根据所述回注潜力指数H确定所述回注目标层是否合适。
根据本发明的一种实施方式，在步骤S8中，所述回注目标层的回注潜力指数H
按照下式确定：
其中A＝0.1365；B＝0.2385；C＝0.6250
H指数是通过AHP层次分析法计算得到的一个无量纲数值，其代表的是目标地层回注潜能指数。
根据本发明的一种实施方式，在步骤S1中，通过获取伽马测井曲线，根据以下公式计算地层泥质含量剖面，
其中：
Vsh-泥质含量，％
GR-测井伽马值
GRsand-纯砂岩伽马值
GRmud-纯泥岩伽马值。
根据本发明的一种实施方式，在步骤S2中所述地层岩性纵向分布按如下标准判断，
若泥质含量(Vsh)≤40，且该层的厚度大于≥4m，则该层为砂岩层；
若泥质含量(Vsh)≤40，且该层的厚度大于<4m，则该层为砂岩夹层；
若泥质含量(Vsh)≥60，且该层的厚度大于≥4m，则该层为泥岩；
若泥质含量(Vsh)≥60，且该层的厚度大于<4m，则该层为泥岩夹层；
若40<泥质含量(Vsh)<60,且该层的厚度大于≥4m，则该层为泥质砂岩层。
若40<泥质含量(Vsh)<60,且该层的厚度大于<4m，则该层为泥质砂岩夹层。
根据本发明的一种实施方式，在步骤S3中，结合回注作业区(陆地或海上)通过岩性纵向分布判断所选砂岩层顶部距离地表(陆地回注)或海底(海上平台回注)的距离是否大于600m，若大于600m则进行下一步选择，否则此目标层不符合回注深度要求。
根据本发明的一种实施方式，在步骤S4中，通过所述地层岩性纵向分布判断所述回注目标层的厚度，若该回注目标层的厚度大于10m，则满足回注要求，否则此目标层不符合回注厚度要求。
根据本发明的一种实施方式，在步骤S5中，通过所述地层岩性纵向分布判断所述回注目标层盖层的厚度，若盖层厚度大于20m，则满足回注盖层的要求，否则该盖层不满足要求。
根据本发明的一种实施方式，在步骤S6中，通过伽马测井值、声波测井值和密度测井值分别计算所述回注目标层的孔隙度、渗透率、破裂压力以及所述盖层的破裂压力剖面。
根据本发明的一种实施方式，在步骤S7中，若所述盖层底部破裂压力大于所述回注目标层顶部破裂压力，且两者的差值大于6.89MPa(即1000psi)，则所述盖层破裂压力满足要求，否则不符合回注要求。
根据本发明的一种实施方式，在步骤S8中，分别将步骤S6获得的目标地层孔隙度、渗透率和破裂压力剖面做归一化处理，然后计算所述回注目标层的归一化孔隙度平均值、归一化渗透率平均值和归一化破裂压力平均值，再计算所述目标层回注潜力指数H；
在步骤S9中，当H≥0.5，则所述回注目标层适合回注；当0.3<H<0.5
，则所述回注目标层应谨慎回注；当H≤0.3，则所述回注目标层不适合回注；
在步骤S3、S4、S5、S7中任一步不符合要求，则所述回注目标层不适合回注，需要重新返回第一步，选取另外目标层进行判断。
本发明基于声波测井、伽马测井和密度测井资料，通过对回注目的层及盖层一系列参数的计算和判断，从而精确筛选回注目的层。该方法可提供完整的、有科学依据的岩屑回注选层流程，可有效解决经验选层的弊端，且计算成本低，筛选方便。
附图说明
图1为本发明一实施例回注选层流程示意图。
具体实施方式
以下将结合附图对本发明的较佳实施例进行详细说明，以便更清楚理解本发明的目的、特点和优点。应理解的是，附图所示的实施例并不是对本发明范围的限制，而只是为了说明本发明技术方案的实质精神。
本发明提供一种基于常规测井资料的岩屑回注选层的方法，这种方法是基于声波测井、伽马测井和密度测井资料，通过对目的层及盖层一系列参数的计算和判断，从而精确筛选回注目的层。该方法可提供完整的、有科学依据的选层流程，可有效解决经验选层的弊端，且计算成本低，筛选方便。
本实施方式一种基于测井资料进行岩屑回注选层的方法，通过声波测井、伽马测井和密度测井资料计算地层参数，根据筛选标准进行回注选层，如图1所示，主要包括如下步骤：
步骤1：根据伽马测井资料计算地层泥质含量剖面；
步骤2：根据岩性分类标准，确定地层岩性纵向分布；
步骤3：选择砂岩层为回注目标层，判断目标层深度；
步骤4：判断目标层厚度；
步骤5：判断盖层厚度；
步骤6：将测井资料代入计算公式，分别计算回注目标层孔隙度、渗透率、破裂压力和盖层破裂压力剖面；
步骤7：判断目标层和盖层破裂压力差；
步骤8：计算目标层归一化孔隙度平均值，归一化渗透率平均值，归一化破裂压力平均值，并代入公式计算回注潜力指数H；
步骤9：根据选层判别标准选取回注层。
下面通过实施方式将更具体地对每步骤进行详细介绍。
根据本发明的一种实施方式，在具体执行上述步骤时，采取了以下实施方式。
步骤1：获取伽马测井曲线，根据公式(1)计算地层泥质含量剖面。
其中：
Vsh-泥质含量，％
GR-测井伽马值
GRsand-纯砂岩伽马值
GRmud-纯泥岩伽马值
步骤2：根据下表标准将地层进行分层，确定地层岩性纵向分布。
| 地层类型 | 泥质含量(Vsh) | 厚度 |
| 砂岩层 | ≤40 | ≥4m |
| 砂岩夹层 | ≤40 | ＜4m |
| 泥岩层 | ≥60 | ≥4m |
| 泥岩夹层 | ≥60 | ＜4m |
| 泥质砂岩层 | 40＜Vsh＜60 | ≥4m |
| 泥质砂岩夹层 | 40＜Vsh＜60 | ＜4m |
步骤3：根据岩性纵向分布选择厚度较大的砂岩层作为回注目标层，首先结合回注作业区(陆地或海上)通过岩性纵向分布判断所选砂岩层顶部距离地表(陆地回注)或海底(海上平台回注)的距离是否大于600m，若深度大于600m可进行下一步选择，否则此目标层不符合回注深度要求，需另选其它层。
步骤4：当目标层深度满足要求后，通过岩性纵向分布判断目标层厚度，若目标地层厚度大于10m，满足回注层要求，否则此目标层不符合回注厚度要求，需另选其它层。
步骤5：当目标层深度、厚度满足要求后，通过岩性纵向分布判断目标层的盖层厚度，若盖层厚度大于20m，满足回注盖层要求，否则盖层不满足要求，需另选其它层。
步骤6：当目标层深度、厚度和盖层厚度都满足要求后，需将伽马测井、声波测井和密度测井带入公式(2-9)，分别计算目标层孔隙度、渗透率、破裂压力和盖层破裂压力剖面。
σv＝∫0 hρrock*g*dh (3)
若
取
当计算结果Swi≤15％，取Swi＝15％
其中：
pf-破裂压力，MPa；
σv-上覆岩层压力，MPa；
σh-最小地应力，MPa；
σH-最大地应力，MPa；
K-渗透率，D；
St-岩石抗拉强度，实验获得，MPa；
α-Biot系数，经验值φ＜α＜1；
μ-岩石静态泊松比，通过实验获得；
φ-孔隙度，％；
ρrock-岩石密度，g/cm3,由密度测井可得；
Pp-孔隙压力，MPa；
Δt-测井声波时差，us/ft；
Δtma-岩石骨架声波时差，砂岩182us/m,石灰岩156us/m，白云岩143us/m，此处取150us/m；
Δtf-泥浆声波时差，盐水608us/m，淡水620us/m；
Δtsh-泥质(泥岩)声波时差，283us/m；
Δtn-正常压实泥岩声波时差,将实际测井的声波取对数后，进行一元线性回归，求取声波的趋势线；
h-井深，m；
Swi-束缚水饱和度，％；
Pn-静液柱压力，MPa；
x-Eaton指数，由区域规律或实钻数据确定，取1.2；
ε1、ε2-构造应力系数，根据实验数据反算；
步骤:7：通过步骤6中计算出的破裂压力剖面判断盖层底部和目标层顶部破裂压力差值，若盖层底部破裂压力大于目标层顶部破裂压力，且差值大于6.89MPa(1000psi)，盖层破裂压力满足要求，否则需另选其他层。
步骤8：当目标层及盖层满足步骤1到步骤7的要求后，分别将步骤6计算得到的目标层孔隙度、渗透率和破裂压力剖面按照公式(10-11)进行归一化处理，然后计算目标层归一化孔隙度平均值、归一化渗透率平均值和归一化破裂压力平均值并代入公式(12)计算回注潜力指数H。
对于孔隙度和渗透率越大越优型指标归一化处理公式：
对于破裂压力越小越优型指标归一化处理公式：
其中Xmax为样本数据的最大值，对应待求参数剖面中最大数据值，Xmin为样本数据的最小值，对应待求参数剖面中最小数据值，Xi为待计算样本数据，对应剖面中每个数据点，Y为归一值。
其中A＝0.1365；B＝0.2385；C＝0.6250
其中，A、B、C是根据下述的AHP层次分析法得到，在公式(12)中A、B、C分别代表孔隙度、渗透率、破裂压力对回注层回注潜力指数的影响权重，无实际物理意义。
需要说明的是，本发明根据AHP层次分析法得到孔隙度、渗透率和破裂压力对于回注潜力影响最大。定义回注潜力指数H等于回注目标层归一化孔隙度平均值、归一化渗透率平均值、归一化破裂压力平均值与其影响权重A、B、C的乘积之和，然后通过H大小判断该层是否适合回注。
需要说明的是在本发明中“剖面”一词的定义是某参数与深度数据点的一一对应值组成的数据集合，表现为某参数值随深度的变化。
需要说明的归一化的目的是将某参数实际数值转换为0-1的标准取值范围内的数值，归一化平均值是指计算该参数归一化后的数据平均值。
步骤9：根据步骤8计算所得H，参照表2确定目标层是否适合进行回注作业。
| 回注建议 | 适合回注 | 谨慎回注 | 不适合回注 |
| H | H≥0.5 | 0.3<H<0.5 | H≤0.3 |
需要说明的是H指数的取值范围并不是随意指定，而是根据多组室内实验效果确定。
其中，公式(12)中系数A、B、C的确定方法可采用以下的AHP层次分析法：
1.构建压裂效果判断矩阵。
通过比较渗透率K(C1)、孔隙度φ(C2)、破裂压力Pf(C3)对回注效果的敏感性大小，确定影响回注效果的判断矩阵为：
| C1 | C2 | C3 | |
| C1 | 1 | 2 | 1/3 |
| C2 | 1/2 | 1 | 1/4 |
| C3 | 3 | 4 | 1 |
2.各参数的权重计算与一致性检验步骤。
(1)将矩阵各行向量进行几何平均：
(2)规范化处理
(3)各影响参数对上一层次的权重：
W＝[W1,W2,…,Wn]T
(4)判断矩阵A的最大特征根：
(5)判断A矩阵的一致性：
当CR<0.10时，即可认为判断矩阵具有符合的一致性，求得权重符合要求。若CR>0.10，则需要调整判断矩阵，直至CR<0.10。其中，一致性指标RI与n之间的关系见下表：
| n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
| RI | 0.00 | 0.00 | 0.58 | 0.90 | 1.12 | 1.24 | 1.32 | 1.41 | 1.45 |
本发明实施方式中：
λmax＝3.0183
一致性符合要求，所以各参数对回注效果的影响权重如下：
| 参数 | 渗透率(K) | 孔隙度(φ) | 破裂压力(P<sub>f</sub>) |
| 权重 | 0.2385 | 0.1365 | 0.6250 |
3.根据多组实验结果对最终回注潜力指数H参数进行范围取值，最终结果如下：
H≥0.5可作为回注层
0.3<H<0.5谨慎回注
H≤0.3不可回注。
需要说明的是，在本文中，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的系统或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制；诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
此外，在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
上述各实施例仅用于说明本发明，其中实施例的各零部件、装置都是可以有所变化的，各实施方式都可根据需要进行组合或删减，附图中并非所有部件都是必要设置，本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下，在其它实施例中实现。因此，本申请将不会被限制于本文所述的这些实施例，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。
Claims (10)
1.一种基于测井资料进行岩屑回注选层的方法，其特征在于，所述方法包括以下步骤：
S1、根据伽马测井资料计算待回注井对应地层的泥质含量剖面；
S2、根据所述泥质含量剖面，确定地层岩性纵向分布；
S3、选择所述地层岩性纵向分布中砂岩所在层为回注目标层，判断该目标层的深度是否符合要求；
S4、判断所述回注目标层的厚度是否符合要求；
S5、判断所述回注目标层盖层的厚度是否符合要求；
S6、根据测井资料，分别计算回注目标地层孔隙度、渗透率、破裂压力和盖层破裂压力剖面；
S7、判断所述回注目标层和所述盖层的破裂压力差是否符合要求；
S8、若前述判断步骤均符合要求，根据所述目标层孔隙度、渗透率和破裂压力剖面分别计算归一化孔隙度平均值，归一化渗透率平均值，归一化破裂压力平均值，并确定所述回注目标层的回注潜力指数H；
S9：根据所述回注潜力指数H确定所述回注目标层是否合适。
2.根据权利要求1所述的基于测井资料进行岩屑回注选层的方法，其特征在于，在步骤S8中，所述回注目标层的回注潜力指数H按照下式确定：
其中A＝0.1365；B＝0.2385；C＝0.6250
3.根据权利要求1或2所述的基于测井资料进行岩屑回注选层的方法，其特征在于，在步骤S1中，通过获取伽马测井曲线，根据以下公式计算地层泥质含量剖面，
其中：
Vsh-泥质含量，％
GR-测井伽马值
GRsand-纯砂岩伽马值
GRmud-纯泥岩伽马值。
4.根据权利要求1至3任一项所述的基于测井资料进行岩屑回注选层的方法，其特征在于，所述地层岩性纵向分布判断标准为：
若泥质含量Vsh≤40，且该层的厚度大于≥4m，则该层为砂岩层；
若泥质含量Vsh≤40，且该层的厚度大于＜4m，则该层为砂岩夹层；
若泥质含量Vsh≥60，且该层的厚度大于≥4m，则该层为泥岩；
若泥质含量Vsh≥60，且该层的厚度大于＜4m，则该层为泥岩夹层；
若40＜泥质含量Vsh＜60，且该层的厚度大于≥4m，则该层为泥质砂岩层；
若40＜泥质含量Vsh＜60，且该层的厚度大于＜4m，则该层为泥质砂岩夹层。
5.根据权利要求1至3任一项所述的基于测井资料进行岩屑回注选层的方法，其特征在于，在步骤S3中，结合回注作业区，通过岩性纵向分布判断所选砂岩层顶部距离地表或海底的距离是否大于600m，若该深度大于600m则进行下一步选择，否则此目标层不符合回注深度要求。
6.根据权利要求1至3任一项所述的基于测井资料进行岩屑回注选层的方法，其特征在于，在步骤S4中，通过所述地层岩性纵向分布判断所述回注目标层的厚度，若该回注目标层的厚度大于10m，则满足回注要求，否则此目标层不符合回注厚度要求。
7.根据权利要求1至3任一项所述的基于测井资料进行岩屑回注选层的方法，其特征在于，在步骤S5中，通过所述地层岩性纵向分布判断所述回注目标层的盖层的厚度，若该盖层的厚度大于20m，则满足回注盖层的要求，否则该盖层不满足要求。
8.根据权利要求1至3任一项所述的基于测井资料进行岩屑回注选层的方法，其特征在于，在步骤S6中，通过伽马测井值、声波测井值和密度测井值分别计算所述回注目标层的孔隙度、渗透率、破裂压力以及所述盖层的破裂压力剖面：
若
取
当计算结果Swi≤15％，取Swi＝15％
其中：
pf-破裂压力，MPa；
σv-上覆岩层压力，MPa；
σh-最小地应力，MPa；
σH-最大地应力，MPa；
K-渗透率，D；
St-岩石抗拉强度，实验获得，MPa；
α-Biot系数，经验值φ＜α＜1；
μ-岩石静态泊松比，通过实验获得；
φ-孔隙度，％；
ρrock-岩石密度，g/cm3，由密度测井可得；
Pp-孔隙压力，MPa；
Δt-测井声波时差，us/ft；
Δtma-岩石骨架声波时差，砂岩182us/m，石灰岩156us/m，白云岩143us/m；
Δtf-泥浆声波时差，盐水608us/m淡水620us/m；
Δtsh-泥质(泥岩)声波时差，283us/m；
Δtn-正常压实泥岩声波时差，将实际测井的声波取对数后，进行一元线性回归，求取声波的趋势线；
h-井深，m；
Swi-束缚水饱和度，％；
Pn-静液柱压力，MPa；
x-Eaton指数，由区域规律或实钻数据确定，取1.2；
ε1、ε2-构造应力系数，根据实验数据反算。
9.根据权利要求8所述的基于测井资料进行岩屑回注选层的方法，其特征在于，在步骤S7中，若所述盖层的底部破裂压力大于所述回注目标层的顶部破裂压力值，且两者的差值大于6.89MPa，则所述盖层的应力满足要求，否则不符合回注要求。
10.根据权利要求9所述的基于测井资料进行岩屑回注选层的方法，其特征在于，在步骤S8中，将步骤S6获得的目标层孔隙度、渗透率和破裂压力值分别做归一化处理，然后计算所述回注目标层的归一化孔隙度平均值、归一化渗透率平均值和归一化破裂压力平均值，再计算所述回注潜力指数H；
在步骤S9中，当H≥0.5，则所述回注目标层适合回注；当0.3＜H＜0.5，则所述回注目标层应谨慎回注；当H≤0.3，则所述回注目标层不适合回注；
在步骤S3、S4、S5、S7中任一步不符合要求，则所述回注目标层不适合回注。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110154307.7A CN113404451B (zh) | 2021-02-04 | 2021-02-04 | 一种基于测井资料进行岩屑回注选层的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110154307.7A CN113404451B (zh) | 2021-02-04 | 2021-02-04 | 一种基于测井资料进行岩屑回注选层的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113404451A true CN113404451A (zh) | 2021-09-17 |
| CN113404451B CN113404451B (zh) | 2022-08-02 |
Family
ID=77675893
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110154307.7A Active CN113404451B (zh) | 2021-02-04 | 2021-02-04 | 一种基于测井资料进行岩屑回注选层的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113404451B (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0010017D0 (en) * | 1999-04-26 | 2000-06-14 | Vetco Gray Inc Abb | Drill cuttings injection system |
| AU2001231621A1 (en) * | 2000-11-28 | 2002-08-15 | Apollo Services Uk Ltd | Apparatus and method for transferring dry oil and gas well drill cuttings |
| WO2009036033A1 (en) * | 2007-09-13 | 2009-03-19 | M-I Llc | Method and system for injecting a slurry downhole |
| CN103132935A (zh) * | 2013-02-23 | 2013-06-05 | 中国石油大学(华东) | 油田钻井钻屑分级真空干燥处理工艺及处理装置 |
| CN103899280A (zh) * | 2014-04-16 | 2014-07-02 | 杰瑞能源服务有限公司 | 一种钻井废弃物的回注系统及回注方法 |
| CN105158802A (zh) * | 2015-08-21 | 2015-12-16 | 西南石油大学 | 坳陷湖盆重力流沉积物测井定量识别方法 |
| CN105888613A (zh) * | 2014-12-23 | 2016-08-24 | 桂淑华 | 钻屑深井注入工艺 |
| CN206816185U (zh) * | 2017-05-10 | 2017-12-29 | 李斌 | 一种钻井套管 |
| CN109162656A (zh) * | 2018-09-07 | 2019-01-08 | 洲际海峡能源科技有限公司 | 一种利用枯竭性油气藏处置固体废弃物的技术操作方法 |
| CN110240373A (zh) * | 2019-06-17 | 2019-09-17 | 菁华绿尔生态环境科技(天津)有限公司 | 一种处理含聚油泥的清洗剂及处理方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6585115B1 (en) * | 2000-11-28 | 2003-07-01 | Baker Hughes Incorporated | Apparatus and method for transferring dry oil and gas well drill cuttings |
-
2021
- 2021-02-04 CN CN202110154307.7A patent/CN113404451B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0010017D0 (en) * | 1999-04-26 | 2000-06-14 | Vetco Gray Inc Abb | Drill cuttings injection system |
| AU2001231621A1 (en) * | 2000-11-28 | 2002-08-15 | Apollo Services Uk Ltd | Apparatus and method for transferring dry oil and gas well drill cuttings |
| WO2009036033A1 (en) * | 2007-09-13 | 2009-03-19 | M-I Llc | Method and system for injecting a slurry downhole |
| CN103132935A (zh) * | 2013-02-23 | 2013-06-05 | 中国石油大学(华东) | 油田钻井钻屑分级真空干燥处理工艺及处理装置 |
| CN103899280A (zh) * | 2014-04-16 | 2014-07-02 | 杰瑞能源服务有限公司 | 一种钻井废弃物的回注系统及回注方法 |
| CN105888613A (zh) * | 2014-12-23 | 2016-08-24 | 桂淑华 | 钻屑深井注入工艺 |
| CN105158802A (zh) * | 2015-08-21 | 2015-12-16 | 西南石油大学 | 坳陷湖盆重力流沉积物测井定量识别方法 |
| CN206816185U (zh) * | 2017-05-10 | 2017-12-29 | 李斌 | 一种钻井套管 |
| CN109162656A (zh) * | 2018-09-07 | 2019-01-08 | 洲际海峡能源科技有限公司 | 一种利用枯竭性油气藏处置固体废弃物的技术操作方法 |
| CN110240373A (zh) * | 2019-06-17 | 2019-09-17 | 菁华绿尔生态环境科技(天津)有限公司 | 一种处理含聚油泥的清洗剂及处理方法 |
Non-Patent Citations (3)
| Title |
|---|
| 周锦程: "基于层次分析法的岩屑回注层筛选决策研究", 《中国石油和化工标准与质量》 * |
| 李斌: "岩屑回注技术及其在我国海洋油气田应用前景展望", 《海洋工程装备与技术》 * |
| 黄智辉: "《地球物理测井资料在分析沉积环境中的应用》", 30 December 1986 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113404451B (zh) | 2022-08-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110412661B (zh) | 细粒岩油气藏甜点段优势段簇的评价方法及装置 | |
| CN109653725B (zh) | 基于沉积微相和岩石相的混积储层水淹程度测井解释方法 | |
| Dreyer et al. | Minipermeameter-based study of permeability trends in channel sand bodies | |
| CN112132454B (zh) | 一种煤层顶板或底板含水层富水性综合评价方法 | |
| CN1749780A (zh) | 油气藏三元分类方法 | |
| CN109187302A (zh) | 一种基于盖层非均质性的盖层有效厚度计算方法 | |
| CN107589469A (zh) | 油水界面的确定方法和装置 | |
| CN113109228A (zh) | 基于煤矿地质钻孔的松散承压含水层渗透系数确定方法 | |
| Glover et al. | Clustering of facies in tight carbonates using machine learning | |
| Melas et al. | Petrophysical characteristics of the Jurassic Smackover formation, jay field, Conecuh embayment, Alabama and Florida | |
| CN113404451B (zh) | 一种基于测井资料进行岩屑回注选层的方法 | |
| CN110688781A (zh) | 一种低渗非均质气藏储层测井解释方法 | |
| CN112228050A (zh) | 一种定量评价致密油储层宏观非均质性的方法及其应用 | |
| CN111101914A (zh) | 水平井压裂段簇优选方法及设备 | |
| WO2022242200A1 (zh) | 一种海相碳酸盐岩生屑灰岩油藏高渗条带识别方法及装置 | |
| CN115905917A (zh) | 综合静动态数据构建海域低渗气层分类曲线的方法 | |
| CN114109374B (zh) | 页岩气储层靶窗位置的确定方法及装置 | |
| CN112663417B (zh) | 一种对岩溶路基段处置范围的筛选方法 | |
| EP3052970A1 (en) | Pore size classification in subterranean formations based on nuclear magnetic resonance (nmr) relaxation distributions | |
| Rafiee et al. | A new cementation factor correlation in carbonate parts of oil fields in south-west Iran | |
| Chen | The influence of tensile strength difference with variable gravel sizes on the hydraulic fracture propagation in the conglomerate reservoir | |
| CN114592848A (zh) | 孔隙度-电阻率-岩性匹配关系法识别低阻油气层的方法 | |
| Bigeldiyev et al. | Dynamic modeling of the gas discharge of a mine in the Karaganda coal basin under high uncertainty using a multiple realization approach | |
| CN116894393B (zh) | 一种多参量信息融合的顶板含水层富水性判别方法 | |
| CN113818859B (zh) | 一种低渗透油藏的极限井距确定方法、判定方法与装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |