CN115045650A - 初始地温测量方法、装置、计算设备及存储介质 - Google Patents

初始地温测量方法、装置、计算设备及存储介质 Download PDF

Info

Publication number
CN115045650A
CN115045650A CN202210738242.5A CN202210738242A CN115045650A CN 115045650 A CN115045650 A CN 115045650A CN 202210738242 A CN202210738242 A CN 202210738242A CN 115045650 A CN115045650 A CN 115045650A
Authority
CN
China
Prior art keywords
target
temperature
initial
target position
locations
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.)
Pending
Application number
CN202210738242.5A
Other languages
English (en)
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.)
China Oilfield Services Ltd
China National Offshore Oil Corp CNOOC
Original Assignee
China Oilfield Services Ltd
China National Offshore Oil Corp CNOOC
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 China Oilfield Services Ltd, China National Offshore Oil Corp CNOOC filed Critical China Oilfield Services Ltd
Priority to CN202210738242.5A priority Critical patent/CN115045650A/zh
Publication of CN115045650A publication Critical patent/CN115045650A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/06Measuring temperature or pressure
    • E21B47/07Temperature
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B33/00Sealing or packing boreholes or wells
    • E21B33/10Sealing or packing boreholes or wells in the borehole
    • E21B33/13Methods or devices for cementing, for plugging holes, crevices or the like
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/10Geothermal energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Computational Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Computing Systems (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Geophysics (AREA)
  • Algebra (AREA)

Abstract

本发明涉及勘测领域技术领域,公开了一种初始地温测量方法、装置、计算设备及存储介质,方法包括:获取井眼的H个目标位置中每一目标位置的至少一个目标温度,H为大于1的正整数,至少一目标位置高于另一目标位置;根据目标温度计算每一目标位置的初始地温;根据H个目标位置对应的地层深度和H个初始地温,拟合初始地温函数。通过上述方式,本发明实施例由于至少一目标位置高于另一目标位置,则可以通过H个目标位置的地层深度参数以及H个初始地温进行数学建模,得到初始地温函数,当用户想要得知任意一个地层深度的初始地温时,则输入对应的地层深度参数,根据该初始地温函数即可获得对应地层深度的初始地温值,提高工作效率。

Description

初始地温测量方法、装置、计算设备及存储介质
技术领域
本发明属于勘测技术领域,具体涉及一种初始地温测量方法、装置、计算设备及存储介质。
背景技术
井筒循环温度是石油钻井和固井行业内非常关注的参数,对固井质量及成本影响很大,而地层的初始温度对于井筒循环温度有非常关键的影响,初始地温的确定往往决定了最终固井循环水泥浆体系的设计,对固井质量的影响很大。
目前的初始地温往往由地质勘探预估得到,其准确性无法验证。
因此,有必要设计一种初始地温测量方法、装置、计算设备及存储介质,以克服上述问题。
发明内容
鉴于上述问题,本发明实施例提供了一种初始地温测量方法、装置、计算设备及存储介质,用于解决现有技术中存在的问题。
根据本发明实施例的第一方面,提供了一种初始地温测量方法,所述方法包括:
获取井眼的H个目标位置中每一所述目标位置的至少一个目标温度,H为大于1的正整数,至少一所述目标位置高于另一所述目标位置;
根据所述目标温度计算每一所述目标位置的初始地温;
根据H个所述目标位置对应的地层深度和H个所述初始地温,拟合初始地温函数。
在一些实施例中,H个所述目标位置包括位于所述井眼内的至少两个检测位置,其中一所述检测位置高于另一所述检测位置,所述获取井眼的H个目标位置中每一所述目标位置的至少一个目标温度,进一步包括:
获取每一所述检测位置在流体静置后不同时刻的至少两个所述目标温度。
在一些实施例中,所述目标位置沿井身深度间隔分布有L个,所述获取井眼的H个目标位置中每一所述目标位置的至少一个目标温度,进一步包括:
获取每一所述目标位置在流体静置后N个时刻的N个目标温度,N为大于1的正整数,L大于N。
在一些实施例中,至少一个所述检测位置位于所述井眼的井底,至少一个所述检测位置高于所述井底。
在一些实施例中,H个所述目标位置包括位于所述井眼内的至少一个检测位置和位于井眼表层的地表位置。
在一些实施例中,H个所述目标位置包括位于所述井眼内的M个检测位置和位于井眼表层的地表位置,
第j个所述检测位置的所述初始地温
Figure BDA0003716678670000021
Figure BDA0003716678670000022
其中,Ti,j为第j个检测位置的第i个目标温度,ΔTj为第j个所述目标位置根据最小二乘法拟合得到误差最小条件下的温差,
Figure BDA0003716678670000023
ψ(ω)=J1(ω)J0(A1/2ω)-βJ0(ω)J1(A1/2ω),φ(ω)=J1(ω)Y0(A1/2ω)-βJ0(ω)Y1(A1/2ω),
Figure BDA0003716678670000024
其中,a1为井眼内流体的热扩散率,a2为地层的热扩散率,λ1为井眼内流体的导热系数,λ2为岩层的导热系数,t为所述目标温度对应的流体静置时间,R1为井眼半径,ω为虚拟复变量,J0为零阶第一类贝塞尔函数,J1为一阶第一类贝塞尔函数,Y0为零阶第二类贝塞尔函数,Y1为一阶第二类贝塞尔函数,M和N均为大于0的正整数。
在一些实施例中,所述初始地温函数
Figure BDA0003716678670000031
其中,Tf0为所述地表位置的所述初始温度,Tf,j为第j个检测位置的初始地温,Zj为第j个所述检测位置对应的地层深度,Z为地层深度的变量。
根据本发明实施例的第二方面,提供了一种初始地温测量装置,所述测量装置包括:
获取模块,用于获取井眼的H个目标位置中每一所述目标位置的至少一一个目标温度,H为大于1的正整数,至少一所述目标位置高于另一所述目标位置;
计算模块,用于根据所述目标温度计算每一所述目标位置的初始地温;
执行模块,用于根据H个所述目标位置对应的地层深度和H个所述初始地温,拟合初始地温函数。
根据本发明实施例的第三方面,提供了一种计算设备,包括:处理器、存储器、通信接口和通信总线,所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信;
所述存储器用于存放至少一可执行指令,所述可执行指令使所述处理器执行如上述任一项所述的初始地温测量方法的操作。
根据本发明实施例的第四方面,提供了一种计算机可读存储介质,所述存储介质中存储有至少一可执行指令,所述可执行指令在运行时执行如上述任一项所述的初始地温测量方法的操作。
本发明实施例由于至少一所述目标位置高于另一所述目标位置,则可以通过H个目标位置的地层深度参数以及H个初始地温进行数学建模,得到初始地温函数,此时,只需要输入对应的地层深度参数,即可得到对应地层深度的初始地温,当用户想要得知任意一个地层深度的初始地温时,则输入对应的地层深度参数,根据该初始地温函数即可获得对应地层深度的初始地温值,用户无需再去测量对应地层深度的温度以计算对应的初始地温,减少工作量,提高工作效率。此外,由于H为大于1的正整数,能够减小目标位置对应目标温度的误差以及初始地温的计算误差,使得初始地温的计算结果较为准确,当H值越大和/或每一目标位置的目标温度数量越多时,越能消除目标位置对应目标温度的误差以及初始地温的计算误差,使得初始地温的计算结果越准确,初始地温函数的拟合结果越好。同时,该方法能够通过某个目标位置短时间内的测量目标温度计算该目标位置的初始地层温度,在保证得到较为精确的初始地温和地温梯度的同时大大减少了等待时间,提高作业效率
上述说明仅是本发明实施例技术方案的概述,为了能够更清楚了解本发明实施例的技术手段,而可依照说明书的内容予以实施,并且为了让本发明实施例的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明的具体实施方式。
附图说明
附图仅用于示出实施方式,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1示出了本发明实施例提供的初始地温测量方法的流程示意图;
图2示出了本发明实施例提供的第j个检测位置的示意图;
图3示出了本发明实施例提供的初始地温测量装置的结构示意图;
图4示出了本发明实施例提供的计算设备的结构示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。
针对现有技术中,初始地温往往由地质勘探预估得到,其准确性无法验证的问题,发明人发现,现有技术中,固井前通常会在静置时段对井内流体的温度进行测量,仅由井底某一处的电测温度演算初始地温,无法得到整个井眼深度对应地层的温度分布梯度,当需要得知其他地层的初始地温时,需要另外测量其他目标位置的电测温度以对应演算该目标位置的初始地温,影响工作效率。此外,在实际测量中,地层温度是一个缓慢恢复的过程,通常需要等待24小时以上才能恢复到初始温度,长时间等待严重影响钻井周期,增加钻井费用,特别是海上钻井。因此需要一种方法,能够根据短时间等待周期所测量的地层温度计算初始地层温度并能精确反应地层温度梯度。
发明人通过分析发现,通过结合井眼的多个不同深度的温度测量位置和对应的多个初始地温进行建模,得到拟合的初始地温函数,从而使得初始地温的预测值更加准确,此外,通过该初始地温函数,可以得出井内的任意深度的初始地温,无需再去测量对应深度的温度以计算对应的初始地温,从而提高工作效率。
本发明实施例中,初始地温测量方法由计算设备执行,例如手机、平板电脑、个人数字助理(PDA)、台式计算机、笔记本电脑等各种计算设备。
具体实施例如下:
图1示出了本发明实施例提供的初始地温测量方法的流程图,该方法包括以下步骤:
步骤110:获取井眼的H个目标位置中每一目标位置的至少一个目标温度,H为大于1的正整数,至少一目标位置高于另一目标位置。
步骤120:根据目标温度计算每一目标位置的初始地温。
步骤130:根据H个目标位置对应的地层深度和H个初始地温,拟合初始地温函数。
其中,步骤110和步骤120中,井眼具有井身深度,H个目标位置中至少两个目标位置沿井身深度间隔设置,即至少一目标位置高于另一目标位置,以便进行函数建模。其中,H为大于1的正整数,换句话说,即是H大于或等于2,目标位置至少设有两个。目标温度至少设有一个,以便通过目标温度计算初始地温,初始地温表示目标位置所在地层深度的地层初始温度,其中,每一目标位置温度测量次数越多,则每一目标位置的目标温度越多,越能消除根据目标温度得出初始地温的计算误差。
在一些实施例中,H个目标位置可以全部位于井眼内部,至少两个目标位置沿井身深度间隔设置,该目标位置设置于井眼内的流体中,即温度传感器进入到流体中进行温度测量,相应的,温度传感器所在位置即为目标位置,该目标温度为流体静置后检测得到的,流体静置是指从井眼内流体停止循环,以防止流体流动影响检测误差,该初始地温无法直接从所测量的目标温度得出,需要根据目标温度间接计算得出,该初始地温为预估值。其中,初始地温可以通过半无限大复合材料圆筒的非稳态导热方程计算得到,也可以通过有限或无限平板、圆球非稳态方程等方式计算得到,在此不做限定,根据需要设置。
其中,可以在流体静置后的不同时刻测量得到每一个位于井眼内的目标位置的至少两个目标温度,以提高初始地温的计算精度,该不同时刻指不同的流体静置时间,例如在流体静置时间为30分钟时测得第一个目标温度,则在流体静置后的第一时刻为30分钟,在流体静置时间为60分钟时测得第二个目标温度,则在流体静置后的第二时刻为60分钟。其中,目标温度可以是流体静置后的0时刻时的流体温度,也可以是流体静置后间隔一段时间的时刻的流体温度,在此不做限定,根据需要设置。
在一些实施例中,也可以是其中一目标位置设于井眼表层,即地表位置,其他目标位置设于井眼内部,该地表位置高于其他位于井眼内部的目标位置。其中,该地表位置为流体的注入口,因此地表位置所得到的目标温度是准确的,地表位置的初始地温即是对应的目标温度,无需预估计算。在一些情况下,井眼设于海域,该井眼上方为海水,则地表位置的初始地温可以直接通过温度传感器在地表检测,所检测到的目标温度即为地表位置的初始地温;或者,在一些情况下,井眼设于陆地上,此时,井眼上方没有海水,则地表位置的初始地温可以由当地天气获取气温,所获取的气温即为地表位置对应的目标温度,或者也可以通过温度传感器在地表检测得到对应的目标温度,在此不做限定。其中,温度传感器可以在地表位置测量某一时刻的一个目标温度,也可以测量不同时刻的目标温度,以使得所得到的地表位置的初始地温更加准确。
在一些实施例中,温度传感器在井眼内同一地层深度的流体静置后不同时刻测量至少两个目标温度,换句话说,在第一时刻,温度传感器位于第一地层深度的第一目标位置检测得出第一目标温度,在第二时刻,温度传感器位于第一地层深度的第二目标位置检测得出第二目标温度,其中,该第一目标位置和第二目标位置在径向上间隔设置,第一时刻和第二时刻间隔数分钟或者数十分钟或者数小时或者其他时间间隔,在此不做限定,根据需要设置。同理,若想获得3个或者更多的目标温度,则对应在同一地层深度的径向分布不同的检测位置以及不同时刻测量得到,在此不做赘述。通过不同的检测位置可以在同一地层深度的径向上分布,其所得到的目标温度计算得出的初始地温会比较准确。
步骤130中,由于目标位置至少设有两个,且至少一目标位置高于另一目标位置,则可以通过H个目标位置的地层深度参数以及H个初始地温进行数学建模,得到初始地温函数,此时,只需要输入对应的地层深度参数,即可得到对应地层深度的初始地温,当用户想要得知任意一个井深的地层深度时,则输入对应的地层深度参数,根据该初始地温函数即可获得对应地层深度的初始地温值,用户无需再去测量对应地层深度的温度以计算对应的初始地温,从而提高工作效率。
在一些实施例中,目标位置包括位于井底的目标位置和位于井表的地表位置。例如,假设井深为2000m,其中井底的地层深度为2000,地表位置的地层深度为0m。位于井底的目标位置的温度可通过温度传感器检测得到,通过计算即可得到对应的初始地温,其中,由于井底有底部边界,在设置温度传感器时能够依据该边界定位,使得温度传感器能够较为准确的定位在该井底,同时井底的地层深度也是可以确定的,相应的,设于井底的温度传感器的位置误差较小,对应的目标温度也可以较为准确,根据该目标温度所得的初始地温的计算结果较为准确,在此基础上,为了使得井底对应地层深度的温度检测结果以及初始地温结果更准确,可以在井底在不同静置时刻下进行多次检测,从而得到井底目标位置的多个目标温度,而后进一步计算井底对应的地层深度的初始地温。地表位置的初始地温也可通过温度传感器检测得到,所检测到的目标温度即为地表位置的初始地温。相应的,则可建立关于初始地温与地层深度的初始地温函数,通常为线性函数。
在一些实施例中,目标位置包括位于井眼内沿井身深度间隔设置的两个目标位置。例如,假设井深为2000m,其中一个目标位置对应的地层深度为1000m,另一个目标位置对应的地层深度为500m,通过温度传感器对应检测1000m以及500m地层深度的温度,从而计算得到两个目标位置对应的初始地温,建立关于初始地温与地层深度的初始地温函数。
在一些实施例中,H大于或等于3,此时,目标位置至少设有3个,在拟合初始地温函数时,通常通过最小二乘法降低拟合误差,即使目标位置并不是设于地表位置和井底位置,也能减小温度传感器检测的温度误差以及对应的初始地温误差,从而使得初始地温函数拟合结果更好,以及使得初始地温的计算结果更准确,其中,初始地温函数可以为线性函数,也可以为非线性函数,如二次多项式拟合函数,其非线性函数为y=a+bx+cx2,或者非线性函数也可根据实际拟合结果形成指数函数,或者其他曲线形式的非线性函数,在此不做限定,根据需要设置。
当H越大时,和/或每一目标位置在不同时刻的测温次数越多时,在拟合初始地温函数时,目标位置的数量越多,以及每一目标位置的目标温度个数越多,相应的,H个目标位置的温度误差越小,从而使得初始地温的计算误差越小,使得初始地温的计算结果越准确,初始地温函数的拟合结果越好。
通过步骤110至步骤130,由于至少一目标位置高于另一目标位置,则可以通过H个目标位置的地层深度参数以及H各初始地温进行数学建模,得到初始地温函数,此时,只需要输入对应的地层深度参数,即可得到对应地层深度的初始地温,当用户想要得知任意一个地层深度的初始地温时,则输入对应的地层深度参数,根据该初始地温函数即可获得对应地层深度的初始地温值,用户无需再去测量对应深度的温度以计算对应的初始地温,减少工作量,提高工作效率。此外,由于H为大于1的正整数,能够减小目标位置对应目标温度的误差以及初始地温的计算误差,使得初始地温的计算结果较为准确,当H值越大和/或每一目标位置的目标温度数量越多时,越能消除目标位置对应目标温度的误差以及初始地温的计算误差,使得初始地温的计算结果越准确,初始地温函数的拟合结果越好。同时,该方法能够通过某个目标位置短时间内的测量目标温度计算该目标位置的初始地层温度,在保证得到较为精确的初始地温和地温梯度的同时大大减少了等待时间,提高作业效率。
在一些实施例中,H个目标位置包括位于井眼内的至少两个检测位置,其中一检测位置高于另一检测位置,步骤110进一步包括:
步骤a01:获取每一检测位置在流体静置后不同时刻的至少两个目标温度。
其中,步骤a01中,针对至少两个目标位置位于井眼内的情况,此时该位于井眼内的目标位置为检测位置,相应的,需要通过温度传感器于该检测位置测量目标温度。其中该至少两个检测位置可以均位于井底与井表之间;或者其中一检测位置位于井底,其他检测位置位于井底与井表之间;当检测位置设有至少三个时,三个检测位置可以均位于井底与井表之间,也可以是其中一检测位置位于井底与井表之间,至少两个检测位置位于井底,也可以是其中一检测位置位于井底,至少两个检测位置位于井底与井表之间。检测位置的设置不做限定,根据需要设置。
每一检测位置在流体静置后不同时刻通过温度传感器检测得到至少两个目标温度,根据该至少两个目标温度计算得出初始地温,以减小初始地温的计算误差,提高初始地温的计算准确度。其中,在一些实施例中,温度传感器为静止的,可以在同一检测位置测量流体静置后不同时刻对应的至少两个目标温度,换句话说,温度传感器静止在同一检测位置,而后测量在流体静置后不同时刻的至少两个目标温度。
在一些实施例中,目标位置沿井身深度间隔分布有L个,步骤110进一步还包括:
步骤a02:获取每一目标位置在流体静置后N个时刻的N个目标温度,N为大于1的正整数,L大于N。
其中,L大于N表示,表示沿井身深度间隔分布的目标位置的设置数量多于每一目标位置的目标温度的设置数量,在一方面,可以在初始地温计算结果较为准确的基础上控制每一目标位置对应的目标温度的测量时长,提高工作效率,另一方面,可以保证关于地层深度的样本量较多,提高初始地温函数的拟合效果。例如,每一目标位置检测得到不同时刻的两个目标温度,相应的沿井身深度,至少设有三个间隔设置的目标位置,从而使得沿井身深度目标位置较多,初始地温函数的拟合效果更好。
在一些实施例中,至少两个检测位置位于井眼的井底,至少一个检测位置高于井底。
其中,至少两个检测位置位于井眼内的检测位置设于井底,可以进一步减小井底地层深度的初始地温的计算误差,使得初始地温函数拟合效果更好。至少一个检测位置高于井底,以根据不同地层深度的检测位置具有至少两个,可以拟合初始地温函数。
在一些实施例中,H个目标位置包括位于井眼内的至少一个检测位置和位于井眼表层的地表位置。
其中,H个目标位置包括位于井眼表层的地表位置,该地表位置为流体的注入口,因此,地表位置的初始地温不受流体的影响,地表位置所得到的目标温度是准确的,地表位置的初始地温即是对应的目标温度,无需预估计算,可以保证在进行初始地温函数拟合不至于偏差过多,使得拟合效果更好。
在一些实施例中,若同时根据地表位置和位于井底的检测位置对应的地层深度和初始地温拟合初始地温函数,由于地表位置和位于井底的检测位置的相关数据均比较准确,能够保证初始地温函数具有良好的拟合效果。
在一些实施例中,如图2所示,H个目标位置包括位于井眼内的M个检测位置和位于井眼表层的地表位置,
第j个检测位置的初始地温
Figure BDA0003716678670000101
其中,Ti,j为第j个检测位置的第i个目标温度,
ΔTj为第j个检测位置根据最小二乘法拟合得到误差最小条件下的温差,
Figure BDA0003716678670000102
ψ(ω)=J1(ω)J0(A1/2ω)-βJ0(ω)J1(A1/2ω),φ(ω)=J1(ω)Y0(A1/2ω)-βJ0(ω)Y1(A1/2ω),
Figure BDA0003716678670000103
其中,a1为井眼内流体的热扩散率,a2为地层的热扩散率,λ1为井眼内流体的导热系数,λ2为岩层的导热系数,t为目标温度对应的流体静置时间,R1为井眼半径,ω为虚拟复变量,J0为零阶第一类贝塞尔函数,J1为一阶第一类贝塞尔函数,Y0为零阶第二类贝塞尔函数,Y1为一阶第二类贝塞尔函数,M和N均为大于0的正整数。
其中,i=1…N,j=1…M,上述计算公式仅可根据第j个检测位置对应的目标温度获得井眼的第j个检测位置处的初始地温,该公式需要通过目标温度对应的流体静置时间计算对应的初始地温,根据相同的算法,可获得其他井深位置的初始地温Tf,j,从而得到M个检测位置的初始地温。
在一些实施例中,初始地温函数
Figure BDA0003716678670000111
其中,Tf0为地表位置的初始温度,Tf,j为第j个检测位置的初始地温,Zj为第j个检测位置对应的地层深度,Z为地层深度。
其中,地层温度通常可认为呈一定规律分布,线性分布是较简单也是行业内普遍接收的一种分布规律。则根据不同深度的地层初始温度,进行一定规律拟合后即可获得初始地温梯度。本发明实施例以线性分布为例介绍具体计算过程:
根据不同时刻、不同检测位置的目标温度,经上述计算后可得到不同检测位置处的初始地温数据对[Tf,j,Zj],j=1…M;Tf,j为j检测位置处的地层初始温度;Zj为第j检测位置处的地层深度。
假设地层温度呈线性分布,则设初始地温函数为:
Tf=Tf0+b·Z
其中Tf0为地层表面温度,通常为已知值,b为拟合系数,Z为地层深度的变量。
测量数据点与拟合直线之间的均方差为:
Figure BDA0003716678670000112
均方差极小值要满足:
Figure BDA0003716678670000121
因此可拟合得到系数为:
Figure BDA0003716678670000122
因此,初始地温函数为:
Figure BDA0003716678670000123
即可得到不同地层深度Z处的初始地温,当用户想要得到某一地层深度处的初始地温时,则对应将该地层深度值输入到初始地温函数,从而得到对应的初始地温,方便用户使用。
此外,通过在时间和空间两个维度对初始地温进行了演算,计算精度更高,初始地温的计算结果更为准确。
如图3所示,为本发明实施例提供的初始地温测量装置200的结构示意图,该测量装置200包括:
获取模块201,用于获取井眼的H个目标位置中每一目标位置的至少一个目标温度,H为大于1的正整数,至少一目标位置高于另一目标位置;
计算模块202,用于根据目标温度计算每一目标位置的初始地温;
执行模块203,用于根据H个目标位置对应的地层深度和H个初始地温,拟合初始地温函数。
在一些实施例中,H个目标位置包括位于井眼内的至少两个检测位置,其中一检测位置高于另一检测位置,获取模块201进一步包括:
第一取得单元,用于获取每一检测位置在流体静置后不同时刻的至少两个目标温度。
在一些实施例中,目标位置沿井身深度间隔分布有L个,获取模块201进一步包括:
第二取得单元,用于获取每一目标位置在流体静置后N个时刻的N个目标温度,N为大于1的正整数,L大于N。
图4是本发明实施例提供的一种计算设备,该计算设备可以包括:处理器(processor)302、通信接口(Communications Interface)304、存储器(memory)306、以及通信总线308。
其中:处理器302、通信接口304、以及存储器306通过通信总线308完成相互间的通信。通信接口304,用于与其它设备比如客户端或其它服务器等的网元通信。处理器302,用于执行程序310,具体可以执行上述初始地温测量方法实施例中的相关步骤。
具体地,程序310可以包括程序代码,该程序代码包括计算机操作指令。
处理器302可能是CPU,或者是特定集成电路ASIC(Application SpecificIntegrated Circuit),或者是被配置成实施本发明实施例的一个或多个集成电路。凑单推荐设备包括的一个或多个处理器,可以是同一类型的处理器,如一个或多个CPU;也可以是不同类型的处理器,如一个或多个CPU以及一个或多个ASIC。
存储器306,用于存放程序310。存储器306可能包含高速RAM存储器,也可能还包括非易失性存储器(non-volatile memory),例如至少一个磁盘存储器。
本发明实施例提供了一种计算机可读存储介质,存储介质中存储有至少一可执行指令,可执行指令在凑单推荐设备上运行时,使得凑单推荐设备执行如上述任意一项的初始地温测量方法的操作。
在此提供的算法或显示不与任何预设计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述,构造这类系统所要求的结构是显而易见的。此外,本发明实施例也不针对任何预设编程语言。应当明白,可以利用各种编程语言实现在此描述的本发明实施例的内容,并且上面对预设语言所做的描述是为了披露本发明实施例的一些实施方式。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实施例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
类似地,应当理解,为了精简本发明实施例并帮助理解各个发明方面中的一个或多个,在上面对本发明的示例性实施例的描述中,本发明实施例的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本发明实施例要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如下面的权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明的单独实施例。
本领域那些技术人员可以理解,可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件,以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外,可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述,本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外,本领域的技术人员能够理解,尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在下面的权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本发明的各个部件实施例可以以硬件实现,或者以在一个或者多个处理器上运行的软件模块实现,或者以它们的组合实现。本领域的技术人员应当理解,可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的一些或者全部部件的一些或者全部功能。本发明实施例还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如,计算机程序和计算机程序产品)。这样的实现本发明实施例的程序可以存储在计算机可读介质上,或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到,或者在载体信号上提供,或者以任何其他形式提供。
应该注意的是上述实施例对本发明实施例进行说明而不是对本发明进行限制,并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中,不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明实施例可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中,这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。上述实施例中的步骤,除有特殊说明外,不应理解为对执行顺序的限定。

Claims (10)

1.一种初始地温测量方法,其特征在于,所述方法包括:
获取井眼的H个目标位置中每一所述目标位置的至少一个目标温度,H为大于1的正整数,至少一所述目标位置高于另一所述目标位置;
根据所述目标温度计算每一所述目标位置的初始地温;
根据H个所述目标位置对应的地层深度和H个所述初始地温,拟合初始地温函数。
2.根据权利要求1所述的初始地温测量方法,其特征在于,H个所述目标位置包括位于所述井眼内的至少两个检测位置,其中一所述检测位置高于另一所述检测位置,所述获取井眼的H个目标位置中每一所述目标位置的至少一个目标温度,进一步包括:
获取每一所述检测位置在流体静置后不同时刻的至少两个所述目标温度。
3.根据权利要求1所述的初始地温测量方法,其特征在于,所述目标位置沿井身深度间隔分布有L个,所述获取井眼的H个目标位置中每一所述目标位置的至少一个目标温度,进一步包括:
获取每一所述目标位置在流体静置后N个时刻的N个目标温度,N为大于1的正整数,L大于N。
4.根据权利要求2~4中任一项所述的初始地温测量方法,其特征在于,至少一个所述检测位置位于所述井眼的井底,至少一个所述检测位置高于所述井底。
5.根据权利要求1所述的初始地温测量方法,其特征在于,H个所述目标位置包括位于所述井眼内的至少一个检测位置和位于井眼表层的地表位置。
6.根据权利要求1所述的初始地温测量方法,其特征在于,H个所述目标位置包括位于所述井眼内的M个检测位置和位于井眼表层的地表位置,
第j个所述检测位置的所述初始地温
Figure FDA0003716678660000021
Figure FDA0003716678660000022
其中,Ti,j为第j个检测位置的第i个目标温度,ΔTj为第j个所述目标位置根据最小二乘法拟合得到误差最小条件下的温差,
Figure FDA0003716678660000023
ψ(ω)=J1(ω)J0(A1/2ω)-βJ0(ω)J1(A1/2ω),φ(ω)=J1(ω)Y0(A1/2ω)-βJ0(ω)Y1(A1/2ω),
Figure FDA0003716678660000024
其中,a1为井眼内流体的热扩散率,a2为地层的热扩散率,λ1为井眼内流体的导热系数,λ2为岩层的导热系数,t为所述目标温度对应的流体静置时间,R1为井眼半径,ω为虚拟复变量,J0为零阶第一类贝塞尔函数,J1为一阶第一类贝塞尔函数,Y0为零阶第二类贝塞尔函数,Y1为一阶第二类贝塞尔函数,M和N均为大于0的正整数。
7.根据权利要求5所述的初始地温测量方法,其特征在于,所述初始地温函数
Figure FDA0003716678660000025
其中,Tf0为所述地表位置的所述初始温度,Tf,j为第j个检测位置的初始地温,Zj为第j个所述检测位置对应的地层深度,Z为地层深度的变量。
8.一种初始地温测量装置,其特征在于,所述测量装置包括:
获取模块,用于获取井眼的H个目标位置中每一所述目标位置的至少一个目标温度,H为大于1的正整数,至少一所述目标位置高于另一所述目标位置;
计算模块,用于根据所述目标温度计算每一所述目标位置的初始地温;
执行模块,用于根据H个所述目标位置对应的地层深度和H个所述初始地温,拟合初始地温函数。
9.一种计算设备,其特征在于,包括:处理器、存储器、通信接口和通信总线,所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信;
所述存储器用于存放至少一可执行指令,所述可执行指令使所述处理器执行如权利要求1-7中任一项所述的初始地温测量方法的操作。
10.一种计算机可读存储介质,其特征在于,所述存储介质中存储有至少一可执行指令,所述可执行指令在运行时执行如权利要求1-7中任一项所述的初始地温测量方法的操作。
CN202210738242.5A 2022-06-27 2022-06-27 初始地温测量方法、装置、计算设备及存储介质 Pending CN115045650A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210738242.5A CN115045650A (zh) 2022-06-27 2022-06-27 初始地温测量方法、装置、计算设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210738242.5A CN115045650A (zh) 2022-06-27 2022-06-27 初始地温测量方法、装置、计算设备及存储介质

Publications (1)

Publication Number Publication Date
CN115045650A true CN115045650A (zh) 2022-09-13

Family

ID=83162596

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210738242.5A Pending CN115045650A (zh) 2022-06-27 2022-06-27 初始地温测量方法、装置、计算设备及存储介质

Country Status (1)

Country Link
CN (1) CN115045650A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117892562A (zh) * 2024-03-15 2024-04-16 山东省地矿工程勘察院(山东省地质矿产勘查开发局八〇一水文地质工程地质大队) 一种地温梯度测量过程中的数据采集优化方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117892562A (zh) * 2024-03-15 2024-04-16 山东省地矿工程勘察院(山东省地质矿产勘查开发局八〇一水文地质工程地质大队) 一种地温梯度测量过程中的数据采集优化方法
CN117892562B (zh) * 2024-03-15 2024-05-17 山东省地矿工程勘察院(山东省地质矿产勘查开发局八〇一水文地质工程地质大队) 一种地温梯度测量过程中的数据采集优化方法

Similar Documents

Publication Publication Date Title
Šimůnek et al. Estimating soil hydraulic parameters from transient flow experiments in a centrifuge using parameter optimization technique
US6182013B1 (en) Methods and apparatus for dynamically estimating the location of an oil-water interface in a petroleum reservoir
US10502859B2 (en) Method and apparatus for determining resistivity of a formation
CA2940810A1 (en) Automatic method for three-dimensional structural interpretation of borehole images acquired in high-angle and horizontal wells
WO2013112827A1 (en) Method of estimating a subterranean formation property
EP2936213A1 (en) Method and apparatus for optimizing deep resistivity measurements with multi-component antennas
Abughaban et al. Advanced trajectory computational model improves calculated borehole positioning, tortuosity and rugosity
CN115045650A (zh) 初始地温测量方法、装置、计算设备及存储介质
WO2021188145A1 (en) Flow meter measurement for drilling rig
AU2014290779A1 (en) System and method for estimating porosity distribution in subterranean reservoirs
CN110192004B (zh) 到矿床边界的距离反演解的像素化
CN110007344B (zh) 一种断溶体储层连通性的地震识别方法及装置
US10935690B2 (en) Methods for inverting electromagnetic logging measurements
CN114756826A (zh) 气井产能预测方法、装置、电子设备和存储介质
US10488553B2 (en) Stress tensor computation using Mindlin formulation
Bikmukhametov et al. Statistical analysis of effect of sensor degradation and heat transfer modeling on multiphase flowrate estimates from a virtual flow meter
CN112711076B (zh) 提取石油钻井中泥浆侵入地层深度的方法和装置
Sawdey et al. Automated inverse computer modeling of borehole flow data in heterogeneous aquifers
Bauer et al. Computation of the equivalent macroscopic permeability tensor of discrete networks with heterogeneous segment length
EP3359777B1 (en) Systems and methods to calibrate individual component measurement
Sun et al. Measuring Groundwater Velocity: Method Based on Groundwater Flow–Induced Cooling
Zhuravchak et al. Study of the Anisotropy Effect in Piecewise Homogeneous Media by Boundary and Near-Boundary Element Methods
US10690801B2 (en) Skin effect correction for focused electrode devices based on analytical model
Sindalovskiy et al. Recovery Tests
CN117951857A (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