CN112699554B - 一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法 - Google Patents

一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法 Download PDF

Info

Publication number
CN112699554B
CN112699554B CN202011602457.1A CN202011602457A CN112699554B CN 112699554 B CN112699554 B CN 112699554B CN 202011602457 A CN202011602457 A CN 202011602457A CN 112699554 B CN112699554 B CN 112699554B
Authority
CN
China
Prior art keywords
horizontal well
fracturing
section
constraint
fracture
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
CN202011602457.1A
Other languages
English (en)
Other versions
CN112699554A (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.)
Xian Shiyou University
Original Assignee
Xian Shiyou 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 Xian Shiyou University filed Critical Xian Shiyou University
Priority to CN202011602457.1A priority Critical patent/CN112699554B/zh
Publication of CN112699554A publication Critical patent/CN112699554A/zh
Application granted granted Critical
Publication of CN112699554B publication Critical patent/CN112699554B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • 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

Abstract

本发明公开一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，包括如下步骤：在压裂过程中，获取地面返排不同示踪剂的浓度采出曲线，进而获得变产数据；利用反卷积约束下的压裂示踪剂流量解释模型将获得的水平井各段内的变产数据转化为定产数据；利用所述定产数据，确定基于压裂示踪约束的单一无限和有限导流垂直裂缝水平井试井模型，利用基于压裂示踪约束的单一无限和有限导流垂直裂缝水平井试井模型获得水平井各段的地层参数；根据所获得的水平井各段的地层相关参数，绘制油藏水平井各段的产出状况图，计算水平井各段贡献率及获取水平井整段的地层信息。本发明克服了常规单一监测手段的局限，实现了多种监测手段的联合约束解释。

Description

一种基于压裂示踪约束的致密油藏水平井压后分段试井分析 方法
技术领域
本发明涉及石油工程学领域，具体涉及的是一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法。
背景技术
我国致密油气资源具有分布广、储量多、潜力大、前景广等特点。经过众多科技人员的研究探索，水平井分段压裂等技术已成为开发此类油气田的有力工具，是提高油气田采收率、油田勘探开发综合效益的重要途径。储层在大规模改造后，准确有效获取各层段压裂裂缝半长、压裂裂缝导流能力、储层渗透率等压裂裂缝参数及储层参数，是科学制定合理生产制度和最大效益开发的重要前提。通过对压裂裂缝的监测方法进行调研，其中包括远离裂缝直接成像技术、接近井眼裂缝监测技术、间接监测技术等。这些方法与手段均可从不同角度进行监测与表征，但大多具有一定的局限性。经过深入的调研发现，压裂示踪监测技术是地层产量评价的重要手段之一，它在估算裂缝导流能力、评价分段压裂的效率等方面具有重要意义，但其不能实现对裂缝的定量表征而且获取参数较为单一；试井能获取更多参数、确定井的生产能力和研究储层参数及储层动态，然而现有的模型均假定水平井筒无限导流，且大多需要关井测压，解释结果具有较多不确定性与多解性。
发明内容
为解决现有技术中存在的问题，本发明的目的在于提供一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，本发明克服了常规单一监测手段的局限，实现了多种监测手段的联合约束解释。
本发明采用的技术方案如下：
一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，包括如下步骤：
S1，在均质箱型致密油藏压裂过程中，将含有不同的示踪剂的压裂液泵送至地层，获取地面返排示踪剂的浓度采出曲线，通过对示踪剂产出标准浓度进行劈分，获得水平井各段的产量变化史，通过水平井各段的产量变化史获得变产数据；
S2，利用反卷积处理方法，获取反卷积约束下的压裂示踪剂流量解释模型，利用反卷积约束下的压裂示踪剂流量解释模型将获得的水平井各段内的变产数据转化为定产数据；
S3，利用所述定产数据，确定基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型及基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型；利用基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型及基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型，获得水平井各段的地层参数；
S4，根据所获得的水平井各段的地层相关参数，绘制油藏水平井各段的产出状况图，计算水平井各段贡献率及获取水平井整段的地层信息。
所述S1中，对示踪剂产出标准浓度进行劈分，获得水平井各段的产量变化史的过程包括如下步骤：
S1.1，将示踪剂浓度标准化；
S1.2，利用标准化的示踪剂浓度计算水平井各段的返排总液量；
S1.3，通过水平井各段的返排总液量得到水平井各段的返排总液量的积分形式；
S1.4，利用水平井各段的返排总液量的积分形式获得水平井各段的产量变化史。
所述水平井各段的产量变化史包括水平井各段示踪剂累加标准浓度的时间变化史、研究段示踪剂标准浓度的时间变化史和地面总返排量的时间变化史。
所述S1.1中，将示踪剂浓度标准化时，水平井第j段第i个取样点的示踪剂标准浓度ρji为：
Figure BDA0002869186180000031
S1.2中，水平井第j段的返排总液量Qj为：
Figure BDA0002869186180000032
S1.3中，水平井第j段的返排总液量Qj的积分形式为：
Figure BDA0002869186180000033
式中，M为压裂分段数与压裂示踪剂种类，C0j为注入第j段的示踪剂初始浓度，Cji为第j段第i个取样点的示踪剂浓度，n为整个示踪监测取样个数，ρ(t)为各段累加标准浓度随时间的变化函数；ρj(t)为第j段标准浓度随时间的变化函数；Q(t)为地面总返排量随时间变化的函数，t为时间。
所述S2中，反卷积约束下的压裂示踪剂流量解释模型如下：
Figure BDA0002869186180000034
式中，Δp(t)为第j个裂缝压降，MPa；σ为时间的对数函数，σ＝lnΔt；z(σ)为反卷积压力响应导数，
Figure BDA0002869186180000035
Δpu为重整压力；t为生产时间，d。
所述基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型P1(t)如下：
Figure BDA0002869186180000036
其中，yeD为储层纵向边界无量纲位置，
Figure BDA0002869186180000037
ye为储层纵向边界位置，m；Lf为参考长度，m；xeD为储层横向边界无量纲位置，
Figure BDA0002869186180000041
xe为储层纵向边界位置，m；yw1D为第1个裂缝的无量纲纵坐标，
Figure BDA0002869186180000042
yw1为第1条裂缝纵坐标，m；xw1D为第1个裂缝的无量纲横坐标，
Figure BDA0002869186180000043
xw1为第1条裂缝横坐标，m；yD为任意位置的纵坐标无量纲位置，
Figure BDA0002869186180000044
xD为任意位置的横坐标无量纲位置，
Figure BDA0002869186180000045
xfD为裂缝无量纲长度，
Figure BDA0002869186180000046
Figure BDA0002869186180000047
u为Laplace变量，上标“～”表示Laplace变换量；K为地层渗透率，μm2；μ；流体粘度，mPa·s；h为储层高度，m；qf为裂缝产出液体流量，m3/d；B为流体体积系数，m3/m3。
所述基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型P2(t)如下：
Figure BDA0002869186180000048
其中，yeD为储层纵向边界无量纲位置，
Figure BDA0002869186180000049
ye为储层纵向边界位置，m；Lf为参考长度，m；xeD为储层横向边界无量纲位置，
Figure BDA00028691861800000410
xe为储层纵向边界位置，m；yw1D为第1个裂缝的无量纲纵坐标，
Figure BDA00028691861800000411
yw1为第1条裂缝纵坐标，m；xw1D为第1个裂缝的无量纲横坐标，
Figure BDA00028691861800000412
xw1为第1条裂缝横坐标，m；yD为任意位置的纵坐标无量纲位置，
Figure BDA00028691861800000413
xD为任意位置的横坐标无量纲位置，
Figure BDA00028691861800000414
xfD为裂缝无量纲长度，
Figure BDA00028691861800000415
Figure BDA00028691861800000416
u为Laplace变量，上标“～”表示Laplace变换量；K为地层渗透率，μm2；μ；流体粘度，mPa·s；h为储层高度，m；qf为裂缝产出液体流量，m3/d；B为流体体积系数，m3/m3。Kfw为裂缝导流系数，μm2·cm；FCD为无量纲导流系数，
Figure BDA0002869186180000051
获得水平井各段的地层参数时，基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型及基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型满足如下约束：
反卷积约束下的压裂示踪剂流量解释模型和实测的压力数据之间的标准偏差最小；
反卷积压力响应导数随时间的对数函数曲线拥有最小的总曲率值；
校正输出时，在相同质量的拟合条件下，选择输出要求最小的输出。
所述水平井各段的地层相关参数包括渗透率、裂缝半长和导流系数。
本发明具有如下有益效果：
本发明在压裂示踪分段流量解释模型的基础上，将均质箱型油藏中的分段压裂水平井试井模型转换为受不同段示踪剂返排浓度约束的一系列无限或有限导流单裂缝的试井模型，并利用分段产量反卷积处理方法，将各压裂段运用定产的压降模型解释图版拟合，即实现了对多段压裂水平井缝间无干扰进行解释，获得各压裂段的裂缝及地层的相关参数。本发明进一步降低了分段压裂水平井试井解释的多解性，提高所得参数的可靠性，可为分段压裂水平井压后评价提供一种新途径，克服了常规单一监测手段的局限，实现了多种监测手段的联合约束解释。本发明寻求压裂示踪与压后试井在联合解释上的协同优势，扬长避短，相互约束控制以减小解释偏差，从而获取水平井裂缝及储层更有效的多种可靠信息，降低了监测成本。
附图说明
图1为本发明基于压裂示踪约束的致密油藏水平井压后分段试井分析方法的流程图；
图2(a)为本发明实施例中压裂第一段水平井试井拟合图；
图2(b)为本发明实施例中压裂第二段水平井试井拟合图；
图3为本发明实施例获得的压裂水平井各段产液量日贡献率示意图。
具体实施方式
下面结合附图和列举实施例的方式对本发明的技术方案进一步做出具体的阐述，但本发明不限于以下所列举的实施例。
本发明应用压裂示踪监测分段流量评价与试井两者的优势，提出了一种基于压裂示踪约束的水平井压后分段试井分析方法。在压裂示踪分段流量解释模型的基础上，利用分段产量反卷积处理技术将均质、箱型油藏中的压裂水平井试井模型转换为受不同段示踪剂返排浓度约束的一系列无限及有限导流单裂缝的试井模型，即实现了对压裂水平井的分段试井，获得各压裂段的裂缝及地层的相关参数，从而获得整个水平井段的地层信息。本发明的具体方案如下：
参照图1，本发明基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，包括以下步骤：
步骤一：在均质箱型致密油藏压裂过程中，将含有不同的示踪剂的压裂液(且各示踪剂间互不影响)泵送至地层，获取地面所返排示踪剂的浓度采出曲线，通过对示踪剂产出标准浓度进行劈分，获得水平井各段的产量变化史，如图1中的中步骤11。
步骤二：利用反卷积处理方法，推导出反卷积约束下的压裂示踪剂流量解释模型，利用反卷积约束下的压裂示踪剂流量解释模型将获得到的水平井各段内的变产问题转化为定产问题，即水平井各段内的变产数据转化为定产数据，如图1中步骤12，为后续进行约束压裂水平井试井解释提供定产基础。
步骤三：将压裂水平井各段的变量转化为定产后，利用所述定产数据，推导出基于压裂示踪约束下的压裂水平井单一垂直裂缝试井解释公式，即利用均质箱型油藏中的基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型及基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型进行解释，在压力传导慢的致密油藏中要求生产初期曲线与理论曲线拟合较好，通过拟合进一步可获得地层各段间渗透率、裂缝半长和导流系数这些地层参数，如图1中的步骤13。
步骤四：整理所获得的水平井各段的地层相关参数，绘制该油藏各段的产出状况图，计算水平井各段贡献率及了解水平井整段的地层信息。
上述步骤一的具体方法如下：
将均质箱型致密油藏水平井进行压裂后，采用封隔器封隔不同的压裂段，将类型不同且互不影响的示踪剂添加到压裂液中，泵送至每个水平井的压裂段中；之后，进行关井，以确保示踪剂能够完全渗入地下，从而使得获得的信息更加接近地层；然后，在井口进行示踪剂采样和监测；最后，对监测到的数据进行整理，将获得各段示踪剂累加标准浓度的时间变化史、研究段示踪剂标准浓度的时间变化史和地面总返排量的时间变化史采用示踪剂标准浓度劈分处理技术获得研究段的产量变化史。
示踪剂标准浓度劈分处理的过程包括如下步骤：
1、示踪剂浓度的标准化
由于不同压裂段示踪监测矿场实施具有差异性，为避免引起计算与分析误差，示踪剂浓度统一采用标准化浓度，则第j段第i个取样点的示踪剂标准浓度为：
Figure BDA0002869186180000071
2、分段返排液量的确定
设第i个取样点时对应的压裂液返排量为Qi，则有：
Figure BDA0002869186180000072
第j段第i个取样点对应的压裂液返排量可用各压裂段示踪剂产出标准浓度劈分，则有：
Figure BDA0002869186180000081
设整个示踪监测取样个数为n，则整个监测过程中第j段压裂液返排总液量为：
Figure BDA0002869186180000082
积分形式的确定
构造函数：各段累加标准浓度随时间的变化函数为ρ(t)；第j段标准浓度随时间的变化函数为ρj(t)；地面总返排量随时间变化的函数为Q(t)。
则式(4)用积分形式表示为：
Figure BDA0002869186180000083
其中，压裂分段数与压裂示踪剂种类均为M，注入第j段的示踪剂初始浓度为C0j，第j段第i个取样点的示踪剂浓度为Cji，对应该取样点返排量为qji，监测过程中第j段压裂液累计返排量为Qj。
所述步骤二中在采用劈分处理方法，获得水平井各段的产量变化史后，在通过反卷积方法推导出基于反卷积约束下的压裂示踪剂流量解释模型，如式(6)所示，通过式(6)将变产数据转化为所需的定产数据。
Figure BDA0002869186180000084
式中，Δp(t):第j个裂缝压降，MPa；σ为时间的对数函数，σ＝lnΔt；z(σ)为反卷积压力响应导数，
Figure BDA0002869186180000085
Δpu为重整压力；t为生产时间，d。
所述步骤三的具体方法如下：
进一步根据致密储层中压力传导慢、水平井处于生产初期时各裂缝间互不干扰的特点，在压裂示踪剂流量解释模型和成熟的压裂水平井试井解释方法的基础上，推导出压裂示踪约束下的压裂水平井单一垂直裂缝试井解释公式，使用试井软件进行拟合，压裂示踪约束下的压裂水平井单一垂直裂缝试井解释公式包括基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型及基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型，具体如下：
基于压裂示踪约束的单一无限导流垂直裂缝压裂水平井试井分析模型如下：
Figure BDA0002869186180000091
基于压裂示踪约束的单一有限导流垂直裂缝压裂水平井试井分析模型：
Figure BDA0002869186180000092
其中，yeD为储层纵向边界无量纲位置，
Figure BDA0002869186180000093
ye为储层纵向边界位置，m；Lf为参考长度，m；xeD为储层横向边界无量纲位置，
Figure BDA0002869186180000094
xe为储层纵向边界位置，m；yw1D为第1个裂缝的无量纲纵坐标，
Figure BDA0002869186180000095
yw1为第1条裂缝纵坐标，m；xw1D为第1个裂缝的无量纲横坐标，
Figure BDA0002869186180000096
xw1为第1条裂缝横坐标，m；yD为任意位置的纵坐标无量纲位置，
Figure BDA0002869186180000097
xD为任意位置的横坐标无量纲位置，
Figure BDA0002869186180000098
xfD为裂缝无量纲长度，
Figure BDA0002869186180000099
Figure BDA00028691861800000910
u为Laplace变量，上标“～”表示Laplace变换量；K为地层渗透率，μm2；μ；流体粘度，mPa·s；h为储层高度，m；qf为裂缝产出液体流量，m3/d；B为流体体积系数，m3/m3。Kfw为裂缝导流系数，μm2·cm；FCD为无量纲导流系数，
Figure BDA0002869186180000101
本发上述方法中，采用saphir试井软件进行拟合，在将压力史和产量史代入软件拟合过程中，需要作出以下限制：①确保卷积模型和测得的压力数据之间的标准偏差最小；②确保在整个拟合过程中反卷积压力响应导数随时间的对数函数曲线拥有最小的总曲率值；③校正输出时，请确保输出的调整尽可能小，即在相同质量的拟合条件下，选择输出要求最小的输出作为解决问题的方法。
实施例
基于本发明的上述方法，针对某均质箱型致密油藏水平井压后分段试井分析，该均质箱型致密油藏多段压裂水平井，水平井长度为800m，井半径为0.0878m，油层厚度为16.15m，孔隙度为0.15，原油黏度为1.6mPa·s，体积系数为1.2，综合压缩系数为MPa-1，共生产4000h，其中关井测试时间600h。该水平井进行了多段压裂，形成了2条裂缝，其裂缝间距为100m。
根据步骤一可获得该油藏各研究段的产量变化史，将得到的产量变化史运用试井软件中的反卷积功能处理成定产问题，再选择箱型均质的油藏形态以及地层条件按步骤三进行拟合，在拟合过程中注意受到限制的约束，油藏生产数据如表1所示，拟合结果图如图2(a)和图2(b)所示。拟合结果显示：第一压裂段在生产初期试井曲线与理论双对数拟合很好，第二压裂段与理论双对数曲线总体拟合很好，保证了该方法的可靠性。
表1
Figure BDA0002869186180000102
Figure BDA0002869186180000111
根据由试井软件获得的拟合结果可得到相关地层参数，如表2所示，进一步可得到该油藏各段产出状况表及各段产液量日贡献率图，如表3和图3所示。表2显示第一压裂段的渗透率为0.07mD、裂缝半长为0.8321m、导流系数Fc为20.22μm2·cm；第二压裂段的渗透率为1.164mD、裂缝半长为67.9704m、导流系数Fc为424.18μm2·cm。表3显示出第一段为次产层，第二段为主产层。
表2
Figure BDA0002869186180000121
表3该油藏各段产出状况表
Figure BDA0002869186180000122
本实例解释分析结果表明：本发明的分析方法分段试井曲线拟合较好，进一步降低了分段压裂水平井试井解释的多解性，提高所得参数的可靠性。
本发明研究成果可为分段压裂水平井压后评价提供一种新途径，克服了常规单一监测手段的局限，实现了多种监测手段的联合约束解释。
以上所述，仅是本发明的部分较佳实施例，任何熟悉本领域的研究人员均可能利用上述技术方案加以修改。因此，依据本发明的技术方案所进行的任何简单修改或等同置换，属于本发明要求保护的范围。

Claims (9)

1.一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，其特征在于，包括如下步骤：
S1，在均质箱型致密油藏压裂过程中，将含有不同的示踪剂的压裂液泵送至地层，获取地面返排示踪剂的浓度采出曲线，通过对示踪剂产出标准浓度进行劈分，获得水平井各段的产量变化史，通过水平井各段的产量变化史获得变产数据；
S2，利用反卷积处理方法，获取反卷积约束下的压裂示踪剂流量解释模型；利用反卷积约束下的压裂示踪剂流量解释模型将获得的水平井各段内的变产数据转化为定产数据；
S3，利用所述定产数据，确定基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型及基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型；利用基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型及基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型，获得水平井各段的地层参数；
S4，根据所获得的水平井各段的地层相关参数，绘制油藏水平井各段的产出状况图，计算水平井各段贡献率及获取水平井整段的地层信息。
2.根据权利要求1所述的一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，其特征在于，S1中，对示踪剂产出标准浓度进行劈分，获得水平井各段的产量变化史的过程包括如下步骤：
S1.1，将示踪剂浓度标准化；
S1.2，利用标准化的示踪剂浓度计算水平井各段的返排总液量；
S1.3，通过水平井各段的返排总液量得到水平井各段的返排总液量的积分形式；
S1.4，利用水平井各段的返排总液量的积分形式获得水平井各段的产量变化史。
3.根据权利要求2所述的一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，其特征在于，水平井各段的产量变化史包括水平井各段示踪剂累加标准浓度的时间变化史、研究段示踪剂标准浓度的时间变化史和地面总返排量的时间变化史。
4.根据权利要求2所述的一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，其特征在于，S1.1中，将示踪剂浓度标准化时，水平井第j段第i个取样点的示踪剂标准浓度ρji为：
Figure FDA0003976044090000021
S1.2中，水平井第j段的返排总液量Qj为：
Figure FDA0003976044090000022
S1.3中，水平井第j段的返排总液量Qj的积分形式为：
Figure FDA0003976044090000023
式中，M为压裂分段数与压裂示踪剂种类，C0j为注入第j段的示踪剂初始浓度，Cji为第j段第i个取样点的示踪剂浓度，n为整个示踪监测取样个数，ρ(t)为各段累加标准浓度随时间的变化函数；ρj(t)为第j段标准浓度随时间的变化函数；Q(t)为地面总返排量随时间变化的函数，t为生产时间。
5.根据权利要求4所述的一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，其特征在于，S2中，反卷积约束下的压裂示踪剂流量解释模型如下：
Figure FDA0003976044090000024
式中，Δp(t)为第j个裂缝压降；σ为时间的对数函数，σ＝lnΔt；z(σ)为反卷积压力响应导数，
Figure FDA0003976044090000025
Δpu为重整压力；t为生产时间。
6.根据权利要求5所述的一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，其特征在于，基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型P1(t)如下：
Figure FDA0003976044090000031
其中，yeD为储层纵向边界无量纲位置，
Figure FDA0003976044090000032
ye为储层纵向边界位置；Lf为参考长度；xeD为储层横向边界无量纲位置，
Figure FDA0003976044090000033
xe为储层横向边界位置yw1D为第1个裂缝的无量纲纵坐标，
Figure FDA0003976044090000034
yw1为第1条裂缝纵坐标；xw1D为第1个裂缝的无量纲横坐标，
Figure FDA0003976044090000035
xw1为第1条裂缝横坐标；yD为任意位置的纵坐标无量纲位置，
Figure FDA0003976044090000036
xD为任意位置的横坐标无量纲位置，
Figure FDA0003976044090000037
xfD为裂缝无量纲长度，
Figure FDA0003976044090000038
Figure FDA0003976044090000039
u为Laplace变量，上标～表示Laplace变换量；K为地层渗透率，单位为μm2；μ为流体粘度，单位为mPa·s；h为储层高度；qf为裂缝产出液体流量；B为流体体积系数。
7.根据权利要求5所述的一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，其特征在于，基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型P2(t)如下：
Figure FDA00039760440900000310
其中，yeD为储层纵向边界无量纲位置，
Figure FDA00039760440900000311
ye为储层纵向边界位置；Lf为参考长度；xeD为储层横向边界无量纲位置，
Figure FDA0003976044090000041
xe为储层横向边界位置；yw1D为第1个裂缝的无量纲纵坐标，
Figure FDA0003976044090000042
yw1为第1条裂缝纵坐标；xw1D为第1个裂缝的无量纲横坐标，
Figure FDA0003976044090000043
xw1为第1条裂缝横坐标；yD为任意位置的纵坐标无量纲位置，
Figure FDA0003976044090000044
xD为任意位置的横坐标无量纲位置，
Figure FDA0003976044090000045
xfD为裂缝无量纲长度，
Figure FDA0003976044090000046
Figure FDA0003976044090000047
u为Laplace变量，上标～表示Laplace变换量；K为地层渗透率；μ为流体粘度；h为储层高度；qf为裂缝产出液体流量；B为流体体积系数；Kfw为裂缝导流系数，单位为μm2·cm；FCD为无量纲导流系数，
Figure FDA0003976044090000048
8.根据权利要求6或7所述的一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，其特征在于，获得水平井各段的地层参数时，基于压裂示踪约束的单一无限导流垂直裂缝水平井试井模型及基于压裂示踪约束的单一有限导流垂直裂缝水平井试井模型满足如下约束：
反卷积约束下的压裂示踪剂流量解释模型和实测的压力数据之间的标准偏差最小；
反卷积压力响应导数随时间的对数函数曲线拥有最小的总曲率值；
校正输出时，在相同质量的拟合条件下，选择输出要求最小的输出。
9.根据权利要求1所述的一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法，其特征在于，所述水平井各段的地层相关参数包括渗透率、裂缝半长和导流系数。
CN202011602457.1A 2020-12-29 2020-12-29 一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法 Active CN112699554B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011602457.1A CN112699554B (zh) 2020-12-29 2020-12-29 一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011602457.1A CN112699554B (zh) 2020-12-29 2020-12-29 一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法

Publications (2)

Publication Number Publication Date
CN112699554A CN112699554A (zh) 2021-04-23
CN112699554B true CN112699554B (zh) 2023-03-14

Family

ID=75511585

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011602457.1A Active CN112699554B (zh) 2020-12-29 2020-12-29 一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法

Country Status (1)

Country Link
CN (1) CN112699554B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115288649B (zh) * 2022-05-10 2024-04-16 西安石油大学 用于煤层气储层的示踪剂体系及煤层气水平井压裂监测方法
CN115993666B (zh) * 2023-03-23 2023-07-07 成都理工大学 一种油基硅包dna示踪剂的制备方法及其应用

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105370259A (zh) * 2014-08-29 2016-03-02 中国石油化工股份有限公司 水平井分段压裂方法
CN109312616A (zh) * 2016-05-03 2019-02-05 斯伦贝谢技术有限公司 用于水力压裂储层的分析的方法和系统
CN110541704A (zh) * 2019-09-10 2019-12-06 大庆亿莱检验检测技术服务有限公司 一种利用示踪剂评价致密油多段压裂井分段产水量的方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100137163A1 (en) * 2006-01-11 2010-06-03 Link Darren R Microfluidic Devices and Methods of Use in The Formation and Control of Nanoreactors
US8244509B2 (en) * 2007-08-01 2012-08-14 Schlumberger Technology Corporation Method for managing production from a hydrocarbon producing reservoir in real-time
CN109726945B (zh) * 2019-03-06 2020-11-13 中国石油大学(北京) 用于致密油藏复杂缝网的试井解释方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105370259A (zh) * 2014-08-29 2016-03-02 中国石油化工股份有限公司 水平井分段压裂方法
CN109312616A (zh) * 2016-05-03 2019-02-05 斯伦贝谢技术有限公司 用于水力压裂储层的分析的方法和系统
CN110541704A (zh) * 2019-09-10 2019-12-06 大庆亿莱检验检测技术服务有限公司 一种利用示踪剂评价致密油多段压裂井分段产水量的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Pressure transient analysis of multistage fracturing horizontal wells with finite fracture conductivity in shale gas reservoirs;Dai Yu et al;《Environmental Earth Sciences》;20161231;全文 *
水平井压后分段测试试井解释方法研究;马威奇;《中国优秀硕士学位论文全文数据库电子期刊 工程科技I辑》;20180215;第2018年卷(第2期);全文 *

Also Published As

Publication number Publication date
CN112699554A (zh) 2021-04-23

Similar Documents

Publication Publication Date Title
CN109594968B (zh) 一种页岩气多段压裂水平井压后裂缝参数评价方法及系统
CN107506948B (zh) 一种基于动态泄流体积的页岩油气综合产量分析方法
US6842700B2 (en) Method and apparatus for effective well and reservoir evaluation without the need for well pressure history
US9745833B2 (en) Optimizing oil recovery and reducing water production in smart wells
RU2274747C2 (ru) Методика оптимизации добычи из многослойных смешанных пластов с использованием данных о динамике изменения дебита смешанных пластов и данных геофизических исследований в эксплуатационных скважинах
US10233749B2 (en) Multi-layer reservoir well drainage region
AU2013201757B2 (en) Wellbore real-time monitoring and analysis of fracture contribution
CN112699554B (zh) 一种基于压裂示踪约束的致密油藏水平井压后分段试井分析方法
EP1319116A1 (en) Evaluation of multilayer reservoirs
US10982516B2 (en) Systems and methods for operating downhole inflow control valves to provide sufficient pump intake pressure
CN116127675A (zh) 一种页岩油水平井体积压裂最大可采储量的预测方法
WO2017040457A2 (en) Coning transient multi-rate test
Almasoodi et al. Drawdown-management and fracture-spacing optimization in the Meramec Formation: numerical-and economics-based approach
US20190010789A1 (en) Method to determine a location for placing a well within a target reservoir
RU2737043C1 (ru) Способ разработки нефтяного пласта многопластового нефтегазоконденсатного месторождения
Clarkson et al. Analysis of multiple flow/buildup tests including a 5-year buildup: case study of an Australian shale gas well
CN111553067B (zh) 页岩油藏压裂水平井不停井试井测试和解释方法、装置
US11686184B2 (en) Opportunistic techniques for production optimization of gas-lifted wells
CN115809536A (zh) 一种页岩气井多段压裂改造的评价方法
CN108180006B (zh) 一种基于体积压裂后地层能量抬升的水平井产能预测方法
Mohammed et al. Combining Statistical Analysis with Simulation to Optimize Unconventional Completions-Upper and Lower Montney Formations, Canada
Bennis et al. Decline Curve Analysis in Vaca Muerta with Choke Size Normalization of Gas Rates
Carpenter Method Integrates Pressure-Transient and Fracture Area To Detect Well Interference
RU2320855C1 (ru) Способ эксплуатации скважины
Rezaei et al. Fracture diagnostic testing

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