CN112268923B - 一种基于测井曲线获取地层热导率的方法 - Google Patents
一种基于测井曲线获取地层热导率的方法 Download PDFInfo
- Publication number
- CN112268923B CN112268923B CN202011141892.9A CN202011141892A CN112268923B CN 112268923 B CN112268923 B CN 112268923B CN 202011141892 A CN202011141892 A CN 202011141892A CN 112268923 B CN112268923 B CN 112268923B
- Authority
- CN
- China
- Prior art keywords
- thermal conductivity
- formation
- content
- mudstone
- logging curve
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/18—Investigating or analyzing materials by the use of thermal means by investigating thermal conductivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/10—Geothermal energy
Landscapes
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Geophysics (AREA)
- Health & Medical Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种基于测井曲线获取地层热导率的方法,属于地层热导率计算领域。一种基于测井曲线获取地层热导率的方法,包括:1)获取自然伽马测井曲线、声波时差测井曲线和含水饱和度测井曲线;2)基于自然伽马测井曲线上的自然伽马值进行计算,得到泥岩含量和砂岩含量;其中,砂岩含量为1‑泥岩含量;基于泥岩含量和砂岩含量进行计算,得到地层固体骨架热导率;基于声波时差测井曲线上的值进行计算,得到地层孔隙度;3)基于地层固体骨架热导率、地层孔隙度和含水饱和度,利用通用的多孔介质热导率计算公式,计算得到地层热导率,进而得到地层热导率‑深度分布曲线。
Description
技术领域
本发明属于地层热导率计算领域,尤其是一种基于测井曲线获取地层热导率的方法。
背景技术
地层热导率反映地层的热传导能力,单位为W/(m.K),其物理意义为单位温度梯度在单位时间内经单位地层导热面积所传递的热量。地层热导率是地源热泵、地热流体开采、油气田开发等地下工程的重要物性参数,决定着开采介质与地层之间的传热性能,与地下工程的设计、施工、运行等各个流程息息相关。然而,地层热导率与孔隙率、含水率、岩石组分等多个因素有关,且随深度波动明显,导致其准确获取较困难。
现有的地层热导率获取主要有两种方法:一种是基于钻孔取样实际测量,推算得全深度的热导率(如CN107727686B、CN105303043B);另一种为通过测量埋管沿线的热响应,反向推算得出地层热导率(如CN106770439A、CN103884738A)。基于钻孔取样的方法,需要在实验室进行激光闪光法、移动点热源TCS法等测量得到热导率,但采集出试样的围压、含水率等已与其地下情况产生较大变化,导致测量结果不准确,且采样只能得到单一深度下的参数。基于埋管沿线热响应进行反算只能得到浅层地层的热导率,但实际地质情况比较复杂,10m之内就会出现地层性质的巨大差异,并且需要传感器测量埋管内各层稳态的流体温度,在中深层埋管中较难实施,很难获取中深层地层热导率。为了对实际地下工程的设备结构参数设计和介质选择提出指导意义,准确获取全深度的地层热导率是亟待发展的技术。
发明内容
本发明的目的在于克服现有技术无法获得中深层地层(1000m以下)热导率在深度上的分布缺点,提供一种基于测井曲线获取地层热导率的方法。
为达到上述目的,本发明采用以下技术方案予以实现:
一种基于测井曲线获取地层热导率的方法,包括以下步骤:
1)获取自然伽马测井曲线、声波时差测井曲线和含水饱和度测井曲线;
2)基于自然伽马测井曲线上的自然伽马值进行计算,得到泥岩含量和砂岩含量,进而得到泥岩含量-深度分布曲线和砂岩含量-深度分布曲线;其中,砂岩含量为1-泥岩含量;
基于泥岩含量和砂岩含量进行计算,得到地层固体骨架热导率,进而得到地层固体骨架热导率-深度分布曲线;
基于声波时差测井曲线上的值进行计算,得到地层孔隙度,进而得到地层孔隙度-深度分布曲线;
3)基于地层固体骨架热导率、地层孔隙度和含水饱和度,利用通用多孔介质热导率计算公式,计算得到地层热导率,进而得到地层热导率-深度分布曲线。
进一步的,步骤2)中,基于自然伽马测井曲线上的值进行计算,得到泥岩含量计算过程为:
式(1)中,GR为自然伽马值,GRmin为自然伽马值的最小值,GRmax为自然伽马值的最大值,IGR为归一化后的伽马值;
式(2)中,Vsh为泥岩含量,GCUR为与地层年龄有关的值。
进一步的,步骤2)中基于泥岩含量和砂岩含量进行计算,得到的地层固体骨架热导率λm的计算过程为:
λm=Vshλsh+Vsaλsa (3)
式中,Vsh为泥岩含量,Vsa为砂岩含量,λsh、为泥岩的导热系数,λsa为砂岩的导热系数。
进一步的,步骤2)中基于声波时差测井曲线上的值进行计算,得到地层孔隙度的计算过程为:
式中,φ为地层孔隙度,Δtp为由声波时差曲线获取的地层声波时差,Δtma、为固体骨架时差,Δtsh为泥岩时差,Δtf为流体时差。
进一步的,步骤3)的具体过程为:
λ0=λm(1-φ)+λwSwφ (5)
式中,λm为地层固体骨架热导率,φ为地层孔隙度,λw、为水的导热系数,Sw为含水饱和度。
与现有技术相比,本发明具有以下有益效果:
本发明的基于测井曲线获取地层热导率的方法,基于测井曲线计算出地层热导率在深度上的分布值,无需额外的实验测试,方便实施;测井曲线是利用声学特性、放射性等地球物理特性,测量地层物理参数的方法,本发明采用了自然伽马测井曲线、声波时差测井曲线、含水饱和度测井曲线分别获取泥岩含量、孔隙度、含水饱和度的全深度分布,再根据参数的换算关系得到地层固体骨架热导率和地层热导率;本发明能够获得中深层地层热导率的全深度分布,改善现有基于取样测试获得特定深度下热导率现状,给地下工程的传热性能准确分析提供重要参数,从而对设备结构参数设计和介质选择具有重要意义。
附图说明
图1为实施例的自然伽马、声波时差、含水饱和度测井曲线;
图2为实施例的IGR和泥岩含量随深度的分布图;
图3为实施例的地层导热系数随深度的分布图;
图4为采用本发明获取的热导率计算地埋管出口温度与实测值的对比图;
图5为采用本发明获取的热导率计算1500m深度处管外壁温度与实测值的对比图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
一种基于测井曲线获取地层热导率的方法,通过钻井完井后的测井参数计算出所有深度的地层热导率,无需额外的实验测试,方便实施。测井曲线是利用声学特性、放射性等地球物理特性,测量地层物理参数的方法,本发明采用了自然伽马测井曲线、声波时差测井曲线、含水饱和度测井曲线分别获取泥岩含量、孔隙度、含水饱和度的全深度分布,再根据参数的换算关系分别得到地层固体骨架热导率和地层热导率。
实施例
获取试验井的测井曲线,从中摘取出自然伽马测井曲线、声波时差测井曲线、含水饱和度测井曲线,得到图1,自然伽马测井曲线是测量地层中自然存在的放射性核素衰变放射出的γ射线强度随地层深度的变化曲线,可以用于划分岩性,计算泥岩含量,计算公式如下:
式(1)中,GR代表自然伽马值,IGR为归一化后的值;
式(2)中,Vsh为泥岩含量,GCUR值与地层的年龄有关,本实施例中的试验井所在地层为新生界第四系,GCUR值取3.7,计算得到的泥岩含量如图2所示,图2有的两条曲线分别为IGR-深度曲线、Vsh-深度曲线。
根据泥岩含量Vsh,利用泥岩和砂岩的导热系数求解地层固体骨架热导率λm:
λm=Vshλsh+Vsaλsa (3)
其中,λsh、为泥岩的导热系数,λsa为砂岩的导热系数。
根据附近地区的岩石取样数据,泥岩和砂岩的导热系数分别取为1.62W/(m.K)和2.52W/(m.K)。
根据声波时差曲线计算地层孔隙度,计算式如下:
式中,φ为地层孔隙度,Δtp为由声波时差曲线读出的地层声波时差,Δtma、Δtsh、Δtf分别为固体骨架时差、泥岩时差和流体时差,根据附近地区的岩石取样数据,Δtf=189μs/ft,,Δtma=49μs/ft,,Δtsh=100μs/ft。
参见图1,图1的中含水饱和度曲线,试验井地下800-1760m存在渗流,该段的含水饱和度为1,根据通用多孔介质热导率计算公式,计算得到地层热导率λ0:
λ0=λm(1-φ)+λwSwφ (5)
式中,λw为水的导热系数,本实施例中取值为0.6W/(m.K),Sw为含水饱和度。
本实施例的地层导热系数分布如图3所示,图3中600米以内的缺少测井曲线,取的是地层导热系数的平均值。
采用上述方法得到的地层热导率,可以计算得出试验井处地埋管的出口温度和某一深度处管外壁温度随时间的变化,其计算值与实验值的对比分别如下图4、图5所示。由对比结果可知计算值与实测值误差较小,说明本发明获得地层热导率的方法准确性较高。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (3)
1.一种基于测井曲线获取地层热导率的方法,其特征在于,包括以下步骤:
1)获取自然伽马测井曲线、声波时差测井曲线和含水饱和度测井曲线;
2)基于自然伽马测井曲线上的自然伽马值进行计算,得到泥岩含量和砂岩含量,进而得到泥岩含量-深度分布曲线和砂岩含量-深度分布曲线;其中,砂岩含量为1-泥岩含量;
基于泥岩含量和砂岩含量进行计算,得到地层固体骨架热导率,进而得到地层固体骨架热导率-深度分布曲线;
基于声波时差测井曲线上的值进行计算,得到地层孔隙度,进而得到地层孔隙度-深度分布曲线;
步骤2)中,基于自然伽马测井曲线上的值进行计算,得到泥岩含量计算过程为:
式(1)中,GR为自然伽马值,GRmin为自然伽马值的最小值,GRmax为自然伽马值的最大值,IGR为归一化后的伽马值;
式(2)中,Vsh为泥岩含量,GCUR为与地层年龄有关的值;
3)基于地层固体骨架热导率、地层孔隙度和含水饱和度,利用通用多孔介质热导率计算公式,计算得到地层热导率,进而得到地层热导率-深度分布曲线;
步骤3)的具体过程为:
λ0=λm(1-φ)+λwSwφ (5)
式中,λm为地层固体骨架热导率,φ为地层孔隙度,λw、为水的导热系数,Sw为含水饱和度。
2.根据权利要求1所述的基于测井曲线获取地层热导率的方法,其特征在于,步骤2)中基于泥岩含量和砂岩含量进行计算,得到的地层固体骨架热导率λm的计算过程为:
λm=Vshλsh+Vsaλsa (3)
式中,Vsh为泥岩含量,Vsa为砂岩含量,λsh`为泥岩的导热系数,λsa为砂岩的导热系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011141892.9A CN112268923B (zh) | 2020-10-22 | 2020-10-22 | 一种基于测井曲线获取地层热导率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011141892.9A CN112268923B (zh) | 2020-10-22 | 2020-10-22 | 一种基于测井曲线获取地层热导率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112268923A CN112268923A (zh) | 2021-01-26 |
CN112268923B true CN112268923B (zh) | 2022-05-20 |
Family
ID=74342220
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011141892.9A Active CN112268923B (zh) | 2020-10-22 | 2020-10-22 | 一种基于测井曲线获取地层热导率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112268923B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114295674B (zh) * | 2022-01-07 | 2022-07-01 | 东北石油大学 | 一种确定沉积盆地中原位热导率的方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105093304A (zh) * | 2014-05-15 | 2015-11-25 | 中国石油化工股份有限公司 | 一种地球物理勘探利用测井曲线自动求取岩性曲线的方法 |
CN104632202B (zh) * | 2014-12-25 | 2017-06-13 | 中国石油天然气股份有限公司 | 确定干粘土三孔隙度测井参数值的方法及装置 |
CN106951660B (zh) * | 2017-04-05 | 2020-08-11 | 中国石油天然气股份有限公司 | 一种海相碎屑岩水平井储层测井解释方法及装置 |
RU2658856C1 (ru) * | 2017-07-14 | 2018-06-25 | Шлюмберже Текнолоджи Б.В. | Способ определения профиля теплопроводности горных пород в скважине |
CN107907910B (zh) * | 2017-11-16 | 2019-12-13 | 山东胜软科技股份有限公司 | 一种不同岩性油藏横波测井确定方法 |
CN109580453A (zh) * | 2018-12-24 | 2019-04-05 | 核工业北京地质研究院 | 基于最优化测井解释计算砂岩型铀矿砂体孔隙度的方法 |
-
2020
- 2020-10-22 CN CN202011141892.9A patent/CN112268923B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112268923A (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9249659B2 (en) | Formation fluid property determination | |
CN109753755B (zh) | 一种确定储层含水饱和度的方法 | |
US5561245A (en) | Method for determining flow regime in multiphase fluid flow in a wellbore | |
JPS60501875A (ja) | 好ましくは多重サ−ミスタプロ−ブを使用する、海底温度勾配による地表下炭化水素の探索 | |
MX2013012178A (es) | Metodos y sistemas para estimar la resistividad y porosidad de la formacion. | |
CN109138975B (zh) | 一种基于时移测井数据的求解相渗特征曲线的新方法 | |
CA2392618C (en) | Improved method and apparatus for predicting the fluid characteristics in a well hole | |
BRPI1001536A2 (pt) | mÉtodo para determinar permeabilidades eficazes de formaÇÕes terrestres | |
CN104514552A (zh) | 一种煤层气储层识别与丰度预测的方法 | |
CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
CN112282743A (zh) | 用于对钻井泥岩地层压力进行预测的方法 | |
CN112268923B (zh) | 一种基于测井曲线获取地层热导率的方法 | |
CN116430474A (zh) | 一种地热资源储量计算方法、***及电子设备 | |
CN108412488B (zh) | 快速确定页岩气储层有机孔隙度的测井方法 | |
US20210293990A1 (en) | Performing Dynamic Time Warping With Null or Missing Data | |
Li et al. | In situ estimation of relative permeability from resistivity measurements | |
BR112021013205A2 (pt) | Método, e, sistema | |
CN114086938B (zh) | 一种非均质性砂岩储层的含气饱和度预测方法 | |
RU2632800C2 (ru) | Метод определения коэффициента текущей нефтенасыщенности разрабатываемого нефтеносного пласта в скважине | |
Akram et al. | A model to predict wireline formation tester sample contamination | |
CN108647417B (zh) | 一种确定页岩气储层含气饱和度的简易方法 | |
CN113818868A (zh) | 一种储层原始温度的计算方法 | |
CN109209357B (zh) | 一种生产测井解释方法 | |
Ohaeri et al. | Evaluation of reservoir connectivity and hydrocarbon resource size in a deep water gas field using multi-well interference tests | |
CN114441402A (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 |