CN113312796A - 考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用 - Google Patents
考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用 Download PDFInfo
- Publication number
- CN113312796A CN113312796A CN202110693391.XA CN202110693391A CN113312796A CN 113312796 A CN113312796 A CN 113312796A CN 202110693391 A CN202110693391 A CN 202110693391A CN 113312796 A CN113312796 A CN 113312796A
- Authority
- CN
- China
- Prior art keywords
- fracture
- conductivity
- reservoir
- porosity
- fractured
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/041—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Electrochemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Sewage (AREA)
Abstract
本发明公开了一种考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用,其构建方法包含:(1)在各向异性介质中,基于欧姆定律建立沿着电场方向上的电导率,从而获得裂缝型储层中任意裂缝倾角下的电导率;(2)将裂缝型储层中沿着裂缝方向的电导率视为基质电导率与裂缝电导率的并联,裂缝型储层中垂直于裂缝发育方向的电导率视为基质电导率与裂缝电导率的串联,并结合所述裂缝型储层中任意裂缝倾角下的电导率,获得考虑裂缝倾角的裂缝型储层电导率模型。本发明弥补了裂缝型储层孔隙度指数研究工作的短缺,针对裂缝型储层,提供一套更完善的孔隙度指数计算模型,对含水饱和度的准确计算更为有利,能够用于储层孔隙度指数评价工作。
Description
技术领域
本发明涉及一种裂缝型储层电导率模型,具体涉及一种考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用。
背景技术
在最新的碳酸盐岩储层孔隙度指数研究工作中,凡涉及到裂缝角度的储层孔隙度指数计算,均是基于等效思想,将渗透率(permeability)从字面意义上直接等效为介电常数(permittivity),在未经理论推导的情况下,直接利用Parsons的渗透率方程(式(1a))等效出不同裂缝倾角下的储层电导率表达式。
Kp=Kx·cos2θ+Ky·sin2θ (1a)
式(1a)中,x、y为最大渗透率与最小渗透率方向;Kp为沿着压力梯度方向上的渗透率;Kx为在x方向上的最大渗透率;Ky为在y方向上的最大渗透率;Vx、Vy分别表示在x、y方向上的速度;θ为压力梯度方向与x方向的夹角。
上述这种等效法导致电导率在x、y、p方向上的含义混淆不清,该电导率公式难以直接应用于实际储层,在裂缝型储层评价工作中,急需一套完整的考虑裂缝倾角的储层电导率模型。
发明内容
本发明的目的是提供一种考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用,解决了现有在裂缝型储层评价中的电导率模型难以直接应用于实际储层的问题,弥补了裂缝型储层孔隙度指数研究工作的短缺,针对裂缝型储层,提供一套更完善的孔隙度指数计算模型,对含水饱和度的准确计算更为有利,能够用于储层孔隙度指数评价工作。
为了达到上述目的,本发明提供了一种考虑裂缝倾角的裂缝型储层电导率模型,该裂缝型储层电导率模型,为:
式中,σθ为裂缝型储层中任意裂缝倾角下的电导率;σf为裂缝电导率;σo为基质电导率;Φf为裂缝孔隙度;θ为裂缝倾角;基于该裂缝型储层电导率模型,依据σw、σo研究裂缝孔隙度Φf、溶洞孔隙度Φv和裂缝倾角θ对裂缝型储层电导率影响。
本发明的另一目的是提供所述的考虑裂缝倾角的裂缝型储层电导率模型的应用,基于所述的考虑裂缝倾角的裂缝型储层电导率模型结合Archie公式构建得到裂缝型储层孔隙度指数模型,为:
式中,mθ为裂缝型储层孔隙度指数;mb为基质孔隙度指数;Φb表示基质孔隙度;Φ为储层孔隙度。
本发明的另一目的是提供所述的考虑裂缝倾角的裂缝型储层电导率模型的构建方法,该方法包含:
(1)在各向异性介质中,基于欧姆定律建立沿着电场方向上的电导率,从而获得裂缝型储层中任意裂缝倾角下的电导率,为:
σθ=σh·cos2θ+σn·sin2θ (6)
式(6)中,σh为裂缝型储层中沿着裂缝方向的电导率;σn为裂缝型储层中垂直于裂缝发育方向的电导率;
(2)将裂缝型储层中沿着裂缝方向的电导率σh视为基质电导率σo与裂缝电导率σf的并联,裂缝型储层中垂直于裂缝发育方向的电导率σn视为基质电导率σo与裂缝电导率σf的串联,并结合所述裂缝型储层中任意裂缝倾角下的电导率,获得如权利要求1所述的考虑裂缝倾角的裂缝型储层电导率模型;
在步骤(2)中,在裂缝型储层中,沿着裂缝方向的电导率σh视为基质电导率σo与裂缝电导率σf的并联,垂直于裂缝发育方向的电导率σn视为基质电导率σo与裂缝电导率σf的串联,得到:
优选地,该方法包含:利用所述的考虑裂缝倾角的裂缝型储层电导率模型结合Archie公式,并依据当储层完全含水时裂缝电导率为地层水电导率,建立孔隙度指数模型,用以研究裂缝型储层的孔隙度指数变化。
由Archie公式得到:
式(9)中,mb为基质孔隙度指数,mθ为裂缝型储层孔隙度指数,mh、mn分别为沿着裂缝方向与垂直裂缝方向的裂缝型储层孔隙度指数,σw为地层水电阻率,Φ为储层孔隙度;Φb表示基质孔隙度;σo为基质电导率;σθ为裂缝型储层中任意裂缝倾角下的电导率;σh为裂缝型储层中沿着裂缝方向的电导率;σn为裂缝型储层中垂直于裂缝发育方向的电导率。
本发明的考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用,解决了现有在裂缝型储层评价中的电导率模型难以直接应用于实际储层的问题,具有以下优点:
相较于现有从渗透率等效为电导率的方法,本发明基于欧姆定律建立的考虑裂缝倾角的裂缝型储层电导率模型,弥补了裂缝型储层孔隙度指数研究工作的短缺,针对裂缝型储层,可提供一套更完善的孔隙度指数计算模型,对含水饱和度的准确计算更为有利。
附图说明
图1为各向异性介质中的电流密度与场强方向示意图。
图2为本发明裂缝型储层孔隙度指数变化图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
一种考虑裂缝倾角的裂缝型储层电导率模型,其构建方法具体如下:
(1)在各向异性介质中,基于欧姆定律建立沿着电场方向上的电导率,从而获得裂缝型储层中任意裂缝倾角下的电导率,具体如下:
欧姆定律的微分形式可表示为:
式(1)中,为电流密度,A/m2;为电场强度,V/m;σ为电导率,S/m。对于均质介质,电流密度方向与电场方向重合,电导率为常量;在各向异性介质中,由于电导率为张量,除了相互垂直的最大、最小电导率方向,电流密度方向与电场强度方向不再重合,电流密度向着电导率更大的方向偏移。
式(2)中,U为电势,V;Jx为x方向上的电流密度,A/m2;Jy为y方向上的电流密度,A/m2;σx为x方向上的电导率,S/m;σy为y方向上的电导率,S/m。
通过图1,可以得出:
联立式(2)、式(3)与式(4),可得到沿着电场方向上的电导率,为:
σE=σx·cos2θ+σy·sin2θ (5)
对于裂缝型储层,沿着裂缝方向的电导率最大,记做σh,h方向即为裂缝发育方向;垂直于裂缝发育方向的电导率最小,记做σn,n方向即为垂直于裂缝的方向。在井眼测量环境,电场方向可视为水平流入地层,特别是对于深双侧向测井,监督电极的设计就是为了控制电场方向保持水平,此时,σE表示沿着水平方向测得电导率,场强方向与h方向的夹角即为裂缝倾角。
由此,可得到任意裂缝倾角下的电导率σθ表达式:
σθ=σh·cos2θ+σn·sin2θ (6)
(2)将裂缝型储层中沿着裂缝方向的电导率视为基质电导率与裂缝电导率的并联,垂直裂缝方向的电导率视为基质电导率与裂缝电导率的串联,并结合裂缝型储层中任意裂缝倾角下的电导率,获得考虑裂缝倾角的裂缝型储层电导率模型,具体如下:
在裂缝型储层中,沿着裂缝方向的电导率σh可视为基质电导率σo与裂缝电导率σf的并联,垂直裂缝方向的电导率σn可视为基质电导率σo与裂缝电导率σf的串联,由此可得:
式(7)中,σo为基质电导率,Φf为裂缝孔隙度。联立式(6)、式(7)可得考虑裂缝倾角的裂缝型储层电导率模型,为:
由式(8)可知,在取定σw、σo后,即可研究裂缝孔隙度Φf、溶洞孔隙度Φv、裂缝倾角θ对裂缝型储层电导率影响规律。
本发明的考虑裂缝倾角的裂缝型储层电导率模型的应用,具体如下:
利用本发明构建的考虑裂缝倾角的裂缝型储层电导率模型结合Archie公式可建立孔隙度指数模型,用以研究裂缝型储层的孔隙度指数变化,指导裂缝型储层评价工作。
由阿尔奇公式可知:
式(9)中,mb为基质孔隙度指数,mθ为裂缝型储层孔隙度指数,mh、mn分别为沿着裂缝方向与垂直裂缝方向的裂缝型储层孔隙度指数;σw为地层水电阻率,Φ为储层孔隙度;Φb表示基质孔隙度。
当储层完全含水时,裂缝里充填地层水,裂缝电导率为地层水电导率,即σf=σw,联立式(8)、式(9),可得裂缝型储层孔隙度指数模型,为:
根据式(10),可得到裂缝型储层孔隙度指数变化图,如图2所示,为本发明裂缝型储层孔隙度指数变化图,由图2可知:
(1)裂缝倾角为90°时,裂缝型储层mθ大于mb;此时,随着Φ增大,mθ减小,逐渐趋近于mb;
(2)当裂缝倾角大于33°时,随着Φ增大,mθ先增大,部分值大于mb,然后逐渐减小,最后趋近于mb;当裂缝倾角小于33°时,随着Φ增大,mθ逐渐增大,最后趋近于mb;
(3)裂缝型储层中,决定mθ与mb大小关系的裂缝临界角为33°;
(4)相同裂缝倾角下,当Φf增大时,储层mθ减小。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (4)
3.如权利要求1所述的考虑裂缝倾角的裂缝型储层电导率模型的构建方法,其特征在于,该方法包含:
(1)在各向异性介质中,基于欧姆定律建立沿着电场方向上的电导率,从而获得裂缝型储层中任意裂缝倾角下的电导率,为:
σθ=σh·cos2θ+σn·sin2θ (6)
式(6)中,σh为裂缝型储层中沿着裂缝方向的电导率;σn为裂缝型储层中垂直于裂缝发育方向的电导率;
(2)将裂缝型储层中沿着裂缝方向的电导率σh视为基质电导率σo与裂缝电导率σf的并联,裂缝型储层中垂直于裂缝发育方向的电导率σn视为基质电导率σo与裂缝电导率σf的串联,并结合所述裂缝型储层中任意裂缝倾角下的电导率,获得如权利要求1所述的考虑裂缝倾角的裂缝型储层电导率模型;
在步骤(2)中,在裂缝型储层中,沿着裂缝方向的电导率σh视为基质电导率σo与裂缝电导率σf的并联,垂直于裂缝发育方向的电导率σn视为基质电导率σo与裂缝电导率σf的串联,得到:
4.如权利要求2中所述的裂缝型储层孔隙度指数模型的构建方法,其特征在于,该方法包含:
利用如权利要求1所述的考虑裂缝倾角的裂缝型储层电导率模型结合Archie公式,并依据当储层完全含水时裂缝电导率为地层水电导率,建立孔隙度指数模型,用以研究裂缝型储层的孔隙度指数变化;
由Archie公式得到:
式(9)中,mb为基质孔隙度指数,mθ为裂缝型储层孔隙度指数,mh、mn分别为沿着裂缝方向与垂直裂缝方向的裂缝型储层孔隙度指数,σw为地层水电阻率,Φ为储层孔隙度;Φb表示基质孔隙度;σo为基质电导率;σθ为裂缝型储层中任意裂缝倾角下的电导率;σh为裂缝型储层中沿着裂缝方向的电导率;σn为裂缝型储层中垂直于裂缝发育方向的电导率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110693391.XA CN113312796B (zh) | 2021-06-22 | 2021-06-22 | 考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110693391.XA CN113312796B (zh) | 2021-06-22 | 2021-06-22 | 考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113312796A true CN113312796A (zh) | 2021-08-27 |
CN113312796B CN113312796B (zh) | 2022-06-03 |
Family
ID=77380298
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110693391.XA Active CN113312796B (zh) | 2021-06-22 | 2021-06-22 | 考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113312796B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114429044A (zh) * | 2022-01-21 | 2022-05-03 | 成都理工大学 | 一种建立缝洞型储层电导率模型的方法及应用 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5095273A (en) * | 1991-03-19 | 1992-03-10 | Mobil Oil Corporation | Method for determining tensor conductivity components of a transversely isotropic core sample of a subterranean formation |
US20110254552A1 (en) * | 2010-04-15 | 2011-10-20 | Peter Wu | Method and apparatus for determining geological structural dip using multiaxial induction measurements |
US20150083404A1 (en) * | 2013-09-23 | 2015-03-26 | Schlumberger Technology Corporation | Determining proppant and fluid distribution |
US20180334901A1 (en) * | 2016-01-20 | 2018-11-22 | Halliburton Energy Services, Inc. | Fracture Interpretation With Resistivity And Sonic Logs In Biaxial Anisotropic Formations |
CN110967771A (zh) * | 2019-12-19 | 2020-04-07 | 西南石油大学 | 一种利用常规测井计算裂缝、溶蚀孔洞孔隙度的方法 |
CN111177950A (zh) * | 2020-01-17 | 2020-05-19 | 成都理工大学 | 阵列感应测井曲线环境影响规律的信息处理系统及方法 |
-
2021
- 2021-06-22 CN CN202110693391.XA patent/CN113312796B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5095273A (en) * | 1991-03-19 | 1992-03-10 | Mobil Oil Corporation | Method for determining tensor conductivity components of a transversely isotropic core sample of a subterranean formation |
US20110254552A1 (en) * | 2010-04-15 | 2011-10-20 | Peter Wu | Method and apparatus for determining geological structural dip using multiaxial induction measurements |
US20150083404A1 (en) * | 2013-09-23 | 2015-03-26 | Schlumberger Technology Corporation | Determining proppant and fluid distribution |
US20180334901A1 (en) * | 2016-01-20 | 2018-11-22 | Halliburton Energy Services, Inc. | Fracture Interpretation With Resistivity And Sonic Logs In Biaxial Anisotropic Formations |
CN110967771A (zh) * | 2019-12-19 | 2020-04-07 | 西南石油大学 | 一种利用常规测井计算裂缝、溶蚀孔洞孔隙度的方法 |
CN111177950A (zh) * | 2020-01-17 | 2020-05-19 | 成都理工大学 | 阵列感应测井曲线环境影响规律的信息处理系统及方法 |
Non-Patent Citations (7)
Title |
---|
AGUILERA, C.G. ET AL: "Effect of fracture dip on petrophysical evaluation of naturally fractured reservoirs", 《JOURNAL OF CANADIAN PETROLEUM TECHNOLOGY》 * |
R.W.PARSONS: "Permeability of Idealized Fractured Rock", 《SOCIETY OF PETROLEUM ENGINEERS JOURNAL》 * |
W.DAVIDKENNEDY等: "电性各向异性介质中含水饱和度的计算", 《测井技术》 * |
任政勇等: "基于电流密度连续性条件的直流电阻率各向异性问题自适应有限元模拟", 《地球物理学报》 * |
孟奇猛: "电测深资料求解灰岩含水层渗透系数的可行性──以吴家营水源地为例", 《云南地质》 * |
孟昆等: "泥页岩储层水平井随钻电磁波电阻率测井响应特性", 《断块油气田》 * |
邓少贵等: "裂缝-孔隙双孔介质的双侧向测井解释方法初探", 《测井技术》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114429044A (zh) * | 2022-01-21 | 2022-05-03 | 成都理工大学 | 一种建立缝洞型储层电导率模型的方法及应用 |
Also Published As
Publication number | Publication date |
---|---|
CN113312796B (zh) | 2022-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nyikos et al. | Electrochemistry at fractal interfaces: the coupling of ac and dc behaviour at irregular electrodes | |
US8159227B2 (en) | Methods for making directional resistivity measurements | |
CN104101905B (zh) | 一种基于岩电参数的储层分类方法 | |
CN113312796B (zh) | 考虑裂缝倾角的裂缝型储层电导率模型及其构建方法和应用 | |
CN110967771B (zh) | 一种利用常规测井计算裂缝、溶蚀孔洞孔隙度的方法 | |
GB693782A (en) | Improvements in or relating to electrical analogues | |
CN108802073A (zh) | 一种基于数字岩心的岩石电性参数获取方法及装置 | |
CN106529150B (zh) | 复合地层盾构隧道拱顶荷载计算方法 | |
CN102767367A (zh) | 高分辨率侧向测井仪及电阻率测量方法 | |
CN109386281B (zh) | 一种获取裂缝性低孔低渗储层测井饱和度的方法 | |
CN113051779B (zh) | 一种三维直流电阻率法数值模拟方法 | |
US9239402B2 (en) | Focused array laterolog tool | |
CN105447762A (zh) | 一种流体替换的低渗透油藏水淹信息的计算方法 | |
CN103852041B (zh) | 一种尾矿库干滩长度在线监测的方法 | |
CN106133732A (zh) | 一种岩石含水饱和度计算方法 | |
CN108759763A (zh) | 一种长距离轨道平行度测量方法及系统 | |
CN113420441B (zh) | 考虑裂缝倾角的缝洞型储层新三孔隙度模型及构建方法 | |
CN111336985A (zh) | 一种采空区铁塔的监测方法、装置、存储介质及电子设备 | |
CN108873073B (zh) | 一种基于网络并行电法的三维跨孔电阻率层析成像方法 | |
CN114047554B (zh) | 大地电阻率模型建模方法、装置、计算机设备和存储介质 | |
CN110471122A (zh) | 一种基于高密度电法温纳装置的三维电阻率测深应用方法 | |
CN105092107A (zh) | 一种高地应力钻井区存在岩心饼化现象时测量地应力方法 | |
US20140266214A1 (en) | Method and system for monitoring subsurface injection processes using a borehole electromagnetic source | |
CN109025981A (zh) | 曲流河侧积砂体对比方法及系统 | |
CN110968963B (zh) | 一种基于电磁理论计算m值的改进三孔隙度模型及方法 |
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 |