CN103827441A - 用于执行井眼压裂作业的系统和方法 - Google Patents

Abstract

提供了一种在穿透地下地层的井眼附近执行油田作业的方法。该方法包括通过在井眼附近产生裂缝来执行压裂作业。裂缝限定井眼附近的水力裂缝网络(HFN)。该方法还包括通过从HFN外推裂缝数据，在井眼附近产生离散裂缝网络(DFN)。DFN包括裂缝分支和基质块，裂缝分支之间具有交叉点。该方法还包括产生通过DFN的泄油区的深度；限定开采参数；以及基于泄油区的深度和开采参数，执行开采作业，以从地下地层开采流体。开采作业可以包括产生通过DFN的流率；基于流率产生DFN在初始时间的压力分布情况；以及基于压力分布情况产生开采速度。

Description

用于执行井眼压裂作业的系统和方法

对先前申请的交叉引用

本申请要求2011年8月4日提交的美国临时申请No.61/574,521和2011年7月28日提交的美国临时申请No.61/574,131的优先权，通过引用将其全部内容合并于此。

技术领域

本公开总体上涉及用于执行井场作业的方法和系统。更具体地说，本公开致力于用于执行压裂作业的方法和系统，例如，调查地下地层以及表征地下地层中的水力裂缝网络。

背景技术

为了便于从油气井中回收油气，可以通过水力压裂这些井周围的地下地层。可以使用水力压裂来在地表下地层中创建裂纹，以使得油或气能够向井移动。通过一个或多个井眼，以高压、高流率将特别设计的流体(这里称为“压裂流体”或“压裂浆”)引入地层来压裂地层。根据地层中的自然应力，水力裂缝可以从井眼沿着两个相反的方向延伸几百英尺。在某些情况下，它们可以形成复杂的裂缝网络。

压裂流体可以承载有支撑剂，支撑剂是一定尺寸的颗粒，其可以与压裂流体混合以帮助提供用于油气从地层/储层到井眼的开采的有效导路。支撑剂可以包括天然的砂粒或碎石；人造或特别设计的支撑剂，例如纤维、涂覆了树脂的砂石、或高强度陶瓷材料如烧结铝矾土。支撑剂异质或同质地聚集到裂缝中，以在地层中“撑”开新的裂纹或孔。支撑剂创建可渗透导路的面，通过该导路，开采流体可以流入井眼。优选地，压裂流体具有高粘度，并因此能够携带有效数量的支撑剂材料。

压裂流体可以由粘性流体实现，该粘性流体有时被称为“填充物”，被以足以在含油气地层中引起和传播裂缝的流率和压力注入到处理井中。“填充物”的注入一直继续到获得具有充足的几何结构的裂缝以使得能够放置支撑剂颗粒。在注入“填充物”之后，压裂流体可以包括压裂流体和支撑剂材料。压裂流体可以是凝胶、油基、水基、盐水、酸、乳剂、泡沫或任何其它类似的流体。压裂流体可以包含若干添加剂，增粘剂、减阻剂、失水添加剂、阻蚀剂等。为了保持支撑剂悬浮在压裂流体中，直到地层的所有间隔都被如期压裂的时间，支撑剂的密度可以接近所用压裂流体的密度。

支撑剂可以包括任何商业可得熔融材料，例如硅土或氧化物。熔融材料可以包括任何商业可得玻璃或高强度陶瓷产品。在放置支撑剂之后，可以将井关闭足以将压力释放到地层中的时间。这使得裂缝闭合，并且对支撑剂颗粒施加闭合应力。关闭时间长度可以从几分钟到几天变化。

当前的水力裂缝监控方法和系统可以映射裂缝在哪里产生以及裂缝的程度。一些微震监控方法和系统可以通过使用建模的行程时间和/或传播路径，将地震到达时间和极化信息映射到三维空间中，来处理地震事件位置。这些方法和系统可以用来推断水力裂缝随时间的传播。

传统水力裂缝模型还可以假定双翼型诱发裂缝。这些双翼裂缝可能不足以表达一些具有先前存在的自然裂缝的非常规储层中诱发的裂缝的复杂属性。已发布的模型可以基于对微震事件分布的监控来映射离散水力裂缝的复杂几何结构。

在一些情况下，模型可以不必因考虑泵送流体的量或者裂缝和注入流体之间以及裂缝之间的机械相互作用而受约束。一些约束模型可以提供对所涉及的机制的基本理解，但为了提供对水力裂缝传播的准确模拟，在数学描述和/或所需计算处理资源和时间方面可能复杂。

非常规地层，例如页岩，正在发展成为油气开采的源。曾经仅将岩石和密封视为源，现在页岩地层被视为具有紧致孔隙度和低渗透率的非常规储层。通过压裂增产创建的水力裂缝的模式可能是复杂的，并且形成裂缝网络，如关联微震事件分布所指示。已开发复杂水力裂缝网络(HFN)以表示所创建的水力裂缝。裂缝模型的例子在美国专利/申请No.6101447、7363162、7788074、20080133186、20100138196以及20100250215中提供。

可以使用页岩地层的水力压裂来从储层进行增产和开采。已经开发了开采模拟，以估计从储层的开采。对于常规储层，已使用各种开采模拟技术。在下述文章中提供了开采模拟的例子：Warren等人.,"The Behavior ofNaturally Fractured Reservoirs,Soc.Pet.Eng.J.,Vol.3(3):pp.245-255(1963)(下文中称为"Warren&Root");Basquet等人,"Gas Flow Simulation inDiscrete Fracture Network Models".2003年2月3-5日在得克萨斯州的休斯顿召开的SPE Reservoir Simulation Symposium(SPE储层模拟专题讨论会)上提交的论文SPE79708(下文中称为"Basquet");Gong等人.,"Detailed Modelingof the Complex Fracture Network of Shale Gas Reservoirs",2011年1月31日在阿曼的马斯喀特召开的SPE Middle East Unconventional Gas Conferenceand Exhibition(SPE中东非常规天然气会展)上提交的SPE论文142705(下文中称为"Gong");Cinco-Ley等人,"Pressure Transient Analysis for NaturallyFractured Reservoirs",1982年9月26日在路易斯安那州新奥尔良召开的Annual Fall Technical Conference and Exhibition(年度秋季技术会展)上提交的SPE论文11026(下文中称为"Cinco-Ley");Xu等人,"Quick Estimate ofInitial Production from Stimulated Reservoirs with Complex Hydraulic FractureNetwork",2011年10月30日在美国科罗拉多州丹佛召开的SPE AnnualTechnical Conference and Exhibition(SPE年度技术会展)上提交的论文SPE146753(下文中称为"Xu2011");以及C.E.Cohen等人，"Production ForecastAfter Hydraulic Fracturing in Naturally Fractured Reservoir:Coupling aComplex Fracturing Simulator and a Semi-Analytical Production Model",2012年2月8日在美国得克萨斯州伍德兰兹召开的SPE Hydraulic FracturingTechnology Conference and Exhibition(SPE水力压裂技术会展)上提交的论文(SPE152541)，通过引用将其全部内容合并于此。然而，储层有可能是非常规的和/或具有自然裂缝，例如具有页岩的储层。

发明内容

在至少一个方面，本公开涉及一种在穿透地下地层的井眼附近执行开采作业的方法。该地下地层附近具有多个裂缝。该方法包括产生通过离散裂缝网络的流率，该离散裂缝网络是由地下地层中的多个裂缝限定的。该离散裂缝网络包括多个裂缝分支和多个基质块，多个裂缝分支之间具有交叉点。该方法还包括基于流率产生离散裂缝网络在初始时间的压力分布情况，以及基于压力分布情况产生开采速度。

在另一个方面，本公开涉及一种在穿透地下地层的井眼附近执行油田作业的方法。该地下地层中具有储层。该方法包括执行压裂作业，压裂作业包括在井眼附近产生裂缝。裂缝限定井眼附近的水力裂缝网络。该方法还包括通过从水力裂缝网络外推裂缝数据，在井眼附近产生离散裂缝网络。离散裂缝网络包括多个裂缝分支和多个基质块，多个裂缝分支之间具有交叉点。该方法还包括产生通过离散裂缝网络的泄油区的深度，限定至少一个开采参数，以及基于泄油区的深度和至少一个开采参数，执行开采作业，以从地下地层开采流体。

最后，在另一方面，本公开涉及一种在穿透地下地层的井眼附近执行油田作业的方法。该方法包括：通过将流体注入地下地层以使得在井眼附近产生裂缝，来对井眼进行增产；测量裂缝，并基于所测量的裂缝来限定水力裂缝网络。

该方法还包括通过从水力裂缝网络外推裂缝数据，在井眼附近产生离散裂缝网络。该离散裂缝网络包括多个裂缝分支和多个基质块，多个裂缝分支之间具有交叉点。该方法还包括产生通过离散裂缝网络的泄油区的深度；限定至少一个开采参数；基于泄油区的深度和开采参数，估计随时间的开采速度；以及基于所估计的开采速度从地下地层开采流体。

提供此发明内容部分是为了介绍下面在详细描述中将进一步描述的构思的选集。此发明内容部分并不意图确定所要求保护的主题的关键或必要特征，也不意图用来帮助限定所要求保护的主题的范围。

附图说明

参考附图来描述用于表征井眼应力的系统和方法的实施例。在所有附图中，使用相同的附图标记来表示相同的特征和组件。

图1.1-1.4是图解井场处的各种油田作业的示意图；

图2.1-2.4是通过图1.1-1.4的作业收集的数据的示意图；

图3是描绘压裂作业的水力压裂现场的示意图；

图4.1和4.2分别是描绘执行油田作业和开采作业的方法的流程图；

图5是从水力裂缝模拟提取的离散裂缝网络（DFN）的开采模拟的示意性图解；

图6是具有多个基质块的图5的DFN的示意性图解；

图7是通过基质块的流动的近似的示意性图解；

图8.1-8.3分别是图解井的开采、累积开采量和压力的示意性图解；

图9是描绘基质块的裂缝的坐标的示意图；

图10是描绘从基质块到DFN的分支的流率的示意图；

图11.1和11.2是描绘高传导率DFN的压力与时间的关系随时间的变化的曲线图；

图12是高传导率DFN的归一化压力和时间延迟随时间的变化的曲线图；

图13是高传导率DFN的累积开采量随时间的变化的曲线图；

图14.1和14.2是低传导率DFN的压力与时间的关系随时间的变化的曲线图；

图15是低传导率DFN的归一化压力和时间延迟随时间的变化的曲线图；

图16是低传导率DFN的累积开采量随时间的变化的曲线图；

图17是使用非常规开采模型(UPM)的低传导率DFN的归一化压力和时间延迟随时间变化的曲线图；

图18是使用UPM的低传导率DFN的累积开采量随时间的变化的曲线图；

图19是压力和时间延迟随时间变化的曲线图的表；

图20是比较使用储层模拟器和UPM模拟的开采随时间变化的曲线图；

图21.1和21.2分别是描绘由储层模拟器和UPM描绘的DFN的示意图；

图22是对于不同的裂缝传导率比较使用储层模拟器和UPM模拟的开采随时间的变化的曲线图；

图23.1和23.2分别是由储层模拟器、UPM和没有延迟的UPM得到的流率和累积开采量随时间的变化的曲线图。

具体实施方式

下面的描述包括实施这里的主题的技术的示例性系统、装置、方法和指令序列。然而，应当理解，所描述的实施例可以在没有这些具体细节的情况下实施。

本公开涉及用于执行压裂作业以估计和/或预测产量的技术。压裂作业涉及裂缝建模，裂缝建模利用椭圆和线网模型来估计产量。

图1.1-1.4图示了可以在井场执行的各种油田作业，图2.1-2.4图示了可以在井场收集的各种信息。图1.1-1.4图示了代表性油田或井场100的简化示意图，该代表性油田或井场100具有地下地层102，地下地层102中包含例如储层104，并且还图示了在井场100执行的各种油田作业。图1.1图示了由勘探工具如地震车106.1执行以测量地下地层的属性的勘探作业。勘探作业可以是用于产生声振动的地震勘探作业。在图1.1中，由源110产生的一种这样的声振动112在地层116中的多个水平层114处反射。可以由位于地球表面的传感器如地震检波器—接收器118接收声振动112，并且地震检波器118产生电输出信号，在图1.1中称为“所接收的数据”120。

响应于代表声振动112的不同参数(例如幅度和/或频率)的所接收的声振动112，地震检波器118可以产生包含有关地下地层的数据的电输出信号。可以提供所接收的数据120作为对地震车106.1的计算机122.1的输入数据，并且响应于输入数据，计算机122.1可以产生地震和微震数据输出124。可以对地震数据输出124进行存储、发送、或根据期望进行进一步处理如数据约减。

图1.2图示了由钻井工具106.2执行的钻井作业，其中钻井工具106.2由钻机128悬挂，并且被推进到地下地层102中，以形成井眼136或其它通道。可以使用泥浆坑130将钻井泥浆经由管线132抽吸到钻井工具中，以使钻井泥浆循环通过钻井工具上到井眼136并返回地表。钻井泥浆可以被过滤，然后返回泥浆坑。可以使用循环系统来存储、控制或过滤流动的钻井泥浆。在这个图示中，钻井工具被推进到地下地层以到达储层104。每个井可以以一个或多个储层为目标。钻井工具可以被适配以使用随钻测井工具测量井下属性。随钻测井工具还可以被适配以如图所示收取岩心样本133，或者被移除以便可以使用其它工具来收取岩心样本。

可以使用地面单元134来与钻井工具和/或场外作业进行通信。地面单元可以与钻井工具通信，以向钻井工具发送命令，并从钻井工具接收数据。地面单元可以具有计算机设备，以接收、存储、处理和/或分析来自作业的数据。地面单元可以收集钻井作业期间产生的数据，并产生可以被存储或发送的数据输出135。计算机设备，例如地面单元中的计算机设备，可以位于井场附近各种位置处和/或位于远距离处。

可以在油田附近放置诸如计量器的传感器(S)，以收集与先前描述的各种作业有关的数据。如图所示，传感器(S)可以放置在钻井工具中一个或多个位置处和/或位于钻机处，以测量钻井参数，如钻压、钻头扭矩、压力、温度、流率、成分、旋转速度和/或其它作业参数。传感器(S)还可以位于循环系统中的一个或多个位置中。

可以由地面单元和/或其它数据收集源来收集由传感器收集的数据，以进行分析或其它处理。可以单独使用或与其它数据结合使用由传感器收集的数据。可以将数据收集在一个或多个数据库和/或就地或场外发送。可以选择性地使用全部或选择的部分数据来对当前和/或其它井眼进行分析和/或预测作业。数据可以是历史数据、实时数据或它们的组合。可以实时使用实时数据，或将其存储以备以后使用。还可以将数据与历史数据或其它输入组合以进行进一步的分析。可以将数据存储在分离的数据库中，或者组合到单个数据库中。

可以使用所收集的数据来执行分析，如建模作业。例如，可以使用地震数据输出来执行地质、地球物理、和/或储层工程分析。可以使用储层、井眼、地面和/或处理后的数据来执行储层、井眼、地质、以及地球物理或其它模拟。来自作业的数据输出可以从传感器直接产生，或在一些预处理或建模之后产生。这些数据输出可以用作其它分析的输入。

数据可以被收集并存储在地面单元134处。一个或多个地面单元可以位于井场或者在远处连接到井场。地面单元可以是单个单元或多个单元的复杂网络，用于执行整个油田必要的数据管理功能。地面单元可以是手动或自动系统。地面单元134可以由用户操作和/或调整。

地面单元可以具有收发器137，以使得地面单元和当前油井的各个部分或其它位置处之间能够进行通信。地面单元134还可以具有或者功能性地连接到一个或多个控制器，以驱动井场100处的机械装置。然后地面单元134可以响应于所接收的数据向油田发送命令信号。地面单元134可以经由收发器接收命令，或者可以自己执行给控制器的命令。可以提供处理器以(本地或远程)分析数据，做出决定和/或驱动控制器。以这种方式，可以基于所收集的数据选择性地调整操作。可以基于该信息来优化部分操作，例如控制钻井、钻压、泵送速率或其它参数。这些调整可以基于计算机协议自动进行，和/或由操作员手动进行。在一些情况下，可以调整井规划以选择最佳操作条件，或者避免问题。

图1.3图示了由钻机128悬挂并进入图1.2的井眼136的电缆测井工具106.3执行的电缆测井作业。电缆测井工具106.3可以被适配以部署到井眼136中，用来产生测井曲线，执行井下测试和/或收集样本。电缆测井工具106.3可以用来提供另一种执行地震勘探作业的方法和设备。图1.3的电缆测井工具106.3可以例如具有爆炸性、放射性、电学或声学能量源144，该能量源144向周围的地下地层102及其中的流体发送电信号和/或从周围的地下地层102及其中的流体接收电信号。

电缆测井工具106.3可以可操作地连接到例如图1.1的地震车106.1的地震检波器118和计算机122.1。电缆测井工具106.3还可以向地面单元134提供数据。地面单元134可以收集在电缆测井作业期间产生的数据，并产生可以被存储或发送的数据输出135。电缆测井工具106.3可以位于井眼中各种深度，以提供勘察结果或与地下地层有关的其它信息。

可以在井场100附近放置诸如计量器的传感器(S)，以收集与先前描述的各种作业有关的数据。如图所示，传感器(S)放置在电缆测井工具106.3中，以测量涉及例如孔隙率、渗透率、流体成分和/或作业的其它参数。

图1.4图示了由从开采单元或“圣诞树”129部署并进入图3中完成的井眼136的开采工具106.4执行的开采作业，用于将流体从井下储层抽吸到地表设施142。流体从储层104通过套管(未示出)中的射孔并进入井眼136中的开采工具106.4，并经由集输管网146到地面设施142。

可以在油田附近放置诸如计量器的传感器(S)，以收集与先前描述的各种作业有关数据。如图所示，传感器(S)可以放置在开采工具106.4或相关设备如“圣诞树”129、集输管网、地面设施和/或开采设施中，以测量流体参数如流体成分、流率、压力、温度和/或开采作业的其它参数。

尽管仅示出了简化的井场配置，但应当理解油田或井场100可以覆盖具有一个或多个井场的陆地、海洋和/或水域的一部分。为了提高采收率或存储例如碳氢化合物、二氧化碳或水，开采也可以包括注入井(未示出)。一个或多个集输设施可以可操作地连接到一个或多个井场，以从井场选择性地收集井下流体。

应当理解，图1.2-1.4图示的工具不但可以测量油田属性而且可以测量非油田作业的属性，例如矿藏、含水层、储藏以及其它地下设施。而且，尽管图示了特定的数据获取工具，但是应当理解，可以使用能够感测诸如地下地层的地震双向行程时间、密度、电阻率、开采速度等和/或其地质信息的参数的各种测量工具(例如电缆测井、随钻测量(MWD)、随钻测井(LWD)、岩心采样等)。可以沿着井眼和/或监测工具在各种位置处放置各种传感器(S)，以收集和/或监测所期望的数据。还可以从场外位置处提供其它数据源。

图1.1-1.4的油田配置图示了井场100和通过这里提供的技术可以使用的各种作业的示例。油田的部分或全部可以在陆地上、水上和/或海上。而且尽管图示了在单个位置处测量单个油田的情况，但是可以以一个或多个油田、一个或多个处理设施、以及一个或多个井场的任何组合来利用储层工程。

图2.1-2.4分别是由图1.1-1.4的工具收集的数据的示例的图形表示。图2.1表示由地震车106.1收取的图1.1的地下地层的地震道202。地震道可以用于提供诸如在一段时间中的双向响应的数据。图2.2图示了由钻井工具106.2收取的岩心样本133。岩心样本可以用于提供诸如沿着岩心的长度方向上岩心样本的密度、孔隙率、渗透率或其它物理属性。可以在变化的压力和温度下对岩心中的流体执行密度和粘度的测试。图2.3图示了由电缆测井工具106.3收取的图1.3的地下地层的测井曲线204。电缆测井可以提供各种深度处地层的电阻率或其它测量结果。图2.4图示了在地面设施142处测量的流过图1.4的地下地层的流体的产量递减曲线或图表206。产量递减曲线可以提供作为时间t的函数的开采速度Q。

图2.1、2.3和2.4的各个图表图示了可以描述或提供有关地层或其中所包含的储层的物理特性的信息的静态测量结果。这些测量结果可以被分析以限定地层的属性，从而确定测量的精确度和/或检查错误。可以将各个测量中每个的图对齐或缩放，以进行属性的比较与核实。

图2.4图示了通过井眼对流体属性的动态测量结果。随着流体流过井眼，对流体属性如流率、压力、成分等进行测量。如下面所描述的，可以对静态和动态测量结果进行分析，并用于产生地下地层的模型，以确定其特性。也可以使用类似的测量结果来测量地层方面随时间的变化。

油田作业

可以在从井眼产生开采之前、期间或之后模拟开采作业。可以使用各种技术来执行对从复杂的破裂储层的开采的模拟。可以使用双重孔隙度模型来解决裂缝和储层(基质)的其余部分之间的属性差异。双重孔隙度可以考虑彼此连接的两个粗网格，一个用于裂缝网络，另一个用于基质。该方法还可以包括对(例如裂缝网络的)属性取平均和简化，以对两种介质之间的交换项进行建模。该方法可以用于例如天然破裂的储层。对于裂缝网络的近井眼效应，例如在通过水力压裂创建的网络的情况下，可以提供附加的分析。在先前通过引用合并于此的Warren和Root的论文中描述了双重孔隙度技术。

另一种方案包括使用既包括裂缝又包括储层的一种介质以及细化的数值网络。可能需要额外的计算时间来进行处理。可以使用例如专门的储层模拟器来提供网格化的灵活性(例如非结构化的网眼产生)。

又一种方案包括使用关于离散裂缝网络(DFN)的双重孔隙度方程。先前通过引用合并于此的Basquet论文中提供了DFN的例子。可以使用其它的方法来模拟从基质到裂缝的流动。在一些情况下，例如具有压缩储层流体(例如气体)的情况下，可以考虑从每个基质块到DFN的开采历史。可以使用方程组中附加的未知数来将基质块网格化。先前通过引用合并于此的Gong论文提供了网格化的例子。还可以提供分析解来模拟流动。可以从连续性方程的拉普拉斯变换推导出解。先前通过引用合并于此的Cinco-Ley和Xu2011论文提供了分析解的例子。

可以考虑瞬变裂缝压力来获得可以使用时间的数值积分的复杂表达。也可以考虑恒定裂缝压力，也可以获得相对于压力为线性的基质和裂缝之间的流率的表达式。可以在例如DFN内的压力的变化可以忽略(例如处于恒定井眼压力)的传导性裂缝中使用这种解。本公开可以采用一种或多种方案来产生分析解。在水力压裂和/或天然破裂储层的情况下，这种解可以延及一定范围的裂缝传导性。

本公开提供针对在天然破裂储层中进行水力压裂的情况下一定范围的裂缝传导性的分析解。这种模拟可以适用于非常规储层，例如页岩气，尽管其也可以适用于其它地下地层。这些非常规储层具有两个主要特征：低岩石渗透率和稠密的天然裂缝网络。可以提供增产方案来解决非常规或可以包括水平井和大规模水力压裂处理以进行开采的其它储层的开采模式的潜在区别。在一些情况下，这些处理引出与天然裂缝相互作用的水力裂缝，并且可以导致将井连接到储层的复杂裂缝网络。

本公开文件公开了一种模拟从储层(例如在创建了复杂的水力裂缝网络之后的非常规(天然破裂)储层)进行的开采的方法。所公开的方法首先将来自非常规裂缝模型(UFM)的结果外推，然后在一定时间期限内以一定精确度范围，用将向用户给出对井在几年内的开采的预测的方法对其进行处理。本申请的方法将针对全范围的裂缝传导率的半分析模型的有效性延伸以在实际情况下进行考虑。可以相对于通过储层模拟(例如从斯伦贝谢技术公司商业可得的ECLIPSETM(www.slb.com))进行的模拟来验证该模拟器，以阐明该算法对于给定范围的裂缝传导率提供精确结果的能力。

本公开文件还公开了模拟从已通过水力压裂来增产的天然破裂储层进行的开采的方法。该方法的部分可以以模拟水力裂缝处理的软件程序实现。本方法可以首先将来自模拟的结果外推，以重新创建经适配的水力裂缝网络，该经适配的水力裂缝网络在网络的交叉点之间具有平均属性，然后估计每个裂缝面前端的等效块深度。最后，可以输入用于开采条件的参数，并且运行开采模拟器。从与裂缝接触的每个基质块的开采使用分析表达，该分析表达可以延及成问题的参数(传导率、渗透率等)的全范围的实际值。这是通过对于每个裂缝面在每个时间步更新初始开采时间，以补偿每个基质块的开采中的延迟，并在某种程度上保持储层流体在放置/开采中的质量平衡，来实现的。这一更新是通过计算这个初始时间的搜索算法来执行的，使得从基质块的每一侧实际开采的质量等于在邻近裂缝中的当前压力条件已在时间方面变得恒定并且已在更新后的初始时间启动的情况下从基质块的每一次开采的质量。本方法可以与通过现有储层模拟器如ECLIPSETM进行的模拟作比较。可以执行针对一定范围的裂缝传导率的结果，并且与储层模拟器进行交叉检验。

图3是图解根据本公开的用于水力压裂地下地层(这里称为“裂缝现场”)的示例性操作设置。裂缝现场300可以位于陆地上或在水环境中，并且包括延伸进入地下地层的处理井301以及延伸进入地下地层并与处理井301偏移的监控井303。监控井303包括以规定距离间隔开的地震检波器接收器305(例如三元件地震检波器)的阵列，如图所示。

在压裂作业期间，压裂流体被从地表311泵送到处理井301中，使得油气储层307中的周围地层被压裂，并形成水力裂缝网络308。这样的压裂产生微震事件310，微震事件310发射压缩波(也被称为“初波”或“P波”)和剪切波(也称为“次波”或“S波”)，压缩波和剪切波通过大地传播，并被监控井303的地震检波器接收器阵列305记录。

可以通过测量P波和S波的到达时间之差来计算到微震事件310的距离。而且，可以使用检查P波的粒子运动的矢端图分析来确定指向事件的方位角。通过使用阵列305的接收器之间的P波和S波到达延迟来约束事件310的深度。这些微震事件310的距离、方位角和深度值可以用于导出由压裂流体随时间而导致的破裂的几何边界或轮廓，例如由高度h、椭圆纵横比e和长轴a限定的椭圆边界，如图3所示。

现场301还包括用于将高压压裂流体供应到处理井301的压裂流体供应和泵送装置(未示出)。压裂流体可以在其中预先混合了支撑剂(以及可能的其它特殊组分)的情况下存储。可替换地，压裂流体可以在没有预先混合支撑剂或其它特殊组分的情况下存储，而支撑剂(和/或其它特殊组分)通过如美国专利No.7,516,793(通过引用将其全文合并于此)中所描述的过程控制系统，以受控的方式，混合到压裂流体中。处理井301还包括：示意性地示出的流率传感器S，用于测量供应给处理井的压裂流体的泵送速度；和井底压力传感器，用于测量处理井301中压裂流体的井底压力。

数据处理系统309链接到监控井303中的阵列305的接收器和处理井301的传感器S(例如流率传感器和井底压力传感器)。数据处理系统309可以与地面单元134合并，和/或与地面单元134一起工作。数据处理系统309执行图4中示出并在此描述的处理。如本领域技术人员将理解，数据处理系统309包括数据处理功能(例如，一个或多个微处理器、关联存储器、以及其它硬件和/或软件)，以实现这里描述的公开内容。

数据处理系统309可以由工作站或位于现场301的其它合适的数据处理系统来实现。可替换地，数据处理系统309可以由分布式数据处理系统实现，其中数据通过通信链接(典型地为卫星链接)(优选实时地)传送到远程位置，以进行这里描述的数据分析。数据分析可以在工作站或其它合适的数据处理系统(例如计算机簇或计算网格)执行。此外，本公开的数据处理功能可以存储在数据存储设备(例如一个或多个光盘或其它可手持非易失性存储设备、或可通过网络访问的服务器)上，并且按照需要加载到合适的数据处理系统上，以如这里所描述的在其上执行。

图4.1是描绘执行油田作业的方法400.1的流程图。该方法包括420执行压裂作业(实际或模拟)，422在井眼附近产生DFN，424产生通过DFN的泄油区的深度，426限定至少一个开采参数，以及428执行开采作业。

图4.2描绘了执行开采作业的方法400.2。这一开采作业可以与图4.1的开采作业428相同。在方法400.2的版本中，开采作业是模拟的。如图4.2所示，方法400.2包括421产生通过离散裂缝网络的流率，423基于流率产生离散裂缝网络的压力分布情况，以及425基于压力分布情况产生开采速度。该方法还可以包括427验证开采速度。该方法可以具有其它特征，并且可以以任何顺序执行。

执行压裂作业420包括在井眼附近产生裂缝，并且在井眼附近限定水力裂缝网络。这个压裂作业可以通过例如如图3所示实际注入流体来执行。还可以使用水力裂缝模拟器来模拟对井的水力压裂。模拟可以包括在井眼附近产生裂缝网络。在美国专利申请No.20100307755中提供了离散裂缝网络技术。来自实际或模拟水力压裂的数据可以用来产生描述所得DFN的数据。

可以由图5所示的计算机产生的图像来视觉呈现水力裂缝模拟530。水力裂缝模拟530包括形成水力裂缝网络536的多个裂缝534。裂缝网络536中描绘了裂缝网络536的特征，例如浆料538、流体540和堤岸542。

产生DFN422包括从水力裂缝网络外推裂缝数据。可以通过裂缝数据的外推来产生DFN。可以从水力裂缝模拟530外推裂缝数据。这个数据可以自动输出以形成开采网络显现532，如箭头533所描绘。图5示出了从水力裂缝模拟530输出到开采网络显现532的数据的示例。开采网络显现532提供了从所测量的裂缝数据创建模拟水力裂缝网络成为等效DFN网络的例子。可以执行该输出以按照可以由开采模型使用的格式创建DFN535。

在图5所示的例子中，DFN535包括分支544和交叉点(或裂缝端部)546。这些裂缝分支544和交叉点546提取水力裂缝模拟530中描绘通过裂缝网络536的流体流动的部分。其余裂缝534被消除。

DFN535的格式对于每个裂缝分支544处的每个属性考虑唯一的平均值。裂缝分支544被限定为连接两个交叉点536的平面。这些交叉点536可以是裂缝交叉点，或裂缝交叉点和端部。每个裂缝分支544处的属性可以是例如分支极点的空间坐标；分支位置处的平均传导率、平均高度、平均储层压力；和/或分支位置处的平均储层渗透率。

可以由本模型来使用由交叉点546和分支544对DFN535的描述，来计算交叉点546处的压力。这一描述也可以使用分支544以便使两者都连接到交叉点546并且计算从相邻基质块的开采。

返回参考图4，产生424通过DFN535的泄油区的深度的步骤可以使用基质块来执行。如图6所示，图5的开采网络显现532被修改为描绘修改后的DFN535'的开采网络显现532，其中在每个裂缝分支544的前端有基质块648，如图6所示。每个基质块具有深度650。

开采网络显现532'提供了在所有分支544的每一侧将被耗尽的基质深度650的产生。修改后的DFN535'可以用于为每个基质块648自动或手动产生泄油区650的深度。这可以以(不与任何储层边界接触的)给定基质块的实际总体积可以被引流的方式完成。

图7示意性地描绘了流体通过基质块的流动。这个附图图解了等价块长度的定义和等价块深度的计算。在所示出的例子中，对于由相同长度的四个裂缝分支544包围的方形基质块648.1，可以假定基质块648的每个四分之一区域752可以被其所接触的裂缝分支544耗尽。图中描绘了基质块648的被耗尽的体积754.1和要被耗尽的等价块深度755。

假定从基质块648.1到裂缝分支544的线性流动(如这里将进一步详细描述)，还可以假定基质块648的这个四分之一区域752具有裂缝分支544的长度L。所以，基质块648.2的这个“四分之一区域”752的深度必须等于块长度L的四分之一(或者说是L/4)，以便要被耗尽的总体积相同。如箭头733所指示，使用线性流动近似，可以确定要被耗尽的基质块648.2的体积754.2的等价块深度L/4。可以使用更复杂的块形状，但是可能会涉及更加复杂的技术。

再次参考图4，限定426一个或多个开采参数的步骤可以通过获得用户输入来执行。用户可以限定一个或多个开采参数，以在模拟中考虑。用户可以基于某种标准或根据需要选择这些开采参数。可以选择的开采参数的例子包括井底压力(BHP)、储层条件下的储层流体粘度、储层条件下的储层流体压缩率、要模拟的开采的持续时间，以及其它。

执行开采作业的步骤428包括基于泄油区的深度和至少一个开采参数从地下地层开采流体。开采作业可以是实际的或模拟的。实际开采作业包括如图1.4所示将流体开采到地面。模拟开采可以使用开采模拟来执行。还可以提供开采结果的显现。这样的显现可以使得用户能够视觉看到开采衰减和累积开采量，而且还能够显现裂缝网络和基质块中的压力场的动力学。图8.1-8.3提供了开采数据随时间(例如140天)变化的显现的例子。

图8.1是描绘开采速度856.1的曲线图800.1。曲线图800.1绘制了每天的开采量(Mscf/d)(y轴)相对于以天为单位的时间t(x轴)的关系。图8.2是描绘累积开采量856.2的曲线图800.2。曲线图800.2描绘了累积开采量P(MMscf)(y轴)相对于时间t(x轴)的关系。图8.3是描绘储层压力(z轴)相对于距离x(m)(x轴)和距离y(m)(y轴)的关系以及裂缝网络858中和基质块848中的压力的三维曲线图800.3。可以提供这些和其它描绘图。可以基于这些开采估计来调整开采作业。

开采作业

开采作业(428和/或400.2)将分三部分来描述。首先，给出在分析及其分析解中使用的方程。第二，提供关于该模型的传导率的效果，以及针对高和低传导率涉及单个裂缝分支的例子。第三，提供模型的有效性和解决问题的方法，例如传导率。

1.分析解

可以使用控制方程和分析解来确定开采速度。将用于多孔介质中可压缩流体的连续性方程应用于基质和裂缝两者。在裂缝网络中，连续性方程可以如下重写：

$\frac{\∂}{{\∂x}_f}\left(\mathrm{\ρ}{Q}_f(x_f,t)\right)=-\mathrm{\ρ}{Q}_{\mathrm{mf}}(x_f,t)---\left(1\right)$

Q_mf是从基质到储层的流率，Q_f是裂缝内的流率，ρ是流体密度，x_f是沿裂缝的轴。假定裂缝渗透率(传导率除以宽度)如此之大，以致于对于开采模拟所考虑的时间尺度(几天到几年)，连续性方程的瞬变项可以忽略。可以假定裂缝内的达西流率。

${\mathrm{\ρQ}}_{\mathrm{mf}}(x_f,t)=-\frac{M}{2\mathrm{RT}}\frac{\∂}{{\∂x}_f}(-\mathrm{HC}\frac{\∂m\left(P_f\right)}{{\∂x}_f})---\left(2\right)$

P_f是裂缝内的压力，C是传导率，T是温度。函数m是拟压力。参见Al-Hussainy等人.,"The Flow of Real Gases Through Porous Media",Journal ofPetroleum Technology,1966,pp.624-36。

$m\left(P_f\right)={2\∫}_{P_{\mathrm{LB}}}^{P_f}\frac{P}{\mathrm{\μ}\left(P\right)Z\left(P\right)}\mathrm{dP}---\left(3\right)$

在基质内，用于可压缩流体的连续性方程取下面的形式。

P_m是基质内的压力，k_m是基质渗透率，c_t是流体压缩率，μ是粘度，Z是体积因子，

是基质的孔隙度。为了简化，方程4可以重写为：

$\frac{{\∂}^{2}m}{{\∂x}_m^2}=\frac{1}{a}\frac{\∂m}{\∂t}---\left(5\right)$

a在方程6中定义。

为了计算Q_mf，可以求解方程(5)，其中x_m是沿垂直于裂缝964及其坐标xf的轴的坐标。图9提供了裂缝964和基质块648中的坐标的图解说明。

可以使用拉普拉斯变换来求方程5的解，例如在Jeannot,Yves."ThansfertThermique",Textbook,Ecole des Mines de Nancy,2009.http://www.thermique55.com/principal/therrnique.pdf;以及Bello,R.O.,"RateTransient Analysis in Shale Gas Reservoirs with Transient Linear Behavior",博士论文,2009中所说明的。这里提供了用于所给出的方法的方程、实现、算法和变量的详细列表。

对于恒定裂缝压力，可以确定基质中的压力分布情况。该模型的第一个假定是可以用下面的实际气体方程来描述气体特性：

$\mathrm{\ρ}=\frac{M}{\mathrm{RT}}\left(\frac{P_m}{Z\left(P_m\right)}\right)$ 基质内                 (7)

$\mathrm{\ρ}=\frac{M}{\mathrm{RT}}\left(\frac{P_f}{Z\left(P_f\right)}\right)$ 裂缝网络内          (8)

用于基质块内部的线性气体流动的基本方程是

其中

$c_t\left(P_m\right)=\frac{1}{P_m}-\frac{1}{Z\left(P_m\right)}\frac{\∂Z\left(P_m\right)}{\∂P_m}---\left(10\right)$

下面对拟压力的定义将简化前面的方程的解

$m\left(P_m\right)=2{\∫}_{P_{\mathrm{LB}}}^{P_m}\frac{P}{\mathrm{\μ}\left(P\right)Z\left(P\right)}\mathrm{dP}---\left(11\right)$

于是方程(11)变为

$\frac{{{\∂m}_{*}}^2}{\∂x_m^2}=\frac{1}{a}\frac{\∂m_{*}}{\∂t}---\left(12\right)$

其中

m_*=m(P_m)-m(P_{m_0})=m(P_m)-m_{m_0}               (14)并且边界条件为

m_*(x_m=L，t)=m*(x_m=-L，t)=m_f-m_{m_0}                 (15)

$\frac{{\∂m}_{*}}{\∂x_m}(0,t)=0---\left(16\right)$

m_*(x_m，t₀)=0                        (17)

对方程(15)进行拉普拉斯变换给出

$\frac{{\∂\mathrm{\θ}}^2}{\∂{x_m}^2}-q^2\mathrm{\θ}=0---\left(18\right)$

其中

$q^2=\frac{s}{a}---\left(19\right)$

对于如下形式的解

θ(x_m,s)=Axosh(qx_m)+Bsinh(qx_m)            (20)方程(16)的拉普拉斯变化是

$\frac{\∂\mathrm{\θ}}{{\∂x}_m}(0,s)=0---\left(21\right)$

其给出B=0，因此

θ(x_m,s)=Acosh(qX_m)                 (22)

现在，假定裂缝网络中的压力几乎恒定。

m_f(x_f，t)=m_f(x_f)

因此，对方程(15)进行拉普拉斯变化给出

$\mathrm{\θ}(x_m=L,s)=\mathrm{\θ}(x_m=-L,s)=\frac{m_f-{m_m}_o}{s}---\left(23\right)$

由此给出

$A=\frac{m_f-m_{m\_o}}{s\cosh \left(\mathrm{qL}\right)}---\left(24\right)$

并且

$\mathrm{\θ}(x_m,s)=\frac{m_f-m_{m\_o}}{s\cosh \left(\mathrm{qL}\right)}\cosh \left({\mathrm{qx}}_m\right)=\frac{m_f-m_{m\_o}}{s}\frac{e^{\mathrm{qx}}+e^{-\mathrm{qx}}}{e^{\mathrm{qL}}(1+e^{-2\mathrm{qx}})}---\left(25\right)$

通过使用

泰勒级数展开，可以提供下面的结果：

$\mathrm{\θ}(x_m,s)=\frac{m_f-m_{m\_o}}{s}(e^{-q(L-x_m)}+e^{-q(L+x_m)}){\mathrm{\Σ}}_{n=0}^{\∞}{(-1)}^n{e}^{-2\mathrm{nqL}}---\left(26\right)$

逆拉普拉斯变换给出

$\begin{array}{c}m\left(P_m\right)-m_{m\_0}=\\ (m_f-m_{m\_o})\underset{n=0}{\overset{\∞}{\mathrm{\Σ}}}{(-1)}^n\left(\mathrm{erfc}\right[\frac{(2n+1)L-x_m}{2\sqrt{a(t-t_0)}}]+\mathrm{erfc}[\frac{(2n+1)L+x_m}{2\sqrt{a(t-t_0)}}\left]\right)\end{array}---\left(27\right)$

然后可以确定恒定裂缝压力下从基质到裂缝的流率。由达西定律给出从基质到裂缝的流率：

$\mathrm{\ρ}{Q}_{\mathrm{mf}}(x_f,t)=-\frac{{\mathrm{MHk}}_m}{2\mathrm{RT}}{\mathrm{\Σ}}_{k=1}^2\frac{\∂m\left(P_m\right)}{{\∂x}_m}{|}_{x_m={-L}_k,L=L_k}---\left(28\right)$

L_k对应于裂缝k侧上的泄油区的最大长度。

$\frac{\∂m\left(P_m\right)}{{\∂x}_m}=\frac{(m_f-m_{m\_o})}{\sqrt{\mathrm{\πa}(t-t_0)}}{\mathrm{\Σ}}_{n=0}^{\∞}{(-1)}^n({e^{-\left(\frac{(2n+1)L-x_m}{2\sqrt{a(t-t_0)}}\right)}}^2-e^{{-\left(\frac{(2n+1)L+x_m}{2\sqrt{a\left({t-t}_0\right)}}\right)}^2})---\left(29\right)$

这给出

$\begin{array}{c}\mathrm{\ρ}{Q}_{\mathrm{mf}}(x_f,t)=\\ (m_f-m_{m\_o})\frac{{\mathrm{MHk}}_m}{2\mathrm{RT}\sqrt{\mathrm{\πa}}}\underset{k=1}{\overset{2}{\mathrm{\Σ}}}\frac{1}{\sqrt{(t-t_{0,k})}}\underset{n=0}{\overset{\∞}{\mathrm{\Σ}}}{(-1)}^n({e^{-\left(\frac{(n+1)L_k}{\sqrt{a(t-t_{0,k})}}\right)}}^2\\ {-e}^{-{\left(\frac{{\mathrm{nL}}_k}{\sqrt{a(t-t_{0,k})}}\right)}^2})\end{array}---\left(30\right)$

可以确定交叉点i和j之间的裂缝分支内的流动。裂缝网络内的流动由下述方程描述：

$\frac{\∂}{{\∂x}_f}\left({\mathrm{\ρQ}}_f(x_f,t)\right)=-\mathrm{\ρ}{Q}_{\mathrm{mf}}(x_f,t)---\left(31\right)$

其中Q_mf是从基质到裂缝的流率(m³/s)，Q_f是来自裂缝的通量(m²/s)。在气体动力可以由所附实际气体方程描述的假设下，该方程变成

$\frac{M}{2\mathrm{RT}}\frac{\∂}{{\∂x}_f}(-\mathrm{HC}\frac{\∂m\left(P_f\right)}{{\∂x}_f})=-\mathrm{\ρ}{Q}_{\mathrm{mf}}(x_f,t)---\left(32\right)$

边界条件如下：

m(P_f(x_f=0))-m_{m_o}=m_f(x_f=0)-m_{m_o}=m_f,i-m_{m_o}=m_f*,i        (33)

m(P_f(x_f=L_f))-m_{m_o}=m_f(x_f=L_f)-m_{m_o}=m_f,j-m_{m_o}=m_f*,j

其中L_f是两个交叉点之间的裂缝的长度。使用下述公式：

m_f*=m_f(x_f)-m_{m_o}                                    (34)

并且将方程(30)引入方程(32)，获得下述方程：

$\frac{\∂}{{\∂x}_f}\left(\frac{{\∂m}_{f*}}{{\∂x}_f}\right)-{\mathrm{\γ}}^2{m}_{f*}=0---\left(35\right)$

其中

${\mathrm{\γ}}^2=\frac{k_m}{C\sqrt{\mathrm{\πa}}}\underset{k=1}{\overset{2}{\mathrm{\Σ}}}\frac{1}{\sqrt{(t-t_{0,k})}}\underset{n=0}{\overset{\∞}{\mathrm{\Σ}}}{(-1)}^n(e^{{-\left(\frac{n{L}_k}{\sqrt{a(t-t_{0,k})}}\right)}^2}-e^{{-\left(\frac{(n+1)L_k}{\sqrt{a(t-t_{0,k})}}\right)}^2})---\left(36\right)$

方程(35)的解具有如下形式

$m_{f*}={\mathrm{Ae}}^{\mathrm{\γ}{x}_f}+B{e}^{-\mathrm{\γ}{x}_f}---\left(37\right)$

并且方程(34)给出

$m_f\left(x_f\right)=m_{m\_o}+\left(\frac{m_{f*,i}{e}^{-\mathrm{\γ}{L}_f}-m_{f*,j}}{e^{-\mathrm{\γ}{L}_f}-e^{\mathrm{\γ}{L}_f}}\right)e^{\mathrm{\γ}{x}_f}+\left(\frac{m_{f*,j}-m_{f*,i}{e}^{\mathrm{\γ}{L}_f}}{e^{-\mathrm{\γ}{L}_f}-e^{\mathrm{\γ}{L}_f}}\right)e^{-\mathrm{\γ}{x}_f---\left(38\right)}$

也可以确定例如分支i,j的交叉点i处的流率。

${\mathrm{\ρQ}}_{f,i,\mathrm{ij}}=-\frac{{\mathrm{MHC}}_{i,j}}{2\mathrm{RT}}\frac{\∂m_f}{{\∂x}_f}{|}_{x_f=0}---\left(39\right)$

当将方程(48)引入方程(39)时

${\mathrm{\ρQ}}_{f,i,\mathrm{ij}}=-\frac{\mathrm{\γ}{\mathrm{MHC}}_{i,j}}{2\mathrm{RT}}\left[\frac{m_{f*,i}(e^{-\mathrm{\γ}{L}_f}+e^{\mathrm{\γ}{L}_f})-2m_{f*,j}}{e^{-\mathrm{\γ}{L}_f}-e^{\mathrm{\γ}{L}_f}}\right]---\left(40\right)$

如果这是来自管子的元件，则该方程变为

$\mathrm{\ρ}{Q}_{f,i,\mathrm{ij}}=-\frac{M{\mathrm{HC}}_{i,j}}{2\mathrm{RT}}\left[\frac{m_{f*,i}-m_{f*,j}}{L_f}\right]---\left(41\right)$

也可以确定裂缝之间的交叉点处的质量平衡

$\underset{j=1}{\overset{N\_\mathrm{ij}}{\mathrm{\Σ}}}\mathrm{\ρ}{Q}_{f,i,\mathrm{ij}}=0---\left(42\right)$

其中N_ij是达到该交叉点的分支的数量。

$\underset{j=1}{\overset{N\_\mathrm{ij}}{\mathrm{\Σ}}}\mathrm{\ρ}{Q}_{f,i,\mathrm{ij}}=\underset{j=1}{\overset{N\_\mathrm{ij}}{\mathrm{\Σ}}}-\frac{{\mathrm{\γ}}_{\mathrm{ij}}{H}_{\mathrm{ij}}{C}_{\mathrm{ij}}M}{2\mathrm{RT}}\left[\frac{(m_{f,i}-m_{m_{o,\mathrm{ij}}})(e^{-\mathrm{\γ}{L}_{f,\mathrm{ij}}}+e^{\mathrm{\γ}{L}_{f,\mathrm{ij}}})-2(m_{f,j}-m_{m_{o,\mathrm{ij}}})}{e^{-\mathrm{\γ}{L}_{f,\mathrm{ij}}}-e^{\mathrm{\γ}{L}_{f,\mathrm{ij}}}}\right]=0---\left(43\right)$

这可以重新整理为如下形式

$\begin{array}{c}{m}_{f,i}\underset{j=1}{\overset{N\_\mathrm{ij}}{\mathrm{\Σ}}}\frac{{\mathrm{\γ}}_{\mathrm{ij}}{H}_{\mathrm{ij}}{C}_{\mathrm{ij}}\left(e^{-\mathrm{\γ}{L}_{f,\mathrm{ij}}+e^{\mathrm{\γ}{L}_{f,\mathrm{ij}}}}\right)}{e^{-\mathrm{\γ}{L}_{f,\mathrm{ij}}}-e^{\mathrm{\γ}{L}_{f,\mathrm{ij}}}}-\underset{j=1}{\overset{N\_\mathrm{ij}}{\mathrm{\Σ}}}{m}_{f,j}\frac{{\mathrm{\γ}}_{\mathrm{ij}}{H}_{\mathrm{ij}}{C}_{\mathrm{ij}}}{e^{-\mathrm{\γ}{L}_{f,\mathrm{ij}}}-e^{\mathrm{\γ}{L}_{f,\mathrm{ij}}}}\\ =\underset{j=1}{\overset{N\_\mathrm{ij}}{\mathrm{\Σ}}}{m}_{m_o,\mathrm{ij}}\frac{{\mathrm{\γ}}_{\mathrm{ij}}{H}_{\mathrm{ij}}{C}_{\mathrm{ij}}(e^{-\mathrm{\γ}{L}_{f,\mathrm{ij}}}+e^{\mathrm{\γ}{L}_{f,\mathrm{ij}}}-2)}{e^{-\mathrm{\γ}{L}_{f,\mathrm{ij}}}-e^{\mathrm{\γ}{L}_{f,\mathrm{ij}}}}\end{array}---\left(44\right)$

还可以更新通过基质块的流体流动的时间函数t_o,k(t)。可以将目标函数F定义为到目前为止从每个基质块的每个面开采的实际质量，和在DFN内的当前压力场已变得恒定并且分析解中考虑的初始时间已变得恒定的情况下将开采的质量，之间的差，并且等于t_o,k(t)。

$\begin{array}{c}{F}_k(t,t_{0,k}\left(t\right))\\ =\underset{0}{\overset{t}{\∫}}\mathrm{\ρ}{Q}_{\mathrm{tot},\mathrm{ij},k}(\mathrm{\τ},P\left(\mathrm{\τ}\right),t_{0,k}\left(\mathrm{\τ}\right))\mathrm{d\τ}\\ -\underset{t_{0,k}\left(t\right)}{\overset{t}{\∫}}{\mathrm{\ρQ}}_{\mathrm{tot},\mathrm{ij},k}(\mathrm{\τ},P\left(t\right),t_{0,k}\left(t\right))\mathrm{d\τ}\end{array}---\left(45\right)$

或者

$F_k(t,t_{0,k}\left(t\right))=M_{\mathrm{tot},k}\left(t\right)-\underset{t_{0,k}\left(t\right)}{\overset{t}{\∫}}\underset{0}{\overset{L_f}{\∫}}\mathrm{\ρ}{Q}_{\mathrm{mf},k}(\mathrm{\τ},P_f(x_f,t),t_{0,k}\left(t\right))\mathrm{d\τd}{x}_f$

可以通过求使得F_o,k(t)=0的t_o,k(t)的值来计算初始时间t_o,k(t)。因此，从分支的k面开采的总质量等于在特定条件(例如，在裂缝中的当前压力条件已变得相同并且从初始开采时间t_o,k(t)起已变得恒定，和/或在初始时间t_o,k(t)之前没有从这一面开采)下到目前为止应当已由同一分支开采的质量。

对t_o,k(t)的值的搜索是通过牛顿迭代算法：

$t_{0,k}^{n+1}\left(t\right)=t_{0,k}^n\left(t\right)-\frac{F_k^n(t,t_{0,k}^n\left(t\right))}{F_k^{\′,n}(t,t_{0,k}^n\left(t\right))}\mathrm{with}{F}_k^{\′,n}(t,t_{0,k}^n\left(t\right))=\frac{\∂F_k^n(t,t_{0,k}^n\left(t\right))}{\∂t_{0,k}^n\left(t\right)}---\left(46\right)$

使用数值梯度来计算导数

其中：

a=变量

c_t=压缩率(Pa^-1)

C=传导率(m².m)

F_o,k=目标函数(m³)

H=裂缝高度(m)

k_m=基质渗透率(m²)

L=泄油区的最大长度(m)

m=实际气体拟压力(Pa/s)

m^*=归一化实际气体拟压力(Pa/s)

m_f=裂缝中的实际气体拟压力(Pa/s)

m_f ^*=裂缝中的归一化实际气体拟压力(Pa/s)

m_m=基质中的实际气体拟压力(Pa/s)

m_m ^*=基质中的归一化实际气体拟压力(Pa/s)

m_{m_0}=基质中的初始实际气体拟压力(Pa/s)

M=摩尔质量(kg/mol)

P_m=基质压力(Pa)

P_m0=初始基质压力(Pa)

P_f=初始裂缝压力(Pa)

P_LB=低基准压力(Pa)

Q_tot=从基质块到裂缝的总流率(m³/s)

Q_mf=从基质块到裂缝的本地流率(m²/s)

Q_f=基质内的流率(m³/s)

R=普适气体常数(J/mol/K)

t=时间(s)

t_o,k=初始开采时间(s)

T=温度(K)

x_f=裂缝中的坐标(m)

x_m=基质中的坐标(m)

Z=体积因子

μ=粘度(Pa.s)

=孔隙度

ρ=储层流体密度(kg/m³)

γ=变量

然后，在t_o,k(t)已知的情况下，可以如下计算压力分布：

$\begin{array}{c}\mathrm{\ρ}{Q}_{\mathrm{mf}}(x_f,t)=(m_f\\ -m_{m\_o})\frac{{\mathrm{MHk}}_m}{2\mathrm{RT}\sqrt{\mathrm{\πa}}}\underset{k=1}{\overset{2}{\mathrm{\Σ}}}\frac{1}{\sqrt{(t-t_{0,k})}}\underset{n=0}{\overset{\∞}{\mathrm{\Σ}}}{(-1)}^n({e^{-\left(\frac{(n+1)L_k}{\sqrt{a(t-t_{0,k})}}\right)}}^2\\ {-e}^{-{\left(\frac{{\mathrm{nL}}_k}{\sqrt{a(t-t_{0,k})}}\right)}^2})\end{array}---\left(47\right)$

这个解在压力方面是线性的，方程2可以被积分并求解。

$m_f\left(x_f\right)=m_{m\_o}+\left(\frac{m_{f*,i}{e}^{-\mathrm{\γ}{L}_f}-m_{f*,j}}{e^{-\mathrm{\γ}{L}_f}-e^{\mathrm{\γ}{L}_f}}\right)e^{\mathrm{\γ}{x}_f}+\left(\frac{m_{f*,j}-m_{f*,i}{e}^{\mathrm{\γ}{L}_f}}{e^{-\mathrm{\γ}{L}_f}-e^{\mathrm{\γ}{L}_f}}\right)e^{-\mathrm{\γ}{x}_f---\left(48\right)}$

${\mathrm{\γ}}^2=\frac{k_m}{C\sqrt{\mathrm{\πa}}}\underset{k=1}{\overset{2}{\mathrm{\Σ}}}\frac{1}{\sqrt{(t-t_{0,k})}}\underset{n=0}{\overset{\∞}{\mathrm{\Σ}}}{(-1)}^n(e^{{-\left(\frac{n{L}_k}{\sqrt{a(t-t_{0,k})}}\right)}^2}-e^{{-\left(\frac{(n+1)L_k}{\sqrt{a(t-t_{0,k})}}\right)}^2})---\left(49\right)$

已知网络内的压力分布，可以利用达西定律计算开采速度。

可以在对时间步没有约束的情况下对各种裂缝网络进行计算。在一些情况下，例如在以恒定BHP进行开采并且传导率高的情况下，来自基质的流动可以基于裂缝内的压力保持恒定的假设。但是实际上，只有裂缝分支的一部分可以具有高传导率。计算的例子在Cipolla,C.L.,Lolon,E.P.,Mayerhofer,M.J.,"Reservoir Modeling and Production Evaluation in shale-Gas Reservoirs",2009年12月7日在卡塔尔多哈召开的International Petroleum TechnologyConference(国际石油技术会议)上提交的SPE论文13185中给出。

分析解的有效性

可以通过对解进行分析以确定其在给定地层中的应用，来验证427分析解。为了研究针对裂缝传导率的不同值的分析模型的有效性，可以分析复杂网络的单个裂缝分支处的压力和开采的演化。这一研究可以包括两组等间距的平行裂缝，如图10所示。该附图描述了将要分析的井眼1074附近的DFN1072中的单个分支1070。DFN1072的基质块1048被描绘为具有从基质块1048到裂缝分支1070的流率1076。

对于在大约0.0001mD的储层中的裂缝网络中的高传导率(有限)(例如，大约2500mD.ft(762mD.m))，BHP几乎瞬时在网络中扩散，并且从那里开始，与初始储层压力和BHP之间的压力差相比，可以忽略DFN内部的压力变化。

图11.1和11.2是描绘对于1天和365天，储层压力P(z轴)与距离x(m)(x轴)和距离y(m)(y轴)之间的关系的三维曲线图1100.1和1100.2。该附图描绘了对于高传导率DFN在两个不同的开采时间的DFN压力和初始储层压力1178。可以提供这些和其它描绘。可以基于开采估计来调整开采作业。

如图12所示，所选裂缝分支(例如图10的分支1070)内的压力在十年的开采期间可以被视为恒定。该附图描绘了在高传导率的情况下，压力(P_m,o-P_f)(左边的y轴)和按天计的初始时间T(右边的y轴)相对于按天计的时间t(例如在三年的开采期间)(x轴)的关系的曲线图1200。针对压力产生的线1280和针对时间延迟产生的线1281几乎是平的。

图13示出了DFN中这一几乎恒定的压力的结果，其中(从基质块到裂缝分支(例如图10中的1048到1070))开采的累积体积朝着由质量平衡限定的最大可采体积(或适当位置处的初始气体体积)收敛。图13是描绘累积开采量P(y轴)相对于时间t(x轴)的关系的曲线图1300，结果是朝着最大可采体积1382延伸的开采曲线1384。该附图描绘了在高传导率的情况下从裂缝分支的累积开采量相对于时间的关系。因为在此例中，我们在考虑可压缩流体所以可以在地表条件下进行体积的测量。这一收敛指示在高传导率DFN的情况下分析解验证了质量平衡。对于低传导率(或有限传导率)DFN(50mD.ft(15.24mD.ft))进行相同分析可以具有不同的结论。

如图14.1和14.2所示，与该问题的压力范围(例如，BHP，初始储层压力等)相比，DFN中的压力可以改变。这个图描绘了对于低传导率DFN，在不同的开采时间，DFN内部的压力。图14.1和14.2是分别描绘对于1天和365天，储层压力P(z轴)相对于距离x(x轴)和距离y(y轴)的关系的三维曲线图1400.1和1400.2。该附图描绘了对于高传导率DFN，在不同的开采时间，DFN内部的压力。还描绘了初始储层压力1478和DFN的压力1435。

如图15所示，在裂缝分支中记录的压力相对于时间的关系上可以看到这一压力变化。如图15所示，在十年的开采期内可以将所选裂缝分支(例如图10的分支1070)中的压力视为恒定。该附图描绘了在低传导率(无限)的情况下，压力(Pm,o-Pf)(左边的y轴)和按天计的时间延迟T(右边的y轴)随按天计的时间t(例如三年开采期内)(x轴)的变化的曲线图1500。对于归一化压力产生的线1580和对于时间延迟产生的线1581几乎是平的。还描绘了边界条件的变化1584。

DFN中压力的这一变化意味着分析解中的恒定压力边界条件的假定可能需要进一步的分析以确认有效性。结果是所计算的来自基质的流率可能被低估，并且质量平衡可能不正确，如图16所示。图16是描绘累积开采量P(y轴)相对于时间t(x轴)的关系的曲线图1600，结果是朝着最大可采体积1682延伸的开采曲线1684。这个附图描绘了在低传导率的情况下从裂缝分支的累积开采量相对于时间的关系。还示出了开采曲线1684和最大可采体积1682之间的误差1686。

裂缝网络中的低扩散率可能导致块的开采的“延迟”，这取决于它离井眼有多远(或者它如何与井眼相连)。这一观察是将分析解的有效性延伸到低传导率裂缝的方法的起点。

分析解的有效性的扩展

可以扩展这一分析解的有效性，以便以使得到目前为止从基质块开采的体积等于在DFN内当前压力条件下根据分析解应当已开采的体积的方式，修改“初始”时间t₀(t_o,k(t)。通过在每一时间步对每个裂缝的每一侧进行这一搜索，迫使分析解满足质量平衡。通过定义目标函数F以使其最小化来开始对于t₀的搜索。

$F_k(t,t_{0,k}\left(t\right))=M_{\mathrm{tot},k}\left(t\right)-\underset{0}{\overset{t}{\∫}}\underset{0}{\overset{L_f}{\∫}}\mathrm{\ρ}{Q}_{\mathrm{mf},k}(\mathrm{\τ},P_f(x_f,t),t_{0,k}\left(t\right))\mathrm{d\τd}{x}_f=0---\left(50\right)$

M_tot是在时间t从裂缝分支的k侧从基质块开采的体积。将其与在该裂缝分支的长度上来自基质的流率从t_0,k到t的积分进行比较。对使得F等于零的t_0,k的搜索，可以使用方程51所描述的牛顿-拉普森迭代算法。

$t_{0,k}^{n+1}\left(t\right)=t_{0,k}^n\left(t\right)-\frac{F_k^n(t,t_{0,k}^n\left(t\right))}{F_k^{\′,n}(t,t_{0,k}^n\left(t\right))}$ ，其中 $F_k^{\′,n}(t,t_{0,k}^n\left(t\right))=\frac{\∂F_k^n(t,t_{0,k}^n\left(t\right))}{\∂t_{0,k}^n\left(t\right)}---\left(51\right)$

通过数值梯度来计算函数F_0,k的导数。如果t_0,k满足其时间边界，那么优化使用两分法。这一优化算法非常有效，因为使用前一时间步的解作为下一循环的初始猜想。从数值的观点看，近似值的计算需要时间积分，这是该模拟中最耗CPU的部分。以变量间最小依赖性对每个分支的每一侧应用该优化算法，使得算法的这一部分成为平行计算的候选。

为了阐明这一方案背后的机制，可以对具有低传导率(或有限传导率)(50mD.ft(15.24mD.ft))的DFN的单个裂缝分支使用上面的分析。在裂缝分支中记录的压力中可以看到这一压力变化，如图17所示。如图17所示，在十年的开采期内可以将所选裂缝分支(例如图10的分支1070)中的压力视为恒定。该附图描绘了在低传导率(无限)的情况下，归一化压力(Pm,o-Pf)(左边的y轴)和按天计的时间延迟T(右边的y轴)随按天计的时间t(例如三年开采期内)(x轴)的变化的曲线图1700。对于归一化压力产生的线1780和对于时间延迟产生的线1781是倾斜的。

图17还示出了裂缝中所计算的压力和用所提出的方法更新的初始时间t_0,k。t_0,k随时间的增长对于维持来自基质的流率和图18所示的累积开采量可能是必要的。图18是描绘累积开采量(y轴)相对于时间(x轴)的关系的曲线图1800，结果是开采曲线1884朝着最大可采体积1882延伸。这个附图描绘了在低传导率的情况下从裂缝分支的累积开采量相对于时间的关系。

这个附图指示该方法减少了质量平衡的误差，因为与图16相比，图18中累积开采量收敛得更接近最大可采体积，由此指示可以使用分析解扩展该方法的有效性。

图19是描绘对于整个裂缝网络，在不同的时间步t₁(1天)、t₂(200天)和t₃(3年)，由该算法计算的压力P的分布和初始时间延迟T的图表1900。该图表包括分别在时间步t₁、t₂和t₃，针对压力的DFN1935.1、1935.2和1935.3以及针对时间延迟的DFN1935.4、1935.5和1935.6。这个附图示出了在开采的不同时间，储层中的压力和初始时间(或“延迟”)。“压力”柱示出了储层块内的压力和裂缝网络内的压力。“初始时间”T柱示出了由该算法计算的每个块的初始时间。

可以使用非常规开采模型(UPM)来执行上面的分析。为了图解UPM的性能，将该模拟与来自商业储层模拟器的模拟进行比较。这一比较是使用两个不同的裂缝几何图形来完成的：单个双翼裂缝网络和“线网”裂缝网络。

在涉及单个双翼裂缝的例子中，水力裂缝是单个对称裂缝，半段的长度为1263ft(384.96m)，裂缝高度为98.4ft(19.99m)。储层的渗透率为0.0001mD，孔隙度为8%，初始储层压力为4000psi(281.29kg/cm)，并且井底压力为1000psi(70.32kg/cm)。在这个例子中，体积因子Z和气体粘度是恒定的，并且分别等于1和0.02cP。

图20是对于在0.005和5000mD.ft(1524mD.m)之间变化的不同裂缝传导率，储层模拟和UPM之间的模拟累积开采量比较。图20表明，与穿孔(格子的中心)的距离越大，初始时间越小。

图20是地表条件下的累积开采量(y轴)相对于时间t(x轴)的关系的曲线图2000。该附图通过与储层模拟器的比较描绘了有效性。所得实线2088.1-2088.7和所得虚线2089.1-2089.7分别示出了在各种位置处，基于储层模拟器和该开采模型的开采量。这一曲线图2000指示，与穿孔的距离越大，BHP扩散到该位置要花费的时间越长。

对于线网裂缝网络，这种情况代表由每个正交方向上13条相同的裂缝构成的复杂裂缝网络，垂直井位于中间。在这个例子中，储层的渗透率为大约0.001mD，孔隙度为8%，初始储层压力为4000psi(281.29kg/cm)，并且井底压力为1000psi(70.32kg/cm)。在这个例子中还示出，体积因子Z和气体粘度是恒定的，并且分别等于1和0.02cP。

图21.1和21.2提供了由各种模拟器执行的DFN的各种显现。该附图显示了在通过商用储层模拟器和UPM进行的模拟之间的比较中所使用的储层和DFN。图21.1示出了由储层模拟器如ECLIPSE^TM描绘的DFN2135.1和2135.2的例子。图21.2示出了使用UPM产生的DFN2135.3。所描绘的每个DFN可以是相同的DFN，导致所示出的不同图像。

图22-24比较在DFN的传导率可以变化的例子中由储层模拟器和UPM产生的结果。图22是针对双翼裂缝，对于在0.082mD.ft(24.99mD.mm)到大约8200mD.ft(2499.36mD.m)之间变化的不同裂缝传导率，储层模拟器和UPM之间的模拟累积开采量的比较。

图22显示，离穿孔(格子的中心)距离越大，初始时间越小。图22是地表条件下的累积开采量(y轴)相对于时间t(x轴)的关系的曲线图2200。该附图描绘了通过与储层模拟器的比较进行的验证。所得实线2288.1-2288.6和所得虚线2289.1-2289.6分别示出了在各种位置处，基于储层模拟器和UPM的开采量。这一曲线图2200指示，与穿孔的距离越大，BHP扩散到该位置要花费的时间越长。

如这里所使用的，没有“延迟”的UPM意味着UPM使用模型的分析部分，具有等于0的恒定初始时间。当裂缝传导率增加时，储层模拟器和没有“延迟”的UPM模拟之间的差异可以减小。

这些比较显示了特别是在低传导率的情况(更新初始时间的算法发挥主要作用)下两种模拟器之间相当良好的一致。为了图解该算法校正初始时间的重要性，图23.1和23.2比较裂缝传导率等于82mD.ft(24.99mD.m)的情况下的模拟结果。

图23.1是描绘地表条件下的流率的曲线图2300.1。相对于时间t(x轴)绘出了开采量(y轴)。所得到的线2390.1-2390.3分别描绘了由储层模拟器、UPM和无延迟的UPM产生的模拟。图23.2是描绘地表条件下的当前开采量的曲线图2300.2。相对于时间t(x轴)绘出了累积开采量P(y轴)。所得到的线2390.4-2390.6分别描绘了由储层模拟器、UPM和无延迟的UPM产生的模拟。这些附图描绘了商用储层模拟器、UPM和无“延迟”UPM之间的流率比较(图23.1)和累积开采量比较(图23.2)。

应当注意的是，在开发任何实际实施例时，必须做出大量特定于实现的决定，以实现开发者的特定目的，例如符合与系统相关的约束和与商业相关的约束，这些约束将对于每个实现而不同。此外，应当理解，这样的开发工作可能是复杂的，并且耗时的，然而对于本领域技术人员在本公开的启示下将只是例行的工作。此外，这里使用/公开的文章还可以包括一些它们引用的组分之外的组分。在发明内容部分及此详细描述中，每个数值应当被理解为由术语“大约”修饰一次(除非已明确地如此修饰)，然后再理解为没有如此修饰，除非上下文中相反指示。同样，在发明内容和此详细描述中，应当理解，所列出或描述为有用、适合等的浓度范围意在将该范围内的任何和每一个浓度，包括端点，考虑为已被陈述。例如，“范围从1到10”应当被理解为指示沿着大约1和大约10之间的连续统的每个可能数值。因此，即使该范围内明确识别了特定数据点，或者该范围内没有明确识别数据点，或者仅提到几个特定项目，应当明白，发明人认为该范围内的任何及所有数据点应当被考虑为已被指定，并且发明人掌握整个范围以及该范围内所有点的知识。

这里做出的陈述仅仅提供与本公开有关的信息，而不构成现有技术，并且可以描述一些阐述所公开的主题的实施例。这里所引用的所有参考文献都通过引用而将其全部内容合并到本申请中。

已经参考一些实施例给出了前面的描述。本公开所属领域的技术人员将理解，可以在不实质性脱离本申请的原理、范围的情况下，实现所描述的结构和作业方法的变更和改变。相应地，前面的描述不应当被解释为仅属于所描述并在附图中示出的精确结构，而是应当被理解为与所附权利要求一致，并支持所附权利要求，权利要求具有其最完整最公平的范围。

尽管上面仅描述了几个示例性实施例，但是本领域技术人员将容易地理解到，在示例性实施例内可能有许多修改而不实质上脱离本发明。相应地，所有这样的修改都意图包括在本公开如在所附权利要求书中所限定的范围之内。在权利要求书中，装置加功能的条款意图覆盖这里描述的执行所阐述的功能的结构，并且不但包括结构等同物，而且包括等同结构。因此，在固定木头部件的场境中，尽管钉子和螺丝钉可能不是结构等同物，因为钉子采用圆柱形表面来将木头部件固定在一起，而螺丝钉采用螺旋表面，但是钉子和螺丝钉可以是等同结构。除了权利要求明确使用表述“装置，用于……”和相关联的功能一起做出的限制，申请人明确表示不希望引用35U.S.C§112第6款来对这里的任何权利要求做出任何限制。

Claims (34)

1.一种在穿透地下地层的井眼附近执行开采作业的方法，该地下地层附近具有多个裂缝，该方法包括：

产生通过离散裂缝网络的流率，所述离散裂缝网络是从由地下地层中的多个裂缝限定的水力裂缝网络外推出来的，所述离散裂缝网络包括多个裂缝分支和多个基质块，所述多个裂缝分支之间具有交叉点；

基于所述流率产生离散裂缝网络在初始时间的压力分布情况；以及

基于所述压力分布情况产生开采速度。

2.根据权利要求1的方法，其中产生流率包括产生从所述多个基质块中的一个到所述多个裂缝分支中的一个的流率。

3.根据权利要求1的方法，其中产生流率包括产生所述多个裂缝中的至少一个内部的流率。

4.根据权利要求1的方法，其中产生流率包括产生从所述多个基质块中的一个到一个裂缝分支的流率。

5.根据权利要求1的方法，其中产生流率包括产生一个裂缝分支内部在所述离散裂缝网络的两个交叉点之间的流率。

6.根据权利要求1的方法，其中产生流率包括产生所述多个裂缝中的一个的内部在所述离散裂缝网络的交叉点处的流率。

7.根据权利要求1的方法，还包括确定所述多个裂缝分支中的两个之间的交叉点处的质量平衡。

8.根据权利要求1的方法，其中产生压力分布情况包括使用达西定律产生所述压力分布情况。

9.根据权利要求1的方法，其中产生压力分布情况不受时间步的约束。

10.根据权利要求1的方法，还包括定义通过基质块的流体流动的时间函数，所述时间函数具有所述初始时间。

11.根据权利要求10的方法，还包括更新所述通过基质块的流体流动的时间函数。

12.根据权利要求1的方法，还包括在多个时间步更新所述开采速度。

13.根据权利要求1的方法，还包括为所述多个裂缝分支更新所述开采速度。

14.根据权利要求1的方法，还包括通过更新所述初始时间来补偿所述多个基质块中的每一个的开采的延迟，使得从所述多个基质块中的每一个开采的实际质量等于在所述多个裂缝分支中的邻近的一个的当前压力条件在从更新后的初始时间开始的时间里已变得恒定的条件下的质量。

15.根据权利要求1的方法，还包括验证所述开采速度。

16.根据权利要求15的方法，其中所述验证包括比较所述开采速度和通过储层模拟器产生的开采速度。

17.根据权利要求15的方法，其中所述验证是针对具有高传导率、低传导率、双翼裂缝、线网裂缝、时间延迟及它们的组合的离散裂缝网络执行的。

18.根据权利要求15的方法，其中所述验证包括修改所述初始时间，使得对于所述多个裂缝分支中的每一个，随着时间流逝而从所述多个基质块开采的体积满足质量平衡。

19.一种在穿透地下地层的井眼附近执行油田作业的方法，该地下地层中具有储层，该方法包括：

执行压裂作业，所述压裂作业包括在井眼附近产生裂缝，所述裂缝限定所述井眼附近的水力裂缝网络；

通过从所述水力裂缝网络外推裂缝数据，在所述井眼附近产生离散裂缝网络，所述离散裂缝网络包括多个裂缝分支和多个基质块，所述多个裂缝分支之间具有交叉点；

产生通过所述离散裂缝网络的泄油区的深度；

限定至少一个开采参数；以及

基于所述泄油区的深度和所述至少一个开采参数，执行开采作业，以从所述地下地层开采流体。

20.根据权利要求19的方法，还包括测量所述井眼附近的井下数据。

21.根据权利要求19的方法，其中所述执行压裂作业包括通过向所述地下地层注入流体来对从所述井眼的开采进行增产。

22.根据权利要求19的方法，其中所述执行压裂作业包括对执行所述压裂作业进行模拟。

23.根据权利要求19的方法，其中所述离散裂缝网络考虑所述多个裂缝分支中的每一个处的至少一个裂缝属性的平均值。

24.根据权利要求23的方法，其中所述至少一个裂缝属性包括裂缝分支末端处的空间坐标；裂缝分支位置处的传导率、平均传导率、高度、平均高度、储层压力、平均储层压力；所述裂缝分支位置处的渗透率、平均储层渗透率；及它们的组合。

25.根据权利要求19的方法，其中产生泄油区的深度包括评价通过所述离散裂缝网络的所述多个基质块的泄油区的深度。

26.根据权利要求19的方法，其中产生泄油区的深度包括基于通过所述多个基质块的线性流动近似来产生用于所述多个基质块中每一个的泄油区的深度。

27.根据权利要求19的方法，其中产生泄油区的深度包括自动地评价在所述多个裂缝分支中的每一个的前端即将被耗尽的所述多个基质块的泄油区的深度，并且对于所述多个基质块中的每一个补偿要耗尽的体积。

28.根据权利要求19的方法，其中所述至少一个开采参数包括井底压力、储层条件下的储层流体粘度、储层条件下的储层流体压缩率、要模拟的开采的持续时间以及它们的组合。

29.根据权利要求19的方法，其中所述执行开采作业包括在井眼中放置管子，并且将流体从储层传输到地表位置。

30.根据权利要求19的方法，其中所述执行开采作业包括通过模拟从所述井眼的流体开采，估计从所述井眼的开采速度。

31.根据权利要求30的方法，其中所述执行开采作业包括显现开采速度。

32.根据权利要求30的方法，还包括基于所估计的开采速度调整所述开采作业的执行。

33.根据权利要求19的方法，其中基于裂缝参数的范围执行开采作业。

34.一种在穿透地下地层的井眼附近执行油田作业的方法，该地下地层中具有储层，该方法包括：

通过将流体注入所述地下地层以使得在所述井眼附近产生裂缝，来对所述井眼进行增产；

测量所述裂缝，并基于所测量的裂缝来限定水力裂缝网络；

通过从所述水力裂缝网络外推裂缝数据，在所述井眼附近产生离散裂缝网络，所述离散裂缝网络包括多个裂缝分支和多个基质块，所述多个裂缝分支之间具有交叉点；

产生通过所述离散裂缝网络的泄油区的深度；

限定至少一个开采参数；

基于所述泄油区的深度和所述至少一个开采参数，估计随时间的开采速度；以及

基于所估计的开采速度从所述地下地层开采流体。

Images