CN110593860B - 一种参数计算方法、装置及设备 - Google Patents

一种参数计算方法、装置及设备 Download PDF

Info

Publication number
CN110593860B
CN110593860B CN201810516078.7A CN201810516078A CN110593860B CN 110593860 B CN110593860 B CN 110593860B CN 201810516078 A CN201810516078 A CN 201810516078A CN 110593860 B CN110593860 B CN 110593860B
Authority
CN
China
Prior art keywords
pressure
oil well
data
karst cave
historical
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
CN201810516078.7A
Other languages
English (en)
Other versions
CN110593860A (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.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
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 University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN201810516078.7A priority Critical patent/CN110593860B/zh
Publication of CN110593860A publication Critical patent/CN110593860A/zh
Application granted granted Critical
Publication of CN110593860B publication Critical patent/CN110593860B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP 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
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Abstract

本发明提供了一种参数计算方法、装置及设备,本发明中获取油井的参数以及所述油井的历史流量数据和历史压力数据,基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据,构建目标试井的数学模型,根据数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图,基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。通过上述方法,可以计算得到溶洞体积,解决了现有技术中没有一种计算溶洞体积的方法的问题。

Description

一种参数计算方法、装置及设备
技术领域
本发明涉及油藏工程技术领域,更具体的说,涉及一种参数计算方法、装置及设备。
背景技术
碳酸盐岩缝洞型油藏基质基本不含油,储集空间以裂缝、溶洞为主,原油在裂缝、溶洞中的流动存在管流和渗流。原油在深大断裂裂缝和溶洞中流动时,垂向流动明显。
溶洞中存储油,如果已知溶洞体积,就可以确定溶洞中存储的油量,便于开采石油。
但是现有技术中,并没有一种计算溶洞体积的方法。
发明内容
有鉴于此,本发明提供一种参数计算方法、装置及设备,以解决现有技术中没有一种计算溶洞体积的方法的问题。
为解决上述技术问题,本发明采用了如下技术方案:
一种参数计算方法,包括:
获取油井的参数、所述油井的历史流量数据和历史压力数据;
基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据;
构建目标试井的数学模型;其中,所述目标试井包括周围介质不同的油井和溶洞;
根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图;
基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。
优选地,所述构建目标试井的数学模型包括:
获取所述油井周围的第一介质的参数和所述溶洞周围的第二介质的参数;其中,所述第一介质与所述第二介质不同;
基于预设油井参数、预设溶洞参数、所述第一介质的参数和所述第二介质的参数,构建所述目标试井的数学模型。
优选地,根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图,包括:
对所述数学模型进行求解,计算得到所述油井的压力值;
基于所述油井的参数,对所述压力值进行数值反演操作,得到所述油井的压力与时间的对应关系;
根据所述油井的压力与时间的对应关系,得到所述油井在不同时刻的理论压力值;
基于所述油井在不同时刻的理论压力值,通过曲线拟合的方法,绘制得到所述理论压力曲线图。
优选地,所述数学模型中包括油井压力和溶洞压力的对应关系;
所述油井压力和溶洞压力的对应关系的生成过程包括:
获取预设流体流动模型;
基于所述预设流体流动模型,计算得到溶洞中流体的流速值;
根据所述流体的流速值,计算得到流体从所述溶洞流入到所述油井中的速度值;
基于所述速度值以及所述预设流体流动模型,计算得到油井压力和溶洞压力的对应关系。
优选地,基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积,包括:
根据所述实际压力数据和所述理论压力曲线图,计算得到波动系数和阻力系数;
基于波动系数、阻力系数、溶洞半径、溶洞深度之间的对应关系,计算得到溶洞半径和溶洞深度;
根据所述溶洞半径和所述溶洞深度,计算得到所述溶洞体积。
优选地,基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据,包括:
从所述历史流量数据中筛选出持续的流量值在第一预设数值内的历史流量数据;
从所述历史压力数据中筛选出与筛选出的历史流量数据相对应的历史压力数据;
将筛选出的历史压力数据作为所述油井的实际压力数据。
优选地,根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图后,还包括:
调整预设输入参数的数值,以使所述理论压力曲线图与所述油井的实际压力数据的点线重合率大于第二预设数值,以及所述理论压力曲线图对应的理论压力导数曲线图与所述油井的实际压力数据对应的实际压力导数数据的点线重合率大于第三预设数值。
一种参数计算装置,包括:
数据获取模块,用于获取油井的参数、所述油井的历史流量数据和历史压力数据;
第一分析模块,用于基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据;
模型构建模块,用于构建目标试井的数学模型;其中,所述目标试井包括周围介质不同的油井和溶洞;
第二分析模块,用于根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图;
体积计算模块,用于基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。
优选地,所述第二分析模块包括:
第一计算子模块,用于对所述数学模型进行求解,计算得到所述油井的压力值;
关系确定子模块,用于基于所述油井的参数,对所述压力值进行数值反演操作,得到所述油井的压力与时间的对应关系;
第二计算子模块,用于根据所述油井的压力与时间的对应关系,得到所述油井在不同时刻的理论压力值;
曲线图绘制子模块,用于基于所述油井在不同时刻的理论压力值,通过曲线拟合的方法,绘制得到所述理论压力曲线图。
优选地,所述数学模型中包括油井压力和溶洞压力的对应关系;
参数计算装置还包括:
模型获取模块,用于获取预设流体流动模型;
流速值计算模块,用于基于所述预设流体流动模型,计算得到溶洞中流体的流速值;
速度值计算模块,用于根据所述流体的流速值,计算得到流体从所述溶洞流入到所述油井中的速度值;
关系确定模块,用于基于所述速度值以及所述预设流体流动模型,计算得到油井压力和溶洞压力的对应关系。
一种电子设备,包括:存储器和处理器;
其中,所述存储器用于存储程序;
处理器调用程序并用于:
获取油井的参数、所述油井的历史流量数据和历史压力数据;
基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据;
构建目标试井的数学模型;其中,所述目标试井包括周围介质不同的油井和溶洞;
根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图;
基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。
相较于现有技术,本发明具有以下有益效果:
本发明提供了一种参数计算方法、装置及设备,本发明中获取油井的参数以及所述油井的历史流量数据和历史压力数据,基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据,构建目标试井的数学模型,根据数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图,基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。通过上述方法,可以计算得到溶洞体积,解决了现有技术中没有一种计算溶洞体积的方法的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明实施例提供的一种双层均匀介质模型示意图;
图2为本发明实施例提供的一种参数计算方法的方法流程图;
图3为本发明实施例提供的另一种参数计算方法的方法流程图;
图4为本发明实施例提供的再一种参数计算方法的方法流程图;
图5为本发明实施例提供的一种流体微元的结构示意图;
图6为本发明实施例提供的一种对数曲线图;
图7为本发明实施例提供的另一种对数曲线图;
图8为本发明实施例提供的一种参数计算装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供了一种参数计算方法,该参数计算方法适用于缝洞型储层单洞-双层均质。双层均匀介质模型示意图可以参照图1。
油井的井筒周围的介质为均匀介质1,溶洞的周围介质为均匀介质2。均匀介质1和2不同,并且溶洞的数量为一个,所以称为单洞-双层均质。
另外,井筒的深度为h1,半径为rw,溶洞的深度为h2,半径为rv
参照图2,参数计算方法可以包括:
S11、获取油井的参数、所述油井的历史流量数据和历史压力数据;
具体的,油井的参数可以参照表1。
表1油井的参数表
Figure BDA0001673190200000061
表1中,给出了油井的参数,油井的参数包括地层厚度、油井半径、孔隙度、流体粘度、体积系数、压缩系数、中部深度、流体密度和流体压缩系数。油井的参数值具体参照表1中的相应内容。
油井的历史流量数据和历史压力数据是油井在抽油的过程中的数据。油井在抽油的过程中,每一时间均对应有一历史流量数据和历史压力数据。
S12、基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据;
其中,实际压力数据是从历史压力数据中筛选出的压力值较稳定的连续数据。
S13、构建目标试井的数学模型;
其中,所述目标试井包括周围介质不同的油井和溶洞。
可选的,在本实施例的基础上,步骤S13可以包括:
1)获取所述油井周围的第一介质的参数和所述溶洞周围的第二介质的参数;其中,所述第一介质与所述第二介质不同;
其中,第一介质为上文中的均匀介质1,第二介质为上文中的均匀介质2。
第一介质和第二介质不同,可以是土壤密度不同、土壤成分不同或土壤硬度不同等。
2)基于预设油井参数、预设溶洞参数、所述第一介质的参数和所述第二介质的参数,构建所述目标试井的数学模型。
具体的,预设油井参数、预设溶洞参数、所述第一介质的参数和所述第二介质的参数是技术人员根据具体的使用场景进行设定的,具体参数可以参照下文中构建的目标试井的数学模型中的参数。
S14、根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图;
理论压力曲线图是根据推理计算得到理论压力值绘制得到。
S15、基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。
本实施例中,获取油井的参数以及所述油井的历史流量数据和历史压力数据,基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据,构建目标试井的数学模型,根据数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图,基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。通过上述方法,可以计算得到溶洞体积,解决了现有技术中没有一种计算溶洞体积的方法的问题。
可选的,在上述任一参数计算方法的实施例的基础上,参照图3,步骤S14可以包括:
S21、对所述数学模型进行求解,计算得到所述油井的压力值;
其中,数学模型为一公式组,且数学模型可以为:
Figure BDA0001673190200000071
其中,pi原始地层压力;pv溶洞压力;pw井底压力;k外部介质渗透率(um2);μ流体粘度;φ外部介质孔隙度;Ct外部介质压缩系数;Q产量;B体积系数;rw井筒半径;rv溶洞半径;Cw井筒储集系数;Cv为溶洞储集常数;Sw井筒的表皮系数;Sv溶洞的表皮系数;h1与井筒相连地层的厚度,也即井筒的深度;h2与溶洞相连地层的厚度,也即溶洞的深度,t表示时间。
设定,
无因次压力
Figure BDA0001673190200000072
无因次时间
Figure BDA0001673190200000081
无因次半径
Figure BDA0001673190200000082
无因次渗透率
Figure BDA0001673190200000083
无因次井筒储集系数
Figure BDA0001673190200000084
无因次溶洞储集系数
Figure BDA0001673190200000085
无因次高度
Figure BDA0001673190200000086
无因次波动系数
Figure BDA0001673190200000087
无因次阻尼系数
Figure BDA0001673190200000088
无因此窜流系数
Figure BDA0001673190200000089
无因此储水系数
Figure BDA00016731902000000810
方程修正系数α
第二类欧拉积分函数
Figure BDA00016731902000000811
其中,下标带w或1的表示井筒的数据,下标带v或2的表示溶洞的数据。
对所述试井模型进行拉普拉斯变换,解得拉普拉斯空间的油井的压力值为:
Figure BDA00016731902000000812
其中,u为拉普拉斯变量。
S22、基于所述油井的参数,对所述压力值进行数值反演操作,得到所述油井的压力与时间的对应关系;
利用Stehfest数值反演技术,由所述拉普拉斯空间的压力值
Figure BDA0001673190200000091
得到所述真实空间的油井的压力pwD(tD),也即得到所述油井的压力与时间的对应关系。
Figure BDA0001673190200000092
其中,N为偶数,一般取值在8-16之间。
S23、根据所述油井的压力与时间的对应关系,得到所述油井在不同时刻的理论压力值;
具体的,根据油井的压力与时间的对应关系,选取一段时间,计算得到油井在不同时刻的压力值,即理论压力值。
S24、基于所述油井在不同时刻的理论压力值,通过曲线拟合的方法,绘制得到所述理论压力曲线图。
具体的,拟合就是把平面上一系列的点,用一条光滑的曲线连接起来。因为这条曲线有无数种可能,从而有各种拟合方法。
具体的,油井在不同时刻的理论压力值为一个一个的点数据。用一条光滑的曲线将油井在不同时刻的理论压力值进行连接,即可得到理论压力曲线图。
本实施例中,能够在计算得到所述油井的压力与时间的对应关系后,计算得到油井在不同时刻的理论压力值,然后根据油井在不同时刻的理论压力值,绘制得到所述理论压力曲线图。
另外,数学模型包括井钻遇大溶洞并且外部地层为上下两层均匀介质地层,数学模型通过联立力学三大守恒定律即质量守恒、动量守恒和能量守恒方程提出了波动和流动耦合的试井模型,该模型可以对缝洞型油藏进行压力恢复数据的试井解释,使用该模型不仅可以得到地层渗透率,而且还能够计算洞体积、阻尼系数、波动系数等储层信息。
再次,本实施例为缝洞型油藏储层参数计算及储量计算提供了配套的解释方法,能够直接解释缝洞型油藏中大溶洞的体积,为油田开发方案的制定提供技术支撑。本发明具有两大优势:(1)模型简单,求解方便,在拉普拉斯空间可以给出解析解,并且解析解不涉及复杂函数的计算,计算速度快;(2)拟合后的解释结果可直接给出溶洞体积的大小,可方便的用于计算单井储量。
可选的,在上述任一实施例的基础上,所述数学模型中包括油井压力和溶洞压力的对应关系。其中,油井压力和溶洞压力的对应关系为数学模型中的最后一个公式。
进一步,参照图4,所述油井压力和溶洞压力的对应关系的生成过程可以包括:
S31、获取预设流体流动模型;
流体从溶洞流入井筒,再由井筒流出地面,这一过程的流体流动要满足的连续性方程、动量守恒和能量守恒方程为:
Figure BDA0001673190200000101
Figure BDA0001673190200000102
Figure BDA0001673190200000103
上述方程中,ρ为流体密度(单位:kg/m3);v为流体流动速度(单位:m/s);x轴为由井筒圆心向下建立的一维坐标轴;p为压力(单位:Pa);f流体受到的摩擦阻力系数;D为井筒的直径(单位:m);pw和pv分别为井筒和溶洞中的压力(单位:Pa);vw为井筒和溶洞连接处流体的速度(单位:m/s),g为重力。
其中,上述公式123即为预设流体流动模型。
S32、基于所述预设流体流动模型,计算得到溶洞中流体的流速值;
参照图5,取井筒中的一个流体微元,由质量守恒可得到:
Figure BDA0001673190200000104
其中,A为流体微元的横截面积。
在高压状况下,流体存在压缩性,井筒也是一个弹性体,其变形由井筒直径、壁厚及井筒材料的杨氏模量决定,展开方程(4)变成:
Figure BDA0001673190200000105
根据流体力学中全导数和偏导数关系:
Figure BDA0001673190200000111
方程(5)可修改为:
Figure BDA0001673190200000112
考虑流体的压缩性,方程(7)中的密度项可表示成压力的函数:
Figure BDA0001673190200000113
式中:G为流体的体积模量(单位:Pa)。假设井筒为弹性变形,对薄壁的圆管,当压力增加dp时,其径向变形dD与dp的关系为:
Figure BDA0001673190200000114
式中:D为井筒直径(单位:m);e为井筒壁厚(单位:m);E为井筒杨氏模量(单位:Pa)。由井筒面积公式:
Figure BDA0001673190200000115
联合(8)和(9),则(7)可变成:
Figure BDA0001673190200000116
定义:
Figure BDA0001673190200000117
式中:C为管道及流体系统中的波速,m/s。
利用全导数公式,方程(11)可变成:
Figure BDA0001673190200000118
可见压力在x方向的传导是以波的形式进行的,波速为C.
如果将溶洞也视作圆柱,连续性方程同(13),但压力传播的波速C可表示为:
Figure BDA0001673190200000119
式中:E为地层的杨氏模量(单位:Pa);φ为孔隙度。
将连续性方程和动量守恒方程联合,可得:
Figure BDA00016731902000001110
考察流体在溶洞中的流动,由于速度v较小,忽略重力和流体摩擦力,并考虑溶洞的存储常数Cv,由于摩擦力为v的二次方项,可略去二阶小量则有:
Figure BDA0001673190200000121
Figure BDA0001673190200000122
方程(16)的解,即为溶洞中的流体流速值:
Figure BDA0001673190200000123
式中:v0为初始时刻的速度(单位:m/s),可由地面产量确定。
S33、根据所述流体的流速值,计算得到流体从所述溶洞流入到所述油井中的速度值;
具体的,由于rv是大溶洞的半径,这个值很大,所以rv 2项更大,同时这一项是一个常数,可以将该项加入到附加压降中,并且我们在后面对方程进行了修正,所以我们可以忽略方程(18)右边第二项,即:
Figure BDA0001673190200000124
由此得到溶洞提供的产量:
Figure BDA0001673190200000125
于是流体从溶洞流入井筒处的速度:
Figure BDA0001673190200000126
S34、基于所述速度值以及所述预设流体流动模型,计算得到油井压力和溶洞压力的对应关系。
根据井筒及溶洞处的能量守恒方程,则有:
Figure BDA0001673190200000127
按照前面的无量纲定义可将方程(22)化为下边的方程*
Figure BDA0001673190200000128
将公式*进行拉普拉斯逆变换,即可得到油井压力和溶洞压力的对应关系。
本实施例中,给出了油井压力和溶洞压力的对应关系的生成过程,进而可以生成油井压力和溶洞压力的对应关系,并构建数学模型。
可选的,在上述任一实施例的基础上,步骤S15可以包括:
1)根据所述实际压力数据和所述理论压力曲线图,计算得到波动系数和阻力系数;
具体的,通过实际压力数据和所述理论压力曲线图的曲线拟合可以得到波动系数β和阻力系数γ的值。
2)基于所述波动系数、所述阻力系数、溶洞半径、溶洞深度之间的对应关系,计算得到溶洞半径和溶洞深度;
根据无量纲定义,β是跟溶洞半径rv和溶洞的深度h2相关,而γ跟溶洞的半径rv相关,现在得到了β和γ的值就可以反推出rv和h2的值。
3)根据所述溶洞半径和所述溶洞深度,计算得到所述溶洞体积。
设定溶洞是一个圆柱体,所以溶洞的体积计算就是
Figure BDA0001673190200000131
h1是井筒的深度,这是一个已知量,跟表2中的中部深度是相对应的(近似相等)。
另外,本实施例中,还能够计算得到井储常数、井筒表皮、渗透率比、溶洞波动系数和溶洞阻尼系数。
具体计算得到的数据参照表2。
表2曲线拟合得到的数据的示意表
Figure BDA0001673190200000132
本实施例中,给出了计算溶洞体积的方法,进而就可以根据本实施例中的方法计算得到溶洞体积,进而就可以推测溶洞中的储油量,以及较准确的井筒的半径。
可选的,在上述任一实施例的基础上,步骤S12可以包括:
1)从所述历史流量数据中筛选出持续的流量值在第一预设数值内的历史流量数据;
具体的,第一预设数值是技术人员根据具体的使用场景进行设定的。
筛选出持续的流量值在第一预设数值内的历史流量数据是为了获取流量稳定的一段历史流量数据,避免筛选出的历史流量数据是在开井和关井时的数据。
2)从所述历史压力数据中筛选出与筛选出的历史流量数据相对应的历史压力数据;
筛选出历史流量数据后,获取筛选出的历史流量数据对应的时间段,然后获取位于该时间段内的历史压力数据即可。
3)将筛选出的历史压力数据作为所述油井的实际压力数据。
进一步,在本实施例的基础上,步骤S14后,还包括:
调整预设输入参数的数值,以使所述理论压力曲线图与所述油井的实际压力数据的点线重合率大于第二预设数值,以及所述理论压力曲线图对应的理论压力导数曲线图与所述油井的实际压力数据对应的实际压力导数数据的点线重合率大于第三预设数值。
其中,预设输入参数可以参照表3。
表3预设输入参数示意表
Figure BDA0001673190200000141
具体的,预设输入参数可以包括无量纲渗透率、无量纲储水系数、无量纲窜流系数、波动系数、阻尼系数、衰减指数和组合参数。
需要说明的是,预设输入参数是技术人员提前设定的。
得到所述油井在不同时刻的理论压力值与理论压力曲线图、实际压力数据后,对预设输入参数进行调整,以使分散的点与曲线图的拟合度较高,点线重合率较高。
其中,第二预设数值和第三预设数值是技术人员根据具体的场景进行设定的。
需要说明的是,不管是理论压力曲线图,还是实际压力数据,均为压力及压力导数无量纲量双对数图,具体参照图6和图7。
图6中,实线为压力无量纲量对数,虚线为压力导数无量纲量对数。
图7中,实线为理论压力曲线图,正方形为油井在不同时刻的实测压力,即为所述油井的实际压力数据,虚线为理论压力导数曲线图,三角为实测导数,即为实际压力导数数据。
可以看出,理论压力曲线图与油井的实际压力数据进行拟合,理论压力导数曲线图与实际压力导数数据拟合。
本实施例中,可以基于所述历史流量数据和所述历史压力数据,分析得到所述油井的实际压力数据。并且,还可以调整曲线的拟合率,使得曲线拟合的更好。
可选的,在上述参数计算方法的实施例的基础上,本发明的另一实施例提供了一种参数计算装置,参照图8,可以包括:
数据获取模块101,用于获取油井的参数、所述油井的历史流量数据和历史压力数据;
第一分析模块102,用于基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据;
模型构建模块103,用于构建目标试井的数学模型;其中,所述目标试井包括周围介质不同的油井和溶洞;
第二分析模块104,用于根据数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图;
体积计算模块105,用于基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。
可选的,在本实施例的基础上,模型构建模块包括:
参数获取子模块,用于获取所述油井周围的第一介质的参数和所述溶洞周围的第二介质的参数;其中,所述第一介质与所述第二介质不同;
模型构建子模块,用于基于预设油井参数、预设溶洞参数、所述第一介质的参数和所述第二介质的参数,构建所述目标试井的数学模型。
本实施例中,获取油井的参数以及所述油井的历史流量数据和历史压力数据,基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据,构建目标试井的数学模型,根据数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图,基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。通过上述方法,可以计算得到溶洞体积,解决了现有技术中没有一种计算溶洞体积的方法的问题。
需要说明的是,本实施例中的各个模块的工作过程,请参照上述实施例中的相应说明,在此不再赘述。
可选的,在上述参数计算装置的实施例的基础上,所述第二分析模块包括:
第一计算子模块,用于对所述数学模型进行求解,计算得到所述油井的压力值;
关系确定子模块,用于基于所述油井的参数,对所述压力值进行数值反演操作,得到所述油井的压力与时间的对应关系;
第二计算子模块,用于根据所述油井的压力与时间的对应关系,得到所述油井在不同时刻的理论压力值;
曲线图绘制子模块,用于基于所述油井在不同时刻的理论压力值,通过曲线拟合的方法,绘制得到所述理论压力曲线图。
本实施例中,能够在计算得到所述油井的压力与时间的对应关系后,计算得到油井在不同时刻的理论压力值,然后根据油井在不同时刻的理论压力值,绘制得到所述理论压力曲线图。
需要说明的是,本实施例中的各个模块和子模块的工作过程,请参照上述实施例中的相应说明,在此不再赘述。
可选的,在上述任一参数计算装置的实施例的基础上,所述数学模型中包括油井压力和溶洞压力的对应关系;
参数计算装置还包括:
模型获取模块,用于获取预设流体流动模型;
流速值计算模块,用于基于所述预设流体流动模型,计算得到溶洞中流体的流速值;
速度值计算模块,用于根据所述流体的流速值,计算得到流体从所述溶洞流入到所述油井中的速度值;
关系确定模块,用于基于所述速度值以及所述预设流体流动模型,计算得到油井压力和溶洞压力的对应关系。
本实施例中,给出了油井压力和溶洞压力的对应关系的生成过程,进而可以生成油井压力和溶洞压力的对应关系,并构建数学模型。
需要说明的是,本实施例中的各个模块的工作过程,请参照上述实施例中的相应说明,在此不再赘述。
可选的,在上述任一参数计算装置的实施例的基础上,体积计算模块用于基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积时,具体用于:
根据所述实际压力数据和所述理论压力曲线图,计算得到波动系数和阻力系数;
基于所述波动系数、所述阻力系数、溶洞半径、溶洞深度之间的对应关系,计算得到溶洞半径和溶洞深度;
根据所述溶洞半径和所述溶洞深度,计算得到所述溶洞体积。
本实施例中,给出了计算溶洞体积的方法,进而就可以根据本实施例中的方法计算得到溶洞体积,进而就可以推测溶洞中的储油量,以及较准确的井筒的半径。
需要说明的是,本实施例中的各个模块的工作过程,请参照上述实施例中的相应说明,在此不再赘述。
可选的,在上述任一参数计算装置的实施例的基础上,第一分析模块用于基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据时,具体用于:
从所述历史流量数据中筛选出持续的流量值在第一预设数值内的历史流量数据;
从所述历史压力数据中筛选出与筛选出的历史流量数据相对应的历史压力数据;
将筛选出的历史压力数据作为所述油井的实际压力数据。进一步,还包括:
参数调整模块,用于第二分析模块根据数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图后,调整预设输入参数的数值,以使所述理论压力曲线图与所述油井的实际压力数据的点线重合率大于第二预设数值,以及所述理论压力曲线图对应的理论压力导数曲线图与所述油井的实际压力数据对应的实际压力导数数据的点线重合率大于第三预设数值。
本实施例中,可以基于所述历史流量数据和所述历史压力数据,分析得到所述油井的实际压力数据。并且,还可以调整曲线的拟合率,使得曲线拟合的更好。
需要说明的是,本实施例中的各个模块的工作过程,请参照上述实施例中的相应说明,在此不再赘述。
可选的,在上述参数计算方法和装置的实施例的基础上,本发明的另一实施例提供了一种电子设备,包括:存储器和处理器;
其中,所述存储器用于存储程序;
处理器调用程序并用于:
获取油井的参数、所述油井的历史流量数据和历史压力数据;
基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据;
构建目标试井的数学模型;其中,所述目标试井包括周围介质不同的油井和溶洞;
根据数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图;
基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。
本实施例中,获取油井的参数以及所述油井的历史流量数据和历史压力数据,基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据,,构建目标试井的数学模型,根据数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图,基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积。通过上述方法,可以计算得到溶洞体积,解决了现有技术中没有一种计算溶洞体积的方法的问题。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (8)

1.一种参数计算方法,其特征在于,包括:
获取油井的参数、所述油井的历史流量数据和历史压力数据;
基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据;
构建目标试井的数学模型;其中,所述目标试井包括周围介质不同的油井和溶洞;
根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图;
基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积;
所述数学模型中包括油井压力和溶洞压力的对应关系;
所述油井压力和溶洞压力的对应关系的生成过程包括:
获取预设流体流动模型;
基于所述预设流体流动模型,计算得到溶洞中流体的流速值;
根据所述流体的流速值,计算得到流体从所述溶洞流入到所述油井中的速度值;
基于所述速度值以及所述预设流体流动模型,计算得到油井压力和溶洞压力的对应关系。
2.根据权利要求 1 所述的参数计算方法,其特征在于,所述构建目标试井的数学模型包括:
获取所述油井周围的第一介质的参数和所述溶洞周围的第二介质的参数;其中,所述第一介质与所述第二介质不同;
基于预设油井参数、预设溶洞参数、所述第一介质的参数和所述第二介质的参数,构建所述目标试井的数学模型。
3.根据权利要求 1 所述的参数计算方法,其特征在于,根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图,包括:
对所述数学模型进行求解,计算得到所述油井的压力值;
基于所述油井的参数,对所述压力值进行数值反演操作,得到所述油井的压力与时间的对应关系;
根据所述油井的压力与时间的对应关系,得到所述油井在不同时刻的理论压力值;
基于所述油井在不同时刻的理论压力值,通过曲线拟合的方法,绘制得到所述理论压力曲线图。
4.根据权利要求 1 所述的参数计算方法,其特征在于,基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积,包括:
根据所述实际压力数据和所述理论压力曲线图,计算得到波动系数和阻力系数;
基于波动系数、阻力系数、溶洞半径、溶洞深度之间的对应关系,计算得到溶洞半径和溶洞深度;
根据所述溶洞半径和所述溶洞深度,计算得到所述溶洞体积。
5.根据权利要求 3 所述的参数计算方法,其特征在于,基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据,包括:
从所述历史流量数据中筛选出持续的流量值在第一预设数值内的历史流量数据;
从所述历史压力数据中筛选出与筛选出的历史流量数据相对应的历史压力数据;
将筛选出的历史压力数据作为所述油井的实际压力数据。
6.根据权利要求 5 所述的参数计算方法,其特征在于,根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图后,还包括:
调整预设输入参数的数值,以使所述理论压力曲线图与所述油井的实际压力数据的点线重合率大于第二预设数值,以及所述理论压力曲线图对应的理论压力导数曲线图与所述油井的实际压力数据对应的实际压力导数数据的点线重合率大于第三预设数值。
7.一种参数计算装置,其特征在于,包括:
数据获取模块,用于获取油井的参数、所述油井的历史流量数据和历史压力数据;
第一分析模块,用于基于所述历史流量数据,从所述历史压力数据中筛选得到所述油井的实际压力数据;模型构建模块,用于构建目标试井的数学模型;其中,所述目标试井包
括周围介质不同的油井和溶洞;
第二分析模块,用于根据所述数学模型以及所述油井的参数,分析得到所述油井的理论压力曲线图;
体积计算模块,用于基于所述实际压力数据和所述理论压力曲线图,计算得到溶洞体积;
所述数学模型中包括油井压力和溶洞压力的对应关系;
参数计算装置还包括:
模型获取模块,用于获取预设流体流动模型;
流速值计算模块,用于基于所述预设流体流动模型,计算得到溶洞中流体的流速值;
速度值计算模块,用于根据所述流体的流速值,计算得到流体从所述溶洞流入到所述油井中的速度值;
关系确定模块,用于基于所述速度值以及所述预设流体流动模型,计算得到油井压力和溶洞压力的对应关系。
8.根据权利要求 7 所述的参数计算装置,其特征在于,所述第二分析模块包括:
第一计算子模块,用于对所述数学模型进行求解,计算得到所述油井的压力值;
关系确定子模块,用于基于所述油井的参数,对所述压力值进行数值反演操作,得到所述油井的压力与时间的对应关系;
第二计算子模块,用于根据所述油井的压力与时间的对应关系,得到所述油井在不同时刻的理论压力值;
曲线图绘制子模块,用于基于所述油井在不同时刻的理论压力值,通过曲线拟合的方法,绘制得到所述理论压力曲线图。
CN201810516078.7A 2018-05-25 2018-05-25 一种参数计算方法、装置及设备 Active CN110593860B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810516078.7A CN110593860B (zh) 2018-05-25 2018-05-25 一种参数计算方法、装置及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810516078.7A CN110593860B (zh) 2018-05-25 2018-05-25 一种参数计算方法、装置及设备

Publications (2)

Publication Number Publication Date
CN110593860A CN110593860A (zh) 2019-12-20
CN110593860B true CN110593860B (zh) 2021-08-13

Family

ID=68848717

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810516078.7A Active CN110593860B (zh) 2018-05-25 2018-05-25 一种参数计算方法、装置及设备

Country Status (1)

Country Link
CN (1) CN110593860B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11231520B2 (en) 2020-05-06 2022-01-25 Saudi Arabian Oil Company Dynamic hydrocarbon well skin modeling and operation
US11692415B2 (en) 2020-06-22 2023-07-04 Saudi Arabian Oil Company Hydrocarbon well stimulation based on skin profiles
CN112267872B (zh) * 2020-10-27 2021-09-14 成都理工大学 利用生产动态数据确定气井钻遇溶洞系统容积大小的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102339339A (zh) * 2010-07-23 2012-02-01 中国石油化工股份有限公司 一种分析缝洞型油藏剩余油分布的方法
CN106246171A (zh) * 2016-09-09 2016-12-21 西南石油大学 部分连通断层边界的无限大双重介质油藏数学建模方法
CN106599449A (zh) * 2016-12-12 2017-04-26 中国石油化工股份有限公司 一种溶洞体积计算的试井解释方法
CN107563027A (zh) * 2017-08-21 2018-01-09 西南石油大学 用于分析多缝洞单元串联结构的试井解释模型与方法
CN107577831A (zh) * 2017-03-17 2018-01-12 西南石油大学 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8620636B2 (en) * 2005-08-25 2013-12-31 Schlumberger Technology Corporation Interpreting well test measurements
CN103452547B (zh) * 2012-08-09 2016-03-23 凌毅立 试井资料中续流数据的分析处理方法和系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102339339A (zh) * 2010-07-23 2012-02-01 中国石油化工股份有限公司 一种分析缝洞型油藏剩余油分布的方法
CN106246171A (zh) * 2016-09-09 2016-12-21 西南石油大学 部分连通断层边界的无限大双重介质油藏数学建模方法
CN106599449A (zh) * 2016-12-12 2017-04-26 中国石油化工股份有限公司 一种溶洞体积计算的试井解释方法
CN107577831A (zh) * 2017-03-17 2018-01-12 西南石油大学 一种缝洞型碳酸盐岩油气藏溶洞规模大小计算方法
CN107563027A (zh) * 2017-08-21 2018-01-09 西南石油大学 用于分析多缝洞单元串联结构的试井解释模型与方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
缝洞型碳酸盐岩油藏试井分析模型研究;邢翠巧 等;《河北工业科技》;20180131;第35卷(第1期);第31-36页 *

Also Published As

Publication number Publication date
CN110593860A (zh) 2019-12-20

Similar Documents

Publication Publication Date Title
CN110593860B (zh) 一种参数计算方法、装置及设备
CN108798634B (zh) 一种缝洞型储层单洞-单层均值试井解释方法
Brinkgreve Selection of soil models and parameters for geotechnical engineering application
Qiu et al. Controlled installation of spudcan foundations on loose sand overlying weak clay
CN107563899B (zh) 油气井产能预测方法及装置
Silva et al. An analytical study of the effect of penetration rate on piezocone tests in clay
CN110096718B (zh) 一种获取碳酸盐岩油藏中溶洞的体积的方法
EP3455654A1 (en) Three-dimensional fracture abundance evaluation of subsurface formation based on geomechanical simulation of mechanical properties thereof
CN105701299A (zh) 一种动态摩阻扭矩计算方法
CN110984973A (zh) 缝洞型碳酸盐岩气藏单井控制储量的确定方法
CN110107277B (zh) 一种获取碳酸盐岩油藏中井钻遇溶洞的体积的方法
CN109359420A (zh) 不同工况下射孔对封隔器冲击压力预测方法及装置
CN114427445A (zh) 缝洞型油藏无限大地层动态产能计算方法和系统
Wang et al. Large deformation finite element analysis of the installation of suction caisson in clay
EP3455459A2 (en) Three-dimensional fracture abundance evaluation of subsurface formations
CN110532579B (zh) 一种参数计算方法、装置及设备
US20190042676A1 (en) Simulating multi-dimensional flow with coupled one-dimensional flow paths
Ferreri et al. Storm sewer pressurization transient–an experimental investigation
CN110096669B (zh) 一种获取碳酸盐岩油藏中裂缝体积的方法
CN108733895A (zh) 套管参数和水泥环参数的确定方法和装置
CN111520133B (zh) 一种确定地层中洞体积的方法
CN109723435A (zh) 确定页岩气藏水平井重复压裂前现地应力场的方法及设备
CN111767670B (zh) 基于井下仰斜单孔放水试验的含水层地质参数反演方法
CN113887110A (zh) 油井产能预测方法、装置、电子设备及计算机存储介质
CN109522579B (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