CN110133721B - 一种水力压裂过程监测方法及系统 - Google Patents

一种水力压裂过程监测方法及系统 Download PDF

Info

Publication number
CN110133721B
CN110133721B CN201910481804.0A CN201910481804A CN110133721B CN 110133721 B CN110133721 B CN 110133721B CN 201910481804 A CN201910481804 A CN 201910481804A CN 110133721 B CN110133721 B CN 110133721B
Authority
CN
China
Prior art keywords
wave field
seismic
field data
data model
data
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.)
Active
Application number
CN201910481804.0A
Other languages
English (en)
Other versions
CN110133721A (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.)
Nanjing Garboro Information Technology Co ltd
Original Assignee
Nanjing Garboro Information Technology Co ltd
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 Nanjing Garboro Information Technology Co ltd filed Critical Nanjing Garboro Information Technology Co ltd
Priority to CN201910481804.0A priority Critical patent/CN110133721B/zh
Publication of CN110133721A publication Critical patent/CN110133721A/zh
Application granted granted Critical
Publication of CN110133721B publication Critical patent/CN110133721B/zh
Active 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/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
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters
    • G01V2210/6248Pore pressure

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开一种水力压裂过程监测方法及系统。方法包括:获取地震总波场数据,所述地震总波场数据包括含有水力压裂裂缝的管波波场数据和噪声干扰数据;对所述地震总波场数据采用信号分离的方法,得到地震波场数据;对所述地震波场数据进行管波波场反演,得到水力压裂的裂缝宽度。采用本发明的方法不仅能够降低施工难度,节省成本,而且能够达到及时、准确监测压裂效果的目的。

Description

一种水力压裂过程监测方法及系统
技术领域
本发明涉及水力压裂领域,特别是涉及一种水力压裂过程监测方法及系统。
背景技术
在水力压裂生产中,压裂效果的及时获取和压裂产生裂缝的分布,是业内尚未有效解决的难题。现阶段,对压裂裂缝的检测方法可以分为非震方法裂缝检测及地震学方法裂缝检测。
在非震方法中常用的有示踪剂方法、电位法及倾斜仪阵列水准观测方法。示踪剂方法将特定化学物理物质随压裂过程注入,并在压裂井或邻井观测示踪剂的分布,该方法监测在时间上滞后较久,观测范围有限;电位法测量压裂前后的电位差,反推引起这种差别的压裂液分布,它受气候、深度限制;倾斜仪阵列水准观测测量压裂前后的地表或压裂段周围的形变,反演地下压裂裂缝体积的位置和形状,反演的多解性较强。
地震学方法施工灵活、精度较高且及时性较好,得到广泛应用,其中常用的观测手段有井下邻近压裂段的观测、检波器阵排列检测及微地震法观测。井下邻近压裂段的观测方法信噪比及观测精度较高,但是邻近井需要较长时间停产观测,设备成本高昂,裂缝检测及时性差;检波器阵排列检测与勘探地震方法相似,所以裂缝监测精度较高,但该方法检波器布设规模很大,费用较高,施工费时,同时后期数据处理时间较长,裂缝校测及时性较差;微地震方法应用单分量检波器地表或邻近井中布设,施工和处理较为简单,成本较低,但是易出现观测数据信噪比低的情况,导致裂缝检测结果可靠性降低。
在最新的研究中,通过在压裂过程中对沿井壁传播的管波进行观测,施工简单,可以较为及时的获得压裂程度及压裂裂缝渗透率等重要参数,得到了较好的应用效果。但是,在压裂的高压环境进行管波的测量对检测设备的要求高,成本较高。
发明内容
本发明的目的是提供一种水力压裂过程监测方法及系统,不仅能够降低施工难度,节省成本,而且能够达到及时、准确监测压裂效果的目的。
为实现上述目的,本发明提供了如下方案:
一种水力压裂过程监测方法,包括:
获取地震总波场数据,所述地震总波场数据包括含有水力压裂裂缝的管波波场数据和噪声干扰数据;
对所述地震总波场数据采用信号分离的方法,得到地震波场数据;
对所述地震波场数据进行管波波场反演,得到水力压裂的裂缝宽度。
可选的,所述获取地震总波场数据,具体包括:
通过MEMS传感器贴井壁测量获取地震总波场数据。
可选的,所述对所述地震总波场数据采用信号分离的方法,得到信号分离后的地震波场数据,具体包括:
对所述地震总波场数据采用中值滤波方法进行处理,得到地震波场数据P(t);
Figure BDA0002084085830000021
其中,pi(x)(x=1,2,3...n)为第i个采样点pi邻域范围的n个采样点,pi为邻域中点,n为采样点的数量。
可选的,所述对所述地震波场数据进行管波波场反演,得到水力压裂裂缝相关的管波信息,具体包括:
将所述地震波场数据P(t)通过傅里叶变换,得到频率域地震波场数据模型P(ω);
获取地震波场拟合数据模型
Figure BDA0002084085830000022
通过所述地震波场拟合数据模型
Figure BDA0002084085830000023
对所述频率域地震波场数据模型P(ω)进行修正,得到最佳匹配地震波场数据模型Pm(ω);
根据所述最佳匹配地震波场数据模型Pm(ω)进行傅里叶反变换,得到时间域的最佳地震波场数据模型
Figure BDA0002084085830000031
根据所述最佳地震波场数据模型
Figure BDA0002084085830000032
得到裂缝导流能力;
根据所述裂缝导流能力,确定所述裂缝宽度。
一种水力压裂过程监测系统,包括:
获取模块,用于获取地震总波场数据,所述地震总波场数据包括含有水力压裂裂缝的管波波场数据和噪声干扰数据;
信号分离模块,用于对所述地震总波场数据采用信号分离的方法,得到地震波场数据;
反演模块,用于对所述地震波场数据进行管波波场反演,得到水力压裂的裂缝宽度。
可选的,所述获取模块,具体包括:
第一获取单元,用于通过MEM传感器贴井壁测量获取地震总波场数据。
可选的,所述信号分离模块,具体包括:
信号分离单元,用于对所述地震总波场数据采用中值滤波方法进行处理,得到地震波场数据P(t);
Figure BDA0002084085830000033
其中,pi(x)(x=1,2,3...n)为第i个采样点pi邻域范围的n个采样点,pi为邻域中点,n为采样点的数量。
可选的,所述反演模块,具体包括:
第一变换单元,用于将所述地震波场数据P(t)通过傅里叶变换,得到频率域地震波场数据模型P(ω);
第二获取单元,用于获取地震波场拟合数据模型
Figure BDA0002084085830000034
修正单元,用于通过所述地震波场拟合数据模型
Figure BDA0002084085830000035
对所述频率域地震波场数据模型P(ω)进行修正,得到最佳匹配地震波场数据模型Pm(ω);
第二变换单元,用于根据所述最佳匹配地震波场数据模型Pm(ω)进行傅里叶反变换,得到时间域的最佳地震波场数据模型
Figure BDA0002084085830000041
裂缝导流能力确定单元,用于根据所述最佳地震波场数据模型
Figure BDA0002084085830000042
确定裂缝导流能力;
裂缝宽度确定单元,用于根据所述裂缝导流能力,确定所述裂缝宽度。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供一种水力压裂过程监测方法,通过获取地震总波场数据;对所述地震总波场数据采用信号分离的方法,得到地震波场数据;对所述地震波场数据进行管波波场反演,得到水力压裂的裂缝宽度。通过本发明的上述步骤能够对采集的地震总波场数据进行信号分离处理,对分离得到的地震波场数据进行反演计算,得到压裂裂缝参数,本发明不仅能够降低施工难度,节省成本,而且能够及时、准确的获取压裂效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明水力压裂过程监测方法流程图;
图2为本发明中值滤波处理前的地震总波场数据图;
图3为本发明中值滤波处理后的地震波场数据图;
图4为本发明水力压裂反演结果图一;
图5为本发明水力压裂反演结果图二;
图6为本发明水力压裂过程监测系统结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种水力压裂过程监测方法及系统,不仅能够降低施工难度,节省成本,而且能够达到及时、准确监测压裂效果的目的。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明水力压裂过程监测方法流程图。如图1所示,一种水力压裂过程监测方法包括:
步骤101:获取地震总波场数据,所述地震总波场数据包括含有水力压裂裂缝的管波波场数据和噪声干扰数据;通过在井壁布设MEMS检波器,对压裂过程中产生的震动进行实时采集与监测,水力压裂过程中,压裂裂缝产生后会产生与其相关的地震波场,以管波的形式向井口传播。通过在井筒壁布设检波器,可以获取包括含有水力压裂裂缝的管波波场数据和噪声干扰数据的地震总波场数据。
步骤102:对所述地震总波场数据采用信号分离的方法,得到地震波场数据;具体包括:
对所述地震总波场数据采用中值滤波方法进行处理,得到地震波场数据P(t);
Figure BDA0002084085830000051
其中,pi(x)(x=1,2,3...n)为第i个采样点pi邻域范围的n个采样点,pi为邻域中点,n为采样点的数量。
通过针对地震总波场中的噪声进行压制,得到水力压裂裂缝相关的管波地震波场,便于压裂裂缝参数反演的精确计算。
以地震总波场中随机噪音压制情况为例进行说明。对原始采集的地震总波场数据通过中值滤波的方法对信号进行随机噪音压制(如公式1和公式2所示)。
中值滤波的过程是把序列或数字图像中一点的值,用该点邻域中各点值的中值来替代。在本次数据采集中,所采集的数据p1p2p3…pn为地震采集序列中n个采样点,pi(x)(x=1,2,3...n)为第i个采样点pi邻域范围的n个采样点,其中pi为邻域中点,将邻域内各个采集点按大小顺序排列起来,如式(1)所示:
pi(1)≤pi(2)≤pi(3)…≤pi(n) (1)
则该序列中值P为
Figure BDA0002084085830000061
P(t)代表采集数据经过中值滤波处理以后得到的新数据,t=1,2,3...n。
将每个点数据取其领域各点值的中值作为该点的值,即可完成中值滤波处理,得到处理后的管波地震波场数据,图2为本发明中值滤波处理前的地震总波场数据图。图3为本发明中值滤波处理后的地震波场数据图。采用中值滤波的方法能够达到管波地震波场数据获取增强的目的,同时保留有效地震数据的震动能量相对强弱关系,以便后续进行定量裂缝参数反演。
步骤103:对所述地震波场数据进行管波波场反演,得到水力压裂的裂缝宽度,具体包括:
将所述地震波场数据P(t)通过傅里叶变换,得到频率域地震波场数据模型P(ω);
获取地震波场拟合数据模型
Figure BDA0002084085830000062
通过所述地震波场拟合数据模型
Figure BDA0002084085830000063
对所述频率域地震波场数据模型P(ω)进行修正,得到最佳匹配地震波场数据模型Pm(ω);
根据所述最佳匹配地震波场数据模型Pm(ω)进行傅里叶反变换,得到时间域的最佳地震波场数据模型
Figure BDA0002084085830000064
根据所述最佳地震波场数据模型
Figure BDA0002084085830000065
得到裂缝导流能力;
根据所述裂缝导流能力,确定所述裂缝宽度。
本发明中,以基于管波波动方程反演技术的裂缝参数反演为例进行说明。
在反演中,定义裂缝导流能力为kw,其中k是压裂段地层渗透率,w为裂缝宽度。通过建立反演模型,充分考虑填充支撑剂的裂缝内流体流动规律,裂缝的弹性柔度以及裂缝和钻孔之间的耦合关系。最后通过不断的参数迭代,将步骤102中的处理数据与计算模型在频率域进行拟合,得到估计的裂缝参数值,通过佳地震波场数据
Figure BDA0002084085830000066
完成反演裂缝导流能力。
在管波经过水力压裂裂缝系统并反射的过程中,其地震波场的传播规律遵循波动方程。由于管波的波长足够大,压裂裂缝在与钻孔的交界处都有相同的压力响应。因此,管波的反射波是压裂裂缝系统的整体反射,而不是单个裂缝的反射特征。在计算过程中,首先将采集得到的并经过中值滤波的实测数据P(t)转换到频率域,即
P(ω)=FFT(P(t)) (3)
其中P(ω)为所采集数据的频域表达式。
模型的拟合数据
Figure BDA0002084085830000071
在频率域的计算公式如下
Figure BDA0002084085830000072
其中,ρT为井孔中流体密度,cT为管波速度,AT为钻孔横截面积,N为裂缝个数,G*=G/(1-v),G为固体剪切模量,v为泊松比,μ为流体动态粘度,
Figure BDA0002084085830000073
为孔隙度,k为渗透性,w为裂缝宽度,H为裂缝高度。Q(t)是井口流量随时间变化的函数,其频率域表达形式为
Figure BDA0002084085830000074
g(ω)为考虑衰减和因果关系保持色散的双向旅行时间因子,其表达式为
Figure BDA0002084085830000075
其中,h是钻孔长度,ω0是管波相速度等于cT时的参考角频率。
由(4)式可看出,本发明给定管波速度cT、井孔长度h、裂缝个数N、裂缝高度H、孔隙度
Figure BDA0002084085830000076
及流体动态粘度μ以后,通过不断改变kw的值,来使得实测数据频谱P(ω)与模型计算的
Figure BDA0002084085830000077
误差最小,从而使得步骤102得到的经过处理的地震数据与模型计算得到的数据匹配最佳,即得到最佳匹配地震波场数据模型Pm(ω)。将得到最佳匹配地震波场数据模型Pm(ω)进行反傅里叶变化得到时间域的最佳地震波场数据模型
Figure BDA0002084085830000078
Figure BDA0002084085830000079
图4为本发明水力压裂反演结果图一;图5为本发明水力压裂反演结果图二。通过图4和图5可知,将得到的最佳地震波场数据模型
Figure BDA0002084085830000081
与经过中值滤波后的地震波场数据P(t)进行对比,模型与实测数据匹配良好。通过反演确定裂缝导流能力kw,在水力压裂渗透率k已知的情况下,可以进一步反演出裂缝宽度w。
本发明能够实现在井口或地表观测,并采集在压裂过程中地震波数据,从而达到降低成本的目的;本发明能够对采集的地震总波场数据进行信号分离处理,通过对分离得到的地震波场数据进行反演计算,得到压裂裂缝参数,及时、准确的获取压裂效果。
图6为本发明水力压裂过程监测系统结构图。如图6所示,一种水力压裂过程监测系统包括:
获取模块601,用于获取地震总波场数据,所述地震总波场数据包括含有水力压裂裂缝的管波波场数据和噪声干扰数据;
信号分离模块602,用于对所述地震总波场数据采用信号分离的方法,得到地震波场数据;
反演模块603,用于对所述地震波场数据进行管波波场反演,得到水力压裂的裂缝宽度。
所述获取模块601,具体包括:
第一获取单元,用于通过MEM传感器贴井壁测量获取地震总波场数据。
所述信号分离模块602,具体包括:
信号分离单元,用于对所述地震总波场数据采用中值滤波方法进行处理,得到地震波场数据P(t);
Figure BDA0002084085830000082
其中,pi(x)(x=1,2,3...n)为第i个采样点pi邻域范围的n个采样点,pi为邻域中点,n为采样点的数量。
所述反演模块603,具体包括:
第一变换单元,用于将所述地震波场数据P(t)通过傅里叶变换,得到频率域地震波场数据模型P(ω);
第二获取单元,用于获取地震波场拟合数据模型
Figure BDA0002084085830000091
修正单元,用于通过所述地震波场拟合数据模型
Figure BDA0002084085830000092
对所述频率域地震波场数据模型P(ω)进行修正,得到最佳匹配地震波场数据模型Pm(ω);
第二变换单元,用于根据所述最佳匹配地震波场数据模型Pm(ω)进行傅里叶反变换,得到时间域的最佳地震波场数据模型
Figure BDA0002084085830000093
裂缝导流能力确定单元,用于根据所述最佳地震波场数据模型
Figure BDA0002084085830000094
确定裂缝导流能力;
裂缝宽度确定单元,用于根据所述裂缝导流能力,确定所述裂缝宽度。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (4)

1.一种水力压裂过程监测方法,其特征在于,包括:
获取地震总波场数据,所述地震总波场数据包括含有水力压裂裂缝的管波波场数据和噪声干扰数据;
对所述地震总波场数据采用信号分离的方法,得到地震波场数据,具体包括:
对所述地震总波场数据采用中值滤波方法进行处理,得到地震波场数据P(t);
Figure FDA0002381472660000011
其中,pi(x)(x=1,2,3...n)为第i个采样点pi邻域范围的n个采样点,pi为邻域中点,n为采样点的数量;
对所述地震波场数据进行管波波场反演,得到水力压裂的裂缝宽度,具体包括:
将所述地震波场数据P(t)通过傅里叶变换,得到频率域地震波场数据模型P(ω);
获取地震波场拟合数据模型
Figure FDA0002381472660000012
通过所述地震波场拟合数据模型
Figure FDA0002381472660000013
对所述频率域地震波场数据模型P(ω)进行修正,得到最佳匹配地震波场数据模型Pm(ω);
根据所述最佳匹配地震波场数据模型Pm(ω)进行傅里叶反变换,得到时间域的最佳地震波场数据模型
Figure FDA0002381472660000014
根据所述最佳地震波场数据模型
Figure FDA0002381472660000015
得到裂缝导流能力;
根据所述裂缝导流能力,确定所述裂缝宽度。
2.根据权利要求1所述的水力压裂过程监测方法,其特征在于,所述获取地震总波场数据,具体包括:
通过MEMS传感器贴井壁测量获取地震总波场数据。
3.一种水力压裂过程监测系统,其特征在于,包括:
获取模块,用于获取地震总波场数据,所述地震总波场数据包括含有水力压裂裂缝的管波波场数据和噪声干扰数据;
信号分离模块,用于对所述地震总波场数据采用信号分离的方法,得到地震波场数据;
反演模块,用于对所述地震波场数据进行管波波场反演,得到水力压裂的裂缝宽度;
所述信号分离模块,具体包括:
信号分离单元,用于对所述地震总波场数据采用中值滤波方法进行处理,得到地震波场数据P(t);
Figure FDA0002381472660000021
其中,pi(x)(x=1,2,3...n)为第i个采样点pi邻域范围的n个采样点,pi为邻域中点,n为采样点的数量;
所述反演模块,具体包括:
第一变换单元,用于将所述地震波场数据P(t)通过傅里叶变换,得到频率域地震波场数据模型P(ω);
第二获取单元,用于获取地震波场拟合数据模型
Figure FDA0002381472660000024
修正单元,用于通过所述地震波场拟合数据模型
Figure FDA0002381472660000025
对所述频率域地震波场数据模型P(ω)进行修正,得到最佳匹配地震波场数据模型Pm(ω);
第二变换单元,用于根据所述最佳匹配地震波场数据模型Pm(ω)进行傅里叶反变换,得到时间域的最佳地震波场数据模型
Figure FDA0002381472660000022
裂缝导流能力确定单元,用于根据所述最佳地震波场数据模型
Figure FDA0002381472660000023
确定裂缝导流能力;
裂缝宽度确定单元,用于根据所述裂缝导流能力,确定所述裂缝宽度。
4.根据权利要求3所述的水力压裂过程监测系统,其特征在于,所述获取模块,具体包括:
第一获取单元,用于通过MEMS传感器贴井壁测量获取地震总波场数据。
CN201910481804.0A 2019-06-04 2019-06-04 一种水力压裂过程监测方法及系统 Active CN110133721B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910481804.0A CN110133721B (zh) 2019-06-04 2019-06-04 一种水力压裂过程监测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910481804.0A CN110133721B (zh) 2019-06-04 2019-06-04 一种水力压裂过程监测方法及系统

Publications (2)

Publication Number Publication Date
CN110133721A CN110133721A (zh) 2019-08-16
CN110133721B true CN110133721B (zh) 2020-06-16

Family

ID=67580040

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910481804.0A Active CN110133721B (zh) 2019-06-04 2019-06-04 一种水力压裂过程监测方法及系统

Country Status (1)

Country Link
CN (1) CN110133721B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113640882B (zh) * 2021-08-10 2024-06-21 南方海洋科学与工程广东省实验室(湛江) 采集脚印噪声去除方法、电子设备及计算机可读存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102998699A (zh) * 2011-09-08 2013-03-27 中国石油天然气集团公司 一种含有套管谐波的垂直地震剖面资料层速度反演方法
CN103244104A (zh) * 2013-05-13 2013-08-14 中国石油大学(华东) 利用频散校正提取套管波衰减率的方法
CN103713315A (zh) * 2012-09-28 2014-04-09 中国石油化工股份有限公司 一种地震各向异性参数全波形反演方法及装置
CN104133247A (zh) * 2014-07-31 2014-11-05 中国石油天然气集团公司 垂直地震剖面数据中套管波的压制方法及装置
CN105528510A (zh) * 2014-09-30 2016-04-27 中国石油化工股份有限公司 水力压裂微地震震源机制的反演方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015134054A1 (en) * 2014-03-05 2015-09-11 Carbo Ceramics Inc. Systems and methods for locating and imaging proppant in an induced fracture

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102998699A (zh) * 2011-09-08 2013-03-27 中国石油天然气集团公司 一种含有套管谐波的垂直地震剖面资料层速度反演方法
CN103713315A (zh) * 2012-09-28 2014-04-09 中国石油化工股份有限公司 一种地震各向异性参数全波形反演方法及装置
CN103244104A (zh) * 2013-05-13 2013-08-14 中国石油大学(华东) 利用频散校正提取套管波衰减率的方法
CN104133247A (zh) * 2014-07-31 2014-11-05 中国石油天然气集团公司 垂直地震剖面数据中套管波的压制方法及装置
CN105528510A (zh) * 2014-09-30 2016-04-27 中国石油化工股份有限公司 水力压裂微地震震源机制的反演方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
水力压裂微地震监测稳定共振频率信号的解释;于辉 等;《物探化探计算技术》;20170115;第39卷(第1期);第90-95页 *

Also Published As

Publication number Publication date
CN110133721A (zh) 2019-08-16

Similar Documents

Publication Publication Date Title
CN109425896B (zh) 白云岩油气储层分布预测方法及装置
CN101086535B (zh) 表征油藏随时间演变的方法和程序
US20180203144A1 (en) Interferometric Microseismic Imaging Methods and Apparatus
KR20200014387A (ko) 지하 구조물의 검출
CN110954943B (zh) 被动源地震频率谐振勘探方法
CN110348135B (zh) 一种随钻声波测井评价地层渗透率的方法
WO2005096023A1 (en) Characterizing properties of a geological formation by coupled acoustic and electromagnetic measurements
CN116378648B (zh) 一种基于随钻声波前视的近钻头地层探测方法及装置
NO318894B1 (no) Fremgangsmate for a estimerer den hydrauliske konduktiviteten til en petrofysisk diskontinuitet i sideveggen til et borehull
CN106526678A (zh) 一种反射声波测井的波场分离方法及装置
Walter et al. Mapping Rainfall-Triggered Slidequakes and Seismic Landslide-Volume Estimation at Heumoes SlopeAll rights reserved. No part of this periodical may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or any information storage and retrieval system, without permission in writing from the publisher
CN102998703B (zh) 基于地表一致性反褶积进行储层预测的方法及设备
CN104181586A (zh) 一种水陆检波器数据海底反射系数反演的方法
CN109541689B (zh) 一种基于反射波能量特征的介质密实度评价方法
CN110133721B (zh) 一种水力压裂过程监测方法及系统
Delay et al. Is subsurface geophysics as seismic and acoustic investigations a rescue to groundwater flow inversion?
Tabjula et al. Empirical correlations for predicting flow rates using distributed acoustic sensor measurements, validated with wellbore and flow loop data sets
JP5517258B2 (ja) 鉛直アレイ地震計を利用したq値測定方法
CN109324344A (zh) 基于纯纵波和拟声波反演的页岩厚度预测方法及系统
CN113900145B (zh) 一种基于微动信号的地下横波速度分布的推测方法和装置
Li et al. Acoustic parameters inversion and sediment properties in the Yellow River reservoir
CN105572742B (zh) 一种确定海水深度的方法和装置
CN112241025B (zh) 一种井震联合地层压力确定方法及系统
CN103424774B (zh) 一种地震油气藏成像方法及装置
CN113296149A (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
CB03 Change of inventor or designer information

Inventor after: Liang Yao

Inventor after: Shu Mengcheng

Inventor after: Zhou Xuhui

Inventor after: Huo Shoudong

Inventor after: Shu Guoxu

Inventor before: Shu Mengcheng

Inventor before: Huo Shoudong

Inventor before: Zhou Xuhui

Inventor before: Liang Yao

Inventor before: Shu Guoxu

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant