CN104597492A - 压裂中速度对微地震旅行时影响强度的确定方法及装置 - Google Patents
压裂中速度对微地震旅行时影响强度的确定方法及装置 Download PDFInfo
- Publication number
- CN104597492A CN104597492A CN201510048367.5A CN201510048367A CN104597492A CN 104597492 A CN104597492 A CN 104597492A CN 201510048367 A CN201510048367 A CN 201510048367A CN 104597492 A CN104597492 A CN 104597492A
- Authority
- CN
- China
- Prior art keywords
- micro
- seismic event
- tour
- whilst
- microearthquake
- 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.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例提供了一种压裂过程中速度对微地震旅行时影响强度的确定方法及装置,其中,该方法包括:根据初始速度模型获取每个微地震事件的三维位置;根据所述每个微地震事件的三维位置,计算所有微地震事件形成的包络面的空间分布,所述包络面内部的空间是压裂过程中地层速度发生变化的地层范围;计算所述每个微地震事件到达检波器的传播路径;根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,所述传播距离表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度。该方案实现了可以确定压裂过程中地层速度变化对每个微地震事件传播旅行时的影响强度。
Description
技术领域
本发明涉及地球物理勘探技术领域,特别涉及一种压裂中速度对微地震旅行时影响强度的确定方法及装置。
背景技术
在非常规油气开采过程中,尤其是页岩气的开发过程中,由于地层低孔低渗特性,使得常规油气开采技术无法获得可观产量,往往需要采用水力压裂方式提高采收率。水力压裂技术通过压裂设备高压将压裂液注入地层,高压流体破坏地层应力平衡状态,使岩石发生破裂。新增加的裂缝能够提供油气传输通道,有效改善地层渗透性和提高采收率。水力压裂过程中,通过检波器监测微地震信号,对微地震信号的进一步处理能够定位裂缝分布,指导下一阶段的压裂操作和制定相应的开发策略。
尽管微地震研究已较为成熟,但目前对微地震的研究都没有考虑压裂过程中的地层速度变化。在微地震研究中,一般利用测井、VSP(Vertical Seismic Profiling,垂直地震剖面)和标定炮数据构建初始速度模型,并基于该初始速度模型研究微地震传播特征。由于测井、VSP和标定炮数据都是在压裂前采集,因此构建的初始速度模型只能反映压裂前的地层信息,无法考虑压裂过程中地层速度的变化。在压裂过程中,高压流体会产生大量裂缝,同时也会改变地层孔隙压力,这些变化都会改变地层速度,进一步影响微地震的传播方位和旅行时。在微地震定位过程中,通过对接收到的微地震信号处理,获得震源的方位和距离信息,进一步确定裂缝所在位置,如果忽略压裂过程的影响,基于压裂前的初始速度模型分析微地震旅行时等信息,使得微地震旅行时等信息会存在误差,由于现有技术中无法评估压裂过程对微地震传播旅行时的影响强度,使得基于不确定可信度的微地震旅行时进行裂缝定位,会直接影响裂缝定位的精度。
发明内容
本发明实施例提供了一种压裂过程中速度对微地震旅行时影响强度的确定方法,以确定压裂过程中速度变化对微地震传播旅行时的影响强度。该方法包括:根据初始速度模型获取每个微地震事件的三维位置;根据所述每个微地震事件的三维位置,计算所有微地震事件形成的包络面的空间分布,所述包络面内部的空间是压裂过程中地层速度发生变化的地层范围;计算所述每个微地震事件到达检波器的传播路径;根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,所述传播距离表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度。
在一个实施例中,根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,包括:根据每个微地震事件的空间分布,在所述包络面中建立与该个微地震事件的空间分布相内切的三维椭圆;计算该个微地震事件的传播路径在所述三维椭圆内的传播距离。
在一个实施例中,计算所述每个微地震事件到达检波器的传播路径,包括:采用射线追踪算法计算所述每个微地震事件到达检波器的传播路径。
在一个实施例中,所述传播距离与压裂过程中地层速度变化对该个微地震事件传播旅行时的影响强度成正比。
在一个实施例中,通过以下公式确定每个微地震事件的三维位置: 其中,(X,Y,Z)为待求的微地震事件的震源Leni=VpVsΔti/(Vp-Vs)
坐标,(Ai,Bi,Ci)为第i个检波器的空间坐标,Len为检波器到微地震事件的震源的距离,Vp为地层纵波速度,Vs为地层横波速度,Δt为地层纵横波旅行时差。
在一个实施例中,通过以下公式确定每个微地震事件的三维位置: 其中,(X,Y,Z)为待求的
微地震事件的震源坐标,(A,B,C)为检波器的空间坐标,Vp为地层纵波速度,Vs为地层横波速度,Δt为地层纵横波旅行时差,Len为检波器到震源的距离,dy/dx为传播方向中x分量与y分量之比,dz/dl为传播z方向分量与水平分量之比,(α1,α2,α3)表示信号传播的方向,β为相应的权值因子。
在一个实施例中,通过以下公式确定每个微地震事件的三维位置:其中,J是误差值,N为用于微地震事件旅行时计算的检波器个数,M为用于微地震事件方位计算的检波器个数,t0i为基于第i个检波器数据计算得到的微地震发生时间,为基于第i个检波器计算得到的微地震方位,为微地震发生时间的平均值,为微地震方位的平均值。
本发明实施例还提供了一种压裂过程中速度对微地震旅行时影响强度的确定装置,以确定压裂过程对微地震传播旅行时的影响强度。该装置包括:获取模块,用于根据初始速度模型获取每个微地震事件的三维位置;包络面计算模块,用于根据所述每个微地震事件的三维位置,计算所有微地震事件形成的包络面的空间分布,所述包络面内部的空间是压裂过程中地层速度发生变化的地层范围;传播路径计算模块,用于计算所述每个微地震事件到达检波器的传播路径;传播距离计算模块,用于根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,所述传播距离表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度。
在一个实施例中,所述传播距离计算模块,包括:椭圆建立单元,用于根据每个微地震事件的空间分布,在所述包络面中建立与该个微地震事件的空间分布相内切的三维椭圆;传播距离计算单元,用于计算该个微地震事件的传播路径在所述三维椭圆内的传播距离。
在一个实施例中,所述传播路径计算模块,具体用于采用射线追踪算法计算所述每个微地震事件到达检波器的传播路径。
在一个实施例中,所述传播距离与压裂过程中地层速度变化对该个微地震事件传播旅行时的影响强度成正比。
在本发明实施例中,通过获取每个微地震事件的三维位置,来计算所有微地震事件形成的包络面的空间分布,在计算每个微地震事件到达检波器的传播路径后,计算每个微地震事件的传播路径在包络面中的传播距离,由于包络面内部的空间是压裂过程中地层速度发生变化的地层范围,因此,该传播距离可以表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度,从而可以实现确定压裂过程中地层速度变化对每个微地震事件传播旅行时的影响强度,进而根据传播距离可以评价每个微地震事件传播旅行时的可信度,进而避免基于不确定可信度的微地震旅行时进行裂缝定位,降低裂缝定位的精度的问题。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是本发明实施例提供的一种压裂过程中速度对微地震旅行时影响强度的方法的流程图;
图2是本发明实施例提供的一种包络面的空间分布示意图;
图3是本发明实施例提供的一种传播距离与微地震传播旅行时偏差的关系示意图;
图4是本发明实施例提供的一种压裂过程中速度对微地震旅行时影响强度的装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
在本发明实施例中,提供了一种压裂过程中速度对微地震旅行时影响强度的方法,如图1所示,该方法包括:
步骤101:根据初始速度模型获取每个微地震事件的三维位置;
步骤102:根据所述每个微地震事件的三维位置,计算所有微地震事件形成的包络面的空间分布,所述包络面内部的空间是压裂过程中地层速度发生变化的地层范围(也称为压裂带);
步骤103:计算所述每个微地震事件到达检波器的传播路径;
步骤104:根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,所述传播距离表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度。
由图1所示的流程可知,在本发明实施例中,通过获取每个微地震事件的三维位置,来计算所有微地震事件形成的包络面的空间分布,在计算每个微地震事件到达检波器的传播路径后,计算每个微地震事件的传播路径在包络面中的传播距离,由于包络面内部的空间是压裂过程中地层速度发生变化的地层范围,因此,该传播距离可以表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度,从而可以实现评估压裂过程中地层速度变化对每个微地震事件传播旅行时的影响强度,进而根据传播距离可以确定每个微地震事件传播旅行时的可信度,进而避免基于不确定可信度的微地震旅行时进行裂缝定位,降低裂缝定位的精度的问题。
具体实施时,在步骤101中可以通过测井数据、VSP或标定炮数据构建初始速度模型,从微地震记录中拾取纵波、横波初至时间以及传播的偏振方向,继而根据初始速度模型利用单点极化法(式(1)),多点联合反演法(式(2))或搜索法(式(3))计算定位每个微地震事件,以获取每个微地震事件的三维位置,具体的,如果采用单点极化法计算定位每个微地震事件,则可以根据以下公式(1)进行计算:
Leni=VpVsΔti/(Vp-Vs)
其中,(X,Y,Z)为待求的微地震事件的震源坐标,(Ai,Bi,Ci)为第i个检波器的空间坐标,Len为检波器到震源的距离,Vp为地层纵波速度,Vs为地层横波速度,Δt为地层纵横波旅行时差。
如果采用多点联合反演法计算定位每个微地震事件,则可以根据以下公式(2)进行计算:
Len=VpVsΔt/(Vp-Vs)
X=A-Len·α1
Y=B-Len·α2
Z=C-Len·α3 (2)
其中,(X,Y,Z)为待求的微地震事件的震源坐标,(A,B,C)为检波器的空间坐标,Vp为地层纵波速度,Vs为地层横波速度,Δt为地层纵横波旅行时差,Len为检波器到震源的距离,dy/dx为传播方向中x分量与y分量之比,dz/dl为传播z方向分量与水平分量之比,(α1,α2,α3)表示信号传播的方向,β为相应的权值因子。
如果采用搜索法计算定位每个微地震事件,则可以根据以下公式(3)进行计算:
其中,J是误差值,最小的J对应最佳的定位结果;N为用于微地震事件旅行时计算的检波器个数,M为用于微地震事件方位计算的检波器个数。t0i为基于第i个检波器数据计算得到的微地震发生时间,为基于第i个检波器计算得到的微地震方位,为微地震发生时间的平均值,为微地震方位的平均值。为保证时间项与方位项的权重相当,λt和为相应的归一化因子。
具体实施时,在步骤103中,采用射线追踪算法计算所述每个微地震事件到达检波器的传播路径,记录微地震信号从产生到传播至检波器整个过程中微地震信号传播所经历的空间位置。
具体实施时,以在压裂过程前为均匀介质,速度为3000m/s的初始速度模型为例,基于该初始速度模型获取的微地震事件及其包络面三维空间分布如图2所示,图2中‘*’为微地震事件,虚线表示包络面,在压裂过程中,计算每个微地震事件的传播路径在包络面中的传播距离与该个微地震事件传播旅行时偏差的关系,该微地震事件传播旅行时偏差是指基于压裂前的初始速度模型和压裂过程中实时速度模型正演结果的差异,如图3所示,图3中横坐标为微地震事件在包络面内的传播距离,单位是米,纵坐标为微地震事件旅行时偏差,单位为毫秒,由图3可知,微地震事件在包络面内的传播距离与该微地震事件旅行时偏差存在正相关关系,也就是微地震事件在包络面内的传播距离与压裂过程中地层速度变化对该个微地震事件传播旅行时的影响强度成正比,传播距离数值越大,表明压裂过程中地层速度变化对该微地震事件传播旅行时的影响强度越大,该数据用于反演时会产生较大的偏差,即该数据的可信度低,因此,微地震事件在包络面内的传播距离可以作为评价微地震数据可信度的依据。
具体的,根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,包括:根据每个微地震事件的空间分布,在所述包络面中建立与该个微地震事件的空间分布相内切的三维椭圆;计算该个微地震事件的传播路径在所述三维椭圆内的传播距离。
基于同一发明构思,本发明实施例中还提供了一种压裂过程中速度对微地震旅行时影响强度的装置,如下面的实施例所述。由于压裂过程中速度对微地震旅行时影响强度的装置解决问题的原理与压裂过程中速度对微地震旅行时影响强度的方法相似,因此压裂过程中速度对微地震旅行时影响强度的装置的实施可以参见压裂过程中速度对微地震旅行时影响强度的方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图4是本发明实施例的压裂过程中速度对微地震旅行时影响强度的装置的一种结构框图,如图4所示,包括:获取模块401、包络面计算模块402、传播路径计算模块403和传播距离计算模块404,下面对该结构进行说明。
获取模块401,用于根据初始速度模型获取每个微地震事件的三维位置;
包络面计算模块402,与获取模块401连接,用于根据所述每个微地震事件的三维位置,计算所有微地震事件形成的包络面的空间分布,所述包络面内部的空间是压裂过程中地层速度发生变化的地层范围;
传播路径计算模块403,与包络面计算模块402连接,用于计算所述每个微地震事件到达检波器的传播路径;
传播距离计算模块404,与传播路径计算模块403连接,用于根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,所述传播距离表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度。
在一个实施例中,所述传播距离计算模块404,包括:椭圆建立单元,用于根据每个微地震事件的空间分布,在所述包络面中建立与该个微地震事件的空间分布相内切的三维椭圆;传播距离计算单元,与椭圆建立单元连接,用于计算该个微地震事件的传播路径在所述三维椭圆内的传播距离。
在一个实施例中,所述传播路径计算模块403,具体用于采用射线追踪算法计算所述每个微地震事件到达检波器的传播路径。
在一个实施例中,所述传播距离与压裂过程中地层速度变化对该个微地震事件传播旅行时的影响强度成正比。
在本发明实施例中,通过获取每个微地震事件的三维位置,来计算所有微地震事件形成的包络面的空间分布,在计算每个微地震事件到达检波器的传播路径后,计算每个微地震事件的传播路径在包络面中的传播距离,由于包络面内部的空间是压裂过程中地层速度发生变化的地层范围,因此,该传播距离可以表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度,从而可以实现确定压裂过程中地层速度变化对每个微地震事件传播旅行时的影响强度,进而根据传播距离可以评价每个微地震事件传播旅行时的可信度,进而避免基于不确定可信度的微地震旅行时进行裂缝定位,降低裂缝定位的精度的问题。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (11)
1.一种压裂过程中速度对微地震旅行时影响强度的确定方法,其特征在于,包括:
根据初始速度模型获取每个微地震事件的三维位置;
根据所述每个微地震事件的三维位置,计算所有微地震事件形成的包络面的空间分布,所述包络面内部的空间是压裂过程中地层速度发生变化的地层范围;
计算所述每个微地震事件到达检波器的传播路径;
根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,所述传播距离表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度。
2.如权利要求1所述的压裂过程中速度对微地震旅行时影响强度的确定方法,其特征在于,根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,包括:
根据每个微地震事件的空间分布,在所述包络面中建立与该个微地震事件的空间分布相内切的三维椭圆;
计算该个微地震事件的传播路径在所述三维椭圆内的传播距离。
3.如权利要求1所述的压裂过程中速度对微地震旅行时影响强度的确定方法,其特征在于,计算所述每个微地震事件到达检波器的传播路径,包括:采用射线追踪算法计算所述每个微地震事件到达检波器的传播路径。
4.如权利要求1所述的压裂过程中速度对微地震旅行时影响强度的确定方法,其特征在于,所述传播距离与压裂过程中地层速度变化对该个微地震事件传播旅行时的影响强度成正比。
5.如权利要求1至4中任一项所述的压裂过程中速度对微地震旅行时影响强度的确定方法,其特征在于,通过以下公式确定每个微地震事件的三维位置:
Leni=VpVsΔti/(Vp-Vs)
其中,(X,Y,Z)为待求的微地震事件的震源坐标,(Ai,Bi,Ci)为第i个检波器的空间坐标,Len为检波器到微地震事件的震源的距离,Vp为地层纵波速度,Vs为地层横波速度,Δt为地层纵横波旅行时差。
6.如权利要求1至4中任一项所述的压裂过程中速度对微地震旅行时影响强度的确定方法,其特征在于,通过以下公式确定每个微地震事件的三维位置:
Len=VpVsΔt/(Vp-Vs)
X=A-Len·α1
Y=B-Len·α2
Z=C-Len·α3
α1=1/β,α2=(dy/dx)/β,
其中,(X,Y,Z)为待求的微地震事件的震源坐标,(A,B,C)为检波器的空间坐标,Vp为地层纵波速度,Vs为地层横波速度,Δt为地层纵横波旅行时差,Len为检波器到震源的距离,dy/dx为传播方向中x分量与y分量之比,dz/dl为传播z方向分量与水平分量之比,(α1,α2,α3)表示信号传播的方向,β为相应的权值因子。
7.如权利要求1至4中任一项所述的压裂过程中速度对微地震旅行时影响强度的确定方法,其特征在于,通过以下公式确定每个微地震事件的三维位置:
其中,J是误差值,N为用于微地震事件旅行时计算的检波器个数,M为用于微地震事件方位计算的检波器个数,t0i为基于第i个检波器数据计算得到的微地震发生时间,为基于第i个检波器计算得到的微地震方位,为微地震发生时间的平均值,为微地震方位的平均值。
8.一种压裂过程中速度对微地震旅行时影响强度的确定装置,其特征在于,包括:
获取模块,用于根据初始速度模型获取每个微地震事件的三维位置;
包络面计算模块,用于根据所述每个微地震事件的三维位置,计算所有微地震事件形成的包络面的空间分布,所述包络面内部的空间是压裂过程中地层速度发生变化的地层范围;
传播路径计算模块,用于计算所述每个微地震事件到达检波器的传播路径;
传播距离计算模块,用于根据所述包络面的空间分布和所述每个微地震事件到达检波器的传播路径,计算所述每个微地震事件的传播路径在所述包络面中的传播距离,所述传播距离表示压裂过程中地层速度变化对该个微地震事件传播旅行时影响的强度。
9.如权利要求8所述的压裂过程中速度对微地震旅行时影响强度的确定装置,其特征在于,所述传播距离计算模块,包括:
椭圆建立单元,用于根据每个微地震事件的空间分布,在所述包络面中建立与该个微地震事件的空间分布相内切的三维椭圆;
传播距离计算单元,用于计算该个微地震事件的传播路径在所述三维椭圆内的传播距离。
10.如权利要求8所述的压裂过程中速度对微地震旅行时影响强度的确定装置,其特征在于,所述传播路径计算模块,具体用于采用射线追踪算法计算所述每个微地震事件到达检波器的传播路径。
11.如权利要求8至10中任一项所述的压裂过程中速度对微地震旅行时影响强度的确定装置,其特征在于,所述传播距离与压裂过程中地层速度变化对该个微地震事件传播旅行时的影响强度成正比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510048367.5A CN104597492A (zh) | 2015-01-30 | 2015-01-30 | 压裂中速度对微地震旅行时影响强度的确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510048367.5A CN104597492A (zh) | 2015-01-30 | 2015-01-30 | 压裂中速度对微地震旅行时影响强度的确定方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104597492A true CN104597492A (zh) | 2015-05-06 |
Family
ID=53123397
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510048367.5A Pending CN104597492A (zh) | 2015-01-30 | 2015-01-30 | 压裂中速度对微地震旅行时影响强度的确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104597492A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105277971A (zh) * | 2015-10-16 | 2016-01-27 | 中国石油天然气集团公司 | 一种微地震监测系统及方法 |
CN110886599A (zh) * | 2018-09-07 | 2020-03-17 | 中国石油化工股份有限公司 | 基于破裂速度的非压裂事件识别方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101173989A (zh) * | 2006-11-03 | 2008-05-07 | 中国石油集团东方地球物理勘探有限责任公司 | 地震采集三维观测系统定量分析方法 |
CN102841379A (zh) * | 2012-09-06 | 2012-12-26 | 中国石油大学(华东) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 |
CN102879820A (zh) * | 2012-09-20 | 2013-01-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于三角网格的三维表层模型构建方法 |
CN103869363A (zh) * | 2014-03-20 | 2014-06-18 | 中国石油天然气集团公司 | 微地震定位方法及装置 |
EP1938281B1 (en) * | 2004-09-10 | 2015-01-21 | ExxonMobil Upstream Research Company | Geologic models of subsurface sedimentray volumes |
-
2015
- 2015-01-30 CN CN201510048367.5A patent/CN104597492A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1938281B1 (en) * | 2004-09-10 | 2015-01-21 | ExxonMobil Upstream Research Company | Geologic models of subsurface sedimentray volumes |
CN101173989A (zh) * | 2006-11-03 | 2008-05-07 | 中国石油集团东方地球物理勘探有限责任公司 | 地震采集三维观测系统定量分析方法 |
CN102841379A (zh) * | 2012-09-06 | 2012-12-26 | 中国石油大学(华东) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 |
CN102879820A (zh) * | 2012-09-20 | 2013-01-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于三角网格的三维表层模型构建方法 |
CN103869363A (zh) * | 2014-03-20 | 2014-06-18 | 中国石油天然气集团公司 | 微地震定位方法及装置 |
Non-Patent Citations (2)
Title |
---|
宋维琪,等: "解域约束下的微地震事件网格搜索法、遗传算法联合反演", 《石油地球物理勘探》 * |
张晓林,等: "水力压裂对速度场及微地震定位的影响", 《地球物理学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105277971A (zh) * | 2015-10-16 | 2016-01-27 | 中国石油天然气集团公司 | 一种微地震监测系统及方法 |
CN110886599A (zh) * | 2018-09-07 | 2020-03-17 | 中国石油化工股份有限公司 | 基于破裂速度的非压裂事件识别方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103869363B (zh) | 微地震定位方法及装置 | |
US10670754B2 (en) | System and method for processing microseismic data | |
US10359529B2 (en) | Singularity spectrum analysis of microseismic data | |
CN104199090B (zh) | 一种地面监测微地震定位的速度模型构建与求解方法 | |
EP3090279B1 (en) | Methods and systems of determining a fault plane of a microseismic event | |
CN105068117B (zh) | 用于裂缝性介质的avo反演方法、装置和设备 | |
CN103105624B (zh) | 基于数据库技术的纵横波时差定位方法 | |
CN104267429A (zh) | 确定地层压力的方法及装置 | |
CN104375188A (zh) | 一种地震波透射衰减补偿方法及装置 | |
CN105277978A (zh) | 一种确定近地表速度模型的方法及装置 | |
CN103245972B (zh) | 一种在二维空间内确定复杂地质构造的方法 | |
CN104360389A (zh) | 一种致密砂岩储层岩石弹性模量计算方法 | |
AU2014407527B2 (en) | Integrating vertical seismic profile data for microseismic anisotropy velocity analysis | |
CN105089615A (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
CN104730574A (zh) | 构建近地表结构模型的方法 | |
CN104166159B (zh) | 四维微地震监测的裂缝形态处理方法和系统 | |
US20160238724A1 (en) | Methods and systems of generating a velocity model | |
CN105626056A (zh) | 一种确定目标区三维地层孔隙压力的方法和装置 | |
CN104749630A (zh) | 构建微地震监测速度模型的方法 | |
CN104597492A (zh) | 压裂中速度对微地震旅行时影响强度的确定方法及装置 | |
CN103984018A (zh) | 多波联合振幅随入射角变化反演方法 | |
CN106353807A (zh) | 裂缝识别方法和装置 | |
US11585212B2 (en) | Anisotropy model guided fracture properties extraction from VSP data | |
CN103628866A (zh) | 获取地层的地层压力系数的方法 | |
CN105301657A (zh) | 一种基于岩石物理意义的曲线校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150506 |