体积压裂水平井特征点拟合反演方法、电子设备及介质技术

公开了一种体积压裂水平井特征点拟合反演方法、电子设备及介质。该方法包括：建立体积压裂水平井分区渗流数学模型；求解体积压裂水平井分区渗流数学模型，进而计算井底压力与井底压力导数；根据体积压裂水平井井底压力导数曲线，确定多个特征点；建立每个特征点与试井解释信息的拟合表达式；调节特征点，使压力导数曲线、压力曲线与对应的实测试井曲线重合，进而计算试井解释信息。本发明专利技术根据体积压裂水平井模型中不同区域对应的渗流特征，进行数学分析，快速有效的对体积压裂水平井的实测试井数据进行拟合，并获得体积压裂水平井的井筒、储层和压裂改造参数，为体积压裂水平井的试井解释和压裂改造效果评价提供依据。

【技术实现步骤摘要】
体积压裂水平井特征点拟合反演方法、电子设备及介质
本专利技术涉及油藏工程试井
，更具体地，涉及一种体积压裂水平井特征点拟合反演方法、电子设备及介质。
技术介绍
水平井体积压裂是低渗致密油气藏工业开发的有效手段，体积压裂可以在X、Y、Z三个方向上实现油气藏的体积改造。该技术不仅能大幅度提高生产井单井产量，同时还能最大限度提高储层动用程度和采收率。随着体积压裂技术的实施，如何用数学方法描述储层改造区压裂缝网，评价压裂效果，确定体积压裂井储层，井筒和压裂改造参数，是目前面临的一些重要问题，本专利技术建立了一种体积压裂水平井试井特征点拟合动态反演方法及装置，能够快速的对体积压裂水平井的试井测试数据进行拟合，反演得到体积压裂水平井的试井解释信息，评价生产井体积压裂的改造效果。现有技术包括：一种体积压裂水平井的非达西渗流数学模型，利用差分方程近似求解，求解速度较慢，并且试井解释方法为常规方法，不能简洁快速的进行拟合获得解释参数；一种致密气藏压裂水平井试井模型，可以用于致密气藏压裂水平井的试井解释，但不能用于油藏体积压裂水平井的试井解释；2区复合油藏不稳定试井模型，该模型可以用于2区复合状油藏水平井的不稳定试井解释，模型中没有考虑体积压裂水平井的水力裂缝和改造封网的影响，不能解决体积压裂水平井试井解释的问题；针对页岩油藏多段压裂水平井，根据实际试井测试资料压力差导数曲线形态，建立了一种同时考虑主裂缝和次裂缝网的体积压裂改造试井模型，但该模型没有考虑井筒远处未经压裂改造区域的流体供给，而是将其视为封闭边界，并且压力和压力导数拟合及井筒、储层参数反演方法为常规的试井解释方法。随着低渗致密油藏的开发，水平井体积压裂成为该类油藏开发的有效手段。目前亟需解决体积压裂井的压力动态分析，储层、裂缝、井筒参数的获取及压裂效果评价等问题。针对以上问题，目前缺少能够有效描述水平井体积压裂改造特征的试井模型及快速动态反演的方法。因此，有必要开发一种体积压裂水平井特征点拟合反演方法、电子设备及介质。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种体积压裂水平井特征点拟合反演方法、装置、电子设备及介质，其能够根据体积压裂水平井模型中不同区域对应的渗流特征，进行数学分析，快速有效的对体积压裂水平井的实测试井数据进行拟合，并获得体积压裂水平井的井筒、储层和压裂改造参数，为体积压裂水平井的试井解释和压裂改造效果评价提供依据。第一方面，本公开实施例提供了一种体积压裂水平井特征点拟合反演方法，包括：建立体积压裂水平井分区渗流数学模型；求解所述体积压裂水平井分区渗流数学模型，进而计算井底压力与井底压力导数；绘制体积压裂水平井井底压力曲线与体积压裂水平井井底压力导数曲线，根据所述体积压裂水平井井底压力导数曲线，确定多个特征点；建立每个特征点与试井解释信息的拟合表达式；调节所述多个特征点，使压力导数曲线、压力曲线与对应的实测试井曲线重合，通过调节后的特征点与所述拟合表达式，计算试井解释信息。优选地，所述体积压裂水平井分区渗流数学模型包括主裂缝区渗流数学模型、改造区渗流数学模型、非改造区渗流数学模型。优选地，所述主裂缝区渗流数学模型为：其中，plfD为主裂缝无因次压力；plsfD为主裂缝壁面处微裂缝无因次压力；xD为x方向无因次距离；yD为y方向无因次距离；FCD为裂缝无因次导流能力。优选地，所述改造区渗流数学模型为：其中，psfD为微裂缝无因次压力；pmD为基质无因次压力；rD为径向无因次半径；tD为无因次时间；CD为无因次井筒储集系数；λ为窜流系数；ω为储容比。优选地，所述非改造区渗流数学模型为：其中，p2D为非改造区无因次压力；γD为无因次渗透率模数；M1,2为内外区流度比，W为内外区储能系数比，e为自然对数。优选地，求解所述体积压裂水平井分区渗流数学模型，进而计算井底压力与井底压力导数包括：分别求解改造区渗流数学模型与非改造区渗流数学模型的通解；计算源函数；将所述主裂缝区渗流数学模型与所述源函数联立，计算井底压力响应；通过Stehfest算法对所述井底压力响应进行Laplace数值反演，获得所述井底压力；针对所述井底压力进行求导，获得所述井底压力导数。优选地，所述源函数为：其中，G为源函数；A、B均为系数；I0为一类零阶的修正贝塞尔函数；K0为二类零阶的修正贝塞尔函数；u为Laplace空间变量。优选地，所述井底压力响应为：其中，为拉普拉斯空间中包含井筒储集效应的无因次井底压力；为拉普拉斯空间中无因次井底压力。作为本公开实施例的一种具体实现方式，第二方面，本公开实施例还提供了一种电子设备，该电子设备包括：存储器，存储有可执行指令；处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的体积压裂水平井特征点拟合反演方法。第三方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的体积压裂水平井特征点拟合反演方法。其有益效果在于：根据体积压裂水平井模型中不同区域对应的渗流特征，进行数学分析，建立体积压裂水平井模型的渗流数学模型。综合考虑了不同区域对应的渗流特征，提高了渗流数学模型建立的准确性，为后续进行试井解释提供了准确的数据基础。通过半解析法获得渗流数学模型的解。根据渗流数学模型的解，可以绘制出体积压裂水平井模型的水平井试井曲线，根据获得的水平井试井曲线进行体积压裂水平井的试井解释。综合考虑不同地质特征对应的渗流特性，使得构建的渗流数学模型更加准确，进一步提高了试井解释结果的准确性。本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术的体积压裂水平井特征点拟合反演方法的步骤的流程图。图2示出了根据本专利技术的一个实施例的类椭圆形改造区体积压裂水平井渗流物理模型的示意图。图3示出了根据本专利技术的一个实施例的体积压裂水平井井底压力曲线与体积压裂水平井井底压力导数曲线的示意图。图4示出了根据本专利技术的一个实施例的体积压裂水平井试井模型压力导数特征曲线的示意图。图5示出了根据本专利技术的一个实施例的不同CDe2S时的压力导数曲线的示意图。图6示出了根据本专利技术的一个实施例的特征点2的坐标值与log(CDe2S)值的关系图。图7示出了根据本专利技术本文档来自技高网...

【技术保护点】
1.一种体积压裂水平井特征点拟合反演方法，其特征在于，包括：/n建立体积压裂水平井分区渗流数学模型；/n求解所述体积压裂水平井分区渗流数学模型，进而计算井底压力与井底压力导数；/n绘制体积压裂水平井井底压力曲线与体积压裂水平井井底压力导数曲线，根据所述体积压裂水平井井底压力导数曲线，确定多个特征点；/n建立每个特征点与试井解释信息的拟合表达式；/n调节所述多个特征点，使压力导数曲线、压力曲线与对应的实测试井曲线重合，通过调节后的特征点与所述拟合表达式，计算试井解释信息。/n

【技术特征摘要】
1.一种体积压裂水平井特征点拟合反演方法，其特征在于，包括：
建立体积压裂水平井分区渗流数学模型；
求解所述体积压裂水平井分区渗流数学模型，进而计算井底压力与井底压力导数；
绘制体积压裂水平井井底压力曲线与体积压裂水平井井底压力导数曲线，根据所述体积压裂水平井井底压力导数曲线，确定多个特征点；
建立每个特征点与试井解释信息的拟合表达式；
调节所述多个特征点，使压力导数曲线、压力曲线与对应的实测试井曲线重合，通过调节后的特征点与所述拟合表达式，计算试井解释信息。


2.根据权利要求1所述的体积压裂水平井特征点拟合反演方法，其中，所述体积压裂水平井分区渗流数学模型包括主裂缝区渗流数学模型、改造区渗流数学模型、非改造区渗流数学模型。


3.根据权利要求2所述的体积压裂水平井特征点拟合反演方法，其中，所述主裂缝区渗流数学模型为：



其中，plfD为主裂缝无因次压力；plsfD为主裂缝壁面处微裂缝无因次压力；xD为x方向无因次距离；yD为y方向无因次距离；FCD为裂缝无因次导流能力。


4.根据权利要求2所述的体积压裂水平井特征点拟合反演方法，其中，所述改造区渗流数学模型为：



其中，psfD为微裂缝无因次压力；pmD为基质无因次压力；rD为径向无因次半径；tD为无因次时间；CD为无因次井筒储集系数；λ为窜流系数；ω为储容比。


5.根据权利要求2所述的体积压裂水平井特征点拟合反演方法，其中，所述非改造区渗流数学模型为：



其中，p2D为非改造区无因次压力；γD...

【专利技术属性】
技术研发人员：王欢，计秉玉，伦增珉，吕成远，赵淑霞，何应付，于洪敏，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京;11