CN112302607A - 一种致密气藏压裂水平井人工裂缝参数解释方法 - Google Patents

一种致密气藏压裂水平井人工裂缝参数解释方法 Download PDF

Info

Publication number
CN112302607A
CN112302607A CN202010643590.5A CN202010643590A CN112302607A CN 112302607 A CN112302607 A CN 112302607A CN 202010643590 A CN202010643590 A CN 202010643590A CN 112302607 A CN112302607 A CN 112302607A
Authority
CN
China
Prior art keywords
inversion
iteration
temperature profile
temperature
horizontal well
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
Application number
CN202010643590.5A
Other languages
English (en)
Other versions
CN112302607B (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 CN202010643590.5A priority Critical patent/CN112302607B/zh
Publication of CN112302607A publication Critical patent/CN112302607A/zh
Application granted granted Critical
Publication of CN112302607B publication Critical patent/CN112302607B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25Methods for stimulating production
    • E21B43/26Methods for stimulating production by forming crevices or fractures
    • 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
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种致密气藏压裂水平井人工裂缝参数解释方法,包括以下步骤:初始赋值人工裂缝参数向量,将初始迭代步数设为k=0,将初始阻尼因子设为ξ0=0.001;计算当前迭代步的井筒温度剖面
Figure RE-DDA0002594020090000011
第k步反演迭代的井筒温度剖面反演误差
Figure RE-DDA0002594020090000012
第k步反演迭代的温度雅阁比矩阵
Figure RE-DDA0002594020090000013
温度对角矩阵Ωk、当前迭代步
Figure RE-DDA0002594020090000014
的增量
Figure RE-DDA0002594020090000015
计算第k+1迭代步的反演裂缝参数向量
Figure RE-DDA0002594020090000016
第k+1步反演迭代的井筒温度剖面反演误差
Figure RE-DDA0002594020090000017
判断
Figure RE-DDA0002594020090000018
是否满足反演迭代终止条件之一,若满足,则停止迭代;否则令k=k+1,继续迭代,直至满足反演迭代终止条件,输出当前
Figure RE-DDA0002594020090000019
为人工裂缝参数反演解释结果。本发明通过对实测的井筒温度剖面数据进行反演,可以对致密气藏压裂水平井各级人工裂缝参数进行定量解释。

Description

一种致密气藏压裂水平井人工裂缝参数解释方法
技术领域
本发明涉及一种致密气藏压裂水平井人工裂缝参数解释方法,属于油气开采领域。
背景技术
致密气藏作为一种重要的天然气资源,我国拥有非常丰富的储量。但由于致密气藏储层低渗致密且孔隙结构复杂,使得我国致密气产量普遍较低。为了提高致密气藏产量和开发效益,目前主要通过压裂水平井进行致密气藏开采,所以,水力压裂形成的人工裂缝对压裂改造有效性至关重要,也直接决定了致密气藏压裂水平井的产能。为了保证致密气藏水平井压裂改造的有效性,必须对压裂改造效果进行准确评价,因此人工裂缝参数(裂缝半长、导流能力等) 定量诊断显得尤为重要。
目前压裂水平井人工裂缝诊断测试技术主要有近井监测技术(如:示踪剂、温生产测井等)和远程监测技术(如:微地震等),近井监测技术主要用于识别人工裂缝流入流体类型,可提供部分的近井带裂缝信息,但不能从中获知人工裂缝的具体几何尺寸参数;远程监测技术的主要优势在于监测水力压裂压力响应的总体范围,对水力压裂改造的总体范围有直观的判断,但不足之处在于难以获知有效支撑裂缝的特征参数,而且远程监测的裂缝延伸范围远比实际的有效人工裂缝控制范围更大,所以现有的常规测试技术很难直接测出压裂水平井每一条人工裂缝的具体参数。
近年来,随着水平井井筒温度剖面测试技术的不断推广应用,尤其是分布式光纤测温(Distributed Temperature Sensing,DTS)技术的飞速发展,目前可以对压裂水平井井全井段的温度剖面进行实时监测,提供实时连续的压裂水平井井筒温度剖面数据。从实测井筒温度剖面上,可直观地识别和定位压裂改造形成的有效人工裂缝,研究表明气藏压裂水平井井筒温度剖面在有效人工裂缝位置处都存在着明显温度降,且人工裂缝位置处的温度降与裂缝流量和人工裂缝参数之间存在一定的正相关关系,因此,通过数学算法建立反演模型,对实测的井筒温度剖面数据进行翻译,找到人工裂缝参数与井筒温度剖面间的对应关系并实现量化评价,就可以对每一条人工裂缝参数进行定量解释。然而,目前国内根据实测井筒温度剖面数据定量解释致密气藏压裂水平井人工裂缝参数的研究基本空白。
因此,非常有必要建立一套致密气藏压裂水平井人工裂缝参数解释模型及方法,用以定量解释致密气藏压裂水平井各级人工裂缝参数,为致密气藏压裂水平井人工裂缝参数定量诊断和压裂改造效果定量评价提供一种新的方法,从而促进我国致密气藏高效经济开发。
发明内容
本发明主要是克服现有技术中的不足之处,提出一种致密气藏压裂水平井人工裂缝参数解释方法。
本发明解决上述技术问题所提供的技术方案是:一种致密气藏压裂水平井人工裂缝参数解释方法,包括以下步骤:
S1、根据目标压裂水平井的实测井筒温度剖面
Figure RE-GDA0002594020080000021
确定井筒温度剖面在各裂缝位置处的温度降分布
Figure RE-GDA0002594020080000022
根据
Figure RE-GDA0002594020080000023
初始赋值人工裂缝参数向量为
Figure RE-GDA0002594020080000024
将初始迭代步数设为k=0,将初始阻尼因子设为ξ0=0.001;
S2、在第k步反演迭代时,将当前反演的人工裂缝参数向量
Figure RE-GDA0002594020080000025
输入井筒温度剖面预测模型计算当前迭代步的井筒温度剖面
Figure RE-GDA0002594020080000026
S3、通过误差函数方程计算第k步反演迭代的井筒温度剖面反演误差
Figure RE-GDA0002594020080000031
S4、计算第k步反演迭代的温度雅阁比矩阵
Figure RE-GDA0002594020080000032
S5、根据温度雅阁比矩阵
Figure RE-GDA0002594020080000033
计算温度对角矩阵Ωk
S6、根据温度对角矩阵Ωk计算当前迭代步
Figure RE-GDA0002594020080000034
的增量
Figure RE-GDA0002594020080000035
S7、根据当前迭代步
Figure RE-GDA0002594020080000036
的增量
Figure RE-GDA0002594020080000037
计算第k+1迭代步的反演裂缝参数向量
Figure RE-GDA0002594020080000038
S8、将第k+1迭代步的反演裂缝参数向量
Figure RE-GDA0002594020080000039
输入井筒温度剖面预测模型计算第k+1迭代步的井筒温度剖面
Figure RE-GDA00025940200800000310
并采用所述的误差函数方程计算第k+1 步反演迭代的井筒温度剖面反演误差
Figure RE-GDA00025940200800000311
S9、再将步骤S3中的第k步反演迭代的井筒温度剖面反演误差
Figure RE-GDA00025940200800000312
与步骤S8中的第k+1步反演迭代的井筒温度剖面反演误差
Figure RE-GDA00025940200800000313
进行比较;若
Figure RE-GDA00025940200800000314
令ξk+1=10ξk,否则,令ξk+1=0.1ξk,然后继续下一步;
S10、判断
Figure RE-GDA00025940200800000315
是否满足反演迭代终止条件之一,若满足,则停止迭代;否则令k=k+1,然后转步骤S2,直至满足反演迭代终止条件,反演任务完成,输出当前
Figure RE-GDA00025940200800000316
为目标压裂水平井的人工裂缝参数反演解释结果。
进一步的技术方案是,所述步骤S3和S8中的误差方程如下:
Figure RE-GDA00025940200800000317
式中,
Figure RE-GDA00025940200800000318
为第k步井筒温度剖面反演误差,
Figure RE-GDA00025940200800000319
为实测的井筒温度剖面,
Figure RE-GDA00025940200800000320
为第k步反演迭代中预测的井筒温度剖面,上标'表示向量转置。
进一步的技术方案是,所述步骤S4中通过如下方程计算:
Figure RE-GDA0002594020080000041
Figure RE-GDA0002594020080000042
式中,
Figure RE-GDA0002594020080000043
表示温度雅阁比矩阵,T为井筒温度剖面计算值,ej为单位向量中的第j分量,N为目标压裂水平井的有效人工裂缝数量,δ表示微小波动变量,i=1,2…N,j=1,2…N。
进一步的技术方案是,所述步骤S5中通过如下方程计算:
Figure RE-GDA0002594020080000044
式中,Ωk表示温度对角矩阵,diag表示矩阵的对角矩阵运算,
Figure RE-GDA0002594020080000045
表示温度雅阁比矩阵。
进一步的技术方案是,所述步骤S6中通过如下方程计算:
Figure RE-GDA0002594020080000046
式中,
Figure RE-GDA0002594020080000047
为当前迭代步
Figure RE-GDA0002594020080000048
的增量,
Figure RE-GDA0002594020080000049
为误差项,ξk表示阻尼因子,Ωk表示温度对角矩阵,
Figure RE-GDA00025940200800000410
表示温度雅阁比矩阵。
进一步的技术方案是,所述步骤S7中通过如下方程计算:
Figure RE-GDA00025940200800000411
式中,
Figure RE-GDA00025940200800000412
为当前迭代步
Figure RE-GDA00025940200800000413
的增量,
Figure RE-GDA00025940200800000414
为第k+1迭代步的反演裂缝参数向量,
Figure RE-GDA00025940200800000415
为当前迭代步。
进一步的技术方案是,所述步骤S10中的反演迭代终止条件包括:
(1)
Figure RE-GDA0002594020080000051
表示井筒温度剖面反演误差满足温度误差精度εT
(2)
Figure RE-GDA0002594020080000052
表示相邻两步反演迭代获得人工裂缝参数向量之间的差异足够小,满则人工裂缝参数迭代误差精度εFrac
本发明具有以下有益效果:
1)本发明通过对实测的井筒温度剖面数据进行反演,可以对致密气藏压裂水平井各级人工裂缝参数进行定量解释;
2)采用常规测试手段很难直接测出压裂水平井每一条人工裂缝的具体参数,本发明提供理了用于定量解释致密气藏压裂水平井每一条人工裂缝参数的解释模型和方法,可以帮助本领域技术人员对压裂改造效果进行准确评价,确保压裂改造的有效性,从而促进我国致密气藏高效经济高发。
附图说明
图1为致密气藏压裂水平井人工裂缝参数反演解释流程示意图;
图2为致密气藏压裂水平井实测井筒温度剖面示意图;
图3为实测井筒温度剖面在各级人工裂缝位置处的温度降分布示意图;
图4为初始化的人工裂缝半长分布示意图;
图5为反演模拟的井筒温度剖面与实测井筒温度剖面拟合示意图;
图6为反演解释出的人工裂缝半长分布示意图。
具体实施方式
下面将结合实施例,对本发明的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明的一种致密气藏压裂水平井人工裂缝参数解释方法,以致密气藏压裂水平井作为目标压裂水平井,以人工裂缝半长作为待解释的人工裂缝参数为例,阐述采用所述的方法进行致密气藏压裂水平井人工参数定量解释的具体步骤:
(1)根据如图2所示的目标压裂水平井的实测井筒温度剖面
Figure RE-GDA0002594020080000061
确定井筒温度剖面在各裂缝位置处的温度降分布
Figure RE-GDA0002594020080000062
如图3,根据
Figure RE-GDA0002594020080000063
初始赋值人工裂缝参数向量为
Figure RE-GDA0002594020080000064
对应的人工裂缝半长分布如图4所示,将初始迭代步数设为 k=0,将初始阻尼因子设为ξ0=0.001;
(2)在第k步反演迭代时,将当前反演的人工裂缝参数向量
Figure RE-GDA0002594020080000065
输入井筒温度剖面预测模型计算当前迭代步的井筒温度剖面
Figure RE-GDA0002594020080000066
(3)通过如下的误差函数方程计算第k步反演迭代的井筒温度剖面反演误差
Figure RE-GDA0002594020080000067
Figure RE-GDA0002594020080000068
式中,
Figure RE-GDA0002594020080000069
为第k步井筒温度剖面反演误差,
Figure RE-GDA00025940200800000610
为实测的井筒温度剖面,
Figure RE-GDA00025940200800000611
为第k步反演迭代中预测的井筒温度剖面,上标'表示向量转置;
(4)通过如下的方程计算第k步反演迭代的温度雅阁比矩阵
Figure RE-GDA00025940200800000612
Figure RE-GDA00025940200800000613
Figure RE-GDA00025940200800000614
式中,
Figure RE-GDA00025940200800000615
表示温度雅阁比矩阵,T为井筒温度剖面计算值,ej为单位向量中的第j分量,N为目标压裂水平井的有效人工裂缝数量,δ表示微小波动变量,i=1,2…N,j=1,2…N;
(5)根据温度雅阁比矩阵
Figure RE-GDA0002594020080000071
和如下的方程计算温度对角矩阵Ωk
Figure RE-GDA0002594020080000072
式中,Ωk表示温度对角矩阵,diag表示矩阵的对角矩阵运算,
Figure RE-GDA0002594020080000073
表示温度雅阁比矩阵;
(6)根据温度对角矩阵Ωk和如下的方程计算当前迭代步
Figure RE-GDA0002594020080000074
的增量
Figure RE-GDA0002594020080000075
Figure RE-GDA0002594020080000076
式中,
Figure RE-GDA0002594020080000077
为当前迭代步
Figure RE-GDA0002594020080000078
的增量,
Figure RE-GDA0002594020080000079
为误差项,ξk表示阻尼因子,Ωk表示温度对角矩阵,
Figure RE-GDA00025940200800000710
表示温度雅阁比矩阵;
(7)根据当前迭代步
Figure RE-GDA00025940200800000711
的增量
Figure RE-GDA00025940200800000712
和如下的方程计算第k+1迭代步的反演裂缝参数向量
Figure RE-GDA00025940200800000713
Figure RE-GDA00025940200800000714
式中,
Figure RE-GDA00025940200800000715
为当前迭代步
Figure RE-GDA00025940200800000716
的增量,
Figure RE-GDA00025940200800000717
为第k+1迭代步的反演裂缝参数向量,
Figure RE-GDA00025940200800000718
为当前迭代步;
(8)将第k+1迭代步的反演裂缝参数向量
Figure RE-GDA00025940200800000719
输入井筒温度剖面预测模型计算第k+1迭代步的井筒温度剖面
Figure RE-GDA00025940200800000720
并采用所述的误差函数方程计算第k+1 步反演迭代的井筒温度剖面反演误差
Figure RE-GDA00025940200800000721
(9)再将步骤(3)中的第k步反演迭代的井筒温度剖面反演误差
Figure RE-GDA00025940200800000722
与步骤S8中的第k+1步反演迭代的井筒温度剖面反演误差
Figure RE-GDA00025940200800000723
进行比较;若
Figure RE-GDA00025940200800000724
令ξk+1=10ξk,否则,令ξk+1=0.1ξk,然后继续下一步;
(10)判断
Figure RE-GDA00025940200800000725
是否满足反演迭代终止条件之一,若满足,表明根据
Figure RE-GDA00025940200800000726
模拟的井筒温度剖面与实测井筒温度剖面相拟合,如图5所示,停止迭代;否则令k=k+1,然后转步骤(2),直至满足反演迭代终止条件,反演任务完成,输出当前
Figure RE-GDA0002594020080000081
为目标压裂水平井的人工裂缝参数反演解释结果,如图6所示;
其中反演迭代终止条件包括:
1)
Figure RE-GDA0002594020080000082
表示井筒温度剖面反演误差满足温度误差精度εT
2)
Figure RE-GDA0002594020080000083
表示相邻两步反演迭代获得人工裂缝参数向量之间的差异足够小,满则人工裂缝参数迭代误差精度εFrac
本发明所述的井筒温度剖面数据可以但不限于通过分布式光纤、拖动式生产测井工具获取。
所述的井筒温度剖面预测模型为一个综合的温度模型,包括:
致密储层渗流模型:
Figure RE-GDA0002594020080000084
致密储层热学模型:
Figure RE-GDA0002594020080000085
人工裂缝渗流模型:
Figure RE-GDA0002594020080000086
人工裂缝热学模型:
Figure RE-GDA0002594020080000091
水平井筒流动模型:
Figure RE-GDA0002594020080000092
水平井筒热学模型:
Figure RE-GDA0002594020080000093
式中:Cg为气体压缩系数,MPa-1;Cp为热容,J/(kg·K);f为井壁摩擦系数,小数;KJT为焦耳汤普逊系数,K/MPa;KT为岩石热导率,J/(m·s·K);kx为储层 x方向的渗透率,mD;ky为储层y方向的渗透率,mD;kz为储层z方向的渗透率,mD;p为储层压力,MPa;qF为裂缝中的流体流速,m/s;qwb为固井段单位体积岩石向井筒传递热量的速率,J/(m3·s);Rinw为井筒内径,m;t表示生产时间,天;T表示温度,K;TI为流体流入温度,K;vI为流入流体的流速, m/s;v为流体的流速,m/s;
Figure RE-GDA0002594020080000094
为孔隙度,小数;μg为气体粘度,mPa·s;ψ为拟压力,MPa2/mP·s;β为热膨胀系数,1/K;γ为井筒打开程度,小数;ρ为流体密度,kg/m3;ρI为流入流体的密度,kg/m3;θ为水平井筒倾角,°;下标F为人工裂缝,下标x表示在x方向,下标y表示在y方向,下标z表示在z方向。
所述的人工裂缝参数包括但不限于人工裂缝半长、人工裂缝导流能力和人工裂缝渗透率。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。

Claims (7)

1.一种致密气藏压裂水平井人工裂缝参数解释方法,其特征在于,包括以下步骤:
S1、根据目标压裂水平井的实测井筒温度剖面
Figure RE-FDA0002594020070000011
确定井筒温度剖面在各裂缝位置处的温度降分布
Figure RE-FDA0002594020070000012
根据
Figure RE-FDA0002594020070000013
初始赋值人工裂缝参数向量为
Figure RE-FDA0002594020070000014
将初始迭代步数设为k=0,将初始阻尼因子设为ξ0=0.001;
S2、在第k步反演迭代时,将当前反演的人工裂缝参数向量
Figure RE-FDA0002594020070000015
输入井筒温度剖面预测模型计算当前迭代步的井筒温度剖面
Figure RE-FDA0002594020070000016
S3、通过误差函数方程计算第k步反演迭代的井筒温度剖面反演误差
Figure RE-FDA0002594020070000017
S4、计算第k步反演迭代的温度雅阁比矩阵
Figure RE-FDA0002594020070000018
S5、根据温度雅阁比矩阵
Figure RE-FDA0002594020070000019
计算温度对角矩阵Ωk
S6、根据温度对角矩阵Ωk计算当前迭代步
Figure RE-FDA00025940200700000110
的增量
Figure RE-FDA00025940200700000111
S7、根据当前迭代步
Figure RE-FDA00025940200700000112
的增量
Figure RE-FDA00025940200700000113
计算第k+1迭代步的反演裂缝参数向量
Figure RE-FDA00025940200700000114
S8、将第k+1迭代步的反演裂缝参数向量
Figure RE-FDA00025940200700000115
输入井筒温度剖面预测模型计算第k+1迭代步的井筒温度剖面
Figure RE-FDA00025940200700000116
并采用所述的误差函数方程计算第k+1步反演迭代的井筒温度剖面反演误差
Figure RE-FDA00025940200700000117
S9、再将步骤S3中的第k步反演迭代的井筒温度剖面反演误差
Figure RE-FDA00025940200700000118
与步骤S8中的第k+1步反演迭代的井筒温度剖面反演误差
Figure RE-FDA00025940200700000119
进行比较;若
Figure RE-FDA00025940200700000120
令ξk+1=10ξk,否则,令ξk+1=0.1ξk,然后继续下一步;
S10、判断
Figure RE-FDA00025940200700000121
是否满足反演迭代终止条件之一,若满足,则停止迭代;否则令k=k+1,然后转步骤S2,直至满足反演迭代终止条件,反演任务完成,输出当前
Figure RE-FDA00025940200700000122
为目标压裂水平井的人工裂缝参数反演解释结果。
2.根据权利要求1所述的一种致密气藏压裂水平井人工裂缝参数解释方法,其特征在于,所述步骤S3和S8中的误差方程如下:
Figure RE-FDA0002594020070000021
式中,
Figure RE-FDA0002594020070000022
为第k步井筒温度剖面反演误差,
Figure RE-FDA0002594020070000023
为实测的井筒温度剖面,
Figure RE-FDA0002594020070000024
为第k步反演迭代中预测的井筒温度剖面,上标'表示向量转置。
3.根据权利要求1所述的一种致密气藏压裂水平井人工裂缝参数解释方法,其特征在于,所述步骤S4中通过如下方程计算:
Figure RE-FDA0002594020070000025
Figure RE-FDA0002594020070000026
式中,
Figure RE-FDA0002594020070000027
表示温度雅阁比矩阵,T为井筒温度剖面计算值,ej为单位向量中的第j分量,N为目标压裂水平井的有效人工裂缝数量,δ表示微小波动变量,i=1,2…N,j=1,2…N。
4.根据权利要求3所述的一种致密气藏压裂水平井人工裂缝参数解释方法,其特征在于,所述步骤S5中通过如下方程计算:
Figure RE-FDA0002594020070000028
式中,Ωk表示温度对角矩阵,diag表示矩阵的对角矩阵运算,
Figure RE-FDA0002594020070000029
表示温度雅阁比矩阵。
5.根据权利要求4所述的一种致密气藏压裂水平井人工裂缝参数解释方法,其特征在于,所述步骤S6中通过如下方程计算:
Figure RE-FDA0002594020070000031
式中,
Figure RE-FDA0002594020070000032
为当前迭代步
Figure RE-FDA0002594020070000033
的增量,
Figure RE-FDA0002594020070000034
为误差项,ξk表示阻尼因子,Ωk表示温度对角矩阵,
Figure RE-FDA0002594020070000035
表示温度雅阁比矩阵。
6.根据权利要求5所述的一种致密气藏压裂水平井人工裂缝参数解释方法,其特征在于,所述步骤S7中通过如下方程计算:
Figure RE-FDA0002594020070000036
式中,
Figure RE-FDA0002594020070000037
为当前迭代步
Figure RE-FDA0002594020070000038
的增量,
Figure RE-FDA0002594020070000039
为第k+1迭代步的反演裂缝参数向量,
Figure RE-FDA00025940200700000310
为当前迭代步。
7.根据权利要求1所述的一种致密气藏压裂水平井人工裂缝参数解释方法,其特征在于,所述步骤S10中的反演迭代终止条件包括:
(1)
Figure RE-FDA00025940200700000311
表示井筒温度剖面反演误差满足温度误差精度εT
(2)
Figure RE-FDA00025940200700000312
表示相邻两步反演迭代获得人工裂缝参数向量之间的差异足够小,满则人工裂缝参数迭代误差精度εFrac
CN202010643590.5A 2020-07-07 2020-07-07 一种致密气藏压裂水平井人工裂缝参数解释方法 Active CN112302607B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010643590.5A CN112302607B (zh) 2020-07-07 2020-07-07 一种致密气藏压裂水平井人工裂缝参数解释方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010643590.5A CN112302607B (zh) 2020-07-07 2020-07-07 一种致密气藏压裂水平井人工裂缝参数解释方法

Publications (2)

Publication Number Publication Date
CN112302607A true CN112302607A (zh) 2021-02-02
CN112302607B CN112302607B (zh) 2021-09-24

Family

ID=74483505

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010643590.5A Active CN112302607B (zh) 2020-07-07 2020-07-07 一种致密气藏压裂水平井人工裂缝参数解释方法

Country Status (1)

Country Link
CN (1) CN112302607B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114282387A (zh) * 2021-12-29 2022-04-05 西南石油大学 一种基于dts的稠油油藏注蒸汽水平井综合评价方法
CN114562254A (zh) * 2022-01-24 2022-05-31 中油奥博(成都)科技有限公司 基于微结构光纤的地下井中流体压力测量***及测量方法
CN114595504A (zh) * 2022-03-09 2022-06-07 西南石油大学 一种页岩气藏压裂水平井产出剖面反演解释方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140358444A1 (en) * 2013-05-31 2014-12-04 Conocophillips Company Method of hydraulic fracture identification using temperature
CN104358554A (zh) * 2014-12-09 2015-02-18 中国石油集团川庆钻探工程有限公司 页岩气组合加砂压裂裂缝导流能力评价方法
CN105003239A (zh) * 2014-04-15 2015-10-28 中国海洋石油总公司 海上压裂水平井压后效果评价方法
CN106368688A (zh) * 2015-07-21 2017-02-01 中国石油化工股份有限公司 研究岩石受热破裂条件和裂缝分布的实验方法
CN108319738A (zh) * 2017-01-18 2018-07-24 中国石油化工股份有限公司 一种页岩气井产量预测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140358444A1 (en) * 2013-05-31 2014-12-04 Conocophillips Company Method of hydraulic fracture identification using temperature
CN105003239A (zh) * 2014-04-15 2015-10-28 中国海洋石油总公司 海上压裂水平井压后效果评价方法
CN104358554A (zh) * 2014-12-09 2015-02-18 中国石油集团川庆钻探工程有限公司 页岩气组合加砂压裂裂缝导流能力评价方法
CN106368688A (zh) * 2015-07-21 2017-02-01 中国石油化工股份有限公司 研究岩石受热破裂条件和裂缝分布的实验方法
CN108319738A (zh) * 2017-01-18 2018-07-24 中国石油化工股份有限公司 一种页岩气井产量预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
朱世琰: "基于分布式光纤温度测试的水平井产出剖面解释理论研究", 《中国优秀博士学位论文全文数据库(电子期刊)工程科技Ⅰ辑》 *
罗红文等: "基于DTS数据反演的低渗气藏压裂水平井产出剖面解释新方法", 《天然气地球科学》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114282387A (zh) * 2021-12-29 2022-04-05 西南石油大学 一种基于dts的稠油油藏注蒸汽水平井综合评价方法
CN114562254A (zh) * 2022-01-24 2022-05-31 中油奥博(成都)科技有限公司 基于微结构光纤的地下井中流体压力测量***及测量方法
CN114595504A (zh) * 2022-03-09 2022-06-07 西南石油大学 一种页岩气藏压裂水平井产出剖面反演解释方法
CN114595504B (zh) * 2022-03-09 2023-03-10 西南石油大学 一种页岩气藏压裂水平井产出剖面反演解释方法

Also Published As

Publication number Publication date
CN112302607B (zh) 2021-09-24

Similar Documents

Publication Publication Date Title
CN112302607B (zh) 一种致密气藏压裂水平井人工裂缝参数解释方法
CN111581786B (zh) 用于分析缝洞串联模式双孔复合储层的试井解释模型的试井解释方法
CN112302606B (zh) 一种低渗气藏压裂水平井产出剖面反演解释方法
AU2002300917B2 (en) Method of predicting formation temperature
CN111927417B (zh) 一种页岩气分段压裂水平井组储量动用状况评价方法
CN110656915B (zh) 一种页岩气多段压裂水平井多工作制度产能预测方法
CN108843303B (zh) 一种基于泥岩蠕变模型的油水井套管损坏预测方法
Zhang et al. Inversion of Downhole temperature measurements in multistage-fracturing stimulation of horizontal wells in unconventional reservoirs
CN114595504B (zh) 一种页岩气藏压裂水平井产出剖面反演解释方法
CN113919240A (zh) 一种基于油气井生产数据的高压气藏储层参数的计算方法
Luo et al. A novel inversion approach for fracture parameters and inflow rates diagnosis in multistage fractured horizontal wells
CN111648764A (zh) 一种多层气藏井下分布式温度监测产出剖面解释评价方法
CN112780257A (zh) 基于分布式光纤传感钻井液漏失监测***及监测方法
CN116128083A (zh) 一种页岩油水平井体积压裂裂缝体积的定量表征方法
CN113109162B (zh) 一种基于热流固耦合的岩石起裂压力计算方法
CN111734394A (zh) 一种确定致密油藏压裂井不定常流井底压力的方法
RU2474687C1 (ru) Способ определения профиля притока флюидов многопластовых залежей
CN106934075B (zh) 钻井液密度确定方法及静态当量密度确定方法
CN110750918A (zh) 一种二氧化碳压裂过程中井筒温度的预测方法
CN115293066B (zh) 考虑地层渗流传热效应的气井温度场计算方法
CN114282387A (zh) 一种基于dts的稠油油藏注蒸汽水平井综合评价方法
CN105003238A (zh) 利用井筒压力温度剖面分析井下蒸汽干度方法
CN113065246B (zh) 一种基于射孔测试联作***的地层参数反演方法
Jiang et al. Analysis of dynamic thermal behaviors for multi-stage hydraulic fracturing treatments in horizontal shale oil and shale gas wells
US20220010672A1 (en) The method of determining a production well flow profile, including determination of hydrodynamic characteristics of reservoir pay zone

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
CB03 Change of inventor or designer information

Inventor after: Luo Hongwen

Inventor after: Li Haitao

Inventor after: Xiang Yuxing

Inventor after: Li Ying

Inventor after: Jiang Beibei

Inventor after: Cui Xiaojiang

Inventor after: Liu Chang

Inventor after: Zhu Xiaoping

Inventor after: Li Baoji

Inventor before: Li Haitao

Inventor before: Luo Hongwen

Inventor before: Li Ying

Inventor before: Jiang Beibei

Inventor before: Cui Xiaojiang

Inventor before: Liu Chang

Inventor before: Zhu Xiaoping

Inventor before: Xiang Yuxing

Inventor before: Li Baoji

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant