CN113958299A - 一种基于温度测量的水平井吸水剖面反演方法 - Google Patents
一种基于温度测量的水平井吸水剖面反演方法 Download PDFInfo
- Publication number
- CN113958299A CN113958299A CN202111239021.5A CN202111239021A CN113958299A CN 113958299 A CN113958299 A CN 113958299A CN 202111239021 A CN202111239021 A CN 202111239021A CN 113958299 A CN113958299 A CN 113958299A
- Authority
- CN
- China
- Prior art keywords
- temperature
- equation
- water injection
- well
- inversion
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 39
- 238000009529 body temperature measurement Methods 0.000 title claims abstract description 21
- 235000020681 well water Nutrition 0.000 title claims abstract description 16
- 239000002349 well water Substances 0.000 title claims abstract description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 109
- 238000002347 injection Methods 0.000 claims abstract description 101
- 239000007924 injection Substances 0.000 claims abstract description 101
- 238000012546 transfer Methods 0.000 claims abstract description 40
- 230000008569 process Effects 0.000 claims abstract description 22
- 239000000243 solution Substances 0.000 claims abstract description 19
- 238000010276 construction Methods 0.000 claims abstract description 10
- 239000013307 optical fiber Substances 0.000 claims abstract description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 27
- 239000012530 fluid Substances 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000004134 energy conservation Methods 0.000 claims description 8
- 239000012466 permeate Substances 0.000 claims description 6
- 238000009795 derivation Methods 0.000 claims description 4
- 229920006395 saturated elastomer Polymers 0.000 claims description 4
- 230000035939 shock Effects 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 4
- 239000011541 reaction mixture Substances 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 description 9
- 238000005259 measurement Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 238000007418 data mining Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000819 phase cycle Methods 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/20—Displacing by water
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
- E21B47/07—Temperature
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing 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
-
- 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
- E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- Theoretical Computer Science (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提出了一种基于温度测量的水平井吸水剖面反演方法,包括:构建吸水剖面反演模型,包括:根据注水过程中井筒与远井地带的传热方程以及能量方程推导出适合反演注水量的反演模型;构建反演误差函数,包括:将反演模型推导出的温度梯度与传热模型半解析解推导出的温度梯度构建成误差函数;测量井壁温度并划分区间,包括:利用分布式光纤传感器测量注水施工结束后注水井的井壁温度,并根据所得的温度数据将水平井段划分为多个区间;计算每个区间的吸水量,包括:将各个坐标点的温度值代入各区间的误差函数中,求出各区间误差函数的最小值,在最小化误差函数的过程中迭代反演出各区间的注水量。
Description
技术领域
本发明涉及油气田开发技术领域,特别涉及一种基于温度测量的水平井吸水剖面反演方法。
背景技术
为保持地层压力,增加注水井的注液量,水平井已经被广泛应用于水驱作业。为了在注水作业中最大限度地提高油气采收率,注入的流体必须尽可能的与油藏中的剩余油气相接触。而注入流体的波及效率主要取决于油藏非均质性,高渗透流动通道、裂缝或不同地质相层序都会导致水驱作业中产生不同的波及效率。注入的流体会优先从相对较高渗透率的孔隙中流过,而无法与较低渗透率孔隙中的油气相接触,最终导致生产井过早见水,预期采收率降低。为了提高水平井注水作业的波及效率以及产液井的最终采收率,监测注水井的注入剖面,实行合理的分层配产配注制度是很有必要的。
目前常用于评价水平井眼吸水剖面的工具为生产测井工具(PLT)。然而,PLT在水平井中的应用仍存在一些技术性问题,例如如何将PLT工具快速准确地移动到井眼的水平段以及在水平段的前进控制问题。而关于吸水井吸水量的预测方法主要有劈分系数法、吸水剖面插值法、数值模拟法以及一些基于数据挖掘的吸水剖面预测方法。劈分系数法、吸水剖面插值法以及数值模拟法主要存在计算结果不准确,测试周期长,成本较高等问题。而基于数据挖掘的吸水剖面预测方法主要存在需要少量吸水剖面资料以及预测精度低的问题。
相较于吸水剖面的监测,井下的温度监测要容易的多,如使用分布式光纤温度传感器(DTS)。注水过程中,井筒与远井端地层之间传热的主要形式为热对流和热传导。注水过程中井筒中温度分布的影响因素主要有地层的非均质性以及注入流体的流量。因此,急需开发一种能够快速准确地反演吸水剖面的方法。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此,本发明的目的在于提出一种基于温度测量的水平井吸水剖面反演方法。
为了实现上述目的,本发明面的实施例提供一种基于温度测量的水平井吸水剖面反演方法,包括如下步骤:
步骤S1,构建吸水剖面反演模型,包括:根据注水过程中井筒与远井地带的传热方程以及能量方程推导出适合反演注水量的反演模型;
步骤S2,构建反演误差函数,包括:将所述反演模型推导出的温度梯度与传热模型半解析解推导出的温度梯度构建成误差函数,以使得所述误差函数值最小的注入量作为反演出的注入量;
步骤S3,测量井壁温度并划分区间,包括:利用分布式光纤传感器测量注水施工结束后注水井的井壁温度,并根据所得的温度数据将水平井段划分为多个区间;
步骤S4,计算每个区间的吸水量,包括:将各个坐标点的温度值代入各区间的误差函数中,求出各区间误差函数的最小值,在最小化误差函数的过程中迭代反演出各区间的注水量,以使得各区间误差函数的值最小时的注水量作为反演所得的注水量。
进一步,在所述步骤S1中,所述根据注水过程中井筒与远井地带的传热方程以及能量方程推导出适合反演注水量的反演模型,包括:建立井筒与远井地带的二维传热模型并求出其半解析解;将半解析解带入无量纲能量方程中得到适合反演注水量的能量方程。即推导出适合反演注水量的反演模型。
进一步,在所述步骤S1中,建立的传热方程、能量方程以及反演模型的推导过程,包括:
在井筒中建立柱坐标,根据能量守恒原理得到井筒与远井端地层间二维非稳态传热微分方程的无因次形式:
式(1)中,ξ=r/rw为无因次半径,rw为井筒半径;τ=t/tλ为无因次时间;tλ=ρc(rw)2/λf,ρc为饱和流体地层的有效体积热容,λf为地层的导热系数;l=λ(r)/λf,为考虑井筒附近导热系数的能的变化而引入的;v=g/g*,其中g*=2πλf/cfl;由能量守恒原理得到注水井井筒中的温度分布为:
式(2)中,Cfl为流体的比热容;G(x)为研井筒的质量流量剖面;λf为地层的导热系数;rw为井眼半径;GradTw为近井地带的径向温度梯度;Tm为井内流体平均温度;
假设均质地层(l(ξ)=1)在1<ξ<ξT(τ)区域中的传热是稳态的,ξT(τ)为热冲击半径,即在r>ξT(τ)的区域内的温度都等于未扰动地层温度。代入边界条件T(ξ=1)=Tw和T(ξ=ξT)=Tf得到传热方程(1)在此区域内的稳态解:
式(3)中,ΔT=Tf-Tw;Tf为远井端未扰动地层温度,℃;Tw为注入流体的温度,℃。
由式(3)得到径向温度梯度为:
将式(4)代入到式(2)的无量纲形式可以得到适合反演注水剖面的方程:
式(5)中,kc为考虑注入井不渗透段(套管)热交换的减少而引入的,在非渗透段kc在0.7-0.95范围内,在渗透段kc=1。
进一步,在所述步骤S2中,构建出的反演误差函数为传热方程半解析解推导出的温度梯度和能量方程推导出的温度梯度差的平方和。
进一步,在所述步骤S2中,所述构建反演误差函数,包括:
在第一个非渗透段,方程(5)为:
G1·Γ1=g*·γ(τ,0)·kc·ΔT1 (6)
在剩下的所有井段(i=2,3,……,n)重方程(5)为:
Gi·Γi=g*·γ(τ,v)·ΔTi (7)
由式(6)和式(7)可得:
构造式(4)和式(8)差的平方和就得到了反演误差函数:
对式(9)进行最小化处理,最小化过程中反演出的Gi值即为所求的注水剖面。
进一步,在所述步骤S3中,划分出的每个区间的两个端点处的温度差不为0。
进一步,在所述步骤S4中,所述最小化误差函数的过程,包括:
(4)将所求的Gi+1的值代入下个井段中依次迭代求出Gi+2,Gi+3,…,Gn的值。最终得到的G1,G2,…Gn,即为注水井各井段的注水剖面。
根据本发明实施例的基于温度测量的水平井吸水剖面反演方法,步骤简单,计算方便,节省成本,能快速反演出水平注水井的注水剖面,提高注水效率,为实现合理分层配产配注提供基础。本发明基于注水施工结束后井筒的温度测量数据反演出水平注水井的吸水剖面,对于油田合理分层配产配注具有重要的意义。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明实施例的基于温度测量的水平井吸水剖面反演方法的流程图;
图2为根据本发明实施例的四种最小化误差函数的流程图;
图3为根据本发明实施例的实施例一中反演出的各井段吸水量与PLT测量值的对比;
图4为根据本发明实施例的实施例一中反演出的各井段相对吸水量与PLT测量值的对比。
具体实施方式
下面详细描述本发明的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
如图1所示,本发明实施例的基于温度测量的水平井吸水剖面反演方法,包括如下步骤:
步骤S1,构建吸水剖面反演模型,包括:根据注水过程中井筒与远井地带的传热方程以及能量方程推导出适合反演注水量的反演模型。
在步骤S1中,根据注水过程中井筒与远井地带的传热方程以及能量方程推导出适合反演注水量的反演模型,包括:建立井筒与远井地带的二维传热模型并求出其半解析解;将半解析解带入无量纲能量方程中得到适合反演注水量的能量方程,即推导出适合反演注水量的反演模型。
需要说明的是,在步骤S1中所建立的传热方程和能量方程均为二维,且推导过程中使用的传热方程的解为半解析解。
具体的,建立的传热方程、能量方程以及反演模型的推导过程,包括:
在井筒中建立柱坐标,根据能量守恒原理得到井筒与远井端地层间二维非稳态传热微分方程的无因次形式:
式(1)中,ξ=r/rw为无因次半径,rw为井筒半径;τ=t/tλ为无因次时间;tλ=ρc(rw)2/λf,ρc为饱和流体地层的有效体积热容,λf为地层的导热系数;l=λ(r)/λf,为考虑井筒附近导热系数的能的变化而引入的;v=g/g*,其中g*=2πλf/cfl;
由能量守恒原理得到注水井井筒中的温度分布为:
式(2)中,Cfl为流体的比热容;G(x)为研井筒的质量流量剖面;λf为地层的导热系数;rw为井眼半径;GradTw为近井地带的径向温度梯度;Tm为井内流体平均温度;tinj为注入时间,h.
假设均质地层(l(ξ)=1)在1<ξ<ξT(τ)区域中的传热是稳态的,ξT(τ)为热冲击半径,即在r>ξT(τ)的区域内的温度都等于未扰动地层温度。代入边界条件T(ξ=1)=Tw和T(ξ=ξT)=Tf得到传热方程(1)在此区域内的稳态解:
式(3)中v=(Gi-Gi+1)/(Δxig*);ΔT=Tf-Tw;Tf为远井端未扰动地层温度,℃;Tw为注入流体的温度,℃。
由式(3)得到径向温度梯度为:
将式(4)代入到式(2)的无量纲形式可以得到适合反演注水剖面的方程:
式(5)中,kc为考虑注入井不渗透段(套管)热交换的减少而引入的,在非渗透段kc在0.7-0.95范围内,在渗透段kc=1。
步骤S2,构建反演误差函数,包括:将反演模型推导出的温度梯度与传热模型半解析解推导出的温度梯度构建成误差函数,以使得误差函数值最小的注入量作为反演出的注入量。
在本发明的实施例中,构建出的反演误差函数为传热方程半解析解推导出的温度梯度和能量方程推导出的温度梯度差的平方和。
具体的,构建反演误差函数,包括:
在第一个非渗透段(套管),方程(5)为:
G1·Γ1=g*·γ(τ,0)·kc·ΔT1 (6)
在剩下的所有井段(i=2,3,……,n)重方程(5)为:
Gi·Γi=g*·γ(τ,v)·ΔTi (7)
式(7)中,Γi为第i井段内沿井筒的温度梯度,Γ=dTm/dx;
由式(6)和式(7)可得:
构造式(4)和式(8)差的平方和就得到了反演误差函数:
对式(9)进行最小化处理,最小化过程中反演出的Gi值即为所求的注水剖面。
步骤S3,测量井壁温度并划分区间,包括:利用分布式光纤传感器(DTS)测量注水施工结束后注水井的井壁温度,并根据所得的温度数据将水平井段划分为多个区间。
在本发明的实施例中,所划分的井段区间长度可以相同也可以不同,但区间的两个端点处的温度差不为0,即只需两个坐标点具有明显的温差即可。
步骤S4,计算每个区间的吸水量,包括:将各个坐标点的温度值代入各区间的误差函数中,求出各区间误差函数的最小值,在最小化误差函数的过程中迭代反演出各区间的注水量,以使得各区间误差函数的值最小时的注水量作为反演所得的注水量。
在本发明的实施例中,在最小化反演误差函数的过程中,利用计算机编程将最小化过程具体化,从而提高反演效率。
具体的,最小化误差函数的过程,包括:
(4)将所求的Gi+1的值代入下个井段中依次迭代求出Gi+2,Gi+3,…,Gn的值。最终得到的G1,G2,…Gn,即为注水井各井段的注水剖面。其反演流程如图2所示。
实施例一:
某油田水平注入井进行注水作业,注入流体为水,注入时间为4小时,井径为7英寸,注入流量Gin=1.7kg/s(150m3/d)。水平段开始处的流体温度为40℃,远井端未扰动地层的温度为70℃。该井水平段处有一个28m长的非渗透带(套管)。地层参数及流体参数如表1所示:
表1某油田地层参数及注入流体参数
针对上述实施例结合基于温度测量的吸水剖面反演方法,具体步骤包括:
步骤一:构建吸水剖面反演模型,根据注水过程中井筒与远井地带的传热方程以及能量方程推导出适合反演注水量的反演模型。
根据能量守恒原理得到井筒与远井端地层间二维非稳态传热微分方程的无因次形式为:
式中,ξ=r/rw为无因次半径,rw为井筒半径;τ=t/tλ为无因次时间;tλ=ρc(rw)2/λf,ρc为饱和流体地层的有效体积热容,λf为地层的导热系数;l=λ(r)/λf,为考虑井筒附近导热系数的能的变化而引入的;v=g/g*,其中g*=2πλf/cfl;
由能量守恒原理得到注水井井筒中的温度分布为:
式中,Cfl为流体的比热容;G(x)为研井筒的质量流量剖面;λf为地层的导热系数;rw为井眼半径;GradTw为近井地带的径向温度梯度。Tm为井内流体平均温度;
假设均质地层(l(ξ)=1)在1<ξ<ξT(τ)区域中的传热是稳态的,ξT(τ)为热冲击半径,即在r>ξT(τ)的区域内的温度都等于未扰动地层温度。代入边界条件T(ξ=1)=Tw和T(ξ=ξT)=Tf得到传热方程(1)在此区域内的稳态解:
式中,ΔT=Tf-Tw;Tf为远井端未扰动地层温度,℃;Tw为注入流体的温度,℃。
由式(3)得到径向温度梯度为:
将式(4)代入到式(2)的无量纲形式可以得到适合反演注水剖面的方程:
式中,kc为考虑注入井不渗透段(套管)热交换的减少而引入的,在非渗透段kc在0.7-0.95范围内,在渗透段kc=1。
步骤二:构建反演误差函数,将能量方程推导出的温度梯度与传热方程推导出的温度梯度构建成误差函数。
在第一个非渗透段(套管),方程(5)为:
G1·Γ1=g*·γ(τ,0)·kc·ΔT1 (6)
在剩下的所有井段(i=2,3,……,n)中方程(5)为:
Gi·Γi=g*·γ(τ,v)·ΔTi (7)
式中,Γi为第i井段内沿井筒的温度梯度,Γ=dTm/dx;。
由式(6)和式(7)可得:
构造式(4)和式(8)差的平方和就得到了反演误差函数:
对式(9)进行最小化处理,最小化过程中反演出的Gi值即为所求的注水剖面。
步骤三:利用分布式光纤传感器(DTS)测量注水施工结束后注水井的井壁温度,并根据所得的温度数据将水平井段划分为多个区间。
本实例中的注水井注水施工结束后水平段井筒中的温度分布由分布式光纤测温系统(DTS)测出,如表2所示:
表2井筒温度测量数据
根据测量的水平段井筒温度值将水平段划分为7个不同长度的井段,如表3所示。
表3水平段分段信息
步骤四:将各个坐标点的温度值代入各区间的误差函数中,在最小化误差函数的过程中迭代反演出各区间的注水量。
由生产测井工具测得各井段的相对吸水量,用以验证本发明所提方法的准确性。各井段的地层参数及相对吸水量如表4所示:
表4 PLT测量的各井段相对吸水量
反演出的吸水量与PLT测量的吸水量的对比图如图3和图4所示。
根据本发明实施例的基于温度测量的水平井吸水剖面反演方法,步骤简单,计算方便,节省成本,能快速反演出水平注水井的注水剖面,提高注水效率,为实现合理分层配产配注提供基础。本发明基于注水施工结束后井筒的温度测量数据反演出水平注水井的吸水剖面,对于油田合理分层配产配注具有重要的意义。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。
Claims (6)
1.一种基于温度测量的水平井吸水剖面反演方法,其特征在于,包括如下步骤:
步骤S1,构建吸水剖面反演模型,包括:根据注水过程中井筒与远井地带的传热方程以及能量方程推导出适合反演注水量的反演模型;
步骤S2,构建反演误差函数,包括:将所述反演模型推导出的温度梯度与传热模型半解析解推导出的温度梯度构建成误差函数,以使得所述误差函数值最小的注入量作为反演出的注入量;
步骤S3,测量井壁温度并划分区间,包括:利用分布式光纤传感器测量注水施工结束后注水井的井壁温度,并根据所得的温度数据将水平井段划分为多个区间;
步骤S4,计算每个区间的吸水量,包括:将各个坐标点的温度值代入各区间的误差函数中,求出各区间误差函数的最小值,在最小化误差函数的过程中迭代反演出各区间的注水量,以使得各区间误差函数的值最小时的注水量作为反演所得的注水量;
其中,所述最小化误差函数的过程,包括:
(4)将所求的Gi+1的值代入下个井段中依次迭代求出Gi+2,Gi+3,…,Gn的值。最终得到的G1,G2,…Gn,即为注水井各井段的注水剖面。
2.如权利要求1所述的基于温度测量的水平井吸水剖面反演方法,其特征在于,在所述步骤S1中,所述根据注水过程中井筒与远井地带的传热方程以及能量方程推导出适合反演注水量的反演模型,包括:建立井筒与远井地带的二维传热模型并求出其半解析解;将半解析解带入无量纲能量方程中得到适合反演注水量的能量方程。即推导出适合反演注水量的反演模型。
3.如权利要求1或2所述的基于温度测量的水平井吸水剖面反演方法,其特征在于,在所述步骤S1中,建立的传热方程、能量方程以及反演模型的推导过程,包括:
在井筒中建立柱坐标,根据能量守恒原理得到井筒与远井端地层间二维非稳态传热微分方程的无因次形式:
式(1)中,ξ=r/rw为无因次半径,rw为井筒半径;τ=t/tλ为无因次时间;tλ=ρc(rw)2/λf,ρc为饱和流体地层的有效体积热容,λf为地层的导热系数;l=λ(r)/λf,为考虑井筒附近导热系数的能的变化而引入的;v=g/g*,其中g*=2πλf/cfl;
由能量守恒原理得到注水井井筒中的温度分布为:
式(2)中,Cfl为流体的比热容;G(x)为研井筒的质量流量剖面;λf为地层的导热系数;rw为井眼半径;GradTw为近井地带的径向温度梯度;Tm为井内流体平均温度;
假设均质地层(l(ξ)=1)在1<ξ<ξT(τ)区域中的传热是稳态的,ξT(τ)为热冲击半径,即在r>ξT(τ)的区域内的温度都等于未扰动地层温度。代入边界条件T(ξ=1)=Tw和T(ξ=ξT)=Tf得到传热方程(1)在此区域内的稳态解:
式(3)中,ΔT=Tf-Tw;Tf为远井端未扰动地层温度,℃;Tw为注入流体的温度,℃。
由式(3)得到径向温度梯度为:
将式(4)代入到式(2)的无量纲形式可以得到适合反演注水剖面的方程:
式(5)中,kc为考虑注入井不渗透段(套管)热交换的减少而引入的,在非渗透段kc在0.7-0.95范围内,在渗透段kc=1。
4.如权利要求1所述的基于温度测量的水平井吸水剖面反演方法,其特征在于,在所述步骤S2中,构建出的反演误差函数为传热方程半解析解推导出的温度梯度和能量方程推导出的温度梯度差的平方和。
6.如权利要求1所述的基于温度测量的水平井吸水剖面反演方法,其特征在于,在所述步骤S3中,划分出的每个区间的两个端点处的温度差不为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111239021.5A CN113958299B (zh) | 2021-10-25 | 2021-10-25 | 一种基于温度测量的水平井吸水剖面反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111239021.5A CN113958299B (zh) | 2021-10-25 | 2021-10-25 | 一种基于温度测量的水平井吸水剖面反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113958299A true CN113958299A (zh) | 2022-01-21 |
CN113958299B CN113958299B (zh) | 2023-12-29 |
Family
ID=79466628
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111239021.5A Active CN113958299B (zh) | 2021-10-25 | 2021-10-25 | 一种基于温度测量的水平井吸水剖面反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113958299B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984886A (zh) * | 2018-07-06 | 2018-12-11 | 中国石油大学(华东) | 一种反演多层油藏井间动态连通性的方法 |
US20190338621A1 (en) * | 2018-05-02 | 2019-11-07 | Conocophillips Company | Production logging inversion based on das/dts |
-
2021
- 2021-10-25 CN CN202111239021.5A patent/CN113958299B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190338621A1 (en) * | 2018-05-02 | 2019-11-07 | Conocophillips Company | Production logging inversion based on das/dts |
CN108984886A (zh) * | 2018-07-06 | 2018-12-11 | 中国石油大学(华东) | 一种反演多层油藏井间动态连通性的方法 |
Non-Patent Citations (2)
Title |
---|
宋红伟等: "AGA和L-M算法联合预测生产油井油水流动剖面", 石油天然气学报, vol. 33, no. 10, pages 99 - 104 * |
朱海涛等: "基于Adam优化算法的水平井流动剖面测温反演方法", 《中国石油大学学报(自然科学版)》, vol. 47, no. 02, pages 99 - 107 * |
Also Published As
Publication number | Publication date |
---|---|
CN113958299B (zh) | 2023-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7725301B2 (en) | System and method for estimating multi-phase fluid rates in a subterranean well | |
US7536905B2 (en) | System and method for determining a flow profile in a deviated injection well | |
CN105160161B (zh) | 井筒内热力参数确定方法及装置 | |
US20140365130A1 (en) | Estimating flow rates from multiple hydrocarbon reservoir layers into a production well | |
CN103867174B (zh) | 一种稠油水平井吸汽状况分析方法及系统 | |
CN115293066B (zh) | 考虑地层渗流传热效应的气井温度场计算方法 | |
US10174612B2 (en) | Method for determining a water intake profile in an injection well | |
CN111400893A (zh) | 一种套管式地埋管换热器流体温度场分析方法 | |
CN102425386B (zh) | 一种符合幂律模式的钻井液流变参数控制方法 | |
CN112302607B (zh) | 一种致密气藏压裂水平井人工裂缝参数解释方法 | |
CN102156145B (zh) | 研究岩体性质的方法和设备 | |
RU2474687C1 (ru) | Способ определения профиля притока флюидов многопластовых залежей | |
CN113958299B (zh) | 一种基于温度测量的水平井吸水剖面反演方法 | |
CN105403590B (zh) | 隔热管导热系数测试方法及装置 | |
RU2569522C1 (ru) | Способ определения давления в скважине | |
Hashish et al. | Injection profiling in horizontal wells using temperature warmback analysis | |
CN114254520B (zh) | 一种超深井油套环空液面高度的确定方法 | |
Wooley et al. | Cementing temperatures for deep-well production liners | |
Murphy | Enhanced interpretation of temperature surveys taken during injection or production | |
CN114282387A (zh) | 一种基于dts的稠油油藏注蒸汽水平井综合评价方法 | |
RU2728116C1 (ru) | Способ взаимной калибровки датчиков температуры скважинного флюида, установленных на перфорационной колонне | |
CN114673494B (zh) | 蒸汽吞吐后储层孔渗预测方法 | |
Mégel et al. | Downhole pressures derived from wellhead measurements during hydraulic experiments | |
RU2704068C1 (ru) | Способ оценки межпластовых внутриколонных перетоков в скважине | |
Lavery et al. | Determining Produced Fluid Properties for Accurate Production Profiling During a Drill Stem Test Using Thermal Imaging Technology. |
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 |