CN109142676A - 一种利用测井渗透率进行地层有效裂缝的识别方法 - Google Patents
一种利用测井渗透率进行地层有效裂缝的识别方法 Download PDFInfo
- Publication number
- CN109142676A CN109142676A CN201810876663.8A CN201810876663A CN109142676A CN 109142676 A CN109142676 A CN 109142676A CN 201810876663 A CN201810876663 A CN 201810876663A CN 109142676 A CN109142676 A CN 109142676A
- Authority
- CN
- China
- Prior art keywords
- fluid
- volume
- drilling
- formula
- stratum
- 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
- 230000035699 permeability Effects 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000012530 fluid Substances 0.000 claims abstract description 130
- 238000005553 drilling Methods 0.000 claims abstract description 87
- 238000005325 percolation Methods 0.000 claims abstract description 22
- 239000000706 filtrate Substances 0.000 claims abstract description 17
- 239000007790 solid phase Substances 0.000 claims abstract description 16
- 239000002245 particle Substances 0.000 claims abstract description 14
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 10
- 239000007791 liquid phase Substances 0.000 claims abstract description 5
- 239000011435 rock Substances 0.000 claims description 39
- 239000011707 mineral Substances 0.000 claims description 22
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 19
- 239000011148 porous material Substances 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 13
- 230000009545 invasion Effects 0.000 claims description 11
- 238000002791 soaking Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 239000002734 clay mineral Substances 0.000 claims description 7
- 230000000704 physical effect Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 4
- 238000001595 flow curve Methods 0.000 claims description 3
- 238000002637 fluid replacement therapy Methods 0.000 claims description 3
- 239000008187 granular material Substances 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000012263 liquid product Substances 0.000 claims description 2
- 239000000047 product Substances 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 6
- 230000008030 elimination Effects 0.000 abstract description 3
- 238000003379 elimination reaction Methods 0.000 abstract description 3
- 235000010755 mineral Nutrition 0.000 description 13
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- TZCXTZWJZNENPQ-UHFFFAOYSA-L barium sulfate Chemical compound [Ba+2].[O-]S([O-])(=O)=O TZCXTZWJZNENPQ-UHFFFAOYSA-L 0.000 description 2
- 229910052601 baryte Inorganic materials 0.000 description 2
- 239000010428 baryte Substances 0.000 description 2
- 230000008595 infiltration Effects 0.000 description 2
- 238000001764 infiltration Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-N carbonic acid Chemical compound OC(O)=O BVKZGUZCCUSVTD-UHFFFAOYSA-N 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000009969 flowable effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 239000011799 hole material Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 239000012071 phase Substances 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 235000002639 sodium chloride Nutrition 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Dispersion Chemistry (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明一种利用测井渗透率进行地层有效裂缝的识别方法,首先根据钻井液中的固相颗粒与液相滤液侵入地层的深度,判断地层的渗透性,并由实际工况得到流体流量Q,结合横截面积A、流体渗滤路径L、压差ΔP和流体粘度μ四项参数,根据达西定律计算出储层的渗透率。该方法有效的回避了无法预判储层为裂缝型储层或孔隙‑裂缝型储层,以及裂缝不同发育程度或产状的经验公式的选取困难,切实可行的解决了无法预判地层中是否存在裂缝所造成的被动问题,有效的实现了产能评估的降本增效,具有极大的应用价值。
Description
技术领域
本发明涉及孔隙裂缝或裂缝性地层的渗透率定量计算,确切地说涉及一种利用测井渗透率进行地层有效裂缝的识别方法,属于地层物性评价范畴。
背景技术
岩层是否具有渗透性以及渗透能力的大小,是油气田能否形成产能和产能大小的重要因素,对油气田开发具有重要意义。目前,渗透率的预测手段通常为实验室岩心测定和测井资料计算求取。实验室测定的岩心渗透率,是基于达西定律,以渗滤流体对岩样进行实验的测量结果。这种测量手段尺度较微观,成本较高且受限于样品采集,需采集的岩心样品数足够多时,实验室的分析结果才有可能接近实际储层的渗透率。因此,在实际的储层评价中通常采用测井手段来估算纵向上连续的储层渗透率。
如公开号为102507412B,公开日为2014年7月2日的中国专利文献公开了一种碳酸盐岩油藏等效模型的裂缝-基质渗透率级差判别方法,包括:等效裂缝单一介质数值模拟的裂缝-基质渗透率级差判别;等效基质单一介质数值模拟的裂缝-基质渗透率级差判别;等效双孔单渗数值模拟的裂缝-基质渗透率级差判别。本发明能够快速实现储层类型多样、裂缝与基质搭配关系复杂、渗透率级差变化大、油气水关系复杂、压力分布变化大、模拟区块大、生产史复杂等各类碳酸盐岩油气藏的历史拟合和方案预测。
由测井手段预测渗透率,是一种岩石物理学的概念,通常基于建立区域内测井参数与储层渗透率(通常由实验室测定)间的经验关系进行计算。目前在孔隙-裂缝性储层系统中,通常将储层的渗透率分为孔隙部分的渗透率和裂缝部分的渗透率来考虑。使用上述方法准确预测渗透率的前提是预知该段的孔隙构成究竟是孔隙还是裂缝,而在实际的计算过程中,这是难以实现的,因此对油气田产能的评估和实际开发都造成了十分被动的影响。或者使用核磁共振特殊测井来计算渗透率,然而这依旧是基于孔隙度的算法,其局限性是当地层孔隙特征不明显的时候很难准确计算。
发明内容
本发明旨在针对上述现有技术所存在的缺陷和不足,提供一种利用测井渗透率进行地层有效裂缝的识别方法,本方法切实可行地解决了无法预判地层中是“孔”还是“缝”所造成的被动问题,极大的优化了油气田产能评估流程,提高了油气田开发效率,对推进储层渗透率预测技术的进程具有革新性的意义。
本发明是通过采用下述技术方案实现的:
一种利用测井渗透率进行地层有效裂缝的识别方法,其特征在于:首先根据钻井液中的固相颗粒与液相滤液侵入地层的深度,判断地层的渗透性,并由实际工况得到流体流量Q,结合横截面积A、流体渗滤路径L、压差ΔP和流体粘度μ四项参数,根据达西定律计算出储层的渗透率:
式中,
K ---- 渗透率,μm2;
Q ---- 流体流量, cm3/s;
A ---- 垂直于流体流动方向岩石横截面积cm2;
L ---- 流体渗滤路径,cm;
---- 压差,Pa;
---- 流体粘度, mPa·s;
根据计算出的渗透率,识别地层的有效裂缝。
所述流体流量Q的计算方法如下:
首先利用体积模型法计算侵入钻井液体积,采用流体替换方法计算钻井液侵入体积,选用PE曲线来进行流体替换差值计算:
(1)
式中,
PE原 ---- 原始状态下理论岩心密度值,b/e;
PE骨架 ---- 各骨架矿物的岩性密度值,b/e;
V总 ---- 孔隙和矿物总百分含量,数值为100;
V骨架 ---- 各骨架矿物对应的体积百分比,%,测井计算求得;
PE黏土 ---- 各黏土矿物对应的体积百分比,b/e;
V黏土 ---- 各黏土矿物对应的体积百分比,%,测井计算求得;
PE储层流体 ---- 原状地层连通孔隙中流体的岩性密度值,b/e;
PE束缚流体 ---- 非连通孔隙中流体的岩性密度值,b/e;
连通 ---- 连通孔隙度,测井计算求得,%;
非连通 ---- 非连通孔隙度,测井计算求得,%;
PE钻井液 ---- 钻井液的PE值,b/e,根据钻井液中所加重矿物含量计算得到,或直接测量;
公式(1)表示当地层未被钻开时,原始状态下的岩性密度理论值;
(2)
式中,
PE测 ---- 实际测量岩性密度值,b/e;
---- 侵入钻井液所占孔隙度,%;
公式(2)表示当地层钻开后,部分连通孔隙和裂缝中流体被钻井液替换后的实际测量的岩性密度值;
(3)
式中,
侵入---- 地层被钻开前后岩性密度变化值Δ侵入,b/e;
式(3)表示由测井仪器实测得到的岩性密度值PE测减去原始地层状态下的岩性密度值PE原,将得到地层被钻开后钻井液部分或完全侵入时的连通孔隙中和裂缝中岩性密度的变化值Δ侵入;将(1)(2)式带入(3)式得到:
(4)
通过变换得到:
(5)
公式(5)代表钻井液侵入连通孔隙和裂缝的等效孔隙度,即侵入钻井液的等效百分体积;
测得岩心物性分析岩样质量M,根据体积公式:
(6)
M ---- 岩心物性分析岩样质量,g;
ρ---- 测井测得该深度岩样的密度值,g/cm3
V ---- 该深度岩样的实际体积,cm3;
而该一深度段某一种矿物的实际体积V实应为:
(7)
V百分含量 ---- 为该矿物这一深度段所占百分体积,%;
V实 ---- 该一深度段某一种矿物的实际体积,g;
设百分体积与实际体积之间关系为α,表达为:
(8)
将公式(6)代入公式(8),可得转换系数:
将侵入钻井液的等效百分体积乘以转换系数,即可得到侵入钻井液等效体积:
(9)
式中,
V钻井液 ---- 侵入钻井液等效体积,cm3;
由于钻井液主要是固相颗粒反应PE值,因此式(9)计算得到的是固相颗粒的等效钻井液体积V钻井液,而实际滤液体积为:
(10)
V滤液 ---- 实际侵入钻井液体积,cm3;
C钻井液 ----钻井液固相颗粒含量浓度,%,(实际采集获得);
把某一深度点总浸泡时间分为三部分:钻进时、完钻后等待测井时,测井时:
(11)
T总 ---- 某一深度所测得资料的总浸泡时间,s;
T进 ---- 某一深度钻进结束时浸泡时间,s;
T等---- 完钻到测井开始时间,s,可采集;
T测 ---- 测井开始到测得该点所耗时间,s,可采集;
钻进时,
(12)
T钻进 ---- 钻进总耗时,s;
T钻遇---- 钻遇致该深度点耗时,s,钻时曲线计算获得;
将时间与深度对应、深度与曲线对应,最后流体流量曲线Q为:
(13)
V路径 ----残留在渗滤路径中的钻井液体积,cm3
Q----通过岩样的流体流量等于进入的流体体积减去残留在路径中的流体体积的差值与流体穿过岩样的总时间的比值。
所述的渗滤路径L选择与实验室岩样相同长度为渗滤路径L。
所述的横截面积A的计算方法为:
(16)
式中:
r ---- 钻头半径,cm。
所述的压差ΔP的计算方法为:
压差的计算分为钻进、等待和测井状态下的压差;
钻进时:
(17)
式中:
P进 ---- 某一深度点钻进时对井壁的压力,pa;
ρ钻井液 ---- 钻井液密度,g/cm3;
h ---- 当前深度点,m;
g ---- 重力加速度,m/s2;
P泵 ---- 泵压,pa;
P地层---- 地层压力,pa。
地层压力需实测或选用区域地层压力;
等待压差和测井状态下压差相等:
(18)
式中:
P测 – 等待和测井时钻井液对井壁的压力,pa;
P地层---- 地层压力,pa。
与现有技术相比,本发明所达到的有益效果如下:
本发明首次提出一种基于达西定律和钻井液侵入体积计算渗透率的方法,切实可行的解决了无法预判地层中是“孔”还是“缝”所造成的被动问题,极大的优化了油气田产能评估流程,提高了油气田开发效率,对推进储层渗透率预测技术的进程具有革新性的意义。
附图说明
下面将结合说明书附图和具体实施方式对本发明作进一步的详细说明,其中:
图1为井筒渗滤路径示意图;
图2为基于达西定律和钻井液侵入体积测井渗透率计算流程示意图。
具体实施方式
作为本发明的最佳实施方式,其公开了一种利用测井渗透率进行地层有效裂缝的识别方法,在裂缝--孔隙性储层中,当地层被钻开时,在井筒压力与地层压力压差作用下钻井液势必会侵入地层。根据达西定律,渗透率的表达式为:
式中,
K ---- 渗透率,μm2;
Q ---- 流体流量, cm3/s;
A ---- 垂直于流体流动方向岩石横截面积cm2;
L ---- 流体渗滤路径,cm;
---- 压差,Pa;
---- 流体粘度, mPa·s。
在实际计算中,流体粘度为钻井液粘度,可直接获取;流体流量Q、地层压差ΔP、横截面积A和渗滤路径L则需通过计算求得。
由于孔隙—裂缝型储层主要的流动通道为裂缝,钻井液中的固相颗粒与液相滤液即可随裂缝侵入地层,又因裂缝—孔隙型储层的有效流动截面远大于孔隙型储层,因此钻井液中的固相颗粒又可继续沿裂缝深入到储层内部,且固相颗粒的侵入深度与液相滤液呈正相关性。再者,岩性密度曲线(PE)受钻井液中固相颗粒浓度影响最大,如钻井液中加有固相颗粒的重晶石,其PE值将远大于地层岩性的标准PE值,因此可基于重晶石的渗入,确定地层渗透性。在此基础之上,结合实际工况得到流体流量Q(cm3/s),结合横截面积A(cm2)、流体渗滤路径L(cm)、压差ΔP(Pa)和流体粘度μ(mPa·s)四项参数,即可得到最终的地层渗透率。
计算流程如下:
(1)流体流量计算(Q)
首先利用体积模型法计算侵入钻井液体积。认为地层由骨架矿物、黏土矿物和孔隙组成。在地层被钻开前,孔隙中都充填的为地层原状流体,在地层被钻开后,连通孔隙和裂缝中的可流动流体被替换或部分替换为钻井液,而非连通孔隙中的流体依旧为原流体,因此采用此流体替换方法即可计算钻井液侵入体积。在此选用对钻井液敏感性相对较高的PE曲线来进行流体替换差值计算。
(1)
式中,
PE原 ---- 原始状态下理论岩心密度值,b/e;
PE骨架 ---- 各骨架矿物的岩性密度值,b/e;
V总 ---- 孔隙和矿物总百分含量,数值为100。
V骨架 ---- 各骨架矿物对应的体积百分比,%,测井计算求得;
PE黏土 ---- 各黏土矿物对应的体积百分比,b/e;
V黏土 ---- 各黏土矿物对应的体积百分比,%,测井计算求得;
PE储层流体 ---- 原状地层连通孔隙中流体的岩性密度值,b/e;
PE束缚流体 ---- 非连通孔隙中流体的岩性密度值,b/e;
连通 ---- 连通孔隙度,测井计算求得,%;
非连通 ---- 非连通孔隙度,测井计算求得,%;
PE钻井液 ---- 钻井液的PE值,b/e,可根据钻井液中所加重矿物含量计算得到,或直接测量。
公式(1)表示当地层未被钻开时,原始状态下的岩性密度理论值。
(2)
式中,
PE测 ---- 实际测量岩性密度值,b/e;
---- 侵入钻井液所占孔隙度,%;
公式(2)表示当地层钻开后,部分连通孔隙和裂缝中流体被钻井液替换后的实际测量的岩性密度值。
(3)
式中,
侵入---- 地层被钻开前后岩性密度变化值Δ侵入,b/e。
式(3)表示由测井仪器实测得到的岩性密度值PE测减去原始地层状态下的岩性密度值PE原,将得到地层被钻开后钻井液部分或完全侵入时的连通孔隙中和裂缝中岩性密度的变化值Δ侵入。将(1)(2)式带入(3)式得到:
(4)
通过变换得到:
(5)
公式(5)代表了钻井液侵入连通孔隙和裂缝的等效孔隙度,即侵入钻井液的等效百分体积。
测得岩心物性分析岩样质量M,根据体积公式:
(6)
M ---- 岩心物性分析岩样质量,g;
ρ---- 测井测得该深度岩样的密度值,g/cm3
V ---- 该深度岩样的实际体积,cm3。
而该一深度段某一种矿物的实际体积V实应为:
(7)
V百分含量 ---- 为该矿物这一深度段所占百分体积,%。
V实 ---- 该一深度段某一种矿物的实际体积,g。
设百分体积与实际体积之间关系为α,可表达为:
(8)
将公式(6)代入公式(8),可得转换系数:
将侵入钻井液的等效百分体积乘以转换系数,即可得到侵入钻井液等效体积:
(9)
式中,
V钻井液 ---- 侵入钻井液等效体积,cm3。
由于钻井液主要是固相颗粒反应PE值,因此式(9)计算得到的是固相颗粒的等效钻井液体积V钻井液,而实际滤液体积为:
(10)
V滤液 ---- 实际侵入钻井液体积,cm3;
C钻井液 ----钻井液固相颗粒含量浓度,%,(实际采集获得)。
因为在地层钻开的瞬间,就会发生钻井液侵入,因此越晚钻开的地层,钻井液浸泡时间越短。且浸泡时间与钻时、完钻到测井开始间隔时间和测井测量速度有关,因此根据实际情况,把某一深度点总浸泡时间分为三部分:钻进时、完钻后等待测井时,测井时:
(11)
T总 ---- 某一深度所测得资料的总浸泡时间,s;
T进 ---- 某一深度钻进结束时浸泡时间,s;
T等---- 完钻到测井开始时间,s,可采集;
T测 ---- 测井开始到测得该点所耗时间,s,可采集。
钻进时,
(12)
T钻进 ---- 钻进总耗时,s;
T钻遇---- 钻遇致该深度点耗时,s,钻时曲线计算获得。
将时间与深度对应、深度与曲线对应,最后流体流量曲线Q为:
(13)
V路径 ----残留在渗滤路径中的钻井液体积,cm3
Q----通过岩样的流体流量(等于进入的流体体积减去残留在路径中的流体体积的差值与流体穿过岩样的总时间的比值)。
(2)渗滤路径(L)
实验室测量岩样渗透率中,渗滤路径长度近似选用岩样长度,在此选择与实验室岩样相同长度为渗滤路径L
图1为井筒渗滤路径示意图,r为钻头半径,△R为渗滤路径。
(3)横截面积(A)计算
由于井筒是向外放射的,所以常规计算的井下钻井液横向通过井壁横截面积不同于实验室测量,而是随着侵入深度的加深,横截面积也增大,即为:
(16)
式中:
r ---- 钻头半径,cm。
(4)压差(ΔP)计算
压差的计算分为钻进、等待和测井状态下的压差。钻进时:
(17)
式中:
P进 ---- 某一深度点钻进时对井壁的压力,pa;
ρ钻井液 ---- 钻井液密度,g/cm3;
h ---- 当前深度点,m;
g ---- 重力加速度,m/s2;
P泵 ---- 泵压,pa;
P地层---- 地层压力,pa。
地层压力需实测或选用区域地层压力。
等待压差和测井状态下压差相等:
(18)
式中:
P测 – 等待和测井时钻井液对井壁的压力,pa;
P地层---- 地层压力,pa。
在通过上述步骤求得流体流量Q、地层压差ΔP、横截面积A和渗滤路径L后,再通过实际测量获取粘度,即可根据达西定律求得储层段的渗透率。
Claims (5)
1.一种利用测井渗透率进行地层有效裂缝的识别方法,其特征在于:首先根据钻井液中的固相颗粒与液相滤液侵入地层的深度,判断地层的渗透性,并由实际工况得到流体流量Q,结合横截面积A、流体渗滤路径L、压差ΔP和流体粘度μ四项参数,根据达西定律计算出储层的渗透率:
式中,
K ---- 渗透率,μm2;
Q ---- 流体流量, cm3/s;
A ---- 垂直于流体流动方向岩石横截面积cm2;
L ---- 流体渗滤路径,cm;
---- 压差,Pa;
---- 流体粘度, mPa·s;
根据计算出的渗透率,识别地层的有效裂缝。
2.根据权利要求1所述的一种利用测井渗透率进行地层有效裂缝的识别方法,其特征在于:所述流体流量Q的计算方法如下:
首先利用体积模型法计算侵入钻井液体积,采用流体替换方法计算钻井液侵入体积,选用PE曲线来进行流体替换差值计算:
(1)
式中,
PE原 ---- 原始状态下理论岩心密度值,b/e;
PE骨架 ---- 各骨架矿物的岩性密度值,b/e;
V总 ---- 孔隙和矿物总百分含量,数值为100;
V骨架 ---- 各骨架矿物对应的体积百分比,%,测井计算求得;
PE黏土 ---- 各黏土矿物对应的体积百分比,b/e;
V黏土 ---- 各黏土矿物对应的体积百分比,%,测井计算求得;
PE储层流体 ---- 原状地层连通孔隙中流体的岩性密度值,b/e;
PE束缚流体 ---- 非连通孔隙中流体的岩性密度值,b/e;
连通 ---- 连通孔隙度,测井计算求得,%;
非连通 ---- 非连通孔隙度,测井计算求得,%;
PE钻井液 ---- 钻井液的PE值,b/e,根据钻井液中所加重矿物含量计算得到,或直接测量;
公式(1)表示当地层未被钻开时,原始状态下的岩性密度理论值;
(2)
式中,
PE测 ---- 实际测量岩性密度值,b/e;
---- 侵入钻井液所占孔隙度,%;
公式(2)表示当地层钻开后,部分连通孔隙和裂缝中流体被钻井液替换后的实际测量的岩性密度值;
(3)
式中,
侵入---- 地层被钻开前后岩性密度变化值Δ侵入,b/e;
式(3)表示由测井仪器实测得到的岩性密度值PE测减去原始地层状态下的岩性密度值PE原,将得到地层被钻开后钻井液部分或完全侵入时的连通孔隙中和裂缝中岩性密度的变化值Δ侵入;将(1)(2)式带入(3)式得到:
(4)
通过变换得到:
(5)
公式(5)代表钻井液侵入连通孔隙和裂缝的等效孔隙度,即侵入钻井液的等效百分体积;
测得岩心物性分析岩样质量M,根据体积公式:
(6)
M ---- 岩心物性分析岩样质量,g;
ρ---- 测井测得该深度岩样的密度值,g/cm3
V ---- 该深度岩样的实际体积,cm3;
而该一深度段某一种矿物的实际体积V实应为:
(7)
V百分含量 ---- 为该矿物这一深度段所占百分体积,%;
V实 ---- 该一深度段某一种矿物的实际体积,g;
设百分体积与实际体积之间关系为α,表达为:
(8)
将公式(6)代入公式(8),可得转换系数:
将侵入钻井液的等效百分体积乘以转换系数,即可得到侵入钻井液等效体积:
(9)
式中,
V钻井液 ---- 侵入钻井液等效体积,cm3;
由于钻井液主要是固相颗粒反应PE值,因此式(9)计算得到的是固相颗粒的等效钻井液体积V钻井液,而实际滤液体积为:
(10)
V滤液 ---- 实际侵入钻井液体积,cm3;
C钻井液 ----钻井液固相颗粒含量浓度,%,(实际采集获得);
把某一深度点总浸泡时间分为三部分:钻进时、完钻后等待测井时,测井时:
(11)
T总 ---- 某一深度所测得资料的总浸泡时间,s;
T进 ---- 某一深度钻进结束时浸泡时间,s;
T等---- 完钻到测井开始时间,s,可采集;
T测 ---- 测井开始到测得该点所耗时间,s,可采集;
钻进时,
(12)
T钻进 ---- 钻进总耗时,s;
T钻遇---- 钻遇致该深度点耗时,s,钻时曲线计算获得;
将时间与深度对应、深度与曲线对应,最后流体流量曲线Q为:
(13)
V路径 ----残留在渗滤路径中的钻井液体积,cm3
Q----通过岩样的流体流量等于进入的流体体积减去残留在路径中的流体体积的差值与流体穿过岩样的总时间的比值。
3.根据权利要求1所述的一种利用测井渗透率进行地层有效裂缝的识别方法,其特征在于:所述的渗滤路径L选择与实验室岩样相同长度为渗滤路径L。
4.根据权利要求1所述的一种利用测井渗透率进行地层有效裂缝的识别方法,其特征在于:所述的横截面积A的计算方法为:
(16)
式中:
r ---- 钻头半径,cm。
5.根据权利要求1所述的一种利用测井渗透率进行地层有效裂缝的识别方法,其特征在于:所述的压差ΔP的计算方法为:
压差的计算分为钻进、等待和测井状态下的压差;
钻进时:
(17)
式中:
P进 ---- 某一深度点钻进时对井壁的压力,pa;
ρ钻井液 ---- 钻井液密度,g/cm3;
h ---- 当前深度点,m;
g ---- 重力加速度,m/s2;
P泵 ---- 泵压,pa;
P地层---- 地层压力,pa;
地层压力实测或选用区域地层压力;
等待压差和测井状态下压差相等:
(18)
式中:
P测 – 等待和测井时钻井液对井壁的压力,pa;
P地层---- 地层压力,pa。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810876663.8A CN109142676B (zh) | 2018-08-03 | 2018-08-03 | 一种利用测井渗透率进行地层有效裂缝的识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810876663.8A CN109142676B (zh) | 2018-08-03 | 2018-08-03 | 一种利用测井渗透率进行地层有效裂缝的识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109142676A true CN109142676A (zh) | 2019-01-04 |
CN109142676B CN109142676B (zh) | 2021-05-18 |
Family
ID=64791555
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810876663.8A Active CN109142676B (zh) | 2018-08-03 | 2018-08-03 | 一种利用测井渗透率进行地层有效裂缝的识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109142676B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112439540A (zh) * | 2019-08-30 | 2021-03-05 | 中国石油天然气股份有限公司 | 一种轻重矿物分离装置及方法 |
CN114239365A (zh) * | 2021-12-21 | 2022-03-25 | 西南石油大学 | 一种天然裂缝钻井液污染后壁面等效渗透率的计算方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749651A (zh) * | 2012-07-12 | 2012-10-24 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种储层裂缝测井标定方法 |
CN104075974A (zh) * | 2014-07-14 | 2014-10-01 | 中国地质大学(北京) | 一种利用低场核磁共振精确测定页岩孔隙度的方法 |
CN105089582A (zh) * | 2015-05-28 | 2015-11-25 | 中国石油天然气股份有限公司 | 基于井下流量控制设备的油藏数值模拟方法及装置 |
CN105350960A (zh) * | 2015-12-07 | 2016-02-24 | 西南石油大学 | 确定低渗透非均质气藏压裂水平井裂缝参数的方法 |
CN107063963A (zh) * | 2016-12-28 | 2017-08-18 | 浙江海洋大学 | 一种致密储层微裂缝扩展及渗流特征的测试装置和方法 |
CN107558999A (zh) * | 2017-07-03 | 2018-01-09 | 恒泰艾普集团股份有限公司 | 一种孤立缝洞体剩余油定量预测方法 |
-
2018
- 2018-08-03 CN CN201810876663.8A patent/CN109142676B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749651A (zh) * | 2012-07-12 | 2012-10-24 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种储层裂缝测井标定方法 |
CN104075974A (zh) * | 2014-07-14 | 2014-10-01 | 中国地质大学(北京) | 一种利用低场核磁共振精确测定页岩孔隙度的方法 |
CN105089582A (zh) * | 2015-05-28 | 2015-11-25 | 中国石油天然气股份有限公司 | 基于井下流量控制设备的油藏数值模拟方法及装置 |
CN105350960A (zh) * | 2015-12-07 | 2016-02-24 | 西南石油大学 | 确定低渗透非均质气藏压裂水平井裂缝参数的方法 |
CN107063963A (zh) * | 2016-12-28 | 2017-08-18 | 浙江海洋大学 | 一种致密储层微裂缝扩展及渗流特征的测试装置和方法 |
CN107558999A (zh) * | 2017-07-03 | 2018-01-09 | 恒泰艾普集团股份有限公司 | 一种孤立缝洞体剩余油定量预测方法 |
Non-Patent Citations (3)
Title |
---|
大庆石油学院开发系编: "《采油工人常用名词解释》", 31 October 1978 * |
李青和等: "测井技术在致密储层裂缝识别中的应用 ", 《内蒙古石油化工》 * |
王全柱: "火成岩储层研究", 《西安石油大学学报(自然科学版)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112439540A (zh) * | 2019-08-30 | 2021-03-05 | 中国石油天然气股份有限公司 | 一种轻重矿物分离装置及方法 |
CN114239365A (zh) * | 2021-12-21 | 2022-03-25 | 西南石油大学 | 一种天然裂缝钻井液污染后壁面等效渗透率的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109142676B (zh) | 2021-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cheng et al. | New insights into spontaneous imbibition in tight oil sandstones with NMR | |
CN106285646B (zh) | 基于多信息融合的钻井漏失层位识别方法 | |
Satter et al. | Reservoir engineering: the fundamentals, simulation, and management of conventional and unconventional recoveries | |
Akbar et al. | A snapshot of carbonate reservoir evaluation | |
EP2457082B1 (en) | Method for evaluating shaped charge perforation test cores using computer tomographic images thereof | |
Ren et al. | Borehole characterization of hydraulic properties and groundwater flow in a crystalline fractured aquifer of a headwater mountain watershed, Laramie Range, Wyoming | |
CN105931125B (zh) | 一种致密油分段多簇体积压裂水平井产量预测方法 | |
RU2513895C1 (ru) | Способ разработки нефтяных залежей | |
CN109324345A (zh) | 叠合盆地油气成藏期岩石孔隙度恢复方法 | |
CN102621586B (zh) | 一种识别地层属性的地层数据处理方法 | |
CN109386285A (zh) | 一种特低渗透油层泥浆侵入程度及其影响的评价方法 | |
CN109142676A (zh) | 一种利用测井渗透率进行地层有效裂缝的识别方法 | |
CN102182444B (zh) | 欠平衡钻井适应性评估方法 | |
CN107237626A (zh) | 一种井钻遇溶洞体积的确定方法 | |
Long et al. | Integrating oil and water geochemistry to assess SRV and DRV in the bakken/three forks hybrid play | |
CN102134994B (zh) | 一种基于基岩油藏油水层电阻率的地层数据处理方法 | |
RU2655310C1 (ru) | Способ определения эффективности гидравлического разрыва пласта скважины | |
Liu et al. | An integrated approach to investigating the stimulation mechanisms of soaking in shale/tight oil reservoirs | |
US10598010B2 (en) | Method for constructing a continuous PVT phase envelope log | |
Moradi | Cost-effective and safe oil production from existing and near-future oil fields | |
CN113177363A (zh) | 一种储层大孔道参数定量表征的油藏工程方法 | |
Wang et al. | Reservoir and development characteristics of the Da'anzhai tight oil in Sichuan Basin, SW China | |
Muslimov | Problems of exploration and development modeling of oil fields | |
RU2486337C1 (ru) | Способ определения продуктивности пласта в процессе бурения скважины | |
CN113803054B (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 |