CN111257196B - 一种基于地层因素的岩石热物理参数的预测方法 - Google Patents

一种基于地层因素的岩石热物理参数的预测方法 Download PDF

Info

Publication number
CN111257196B
CN111257196B CN202010114305.0A CN202010114305A CN111257196B CN 111257196 B CN111257196 B CN 111257196B CN 202010114305 A CN202010114305 A CN 202010114305A CN 111257196 B CN111257196 B CN 111257196B
Authority
CN
China
Prior art keywords
pressure
temperature
thermal
wave velocity
normal
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
Application number
CN202010114305.0A
Other languages
English (en)
Other versions
CN111257196A (zh
Inventor
熊健
林海宇
万有维
裴浩辰
丁怀硕
戎成干
廖文海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202010114305.0A priority Critical patent/CN111257196B/zh
Publication of CN111257196A publication Critical patent/CN111257196A/zh
Application granted granted Critical
Publication of CN111257196B publication Critical patent/CN111257196B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/28Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/44Sample treatment involving radiation, e.g. heat
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00Investigating or analyzing materials by the use of thermal means
    • G01N25/20Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/041Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/07Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Acoustics & Sound (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Dispersion Chemistry (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明公开了一种基于地层因素的岩石热物理参数的预测方法,包括以下步骤:采集研究区域不同岩性的岩石并制备试样;测量试样的密度、孔隙度、纵波速度、横波速度、电阻率;在常温常压、温度T压力P条件下,分别测试热导率,并计算热扩散系数和比热容;将上述常温常压下获得的密度、孔隙度、纵波速度、横波速度、电阻率、热导率、热扩散系数和比热容输入BP神经网络模型中进行训练;利用常温常压对温度T压力P条件下岩石热物理参数进行校正得到地层热物理参数剖面。采用该方法预测的热物理参数其应用范围广。

Description

一种基于地层因素的岩石热物理参数的预测方法
技术领域
本发明涉及岩石热物理参数的预测方法领域,更具体的说是涉及一种基于地层因素的岩石热物理参数的预测方法。
背景技术
岩石的热物理性质在地热、岩石圈热结构、岩土工程等领域有着重要的研究与应用,在油气领域也备受关注。通常是通过钻取深部岩心室内测试来获取热物理参数,但是这种传统的方式会受到取样数量、取样条件、时间等因素的限制,不能较为准确地代表地下岩石真实的热物理性质。Evans et al.(1977)建立了沉积岩中热导率与孔隙度、波速和密度的关系,不过这种关系只适用于特定的地区和岩性;Merkl et al.(1976)对石灰岩地层的测井资料分析后获取了岩石主要的矿物成分及含量,通过每种矿物组成的标准热导率,计算得到了岩石骨架的热导率,不过这种方法需要有详细的测井资料来得到岩石矿物组分;Vacquier et al.(1988)、Hartmann et al.(2005)利用密度、声波等测井曲线预测了地层的热导率;欧新功等(2006)利用纵波速度预测了地层的热导率;Esteban et al.(2015)、马峰等(2019)基于室内试验,用纵波速度和孔隙度预测了岩心热导率,这些预测方法都忽略了温度和压力对对岩石热物理参数的影响,造成预测结果的局限性。
发明内容
本发明为了解决上述技术问题,研发一种基于地层因素的岩石热物理参数的预测方法,其对不同岩性岩样在不同温度压力条件下测试的热物理参数,基于密度、孔隙度、纵波速度、横波速度、电阻率,利用BP神经网络模型实现热物理参数的预测,其应用范围广。
本发明通过下述技术方案实现:
一种基于地层因素的岩石热物理参数的预测方法,包括以下步骤:
A、采集研究区域不同岩性的岩石并制备试样;
B、测量试样的密度、孔隙度、纵波速度、横波速度、电阻率;
C、在常温常压、温度T压力P条件下,分别测试热导率,并计算热扩散系数和比热容;
D、将上述常温常压下获得的密度、孔隙度、纵波速度、横波速度、电阻率、热导率、热扩散系数和比热容输入BP神经网络模型中进行训练;所述BP神经网络模型以密度、孔隙度、纵波速度、横波速度、电阻率为输入参数,并以热导率、热扩散系数和比热容为输出参数;
E、利用常温常压对温度T压力P条件下岩石热物理参数进行校正得到地层热物理参数剖面。
本发明与现有技术相比,具有如下的优点和有益效果:
1、本方法对不同岩性岩样在不同温度压力条件下测试的热物理参数,测量密度、孔隙度、纵波速度、横波速度、电阻率,利用常温常压下密度、孔隙度、纵波速度、横波速度、电阻率、热物理参数训练BP神经网络模型,利用BP神经网络模型实现热物理参数的预测,并以常温常压数据校正高温高压的热物理参数,得到地层热物理参数剖面,其应用范围广。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。
图1为本发明的原理图
图2为本发明建立的BP神经网络模型结构。
图3为本发明BP神经网络模型训练界面。
图4为本发明的BP神经网络模型的训练误差图。
图5为本发明热导率计算值和实测值对比图。
图6为本发明比热容计算值和实测值对比图。
图7为本发明热扩散系数计算值和实测值对比图。
图8为实施例2中计算得到地层热物理参数剖面图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
如图1所示的一种基于地层因素的岩石热物理参数的预测方法,包括以下步骤:
A、采集研究区域不同岩性的岩石并制备试样;
B、测量试样的密度、孔隙度、纵波速度、横波速度、电阻率;
C、在常温常压、温度T压力P条件下,分别测试热导率,并计算热扩散系数和比热容;
D、将上述常温常压下获得的密度、孔隙度、纵波速度、横波速度、电阻率、热导率、热扩散系数和比热容输入BP神经网络模型中进行训练;所述BP神经网络模型以密度、孔隙度、纵波速度、横波速度、电阻率为输入参数,并以热导率、热扩散系数和比热容为输出参数;
E、利用常温常压对温度T压力P条件下岩石热物理参数进行校正得到地层热物理参数剖面。
实施例2
基于上述实施例的原理,本实施例公开一具体实施方法并结合实验数据对本方案的创造性进行论证。
一种基于地层因素的岩石热物理参数的预测方法,包括以下步骤:
A、收集研究区块地质、测井资料和不同岩性的岩样,制备试样。根据SY/T 5336—1996《岩心常规分析方法》以及GB/T50266—99《工程岩体试验方法标准》规定的标准,将岩心样品制成30-60块直径为2.5cm、高为5.0cm的圆柱体岩样,岩样高度、直径误差均不超过0.3mm;两端面的不平行度,最大不超过0.05mm;端面垂直于试件轴,最大偏差不超过0.25°。
B、烘干试样,测量试样的几何尺寸和重量,获得试样的密度;使用气测孔渗仪SY/T6385-2016《覆压下岩石孔隙度和渗透率测定方法》获取试样的孔隙度;采用超声波透射法,按照SY/T 6351—1998《岩石声波特性的实验室测定》的要求对岩样进行室内声波测试,得到试样的纵波波速和横波速度;按照SYT 5385-2007《岩石电阻率参数实验室测量及计算方法》的要求,测试试样的电阻率。
C、在条件1-5的室内,进行岩石热学参数测试,获取热导率参数λ,并通过式1计算得到扩散系数α,利用公式2可以计算得到比热容Cp。其中条件1-5的温度分别为20℃、50℃、100℃、150℃和200℃,压力分别为0.1MPa、20MPa、50MPa、100MPa,其中,20℃、0.1MPa构成常温常压环境,对应的,50℃、100℃、150℃、200℃及其对应压力构成高温高压环境。
Figure GDA0002438895810000041
λ=αCpρ (2)
式中:λ为热导率,单位为W/(m·K);α为热扩散系数,单位为mm2/s;Cp为比热容,单位为kJ/kg·K;ρ为岩样密度,单位为g/cm3
D、将上述常温常压下获得的密度、孔隙度、纵波速度、横波速度、电阻率、热导率、热扩散系数和比热容输入BP神经网络模型中进行训练。
其中,BP神经网络模型以密度、孔隙度、纵波速度、横波速度、电阻率为输入参数,并以热导率、热扩散系数和比热容为输出参数,即该模型输入参数为5维输入,输出参数为3维输出,据BP神经网络理论可知,三层BP神经网络就可以完成任意的n维到m维的映射,考虑本发明将热导率、热扩散系数、比热容同时作为输出,将选择双隐含层结构,其中,根据输入与输出节点数确定隐层1的神经元个数为11个,隐层2的神经元个数为3个。本发明所建立的BP神经网络模型结构如图2所示。按照图2所示模型结构模型图,构建输入层为5维、输出层为3维、隐含层1神经元数目为11个,隐含层2神经元数目为3个的四层BP神经网络。隐含层和输出层的传递函数分别为tansig函数、logsig函数,网络训练函数为trainlm函数,网络权值学习函数为learngdm函数,性能函数为mse函数。
训练过程中,最大迭代次数为1000,网络训练误差为0.001。将常温常压环境实验数据分为两部分,一部分用于训练BP神经网络模型,另外一部分用于验证模型的精度。将用于训练的实验数据对BP网络模型进行训练。发明人按此方法搭建的BP神经网络模型训练界面、训练误差及训练状态分别如图3、图4所示。
利用构建的基于BP神经网络的岩石热物理参数模型计算岩石的热物理参数。岩石热物理参数计算值和实测值对比如图5~8所示。从对比图可以看出,除个别点外,由模型计算出来的热物理参数与实测热物理参数较为吻合,热导率计算值和实测值平均相对误差为7.300%,比热容计算值和实测值平均相对误差为3.788%,热扩散系数计算值和实测值平均相对误差为11.797%。
E、由于模型是在常温常压的实验条件下求取的,而压力和温度对热物理参数有着一定的影响,所以需要以常温常压条件下计算的到的热物理参数校正到相应的环境条件下的真实值,即利用常温常压对温度T压力P条件下岩石热物理参数进行校正得到地层热物理参数剖面。
其中,热导率的校正公式3如:
λ(p,T)=g1(x)λ(p0,T0) (3)
式中,λ(p,T)为在压力为pMPa,温度为T℃下的热导率;λ(p0,T0)为常温常压下的热导率;g1(x)为实验中不同压力和温度下测量的热导率值与常温常压下测量的热导率比值,由实验得到的g1(x),其中,
g1(x)=4.2339×10-4p-4.9791×10-4(T-273)+1.015。
扩散系数的校正公式4如:
α(p,T)=g2(x)α(p0,T0) (4)
式中,α(p,T)为在压力为pMPa,温度为T℃下的热扩散系数;α(p0,T0)为常温常压下的热扩散系数;g2(x)为实验中不同压力和温度下测量的热扩散系数与常温常压下测量的热扩散系数的比值,由实验得到的g2(x),其中,
g2(x)=4.0243×10-4p-4.8159×10-4(T-273)+1.012。
比热容的校正公式5如:
Cp(p,T)=g3(x)Cp(p0,T0) (5)
式中,Cp(p,T)为在压力为pMPa,温度为T℃下的比热容值;Cp(p0,T0)为常温常压下的比热容值;g3(x)为实验中不同压力和温度下测量的比热容与常温常压下测量的比热容比值,由实验得到的g3(x),其中,
g3(x)=3.9675×10-4p-4.6658×10-4(T-273)+1.014。
计算得到地层热物理参数剖面,如图8所示。可以看出,在碳酸盐岩地层段预测热导率曲线明显增大,在孔隙度较大的地层热导率预测值明显减小,这些与常规认识结论较为符合。另外,可以看出热导率预测曲线与扩散系数曲线有着较为相似的变化趋势。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种考虑地层因素的岩石热物理参数的预测方法,其特征在于,包括以下步骤:
A、采集研究区域不同岩性的岩石并制备试样;
B、测量试样的密度、孔隙度、纵波速度、横波速度、电阻率;
C、在常温常压、温度T压力P条件下,分别测试热导率,并计算热扩散系数和比热容;
D、将上述常温常压下获得的密度、孔隙度、纵波速度、横波速度、电阻率、热导率、热扩散系数和比热容输入BP神经网络模型中进行训练;所述BP神经网络模型以密度、孔隙度、纵波速度、横波速度、电阻率为输入参数,并以热导率、热扩散系数和比热容为输出参数;
E、利用常温常压对温度T压力P条件下岩石热物理参数进行校正得到地层热物理参数剖面;
对温度T压力P条件下岩石热物理参数进行校正的具体方法为:
λ(p,T)=g1(x)λ(p0,T0),
α(p,T)=g2(x)α(p0,T0),
Cp(p,T)=g3(x)Cp(p0,T0),
式中,λ(p,T)为在压力为pMPa,温度为T℃下的热导率;λ(p0,T0)为常温常压下的热导率;g1(x)为实验中不同压力和温度下测量的热导率值与常温常压下测量的热导率比值;α(p,T)为在压力为pMPa,温度为T℃下的热扩散系数;α(p0,T0)为常温常压下的热扩散系数;g2(x)为实验中不同压力和温度下测量的热扩散系数与常温常压下测量的热扩散系数的比值;Cp(p,T)为在压力为pMPa,温度为T℃下的比热容值;Cp(p0,T0)为常温常压下的比热容值;g3(x)为实验中不同压力和温度下测量的比热容与常温常压下测量的比热容比值;
其中,g1(x)=4.2339×10-4p-4.9791×10-4(T-273)+1.015,g2(x)=4.0243×10-4p-4.8159×10-4(T-273)+1.012,g3(x)=3.9675×10-4p-4.6658×10-4(T-273)+1.014。
2.根据权利要求1所述的一种考虑地层因素的岩石热物理参数的预测方法,其特征在于,所述温度T压力P均高于常温常压。
3.根据权利要求1所述的一种考虑地层因素的岩石热物理参数的预测方法,其特征在于,所述BP神经网络模型为双隐含层结构。
4.根据权利要求3所述的一种考虑地层因素的岩石热物理参数的预测方法,其特征在于,所述BP神经网络模型包括依次的输入层、隐含层1、隐含层2和输出层,所述隐含层1的神经元数目为11个,所述隐含层2的神经元数目为3个。
CN202010114305.0A 2020-02-24 2020-02-24 一种基于地层因素的岩石热物理参数的预测方法 Active CN111257196B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010114305.0A CN111257196B (zh) 2020-02-24 2020-02-24 一种基于地层因素的岩石热物理参数的预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010114305.0A CN111257196B (zh) 2020-02-24 2020-02-24 一种基于地层因素的岩石热物理参数的预测方法

Publications (2)

Publication Number Publication Date
CN111257196A CN111257196A (zh) 2020-06-09
CN111257196B true CN111257196B (zh) 2020-12-15

Family

ID=70951249

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010114305.0A Active CN111257196B (zh) 2020-02-24 2020-02-24 一种基于地层因素的岩石热物理参数的预测方法

Country Status (1)

Country Link
CN (1) CN111257196B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111948246B (zh) * 2020-08-25 2023-09-29 中国矿业大学 一种利用矿物组分计算砂岩热导率的方法
CN111948247B (zh) * 2020-08-25 2023-03-14 中国矿业大学 一种利用矿物含量分计算泥岩热导率的方法
CN114943186A (zh) * 2022-07-19 2022-08-26 数皮科技(湖北)有限公司 基于全岩化学数据的花岗岩热导率极限提升梯度预测方法
CN117740652B (zh) * 2024-02-19 2024-05-10 中国地质大学(武汉) 一种植生型多孔混凝土透砂系数快速测定方法及系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884738B (zh) * 2014-04-04 2016-09-07 中国科学技术大学 地热单井地层热物性分布评估方法
RU2636821C1 (ru) * 2016-05-27 2017-11-28 Шлюмберже Текнолоджи Б.В. Способ определения механических свойств породы пласта-коллектора
CN206114568U (zh) * 2016-09-21 2017-04-19 中国地质大学(武汉) 一种高温高压下岩石热物性参数测试系统
CN108562610B (zh) * 2018-03-13 2021-11-02 中国石油天然气股份有限公司 一种确定岩石热导率的方法及系统

Also Published As

Publication number Publication date
CN111257196A (zh) 2020-06-09

Similar Documents

Publication Publication Date Title
CN111257196B (zh) 一种基于地层因素的岩石热物理参数的预测方法
CN106761677B (zh) 页岩气水平井单井产能的测井预测方法
CN103323352B (zh) 天然气水合物沉积物动三轴力学-声学-电学同步测试的实验装置及方法
CN110348135B (zh) 一种随钻声波测井评价地层渗透率的方法
CN112255688B (zh) 一种基于岩石物理理论的三维地震反演地层压力的方法
CN110824556B (zh) 一种非常规致密砂岩储层的岩石物理模型建立方法及应用
CN109884109B (zh) 一种页岩孔隙水吸附量、游离量及其分布的评价方法
CN103132992A (zh) 一种岩石可钻性各向异性的评价方法及系统
CN113640119B (zh) 一种确定应力相关岩石动态Biot系数的方法
CN109583113B (zh) 一种岩石地层压实系数和有效孔隙体积压缩系数计算方法
Qu et al. Controls on matrix permeability of shale samples from Longmaxi and Niutitang formations, China
CN110596757A (zh) 一种页岩地层纵波及横波速度的校正方法
CN113392567B (zh) 一种双重孔隙煤体的煤层气预测方法和系统
CN109613624A (zh) 一种储层岩石声电性质联合模拟方法
CN108412488B (zh) 快速确定页岩气储层有机孔隙度的测井方法
CN111077174A (zh) 一种页岩储层游离气和吸附气含量计算方法
Wang et al. A modified pulse‐decay approach to simultaneously measure permeability and porosity of tight rocks
CN112505766B (zh) 一种评价井外不同方位下裂缝发育程度的方法
CN111089904B (zh) 一种考虑地层特性的室内测量波速围压、频散校正方法
CN110909311A (zh) 一种计算薄煤层含气量的方法
CN116658157B (zh) 致密砂岩气藏地层压力预测方法及系统
CN116025345B (zh) 一种横观各向同性地层静态有效应力系数的确定方法
CN114895365B (zh) 高温高压岩心电阻率标定的深部地温场预测方法及装置
CN114184764B (zh) 致密碳酸盐岩储层岩石力学层划分方法及系统
CN113935215B (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