CN111624661B - 基于微震散射波品质因子的压裂效果评价方法及装置 - Google Patents

基于微震散射波品质因子的压裂效果评价方法及装置 Download PDF

Info

Publication number
CN111624661B
CN111624661B CN202010361714.0A CN202010361714A CN111624661B CN 111624661 B CN111624661 B CN 111624661B CN 202010361714 A CN202010361714 A CN 202010361714A CN 111624661 B CN111624661 B CN 111624661B
Authority
CN
China
Prior art keywords
microseismic
wave signal
fracturing
scattered wave
scattered
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
CN202010361714.0A
Other languages
English (en)
Other versions
CN111624661A (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.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN202010361714.0A priority Critical patent/CN111624661B/zh
Publication of CN111624661A publication Critical patent/CN111624661A/zh
Application granted granted Critical
Publication of CN111624661B publication Critical patent/CN111624661B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • 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
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明实施例提供一种基于微震散射波品质因子的压裂效果评价方法及装置,该方法包括:在压裂过程中,通过微震监测系统采集压裂区原始微震信号;根据地震波走时区分直达波和散射波;获取微震信号中的直达波和散射波的质心频率,得到压裂区品质因子变化;利用散射波逆时偏移技术对散射体进行定位,确定裂缝发育区域;根据品质因子的变化和裂缝发育区域对压裂区压裂效果进行评价。本发明能够准确地对有效压裂区域进行确定,提高压裂效果评价的准确性。

Description

基于微震散射波品质因子的压裂效果评价方法及装置
技术领域
本发明实施例涉及油气物理探测技术领域,尤其涉及一种基于微震散射波品质因子的压裂效果评价方法及装置。
背景技术
随着非常规油气的开发,尤其是页岩油气的快速发展,水力压裂技术逐渐成为提高非常规油气采收率的重要技术手段之一。在实际生产中,不能仅仅使用水力压裂手就直接对储层进行开采,还要对水力压裂的效果进行评价。
目前,现有技术中主要利用压裂过程中的微震信号中的直达波进行压裂效果的评价。即在压裂过程中通过采集微震信号,对采集的信号进行处理和解释,利用微震信号中的直达波信号获取裂缝参数信息,再利用裂缝参数信评价水力压裂的压裂效果。
但是,在实际地层中,有些裂缝压裂不足,这些裂缝对油井产能无贡献,现有技术的利用裂缝参数评价水力压裂的压裂效果时,不能准确地分辨对油井产能有贡献的裂缝与对油井产能无贡献的裂缝,易造成微震监测得到改造体积大于实际有效的改造体积,影响压裂效果评价的准确性。
发明内容
本发明实施例提供一种基于微震散射波品质因子的压裂效果评价方法及装置,能够准确地对压裂区有效压裂区域进行确定,提高压裂效果评价的准确性。
第一方面,本发明实施例提供一种基于微震散射波品质因子的压裂效果评价方法,包括:
在压裂过程中,通过微震监测系统采集压裂区的原始微震信号;
根据直达波走时区分所述原始微震信号中的直达波信号和散射波信号;
获取微震信号中直达波信号的质心频率和散射波信号的质心频率;
根据所述直达波信号的质心频率和所述散射波信号的质心频率,得到压裂区的品质因子变化;
利用散射波逆时偏移技术对散射体进行定位,确定裂缝发育区域;
根据所述品质因子的变化和所述裂缝发育区域对所述压裂区的压裂效果进行评价。
第二方面,本发明实施例提供一种基于微震散射波品质因子的压裂效果评价装置,包括:
采集单元,用于在压裂过程中,通过微震监测系统采集压裂区的原始微震信号;
处理单元,用于根据地震波走时区分所述原始微震信号中的直达波和散射波信号;获取微震信号中直达波和散射波信号的质心频率;根据所述直达波信号的质心频率和所述散射波信号的质心频率,得到压裂区的品质因子变化;利用散射波逆时偏移技术对散射体进行定位,确定裂缝发育区域;
评价单元,用于根据所述品质因子的变化和所述裂缝发育区域对所述压裂区的压裂效果进行评价。
第三方面,本发明实施例提供一种基于微震散射波品质因子的压裂效果评价设备,包括:至少一个处理器和存储器;
所述存储器存储计算机执行指令;
所述至少一个处理器执行所述存储器存储的计算机执行指令,使得所述至少一个处理器执行如上第一方面以及第一方面各种可能的设计所述的基于微震散射波品质因子的压裂效果评价方法。
第四方面,本发明实施例提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现如上第一方面以及第一方面各种可能的设计所述的基于微震散射波品质因子的压裂效果评价方法。
本实施例提供的基于微震散射波品质因子的压裂效果评价方法及装置,根据直达波走时区分所述原始微震信号中的直达波和散射波信号;获取微震信号中直达波和散射波信号的质心频率;根据所述直达波信号的质心频率和所述散射波信号的质心频率,得到压裂区的品质因子变化;对散射波进行逆时成像,定位裂缝发育充足的压裂区,进而对压裂区的压裂效果进行评价,能够准确评价压裂区的压裂效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一实施例提供的水力压裂过程中微震监测的场景示意图;
图2为本发明实施例提供的基于微震散射波品质因子的压裂效果评价方法的流程示意图;
图3为本发明实施例提供的利用地震波走时区分地震记录中直达波信号和散射波信号示意图;
图4为本发明实施例提供的地震信号衰减前后质心位置变化示意图;
图5为本发明另一实施例提供的基于微震散射波品质因子的压裂效果评价方法的流程示意图;
图6为本发明实施例提供的基于微震散射波品质因子的压裂效果评价装置的结构框图;
图7为本发明实施例提供的基于微震散射波品质因子的压裂效果评价设备的硬件结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面首先对本发明涉及的名词进行解释。
微震:是指由岩石破裂或流体扰动产生微小的震动。广义上的微震可以分为两大类:工程生产上的微震(microseism)和自然产生的微地震(microearthquake)。前者是由人为生产施工导致岩石破裂产生的。如油田压裂微震监测技术,矿山安全生产微震监测预警系统、水电站大坝微震监测系统等等。后者是由天然应力场的变化引起岩石破裂或者岩浆、雨水等流体扰动引起孔压变化产生的。如利用大地震后微震研究大地震产生机理,利用微震数据研究断裂走向和长度,利用微震数据反演火山区速度和品质因子Q等等。
水力压裂:水力压裂在非常规油气开采应用中是一种重要的工程方法,该方法能够改造中低渗透油气层以使油气井增产,可以使非常规油气从低渗透性的致密储层中得到高效回收。水力压裂就是利用井泵把大排量压裂液输送到地层,当注入压裂液到一定程度储层将受到高压作用产生裂缝,再持续向储层注入压裂液,导致裂缝持续扩张。为使裂缝不被挤压闭合,接着注入带有支撑剂的压裂液,一方面可以使裂缝持续存在不会由于压力闭合,另一方面还可以使裂缝继续延伸。最后回收压裂液,支撑剂保留在裂缝中留下多条裂缝,使储层中有了新的流体通道,压裂之后的油气井的产量会大幅增长。因此,在油气勘探和开发中评价压裂效果具有重要意义。
品质因子:波在岩石中传播会发生衰减,这种衰减性质通常用岩石的品质因子Q来表征。本申请中引起品质因子变化的原因是裂缝中异常介质对地震波的吸收能力增强。
散射波:散射波从物理学角度来说,如果界面凹凸不平,在凹(或凸)部分的尺度相对于波长很大时,发生波的反射;相对于波长较小或可比时则发生散射,形成散射波。本申请中所述散射波是裂缝中介质异常产生的散射波。
目前,现有技术中在实际地层中,有些裂缝压裂不足,这些裂缝对油井产能无贡献,现有技术的利用裂缝参数信评价水力压裂的压裂效果时,不能准确地分辨对油井产能有贡献的裂缝与对油井产能无贡献的裂缝,易造成微震监测得到改造体积大于实际有效的改造体积,影响压裂效果评价的准确性。
为了准确地对有效压裂区域的进行确定,本发明实施例提出了如下解决思路:通过切除原始微震信号中的直达波,得到散射波信号,并根据直达波信号的质心频率以及所述散射波信号的质心频率,得到压裂区的品质因子,再结合地层背景品质因子,得到压裂区的品质因子变化,定位裂缝发育充足的压裂区,进而对压裂区的压裂效果进行评价,能够准确评价压裂区的压裂效果。
参考图1,图1为本发明一实施例提供的水力压裂过程中微震监测的场景示意图。压裂井为对储层进行水力压裂的井,微震(或微地震)为压裂过程形成裂缝过程中产生的震动。监测井中布置有微震监测系统,微震监测系统包括传感器(检波器)和台站,其中传感器用于采集微震信号,台站用于收集监测井的微震信号。监测井可以是多个,对应多组传感器和台站,多组台站采集的微震信号汇总至处理设备。其中处理设备可以是计算机或服务器。
参考图2,图2为本发明实施例提供的基于微震散射波品质因子的压裂效果评价方法的流程示意图。本方法可以应用于图1中的处理设备,包括:
S201:在压裂过程中,通过微震监测系统采集压裂区的原始微震信号。
在本实施例中,通过水力压裂技术对压裂区的储层进行改造。在压裂过程中,产生裂缝的过程中,有微震发生。微震通过地层传到监测井,在监测井中的微震系统的传感器采集原始微震信号。
S202:根据地震波走时区分所述原始微震信号中的直达波信号和散射波信号。
本实施例中,利用地震记录中检波器记录到直达波和散射波的走时差区分直达波和散射波。如图3所示。
S203:获取微震信号中直达波信号的质心频率以及所述散射波信号的质心频率。
具体地,将直达波信号的振幅谱表示为:
Figure BDA0002475301060000051
式中,Z(f)为直达波信号的振幅谱,A为控制振幅大小的常数,k为控制振幅对称性的常数,f为直达波信号频率,f0为控制带宽的带宽因子;
基于直达波信号的能量谱Z(f)2得到直达波信号的质心频率表示为:
Figure BDA0002475301060000052
将散射波信号振幅谱表示为:
Figure BDA0002475301060000061
式中,S(f)是散射波信号的振幅谱,Δt是直达波和散射波信号的旅行时差,Q为品质因子;
基于散射波信号的能量谱S(f)2得到散射波信号的质心频率表示为:
Figure BDA0002475301060000062
S204:根据所述直达波信号的质心频率和所述散射波信号的质心频率,得到压裂区的品质因子变化。
在本实施例中,将所述直达波信号的质心频率减去所述散射波信号的质心频率,得到质心频移量,公式如下:
Figure BDA0002475301060000063
根据所述质心频移量Δf,得到压裂区的品质因子Q:
Figure BDA0002475301060000064
式中,Δf为质心频移量。
在本实施例中,步骤S202和203中,地震信号在地层介质中传播时,高频成分往往比低频成分更易被地层吸收,会呈现接收信号趋于低频的现象。观察地震信号波形可以看出波形变宽、主频明显下降,因此根据这一现象提出了质心频率偏移法,该方法是通过计算地震波在传播时子波频谱的质心频率发生偏移的偏移量与地层品质因子Q之间的变化关系来计算地层品质因子Q值的。其中,图4为本发明实施例提供的地震信号衰减前后质心位置变化示意图。
S205:利用散射波逆时偏移技术对散射体进行定位,确定裂缝发育区域。
具体地,散射波逆时偏移技术包括:对每个微震事件,将每个监测点记录的散射波信号在时间轴上反转,并作为新的震源向外散发地震波;将所有监测点向外散发的地震波场进行互相关操作,并根据预设的成像条件,定位散射体位置,即裂缝发育区域。
S206:根据所述品质因子的变化和所述裂缝发育区域对所述压裂区的压裂效果进行评价。
具体地,根据所述品质因子的变化确定发生充分压裂的裂缝区域;将所述充分压裂的裂缝区域中,与储层位置距离满足预设条件的裂缝充分发育区域确定为有效压裂区域。
在本实施例中,每条裂缝的位置与储层的位置关系为每条裂缝到储层的距离。距离储层较远的裂缝对油井的产能无贡献,距离储层较近的裂缝对油井的产能有贡献。
从上述实施例可知,本实施例根据直达波信号的质心频率以及所述散射波信号的质心频率,得到压裂区的品质因子变化,定位裂缝发育的压裂区,进而对压裂区的压裂效果进行评价,能够准确评价压裂区的压裂效果。
参考图5,图5为本发明另一实施例提供的基于微震散射波品质因子的压裂效果评价方法的流程示意图。对上述实施例的步骤S205利用散射波逆时偏移技术对散射体进行定位,确定裂缝发育区域,的过程具体详述如下:
S2051:对每个微震事件,将每个监测点记录的散射波信号在时间轴上反转,并作为新的震源向外散发地震波。
在本实施例中,公式如下:
R(x,z,t)=D(x,τ0-t)*G(x,z,t) (7)
式中R(x,z,t)为反传波场,D(x,τ0-t)为地震记录,τ0为所有监测点记录到该微震信号的最迟时间,G(x,z,t)是格林函数,*表示褶积。
S2052:将所有监测点向外散发的地震波场进行互相关操作,并根据预设的成像条件,获取散射体位置,即裂缝发育区域。
在本实施例中,预设成像条件公式为:
Figure BDA0002475301060000071
式中,t表示时间;R表示空间坐标(x,z)处的反向传播检波器波场;I(x,z)表示(x,z)处的成像计算结果;n为检波器的个数。
在本实施例中,根据所述品质因子的变化确定发生充分压裂的裂缝区域;将所述充分压裂的裂缝区域中,与储层位置较近的裂缝发育区域确定为有效压裂区域。
从上述实施例可知,本实施例根据直达波信号的质心频率以及所述散射波信号的质心频率,得到压裂区的品质因子变化,对散射波逆时定位裂缝发育的压裂区:再根据裂缝充分发育区域与储层的位置关系,最终确定对油井产能有贡献的裂缝区域,准确定位有效压裂区。
对应于上文实施例的基于微震散射波品质因子的压裂效果评价方法,图6为本发明实施例提供的基于微震散射波品质因子的压裂效果评价装置的结构框图。为了便于说明,仅示出了与本发明实施例相关的部分。
参照图6,所述基于微震散射波品质因子的压裂效果评价装置60包括:采集单元601、处理单元602和评价单元603。
采集单元601,用于在压裂过程中,通过微震监测系统采集压裂区的原始微震信号;
处理单元602,用于根据地震波走时区分所述原始微震信号中的直达波和散射波信号;获取微震信号中直达波和散射波信号的质心频率;根据所述直达波信号的质心频率和所述散射波信号的质心频率,得到压裂区的品质因子变化;利用散射波逆时偏移技术对散射体进行定位,确定裂缝发育区域;
评价单元603,用于根据所述品质因子的变化和所述裂缝发育区域对所述压裂区的压裂效果进行评价。
本实施例提供的装置,可用于执行上述方法实施例的技术方案,其实现原理和技术效果类似,本实施例此处不再赘述。
在一种可能的设计中,所述处理单元602,还用于对每个微震事件,将每个监测点记录的散射波信号在时间轴上反转,并作为新的震源向外散发地震波,公式如下:
R(x,z,t)=D(x,τ0-t)*G(x,z,t)
式中R(x,z,t)为反传波场,D(x,τ0-t)为地震记录,τ0为所有监测点记录到该微震信号的最迟时间,G(x,z,t)是格林函数,*表示褶积。
将所有监测点向外散发的地震波场进行互相关操作,并根据预设的成像条件,定位散射体位置,即裂缝发育区域,公式如下:
Figure BDA0002475301060000081
式中,t表示时间;R表示空间坐标(x,z)处的反向传播检波器波场;I(x,z)表示(x,z)处的成像计算结果;n为检波器的个数。
在一种可能的设计中,所述处理单元602,用于将直达波信号的振幅谱表示为:
Figure BDA0002475301060000091
式中,Z(f)为直达波信号的振幅谱,A为控制振幅大小的常数,k为控制振幅对称性的常数,f为直达波信号频率,f0为控制带宽的带宽因子;
基于直达波信号的能量谱Z(f)2得到直达波信号的质心频率表示为:
Figure BDA0002475301060000092
将散射波信号振幅谱表示为:
Figure BDA0002475301060000093
式中,S(f)是散射波信号的振幅谱,Δt是直达波和散射波信号的旅行时差,Q为品质因子;
基于散射波信号的能量谱S(f)2得到散射波信号的质心频率表示为:
Figure BDA0002475301060000094
参考图7,图7为本发明实施例提供的基于微震散射波品质因子的压裂效果评价设备的硬件结构示意图。如图7所示,本实施例的基于微震散射波品质因子的压裂效果评价设备70包括:处理器701以及存储器702;其中,存储器702,用于存储计算机执行指令;
处理器701,用于执行存储器存储的计算机执行指令,以实现上述实施例中终端设备所执行的各个步骤。具体可以参见前述方法实施例中的相关描述。
可选地,存储器702既可以是独立的,也可以跟处理器701集成在一起。
当存储器702独立设置时,该基于微震散射波品质因子的压裂效果评价设备还包括总线703,用于连接所述存储器702和处理器701。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现如上所述的基于微震散射波品质因子的压裂效果评价方法。
在本发明所提供的几个实施例中,应该理解到,所揭露的设备和方法,可以通过其它的方式实现。例如,以上所描述的设备实施例仅仅是示意性的,例如,所述模块的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个模块可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或模块的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能模块可以集成在一个处理单元中,也可以是各个模块单独物理存在,也可以两个或两个以上模块集成在一个单元中。上述模块成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。
上述以软件功能模块的形式实现的集成的模块,可以存储在一个计算机可读取存储介质中。上述软件功能模块存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)或处理器(英文:processor)执行本申请各个实施例所述方法的部分步骤。
应理解,上述处理器可以是中央处理单元(Central Processing Unit,简称CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,简称DSP)、专用集成电路(Application Specific Integrated Circuit,简称ASIC)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合公开所公开的方法的步骤可以直接体现为硬件处理器执行完成,或者用处理器中的硬件及软件模块组合执行完成。
存储器可能包含高速RAM存储器,也可能还包括非易失性存储NVM,例如至少一个磁盘存储器,还可以为U盘、移动硬盘、只读存储器、磁盘或光盘等。
总线可以是工业标准体系结构(Industry Standard Architecture,简称ISA)总线、外部设备互连(Peripheral Component,简称PCI)总线或扩展工业标准体系结构(Extended Industry Standard Architecture,简称EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示,本申请附图中的总线并不限定仅有一根总线或一种类型的总线。
上述存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现,如静态随机存取存储器(SRAM),电可擦除可编程只读存储器(EEPROM),可擦除可编程只读存储器(EPROM),可编程只读存储器(PROM),只读存储器(ROM),磁存储器,快闪存储器,磁盘或光盘。存储介质可以是通用或专用计算机能够存取的任何可用介质。
一种示例性的存储介质耦合至处理器,从而使处理器能够从该存储介质读取信息,且可向该存储介质写入信息。当然,存储介质也可以是处理器的组成部分。处理器和存储介质可以位于专用集成电路(Application Specific Integrated Circuits,简称ASIC)中。当然,处理器和存储介质也可以作为分立组件存在于电子设备或主控设备中。
本领域普通技术人员可以理解:实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时,执行包括上述各方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (4)

1.一种基于微震散射波品质因子的压裂效果评价方法,其特征在于,包括:
在压裂过程中,通过微震监测系统采集压裂区的原始微震信号;
根据地震波走时区分所述原始微震信号中的直达波信号和散射波信号;
将直达波信号的振幅谱表示为:
Figure FDA0002894313580000011
式中,Z(f)为直达波信号的振幅谱,A为控制振幅大小的常数,k为控制振幅对称性的常数,f为直达波信号频率,f0为控制带宽的带宽因子;
基于直达波信号的能量谱Z(f)2得到直达波信号的质心频率表示为:
Figure FDA0002894313580000012
将散射波信号振幅谱表示为:
Figure FDA0002894313580000013
式中,S(f)是散射波信号的振幅谱,Δt是直达波和散射波信号的旅行时差,Q为品质因子;
基于散射波信号的能量谱S(f)2得到散射波信号的质心频率表示为:
Figure FDA0002894313580000014
将所述直达波信号的质心频率减去所述散射波信号的质心频率,得到质心频移量,公式如下:
Figure FDA0002894313580000015
根据所述质心频移量Δf,得到压裂区的品质因子Q:
Figure FDA0002894313580000016
式中,Δf为质心频移量;
对每个微震事件,将每个监测点记录的散射波信号在时间轴上反转,并作为新的震源向外散发地震波,公式如下:
R(x,z,t)=D(x,τ0-t)*G(x,z,t)
式中R(x,z,t)为反传波场,D(x,τ0-t)为地震记录,τ0为所有监测点记录到该微震信号的最迟时间,G(x,z,t)是格林函数,*表示褶积;
将所有监测点向外散发的地震波场进行互相关操作,并根据预设的成像条件,定位散射体位置,即裂缝发育区域,公式如下:
Figure FDA0002894313580000021
式中,t表示时间;R表示空间坐标(x,z)处的反向传播检波器波场;I(x,z)表示(x,z)处的成像计算结果;n为检波器的个数;
根据所述品质因子的变化确定发生充分压裂的裂缝区域;
将所述充分压裂的裂缝区域中,与储层位置距离满足预设条件的裂缝充分发育区域确定为有效压裂区域。
2.一种基于微震散射波品质因子的压裂效果评价装置,其特征在于,
采集单元,用于在压裂过程中,通过微震监测系统采集压裂区的原始微震信号;
处理单元,用于根据地震波走时区分所述原始微震信号中的直达波信号和散射波信号;获取微震信号中直达波信号的质心频率和散射波信号的质心频率;根据所述直达波信号的质心频率和所述散射波信号的质心频率,得到压裂区的品质因子变化;利用散射波逆时偏移技术对散射体进行定位,确定裂缝发育区域;
评价单元,用于根据所述品质因子的变化和所述裂缝发育区域对所述压裂区的压裂效果进行评价;
所述处理单元,用于将直达波信号的振幅谱表示为:
Figure FDA0002894313580000022
式中,Z(f)为直达波信号的振幅谱,A为控制振幅大小的常数,k为控制振幅对称性的常数,f为直达波信号频率,f0为控制带宽的带宽因子;
基于直达波信号的能量谱Z(f)2得到直达波信号的质心频率表示为:
Figure FDA0002894313580000031
将散射波信号振幅谱表示为:
Figure FDA0002894313580000032
式中,S(f)是散射波信号的振幅谱,Δt是直达波和散射波信号的旅行时差,Q为品质因子;
基于散射波信号的能量谱S(f)2得到散射波信号的质心频率表示为:
Figure FDA0002894313580000033
所述处理单元,还用于对每个微震事件,将每个监测点记录的散射波信号在时间轴上反转,并作为新的震源向外散发地震波,公式如下:
R(x,z,t)=D(x,τ0-t)*G(x,z,t)
式中R(x,z,t)为反传波场,D(x,τ0-t)为地震记录,τ0为所有监测点记录到该微震信号的最迟时间,G(x,z,t)是格林函数,*表示褶积;
将所有监测点向外散发的地震波场进行互相关操作,并根据预设的成像条件,定位散射体位置,即裂缝发育区域,公式如下:
Figure FDA0002894313580000034
式中,t表示时间;R表示空间坐标(x,z)处的反向传播检波器波场;I(x,z)表示(x,z)处的成像计算结果;n为检波器的个数。
3.一种基于微震散射波品质因子的压裂效果评价设备,其特征在于,包括:至少一个处理器和存储器;
所述存储器存储计算机执行指令;
所述至少一个处理器执行所述存储器存储的计算机执行指令,使得所述至少一个处理器执行如权利要求1所述的基于微震散射波品质因子的压裂效果评价方法。
4.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有计算机执行指令,当处理器执行所述计算机执行指令时,实现如权利要求1所述的基于微震散射波品质因子的压裂效果评价方法。
CN202010361714.0A 2020-04-30 2020-04-30 基于微震散射波品质因子的压裂效果评价方法及装置 Expired - Fee Related CN111624661B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010361714.0A CN111624661B (zh) 2020-04-30 2020-04-30 基于微震散射波品质因子的压裂效果评价方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010361714.0A CN111624661B (zh) 2020-04-30 2020-04-30 基于微震散射波品质因子的压裂效果评价方法及装置

Publications (2)

Publication Number Publication Date
CN111624661A CN111624661A (zh) 2020-09-04
CN111624661B true CN111624661B (zh) 2021-04-20

Family

ID=72258896

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010361714.0A Expired - Fee Related CN111624661B (zh) 2020-04-30 2020-04-30 基于微震散射波品质因子的压裂效果评价方法及装置

Country Status (1)

Country Link
CN (1) CN111624661B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112859163B (zh) * 2021-01-07 2022-08-30 中国石油大学(北京) 利用散射波确定压裂区地层品质因子变化的方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4354381A (en) * 1980-12-30 1982-10-19 Mobil Oil Corporation Sonic attenuating method for distinguishing between oil, water and gas
CN105652321A (zh) * 2015-12-30 2016-06-08 中国石油大学(华东) 一种粘声各向异性最小二乘逆时偏移成像方法
CN106646616A (zh) * 2017-01-19 2017-05-10 中国石油大学(北京) 一种补偿地层吸收衰减的方法及装置
CN107290722A (zh) * 2017-06-29 2017-10-24 中国石油大学(北京) 微震源的定位方法和装置
CN109490965A (zh) * 2018-10-15 2019-03-19 长江大学 一种定量评价地层非均匀性的方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4354381A (en) * 1980-12-30 1982-10-19 Mobil Oil Corporation Sonic attenuating method for distinguishing between oil, water and gas
CN105652321A (zh) * 2015-12-30 2016-06-08 中国石油大学(华东) 一种粘声各向异性最小二乘逆时偏移成像方法
CN106646616A (zh) * 2017-01-19 2017-05-10 中国石油大学(北京) 一种补偿地层吸收衰减的方法及装置
CN107290722A (zh) * 2017-06-29 2017-10-24 中国石油大学(北京) 微震源的定位方法和装置
CN109490965A (zh) * 2018-10-15 2019-03-19 长江大学 一种定量评价地层非均匀性的方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
含水合物地层地震波衰减特性及衰减系数估计方法研究;李传辉;《中国博士学位论文全文数据库 基础科学辑》;20151115;第71-87页 *
水力压裂裂缝微地震逆时偏移成像;林叶等;《中国地球科学联合学术年会2015》;20151231;第1534-1535页 *

Also Published As

Publication number Publication date
CN111624661A (zh) 2020-09-04

Similar Documents

Publication Publication Date Title
US9506339B2 (en) Active seismic monitoring of fracturing operations and determining characteristics of a subterranean body using pressure data and seismic data
US20160349389A1 (en) Method for developing a geomechanical model based on seismic data, well logs and sem analysis of horizontal and vertical drill cuttings
CN110529087B (zh) 一种评价地层水力压裂效果的方法及装置
US11467308B2 (en) Fibro: a fiber optic data processing software for unconventional reservoirs
MXPA06012109A (es) Mapeo de fractura microsismica que utiliza mediciones de sincronizacion de fuente sismica para calibracion de velocidad.
CN102089497B (zh) 使用被动地震发射监视生产油田中沿着流动通路移动的储层流体
JP2003516516A (ja) 地質流体貯留の4d透水性解析のための方法
US9810803B2 (en) Method for subsurface mapping using seismic emissions
CN101470211A (zh) 利用微震事件生成的s尾波来识别q因子
AU2016425662A1 (en) Multivariate analysis of seismic data, microseismic data, and petrophysical properties in fracture modeling
SA519410848B1 (ar) تقييم ثقوب حفرة مغلفة في مكامن رمل غاز مضغوطة مع تسجيل موجة ستونلي
US20080159074A1 (en) System and method for quality control of noisy data
CN111624661B (zh) 基于微震散射波品质因子的压裂效果评价方法及装置
CN111580159B (zh) 基于微震品质因子的压裂效果评价方法及装置
US11209559B2 (en) Method and system for analyzing a borehole using passive acoustic logging
CA2898192C (en) Method for determining fracture network volume using passive seismic signals
US20080162052A1 (en) System and Method for Filtering Seismic Noise Data to Analyze Seismic Events
CN111025393B (zh) 针对含薄煤层地层的储层预测方法、装置、设备及介质
CA2897052C (en) Method for determining fracture plane orientation using passive seismic signals
CN111427085B (zh) 一种碳酸盐岩层间岩溶储层预测方法及装置
CN114076990B (zh) 油页岩反射能量确定方法、系统、存储介质及电子设备
CN115184990B (zh) 一种微地震监测观测方法
CN112965109B (zh) 一种强反射屏蔽下的储层的预测方法及装置
CN117967335A (zh) Tbm卡机预警方法、装置、电子设备及存储介质
CN117784232A (zh) 一种多尺度天然裂缝刚度系数的计算方法及装置

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210420

CF01 Termination of patent right due to non-payment of annual fee