CN103548032B - 使用微震事件数据的水力裂缝表征的系统和方法 - Google Patents

使用微震事件数据的水力裂缝表征的系统和方法 Download PDF

Info

Publication number
CN103548032B
CN103548032B CN201180070143.4A CN201180070143A CN103548032B CN 103548032 B CN103548032 B CN 103548032B CN 201180070143 A CN201180070143 A CN 201180070143A CN 103548032 B CN103548032 B CN 103548032B
Authority
CN
China
Prior art keywords
plane
towards
data point
group
microseism
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.)
Expired - Fee Related
Application number
CN201180070143.4A
Other languages
English (en)
Other versions
CN103548032A (zh
Inventor
马建福
肯尼斯·E·威廉姆斯
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.)
Landmark Graphics Corp
Original Assignee
Landmark Graphics Corp
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 Landmark Graphics Corp filed Critical Landmark Graphics Corp
Publication of CN103548032A publication Critical patent/CN103548032A/zh
Application granted granted Critical
Publication of CN103548032B publication Critical patent/CN103548032B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/288Event detection in seismic signals, e.g. microseismics
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/10Aspects of acoustic signal generation or detection
    • G01V2210/12Signal generation
    • G01V2210/123Passive source, e.g. microseismics
    • G01V2210/1234Hydrocarbon reservoir, e.g. spontaneous or induced fracturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • User Interface Of Digital Computer (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Complex Calculations (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
  • Lubricants (AREA)

Abstract

使用微震事件数据来识别地下裂缝的朝向、间距以及倾角的水力裂缝表征的系统和方法。

Description

使用微震事件数据的水力裂缝表征的系统和方法
相关申请的交叉引用
不适用。
有关联邦资助研究的声明
不适用。
技术领域
本发明大体涉及使用微震(microseismic)事件数据的水力裂缝(fracture)表征的系统和方法。更具体地,本发明涉及使用微震事件数据来识别地下裂缝的朝向、间距以及倾角的水力裂缝表征。
背景技术
微震事件数据有时被获取作为为了开发油气储量目的而被钻探的井的水力压裂处理的一部分。该数据是从与由水力压裂诱导的裂缝相关的低振幅(微)地震事件获取的。为了提高油气储层的产能,水力压裂的目标是要通过将高压流体和支撑剂注入岩石基质中而在地下诱导人工裂缝。
微震事件位置通常被实时地监控,并且当其发生时,可以验证显示成三维(3D)视图的事件的位置。其也可适用于在水力压裂处理的结论得出之后进行分析,并因此可用来与该地区中的其它井的结果进行比较。微震事件通常沿着地下裂缝或在地下裂缝附近发生,该地下裂缝可以是诱导裂缝或通过水力压裂处理重新开启的先前存在的天然裂缝。当应力朝向与现在的不同时,裂缝的朝向会受到现今的应力型式以及过去不同时期产生的裂缝系统的存在的强烈影响。
被探测和分析的每一个独立各异的微震事件是井下裂缝的结果,该井下裂缝具有能够从倾斜仪(tiltmeter)或地震传感器数据提取的朝向、震级(magnitude)、位置以及其它属性。例如,该裂缝可以利用诸如长度、宽度、高度以及压力等其它参数来表征。存在与每一个微震事件相关的位置不确定性。沿x-y方向的这一不确定性与沿垂直(z)深度域的不同。每一个事件的位置不确定性可以由长椭球(prolate spheroid)来表示。
在一些情况下,微震事件存在明显的朝向和间距,它们遵循在裂缝分析的数学描述中经常使用的典型双翼裂缝概念。在其它情况下,表示包含所有微震数据点的3D体积的密集数据云是诱导裂缝或重新激活裂缝的复杂裂缝型态的证据。在这些情况下,微震数据的分析变得非常具主观性和解释性。即使在这些情况下,在数据云内也存在可能代表存在于地下的裂缝型态的型态。
当今的应力场可能与原始裂缝创建时的不同。诱导的水力裂缝的现今朝向受到地下应力状态的强烈影响。在垂直应力与两个水平应力之间总是存在一定程度的应力各向异性。各向异性越大,由水力压裂模拟诱导的裂缝越平坦且将越拟合传统的双翼模型。岩石的渗透率越大,裂缝将越平坦。应力型式越各向同性,裂缝越容易能够由于岩石中的不连续性而偏转,且越容易能够创建复杂的裂缝网络。
当前,存在一些用于尝试和识别诱导裂缝和天然裂缝的朝向、倾角以及间距的裂缝表征技术。
在一项技术中,识别微震数据点的整个数据云,以构建模拟储层体积(stimulated reservoir volume)SRV或估计模拟体积(estimated stimulatedvolume)ESV。该信息被推断以作为对通过流体和支撑剂模拟的岩石的量的测量。然而,被泵入地层的仅一小部分能量在地表被接收作为可探测的微震事件。事实上,存在许多位于仪器的探测范围以下的微震事件。因此,数据云维数的使用只是一个引导,并且这种裂缝表征技术的结果可能并不是有效的。基于微震事件数据的型态解释裂缝型态是另一项技术,该项技术在M.K.Fisher等人的“使用微震裂缝映射的巴耐特页岩中的优化水平完井技术(Optimizing Horizontal Completion Technologies in the Barnett Shale UsingMicroseismic Fracture Mapping)”中得到描述。地表地震数据(主动和被动)分析是用于裂缝表征的又一项技术。在一些情况下,与地下裂缝的存在有关的速度的横向变化可能存在差异。沿一些方向的该速度比沿其它方向的快。转换波(例如,P到S)也已被用作裂缝表征技术。
还存在几种不同的裂缝表征技术,这些技术能够将微震事件数据与地下的模型数学地关联起来,并产生例如离散裂缝网络(discrete fracture network,DFN)或一组大概的裂缝表征,该技术例如为在第2010/0307755和0211/0029291号美国专利申请公开中描述的技术。
在M.J威廉姆斯等人的“来自微震事件位置的主平面的定量解释及在产能预测中的应用(Quantitative Interpretation of Maj or Planes From MicroseismicEvent Locations with Application in Production Prediction)”中提出了自动裂缝表征技术,以对数据云进行数值分析,并识别来自微震数据点的主平面,用于储层建模。这一技术着手处理从概率框架识别裂缝趋势的问题,然而,其仅识别单个裂缝朝向。
发明内容
因此,通过提供使用微震事件数据来识别地下裂缝的朝向、间距以及倾角的水力裂缝表征的系统和方法,本发明满足上述需要并克服了现有技术中的一个或多个缺陷。
在一个实施例中,本发明包括一种使用微震数据点的水力裂缝表征的方法,包括:i)使用计算机处理器识别与一组或多组平面及其朝向不确定性有关的平均微震数据点的主要走向朝向;ii)考虑到每一个微震数据点与每一组平面内每一个平面之间的接触加权来确定每一组平面的倾角;iii)确定每一组平面内每一个平面的位置;以及iv)对于具有共同的主要走向朝向的每一组平面执行统计分析,以确定与各组地下裂缝的朝向、间距以及倾角有关的不确定性。
在另一个实施例中,本发明包括一种有形地承载计算机可执行指令的非瞬态程序载体装置,用于使用微震数据点的水力裂缝表征,该指令可被执行以实施:i)识别与一组或多组平面及其朝向不确定性有关的平均微震数据点的主要走向朝向;ii)考虑到每一个微震数据点与每一组平面内每一个平面之间的接触加权来确定每一组平面的倾角;iii)确定每一组平面内每一个平面的位置;以及iv)对于具有共同的主要走向朝向的每一组平面执行统计分析,以确定与各组地下裂缝的朝向、间距以及倾角有关的不确定性。
通过各个实施例及相关附图的以下说明,本发明的其它方面、优点以及实施例将变得对本领域技术人员显而易见。
附图说明
下文参考附图描述本发明,其中类似的元件引用类似的附图标记,并且其中:
图1A为示出用于实施本发明的方法的一个实施例的流程图。
图1B为图1A中示出的方法的继续。
图2A示出用于图1A的步骤108中每一个朝向的平均微震数据点的显示的一个格式。
图2B示出用于图1A的步骤108中每一个朝向的平均微震数据点的显示的另一个格式。
图2C示出用于图1A的步骤108中每一个朝向的平均微震数据点的显示的又一个格式。
图3A示出图2A中的格式用于在使用蒙特卡洛方法的标准统计分析完成之后的显示。
图3B示出图2B中的格式用于在使用蒙特卡洛方法的标准统计分析完成之后的显示。
图3C示出图2C中的格式用于在使用蒙特卡洛方法的标准统计分析完成之后的显示。
图4示出代表图1B的步骤113中一组地下裂缝的一些垂直倾角平面的最佳平面朝向和间距的示例性二维地图。
图5为示出图4中所示三个族群平面的统计的表格。
图6为示出用于实施本发明的系统的一个实施例的方框图。
具体实施方式
虽然具体地描述本发明的主题,然而,说明书本身不意欲限制本发明的范围。因此,该主题也可以结合其它技术以其它方式具体实施,以包括不同的步骤或与这里所描述的相似的步骤的组合。此外,虽然术语“步骤”在这里可以用于描述所采用方法的不同要素,但该术语不应解释为隐含这里公开的不同步骤之中或之间的任何特定顺序,除非说明书清楚地限制到特定顺序。虽然以下说明指的是石油和天然气产业,然而本发明的系统和方法不限于此,并且也可以被应用于其它产业以得到类似的结果。
由水力压裂诱导的单独的微震事件产生可以用来识别存在于地下的水压诱导裂缝和天然裂缝的朝向、间距以及倾角(dip)的微震事件数据。平面被限定并被引入三维体积中。限定一平面需要最少三个微震事件交叉点。可以计算每一个平面的交叉点及各自的长椭球,其表示从微震事件数据的每一个识别的中心点沿x、y以及z方向的位置不确定性。以系统的方式经由输入数据旋转这些平面,使得微震事件数据的朝向、倾角以及间距能够识别拟合输入数据的最有可能的一组平面。任何裂缝可以与任何的朝向相关,并且可以处于两个或多个族群的裂缝的交叉点处。
随沿着平面存在的事件的数量增加,通过微震事件数据识别地下实际裂缝平面的可能性增加。由于多个地质事件可能在任何给定地区产生了几组重叠的裂缝网络,因而不假设先前存在的裂缝网络的数量。以下方法解决了如何识别每一个单独的裂缝和裂缝族群的最有可能的朝向、间距以及倾角,该方法也可以用来识别存在于微震事件数据中的任何相关性。
方法说明
现在参照图1,其为示出用于实施本发明的方法100的一个实施例的流程图。
在步骤101中,使用参考图6描述的客户端接口/视频接口输入微震事件数据。从传感器收集微震事件数据作为天然诱导的或由水力压裂诱导的微震事件的结果。因此,微震事件数据可以包括可能表示与一个或多个井有关的独立水力压裂事件的多个数据集。微震事件数据包括微震数据点和与每一个微震数据点的x、y、z位置有关的信息以及与每一个x、y、z位置相关的位置不确定性。该微震事件数据还可以包括与每一个微震数据点有关的矩震级(moment magnitude)和震矩张量(moment tensor)。对于每一个微震数据点的位置不确定性,该位置不确定性可以包括每一个微震数据点的单独的x、y、z不确定值或对可以被应用于所有微震数据点的这些不确定值的估计。可选地,与传感器的位置及其质量和速度有关的其它数据可以被包括作为微震事件数据的一部分。该微震事件数据可以来自任何可用的数据源,该数据源包括但不限于一个或多个监控井中的倾斜仪、地表传感器或井下传感器。此外,该微震事件数据可以通过使用钻孔数据加以补充,该钻孔数据包括但不限于岩芯、偶极(dipole)、声波测井或图像测井。例如,朝向的偶极声波测井能够表明垂直方向与最大和最小水平方向之间的应力差异。在以下步骤中,为了方法100的目的,可以包括和考虑其它数据。
在步骤102中,使用本领域熟知的技术和来自步骤101的微震事件数据识别所有微震数据点的位置及其各自的位置不确定性。
在步骤103中,为平面设定朝向(例如,走向角和倾角)。例如,该走向角(strike angle)和倾角(dip angle)可以被分别设定(初始化)为0°和0°、或360°和0°、或0°和90°、或360°和90°。
在步骤104中,创建一平面,该平面具有法向量cosθ),并穿过用于确定该平面的朝向(走向、倾角)的所有微震数据点的平均位置。还使用本领域熟知的技术计算每一个微震数据点与平面之间的垂直距离。Psi表示走向角,theta(θ)表示倾角。可以通过 d = ( x ‾ - x ) sin θ cos φ + ( y ‾ - y ) sin θ sin φ + ( z ‾ - z ) cos θ 计算每一个微震数据点与平面之间的垂直距离。因此,每一个微震数据点与平面之间的垂直距离具有在平面的一侧为正且在平面的另一侧为负的符号。每一个微震数据点与符号为(±)的垂直距离相关。因此,不同数量的微震数据点将与位于每一个带符号的距离的可能平面相关。这些微震数据点将归类成组,并且可以与可能的裂缝平面相关。
在步骤105中,使用本领域熟知的技术创建与步骤104中创建的平面平行的一组平面,其中每一个新的平面都与由每一个微震数据点的位置不确定性表示的体积相交。
在步骤106中,计数在步骤105中创建的平面的数量以及接触这些平面的微震数据点的数量。还使用本领域熟知的技术来计算用于确定平面朝向的接触每一个平面的平均微震数据点。
在步骤107中,将步骤104中创建的平面的走向角或倾角增量1°,并且重复步骤104-106直到走向角等于0°或360°且倾角等于0°或90°。在重复步骤104-106直到走向角等于0°或360°且倾角等于0°或90°之后,方法100进行到步骤108。
在步骤108中,在每一个朝向上显示步骤106中计算出的平均微震数据点,以使用本领域熟知的技术识别涉及一组平面的平均微震数据点的主要走向朝向及其朝向不确定性。朝向不确定性也被显示用于分析。以此方式,详尽地处理、分析以及显示步骤104-107中不同组合的走向角和倾角的结果。虽然可以在优先基础上考虑很多格式以在本步骤中显示该数据,然而方法100不限于特定类型的格式,并且图2A-图2C和图3A-图3C示出可以被显示的不同格式。
在图2A中,示出了用于步骤108中每一个朝向的平均微震数据点的显示的一个格式。图2A中的格式是每一个朝向的平均微震数据点的走向和倾角的示例性绘图。每一个绘制的x-y点处的微震数据点被显示为灰色阴影。识别出主要走向朝向为112°、132°以及178°的三个族群的裂缝平面。另外,与优先倾角在74°的范围不那么垂直的132°和178°族群相比,112°的裂缝平面族群具有大约84°的优先倾角。图2B示出用于每一个朝向的平均微震数据点的显示的另一个格式。图2B中的格式是每一个朝向的平均微震数据点与走向相比的直方图。图2C示出用于每一个朝向的平均微震数据点的显示的又一个格式。图2C中的格式是每一个朝向的平均微震数据点的辐射测绘(radial plot)。
在图3A中,图2A中的格式用来示出在使用蒙特卡洛方法(Monte Carlomethod)的标准统计分析完成之后显示的相同数据的分析。换句话说,显示了涉及一组平面及其朝向不确定性的平均微震数据点的朝向不确定性或主要走向朝向。以此方式,能够使用蒙特卡洛统计分析来验证涉及一组平面及其朝向不确定性的平均微震数据点的主要走向朝向。沿x、y、z方向的每一个微震事件的不确定性遵循高斯分布。随机过程通过x+σxN(0,1)被应用于每一个微震事件的x-位置,其中σx是沿x-方向的不确定性的半径,N(0,1)是正态分布。类似过程被应用于y和z位置。在所有事件被这一过程扰乱之后,重复步骤104-107以收集平均微震数据点并将其绘制成走向和倾角空间。在对多个实现取样之后,地图中的同相分量将突出,异相分量将除去。因此,可以更加清楚地识别主要走向朝向。在图3B中,图2B中的格式用来示出在使用蒙特卡洛方法的标准统计分析完成之后显示的相同数据的分析。在图3C中,图2C中的格式用来示出在使用蒙特卡洛方法的标准统计分析完成之后显示的相同数据的分析。
在步骤109中,使用来自步骤101的微震事件数据在每一个主要走向朝向为分别每一个的倾角重复步骤104-106。以此方式,考虑每一个微震事件的位置不确定性和震级,可以执行处于主要走向朝向的每一个平面的详细分析。
在步骤110中,使用本领域熟知的技术识别来自所创建的每一组平面且具有最大数量的微震数据点交叉点的平面及其位置。该平面表示相对于来自步骤101的微震数据点具有最佳拟合的平面。
在步骤111中,步骤110中识别的平面与交叉微震数据点之间的接触加权基于下列各项来计算:i)沿x、y以及z方向的每一个微震数据点的位置不确定性;ii)从每一个微震数据点的中心到平面的垂直距离;iii)表示每一个单独的微震事件的矩震级的数值加权值;以及iv)微震数据点与平面之间的相关性。每一个微震数据点表示在中心具有最高程度的位置确定性以及在边缘具有最低程度的位置确定性的体积。矩震级较高的微震数据点具有的位置确定性程度比矩震级较低的(其可能几乎探测不到)微震数据点大。因此,低矩震级表示其位置最不确定的微震数据点。因此,在其中心位置测量的微震数据点的高矩震级可以在.5到1之间。在所有情况下,0.05的值被指定为由沿x、y以及z方向的不确定性限定的扁椭球(oblate spheroid)的外边缘处的接触加权。该值从中心位置到边缘线性地减小,并表示位置不确定性。因此,指定给每一个微震数据点的接触加权的数值是下列各项的结果:i)与如由应用于矩震级的标度(scaling)表示的每一个微震事件的识别置信度对应的值;ii)沿x、y以及z方向的每一个微震数据点的位置不确定性;以及iii)从微震数据点的中心位置到交叉平面的垂直距离。使用平面与微震数据点之间的距离的倒数,可以进一步分析微震数据点与平面之间的统计相关性。如果微震数据点存在于一个裂缝平面中,则这些距离的和应该为最小。
在步骤112中,移除接触步骤110中识别的平面的微震数据点,并且重复步骤109-112直到在步骤110中不存在接触三个或更多个微震数据点的平面。以此方式,记录该平面的参数。在单独地识别接触三个或更多个微震数据点的所有平面之后,方法100进行到步骤113。
在步骤113中,为每一组具有共同的主要走向朝向的平面执行统计分析,该组平面表示使用步骤111中计算出的接触加权的一组地下裂缝,以确定与每一组地下裂缝的朝向、间距以及倾角有关的不确定性。该分析可以包括例如为每一组具有共同的主要走向朝向的平面计算出的基本统计量(走向、倾角以及长度的平均值、中值以及标准偏差)。
在步骤114中,显示步骤113的结果。响应于水力压裂期间的微震事件,与由每一个平面表示的所有潜在地下裂缝的朝向、间距以及倾角有关的不确定性可以被显示给随后可以基于可能发生裂缝的地区的地质历史来接受或拒绝分析结果的分析员。当应力型式可能与现在的不同时,则可能存在早期地质时期时产生的裂缝。因此,分析员可以考虑该地区的地质历史知识来识别任何可能的走向、倾角以及间距的组合,该可能的组合鉴于微震事件数据并非水力裂缝的完整说明(尤其是在非常小的微震事件可能低于能够检测这种微震事件的仪器的分辨率时)的事实而可能无效。此外,通过方法100识别的具有最大数学拟合程度的该组平面可能在地质上并不合理,而其它不同的拟合对于整体分析和裂缝估计来说可能更好。
在图4中,示出代表步骤113中一组地下裂缝的一些垂直倾角平面的最佳平面朝向和间距的示例性二维地图。由地图表示的数据集包括284个微震事件。选择三个主要走向朝向(290°、312°、358°)以示出步骤109-115期间的分析。为了简化,没有选择补充的朝向。总共有31个平面被识别出。7个平面属于290°族群中,10个平面属于312°族群中以及14个平面属于358°族群中。这些平面包括221个微震事件,是在这一数据集中观察的总微震事件的78%。为了使这些平面可视化,所有微震事件(灰点)投射在垂直地图视图中的东北方向平面(x,y平面)上。仅示出所识别的平面(灰线)中的12个,这些平面具有大于85°的倾角。在290°族群中具有三个匹配微震事件的平面,在312°族群中有四个匹配微震事件的平面以及在358°族群中有五个匹配微震事件的平面。图5中的表格示出图4中所示三个族群平面的统计。利用每一个的标准偏差(平均值)来显示每一个族群的朝向(走向)、倾角以及间距。
因此,方法100允许更加准确地识别各个族群的裂缝趋势。在模拟的地区中可以存在多个族群的裂缝,每一个族群的裂缝具有独立的朝向、倾角以及间距。在已经识别出复杂裂缝网络的地区中尤其如此。因此,方法100可以用于:
1)输入到复杂的裂缝模型,在该裂缝模型中对裂缝族群的说明以及这些族群在单个井内和从地区内的一个井到另一个井如何变化进行建模;
2)输入到机械地球模型中,该机械地球模型包含应力、岩石参数、分层,以及与地下的结构、层位(stratigraphy)和性质有关的其它信息;
3)输入到井性能和储层模型中,该井性能和储层模型用于生产率的单个井和多个井分析、排水分析以及最终回收率分析;
4)实时分析,该实时分析使用:i)来自后续模拟中第一个模拟的识别的裂缝型态以查看现有的识别裂缝是否继续重新开启以及新的裂缝趋势是否被井中的压裂所拦截;以及ii)来自附近其它井中识别的裂缝族群的数据,以确定是否有相同的裂缝趋势在当前模拟中被诱导;以及
5)分析微震事件数据以确定单独的微震事件的矩震级和震矩张量与所识别的裂缝平面朝向、间距以及倾角族群之间的任何关系。
系统说明
本发明可以经由一般称为由计算机执行的软件应用或应用程序的计算机可执行指令程序(例如,程序模块)来实施。该软件例如可以包括执行特定任务或实施特定抽象数据类型的例程、程序、对象、组件以及数据结构。所述软件形成接口,以允许计算机根据输入源作出回应。该软件也可以与其它代码段配合,以响应于结合所接收的数据源而接收的数据启动各种任务。该软件可以被存储和/或承载在任何类型的存储介质上,例如CD-ROM、磁盘、磁泡存储器(bubble memory)以及半导体存储器(例如,各种类型的RAM或ROM)。此外,该软件及其结果可以通过各种载体介质(例如光学纤维、金属导线)来发送和/或经由任何类型的网络(例如互联网)来发送。
而且,本领域技术人员将意识到本发明可以利用包括手持装置、多处理器系统、基于微处理器或可编程消费电子产品、迷你计算机、大型计算机等各种计算机系统配置来实践。与本发明一起使用的任何数量的计算机系统和计算机网络是可以接受的。本发明可以在分布式计算环境中实践,其中通过经由通信网络链接的远程处理装置来执行任务。在分布式计算环境中,程序模块可以位于包括存储装置的本地计算机存储介质和远程计算机存储介质两者中。因此,本发明可以在计算机系统或其它处理系统中与各种硬件、软件或其组合连接实施。
现在参照图6,方框图示出了用于在计算机上实施本发明的系统的一个实施例。该系统包括有时称为计算系统的计算单元,该计算单元包含存储器、应用程序、客户端接口、视频接口以及处理单元。该计算单元仅是适合的计算环境的一个示例,并且不意欲建议对本发明的使用或功能性的范围的任何限制。
该存储器主要存储应用程序,该应用程序也可以被描述为包含计算机可执行指令的程序模块,其通过用于实施这里描述并如图1-图5所示的本发明的计算单元来执行。因此,该存储器包括裂缝表征模块,其实现参考图1A-图1B所示和描述的方法。虽然复杂的裂缝模型可以被用以利用裂缝表征模块的结果,然而可以代替使用其它接口应用,或者裂缝表征模块可以被用作独立应用。
虽然计算单元被显示为具有通用存储器,然而计算单元典型地包括各种计算机可读介质。举例来说明,并非限制,计算机可读介质可以包括计算机存储介质。计算系统存储器可以包括呈易失性和/或非易失性存储器形式的计算机存储介质,例如只读存储器(ROM)和随机存取存储器(RAM)。基本输入/输出系统(BIOS)被典型地存储在ROM中,所述输入/输出系统包含例如在启动期间有助于计算单元内的元件之间传递信息的基本例程。RAM典型地包含可立即访问和/或目前被处理单元操作的数据和/或程序模块。举例来说明,并非限制,计算单元包括操作系统、应用程序、其它程序模块以及程序数据。
存储器中示出的组件也可以包括在其它可移动/不可移动、易失性/非易失性计算机存储介质中,或者它们可以经由应用程序接口(“API”)在计算单元中实施,所述应用程序接口可以驻留在经由计算机系统或网络连接的单独的计算单元上。仅举例而言,硬盘驱动器可以从不可移动的非易失性磁性介质读取或者向其写入,磁盘驱动器可以从可移动非易失性磁盘读取或者向其写入,以及光盘驱动器可以从可移动非易失性光盘(例如CD ROM或其它光介质)读取或者向其写入。能够在示例性操作环境中使用的其它可移动/不可移动、易失性/非易失性计算机存储介质可以包括但不限于盒式磁带机、快闪存储卡、数字多用光盘、数字视频磁带、固态RAM、固态ROM等。上文讨论的驱动器及其相关的计算机存储介质提供计算机可读指令、数据结构、程序模块以及用于计算单元的其它数据的存储。
客户可以经由客户端接口将命令和信息输入计算单元,该客户端接口可以是诸如键盘和指点设备(通常被称为鼠标、轨迹球或触摸垫)等输入装置。输入装置可以包括麦克风、操纵杆、卫星天线(satellite dish)、扫描仪等。这些及其它输入装置常经由系统总线连接至处理单元,但是可以通过其它接口和总线结构连接,例如并行端口或通用串行总线(USB)。
监视器或其它类型的显示装置可以经由接口(例如,视频接口)被连接至系统总线。图形用户接口(“GUI”)也可以使用视频接口来从客户端接口接收指令,并且将指令发送到处理单元。除了监视器之外,计算机也可以包括其它外围输出装置(例如,扬声器和打印机),该外围输出装置可以经由输出外围接口连接。
虽然计算单元的许多其它内部组件未示出,然而本领域普通技术人员将理解这种组件及其互连是熟知的。
虽然结合目前优选实施例描述了本发明,但本领域技术人员将理解其并不意欲将本发明限制于那些实施例。因此,在不脱离由所附权利要求及其等同项限定的本发明的精神和范围的情况下,预期可对公开的实施例作出各种替代实施例和变型。

Claims (20)

1.一种使用微震数据点的水力裂缝表征的方法,包括:
使用计算机处理器识别与一组或多组平面及其朝向不确定性有关的平均微震数据点的一个或多个走向朝向,在所述一个或多个走向朝向中各个走向朝向处的平均微震数据点数目大于其他各个走向朝向处的平均微震数据点数目;
考虑到每一个微震数据点与每一组平面内每一个平面之间的接触加权来确定每一组平面的倾角;
确定每一组平面内每一个平面的位置;以及
对于具有共同的所述识别出的走向朝向的每一组平面执行统计分析,以确定与各组地下裂缝的朝向、间距以及倾角有关的不确定性。
2.根据权利要求1所述的方法,还包括:识别所有微震数据点的位置。
3.根据权利要求2所述的方法,还包括:为平面设定初始朝向,所述初始朝向包括走向角和倾角。
4.根据权利要求3所述的方法,其中,与一组或多组平面及其朝向不确定性有关的所述平均微震数据点的所述一个或多个走向朝向通过以下步骤来识别:
创建平面,该平面具有法向量并穿过用于所述平面朝向的所有微震数据点的平均位置;
计算每一个微震数据点与所述平面之间的垂直距离;
创建与所述平面平行的一组平面,其中该组平面中的每一个平面与由每一个微震数据点的位置不确定性表示的体积相交;
对该组平面中的平面的数量以及与该组平面中的所述平面接触的微震数据点的数量计数;
对于该组平面中所述平面的每一个朝向计算与该组平面中每一个平面接触的所述平均微震数据点;
将所述平面的所述走向角或所述倾角增量1度;以及
重复之前的每一个步骤直到所述走向角等于0度或360度且所述倾角等于0度或90度。
5.根据权利要求4所述的方法,其中,每一组平面的所述倾角以及每一组平面内每一个平面的所述位置通过以下步骤来确定:
分别对于每一个倾角在每一个所述识别出的走向朝向上重复权利要求4中的每一个步骤;
从创建的每一组平面中识别出具有最大数量的微震数据点交叉点的平面及其位置;
计算具有所述最大数量的微震数据点交叉点的所述平面与所述微震数据点之间的接触加权;
移除与具有所述最大数量的微震数据点交叉点的所述平面接触的所述微震数据点;以及
重复之前的每一个步骤直到不存在接触三个或更多个微震数据点的具有所述最大数量的微震数据点交叉点的平面。
6.根据权利要求1所述的方法,其中,所述一组或多组平面的所述朝向不确定性通过执行蒙特卡洛统计分析来确定。
7.根据权利要求5所述的方法,其中,考虑到每一组平面和与每一组平面有关的所述微震数据点的所述接触加权,来计算与分别每一组地下裂缝的所述朝向、间距以及倾角有关的所述不确定性。
8.根据权利要求5所述的方法,其中,所述接触加权基于每一个微震数据点的所述位置不确定性、每一个微震数据点与所述平面之间的所述垂直距离、以及表示与每一个微震数据点相关的矩震级的数值加权值。
9.根据权利要求5所述的方法,其中,所述接触加权被分别指定给每一个微震数据点。
10.根据权利要求1所述的方法,还包括:显示所述统计分析。
11.一种用于使用微震数据点的水力裂缝表征的装置,包括:
用于识别与一组或多组平面及其朝向不确定性有关的平均微震数据点的一个或多个走向朝向的模块,在所述一个或多个走向朝向中各个走向朝向处的平均微震数据点数目大于其他各个走向朝向处的平均微震数据点数目;
用于考虑到每一个微震数据点与每一组平面内每一个平面之间的接触加权来确定每一组平面的倾角的模块;
用于确定每一组平面内每一个平面的位置的模块;以及
用于对于具有共同的所述识别出的走向朝向的每一组平面执行统计分析,以确定与各组地下裂缝的朝向、间距以及倾角有关的不确定性的模块。
12.根据权利要求11所述的装置,还包括:用于识别所有微震数据点的位置的模块。
13.根据权利要求12所述的装置,还包括:用于为平面设定初始朝向的模块,所述初始朝向包括走向角和倾角。
14.根据权利要求13所述的装置,其中,与一组或多组平面及其朝向不确定性有关的所述平均微震数据点的所述一个或多个走向朝向通过以下步骤来识别:
创建平面,该平面具有法向量并穿过用于所述平面朝向的所有微震数据点的平均位置;
计算每一个微震数据点与所述平面之间的垂直距离;
创建与所述平面平行的一组平面,其中该组平面中的每一个平面与由每一个微震数据点的位置不确定性表示的体积相交;
对该组平面中的平面的数量以及与该组平面中的所述平面接触的微震数据点的数量计数;
对于该组平面中所述平面的每一个朝向计算与该组平面中的每一个平面接触的所述平均微震数据点;
将所述平面的所述走向角或所述倾角增量1度;以及
重复之前的每一个步骤直到所述走向角等于0度或360度且所述倾角等于0度或90度。
15.根据权利要求14所述的装置,其中,每一组平面的所述倾角以及每一组平面内每一个平面的所述位置通过以下步骤来确定:
对于每一个倾角分别在每一个所述识别出的走向朝向上重复权利要求14中的每一个步骤;
从创建的每一组平面中识别出具有最大数量的微震数据点交叉点的平面及其位置;
计算具有所述最大数量的微震数据点交叉点的所述平面与所述微震数据点之间的接触加权;
移除与具有所述最大数量的微震数据点交叉点的所述平面接触的所述微震数据点;以及
重复之前的每一个步骤直到不存在接触三个或更多个微震数据点的具有所述最大数量的微震数据点交叉点的平面。
16.根据权利要求11所述的装置,其中,所述一组或多组平面的所述朝向不确定性通过执行蒙特卡洛统计分析来确定。
17.根据权利要求15所述的装置,其中,考虑到每一组平面和与每一组平面有关的所述微震数据点的所述接触加权,来计算与分别每一组地下裂缝的所述朝向、间距以及倾角有关的所述不确定性。
18.根据权利要求15所述的装置,其中,所述接触加权基于每一个微震数据点的所述位置不确定性、每一个微震数据点与所述平面之间的所述垂直距离、以及表示与每一个微震数据点相关的矩震级的数值加权值。
19.根据权利要求15所述的装置,其中,所述接触加权被分别指定给每一个微震数据点。
20.根据权利要求11所述的装置,还包括:用于显示所述统计分析的模块。
CN201180070143.4A 2011-04-15 2011-04-15 使用微震事件数据的水力裂缝表征的系统和方法 Expired - Fee Related CN103548032B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US2011/032741 WO2012141720A1 (en) 2011-04-15 2011-04-15 Systems and methods for hydraulic fracture characterization using microseismic event data

Publications (2)

Publication Number Publication Date
CN103548032A CN103548032A (zh) 2014-01-29
CN103548032B true CN103548032B (zh) 2016-10-05

Family

ID=47009618

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180070143.4A Expired - Fee Related CN103548032B (zh) 2011-04-15 2011-04-15 使用微震事件数据的水力裂缝表征的系统和方法

Country Status (9)

Country Link
US (1) US9009010B2 (zh)
EP (1) EP2697734B1 (zh)
CN (1) CN103548032B (zh)
AU (1) AU2011365515B2 (zh)
BR (1) BR112013026513A2 (zh)
CA (1) CA2831251A1 (zh)
EA (1) EA025073B1 (zh)
MX (1) MX2013012009A (zh)
WO (1) WO2012141720A1 (zh)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2743611C (en) * 2011-06-15 2017-03-14 Engineering Seismology Group Canada Inc. Methods and systems for monitoring and modeling hydraulic fracturing of a reservoir field
US9417348B2 (en) 2012-10-05 2016-08-16 Halliburton Energy Services, Inc. Updating microseismic histogram data
US9377547B2 (en) * 2012-10-05 2016-06-28 Halliburton Energy Services, Inc. Analyzing fracture stratigraphy
WO2014110542A1 (en) * 2013-01-14 2014-07-17 Westerngeco Seismic Holdings Limited Method of analyzing seismic data
US8960280B2 (en) 2013-01-23 2015-02-24 Microseismic, Inc. Method for determining fracture plane orientation using passive seismic signals
US9612359B2 (en) * 2013-06-12 2017-04-04 Baker Hughes Incorporated Generation of fracture networks using seismic data
US9551208B2 (en) * 2013-08-26 2017-01-24 Halliburton Energy Services, Inc. Identifying uncertainty associated with a stimulated reservoir volume (SRV) calculation
US9903189B2 (en) 2013-08-26 2018-02-27 Halliburton Energy Services, Inc. Real-time stimulated reservoir volume calculation
US9523275B2 (en) 2013-08-26 2016-12-20 Halliburton Energy Services, Inc. Identifying an axis of a stimulated reservoir volume for a stimulation treatment of a subterranean region
US9529103B2 (en) 2013-08-26 2016-12-27 Halliburton Energy Services, Inc. Identifying overlapping stimulated reservoir volumes for a multi-stage injection treatment
US9529104B2 (en) 2013-08-26 2016-12-27 Halliburton Energy Services, Inc. Indentifying a stimulated reservoir volume from microseismic data
US10006271B2 (en) 2013-09-26 2018-06-26 Harris Corporation Method for hydrocarbon recovery with a fractal pattern and related apparatus
FR3012623B1 (fr) * 2013-10-31 2017-02-17 Cgg Services Sa Systeme et procede pour analyser des evenements microsismiques en utilisant des groupes
CA2943189C (en) * 2014-04-30 2019-11-12 Halliburton Energy Services, Inc. Bin constraints for generating a histogram of microseismic data
WO2015178885A1 (en) * 2014-05-19 2015-11-26 Halliburton Energy Services, Inc. Identifying an error bound of a stimulated reservoir volume of a subterranean region
CA2945467C (en) 2014-05-23 2019-12-17 Halliburton Energy Services, Inc. Enhancing reservoir characterization using real-time srv and fracture evolution parameters
WO2016039773A1 (en) 2014-09-12 2016-03-17 Halliburton Energy Services, Inc. Analysis of microseismic supported stimulated reservoir volumes
US10527744B2 (en) 2014-10-13 2020-01-07 Halliburton Energy Services, Inc. Data-driven estimation of stimulated reservoir volume
US10197704B2 (en) 2014-12-19 2019-02-05 Baker Hughes, A Ge Company, Llc Corrective scaling of interpreted fractures based on the microseismic detection range bias correction
CN105137488B (zh) * 2015-09-30 2018-04-17 河北煤炭科学研究院 基于人工注液的地下水流场描述方法
EP3374596B1 (en) * 2015-11-10 2023-06-07 Landmark Graphics Corporation Fracture network triangle mesh adjustment
US10267132B2 (en) * 2015-12-21 2019-04-23 Baker Hughes, A Ge Company, Llc Eliminating discrete fracture network calculations by rigorous mathematics
CN105629309B (zh) * 2015-12-30 2017-12-12 河北煤炭科学研究院 钻孔注浆浆液空间扩散范围及路径描述方法
CN111175816B (zh) * 2020-01-06 2022-04-15 中国石油化工股份有限公司 油藏改造实时构建微地震裂缝网络的方法及装置
CN112502700B (zh) * 2020-12-03 2022-03-01 西南石油大学 一种深层复杂构造页岩非均匀应力场分布的计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101393267A (zh) * 2004-03-16 2009-03-25 尖顶技术公司 用于分析倾斜数据和微震数据的方法
CN101421640A (zh) * 2006-02-09 2009-04-29 普拉德研究及开发有限公司 用于预测井位的碳氢化合物产量的方法和设备
CN201251564Y (zh) * 2008-05-12 2009-06-03 西安西科测控设备有限责任公司 一种实时监测矿井顶板岩层或混凝土结构稳定性的装置

Family Cites Families (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4698759A (en) 1985-03-28 1987-10-06 The United States Of America As Represented By The United States Department Of Energy Process for structural geologic analysis of topography and point data
FR2725814B1 (fr) 1994-10-18 1997-01-24 Inst Francais Du Petrole Methode pour cartographier par interpolation, un reseau de lignes, notamment la configuration de failles geologiques
FR2725794B1 (fr) 1994-10-18 1997-01-24 Inst Francais Du Petrole Methode pour modeliser la distribution spatiale d'objets geometriques dans un milieu, tels que des failles dans une formation geologique
US5996726A (en) * 1998-01-29 1999-12-07 Gas Research Institute System and method for determining the distribution and orientation of natural fractures
JP4128489B2 (ja) * 2003-06-13 2008-07-30 富士通株式会社 ショック検出装置、ディスク装置、ショック検出方法、ショック検出プログラム
US6985816B2 (en) * 2003-09-15 2006-01-10 Pinnacle Technologies, Inc. Methods and systems for determining the orientation of natural fractures
KR100568285B1 (ko) * 2003-12-13 2006-04-05 삼성전기주식회사 틸트 보상형 전자 나침반의 복각 탐색 방법
US7391675B2 (en) 2004-09-17 2008-06-24 Schlumberger Technology Corporation Microseismic event detection and location by continuous map migration
US7460436B2 (en) 2005-12-05 2008-12-02 The Board Of Trustees Of The Leland Stanford Junior University Apparatus and method for hydraulic fracture imaging by joint inversion of deformation and seismicity
US7254091B1 (en) 2006-06-08 2007-08-07 Bhp Billiton Innovation Pty Ltd. Method for estimating and/or reducing uncertainty in reservoir models of potential petroleum reservoirs
GB2439571B (en) 2006-06-28 2008-11-12 Schlumberger Holdings Method for updating a model of the earth using microseismic measurements
US7746725B2 (en) 2006-12-04 2010-06-29 Schlumberger Technology Corporation Fracture clusters identification
US7777606B2 (en) 2007-01-09 2010-08-17 Westerngeco L.L.C. Fracture cluster mapping
US8494777B2 (en) 2008-04-09 2013-07-23 Schlumberger Technology Corporation Continuous microseismic mapping for real-time 3D event detection and location
GB2466438B (en) 2008-12-17 2011-04-06 Schlumberger Holdings Analysis of fracture networks
US8908473B2 (en) 2008-12-23 2014-12-09 Schlumberger Technology Corporation Method of subsurface imaging using microseismic data
US20100238765A1 (en) 2009-03-20 2010-09-23 Grechka Vladimir Reservoir characterization from multicomponent microseismic data
US20100256964A1 (en) 2009-04-07 2010-10-07 Schlumberger Technology Corporation System and technique to quantify a fracture system
WO2010116236A2 (en) 2009-04-08 2010-10-14 Schlumberger Technology B.V. Methods and systems for microseismic mapping
US8301427B2 (en) 2009-06-05 2012-10-30 Schlumberger Technology Corporation Fracture network characterization method
US8498852B2 (en) 2009-06-05 2013-07-30 Schlumberger Tehcnology Corporation Method and apparatus for efficient real-time characterization of hydraulic fractures and fracturing optimization based thereon
US20110029291A1 (en) * 2009-07-31 2011-02-03 Xiaowei Weng Method for fracture surface extraction from microseismic events cloud
US20110029293A1 (en) * 2009-08-03 2011-02-03 Susan Petty Method For Modeling Fracture Network, And Fracture Network Growth During Stimulation In Subsurface Formations
US7978563B2 (en) 2009-08-18 2011-07-12 Microseismic, Inc. Method for passive seismic emission tomography including polarization correction for source mechanism
US9410421B2 (en) * 2009-12-21 2016-08-09 Schlumberger Technology Corporation System and method for microseismic analysis

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101393267A (zh) * 2004-03-16 2009-03-25 尖顶技术公司 用于分析倾斜数据和微震数据的方法
CN101421640A (zh) * 2006-02-09 2009-04-29 普拉德研究及开发有限公司 用于预测井位的碳氢化合物产量的方法和设备
CN201251564Y (zh) * 2008-05-12 2009-06-03 西安西科测控设备有限责任公司 一种实时监测矿井顶板岩层或混凝土结构稳定性的装置

Also Published As

Publication number Publication date
US9009010B2 (en) 2015-04-14
EA025073B1 (ru) 2016-11-30
EP2697734A1 (en) 2014-02-19
EP2697734B1 (en) 2017-07-05
US20130304437A1 (en) 2013-11-14
EA201391505A1 (ru) 2014-03-31
MX2013012009A (es) 2014-08-21
AU2011365515A1 (en) 2013-10-24
EP2697734A4 (en) 2014-12-31
WO2012141720A1 (en) 2012-10-18
CN103548032A (zh) 2014-01-29
AU2011365515B2 (en) 2014-07-03
CA2831251A1 (en) 2012-10-18
BR112013026513A2 (pt) 2016-12-27

Similar Documents

Publication Publication Date Title
CN103548032B (zh) 使用微震事件数据的水力裂缝表征的系统和方法
US9268050B2 (en) Determining a confidence value for a fracture plane
US20190242233A1 (en) Stimulation Using Fiber-Derived Information and Fracturing Modeling
US11156740B2 (en) Electrofacies determination
CA2943189C (en) Bin constraints for generating a histogram of microseismic data
Rezaei et al. Applications of Machine Learning for Estimating the Stimulated Reservoir Volume (SRV)
Adenan Machine Learning Assisted Framework for Advanced Subsurface Fracture Mapping and Well Interference Quantification
WO2024064009A1 (en) Machine learning training for characterizing water injection and seismic prediction

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20161005

Termination date: 20170415