CN101553742B - 使用图像和声波测井图的组合区分天然裂缝导致的声波各向异性和应力导致的声波各向异性 - Google Patents
使用图像和声波测井图的组合区分天然裂缝导致的声波各向异性和应力导致的声波各向异性 Download PDFInfo
- Publication number
- CN101553742B CN101553742B CN200780033949.XA CN200780033949A CN101553742B CN 101553742 B CN101553742 B CN 101553742B CN 200780033949 A CN200780033949 A CN 200780033949A CN 101553742 B CN101553742 B CN 101553742B
- Authority
- CN
- China
- Prior art keywords
- modulus
- dipole
- transverse wave
- slowness
- shear
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
使用图像和声波测井图的组合区分裂缝导致的声波各向异性和应力导致的声波各向异性。通过已知的技术得到井眼图像和声波测井图。在各向异性地层中的来自单极P波和S波的声波数据和交叉偶极横波声波数据的分析用于估算至少一个压缩模量和两个剪切模量以及偶极快横波方向。图像测井图的裂缝分析能够确定裂缝的类型和几何特性。来自图像测井图的地质学分析和地质力学分析提供天然裂缝和应力导致的裂缝的先前区分。基于由图像测井图解释的裂缝特性的知识的天然裂缝导致的声波各向异性和应力导致的声波各向异性的正演定量模型允许计算快横波方位和快慢横波慢度之间的差。然后优化预测的与观测的声波测量值(即,快横波方位和快横波慢度)之间的不吻合值,以区分具有被张开的天然裂缝、闭合的天然裂缝和由不相等的主应力效应导致的裂缝所影响的弹性介质的深度带。
Description
技术领域
本发明总体涉及井眼地质学和地球物理学领域,本发明更具体地涉及具有弹性介质的深度带的确定和估算,该弹性介质通过天然裂缝或者应力导致效应而被改变。
背景技术
通过将井眼钻到含油气层、完井、然后采收产物而将石油和天然气从地下岩层提取。使用各种传感器来提高井眼的建造和已完井的产率。例如,使用电缆和随钻声波测井仪以利用纵波和横波速度测量值测量井眼周围的地层的动态弹性特性。当地层的弹性特性是各向异性时,根据传播方向和偏振方向,可以测量若干个速度,并使用该速度部分或全部地表征各向异性弹性张量。各种条件可以引起各向异性,包括但不限于岩石固有特性、裂缝、和不相等的主应力。后者条件对于井眼的稳定性、最佳压裂、完井设计和其它地球和石油物理应用来说有一些指示意义。因此正确的识别各向异性的原因是重要的。
用于识别引起各向异性的原因的一些技术是公知的。在各向异性地层中的单极P波和S波,单极斯通利波和交叉偶极横波声波数据可以用来评估一个压缩模量和三个剪切模量[2005年75th Ann.Internat.Mtg.Soc.of Expl.Geophys.中Sinha,B.等人的Radial profiling of three formation shear moduli;2004年3月30日授予Sinha,B等人的题目为”Determination of anisotropicmoduli ofearth formations”的美国专利6,714,480,这些文献所公开的全部内容通过引用在此全文并入]。具有垂直对称轴线的正交地层由三个剪切模量c44、c55、和c66表征。在垂直井眼中,能够通过交叉偶极波形的方位各向异性分析直接评估两个垂直剪切模量(c44和c55)。快横波方位能够使用诸如Alford旋转的方法来计算,而快横波慢度(slowness:或者时差)和慢横波慢度可以由交叉偶极频散的零频率极限来估算[1986年56th Ann.Internat.Mtg.,Soc.of Expl.Geophys.中Alford,R.M.的Shear data in the presence ofazimuthal anisotropy;1994年64th Ann.Internat.Mtg,Soc.of Expl.Geophys.中Esmersoy,C.等人的Dipole shear anisotropy logging;2004年75th Ann.Internat.Mtg.Soc.of Expl.Geophys.中Sinha,B.等人的Radial profiling of three formation shear moduli;1993年5月25日授予Esmersoy,C.的题目为“Method and Apparatus for Determining Properties of Anisotropic ElasticMedia”的美国专利5,214,613;1998年9月15日授予Esmersoy,C.的题目为“Dipole Shear Anisotropy Logging”的美国专利5,808,963,或者对于可选的方法参见2004年4月6日授予Sinha,B.等人的题目为“Determination of dipoleshear anisotropy of earth formations”的美国专利申请6,718,266;1999年Geophysics中TangX等人的Simultaneous inversion of formation shear-wave anisotropy parameters from cross-dipole acoustic-array waveform data,这些文献所公开的全部内容通过引用在此全文并入]。假设施加校正以消除任何近井眼变化和测井仪影响,第三剪切模量c66能够由斯通利波数据来估算[Geophysics,1993,58,1091-1098中Norris,AN.等人的Weak elastic anisotropy and the tube wave;2004年3月30日授予Sinha,B.的题目为“Determination of anisotropic moduli of earth formations”的美国专利6,714,480,这些文献所公开的全部内容通过引用在此全文并入]。然后使用偶极频散曲线识别弹性特性的各向异性的原因:(i)使用偶极曲线的特征交叉识别应力导致的效应(由于远场不相等的主应力和井眼周围近场压力集中)[1996年的Geophysics中Sinha,B.K.等人的Stress-induced azimuthal anisotropy in borehole flexural waves;1998年的Geophysics中Winkler,K.W.等人的Effects of borehole stress concentrations on dipole anisotropy measurements;2000年的Geophysics中Sinha,B.K.等人的Dipole dispersion crossover and sonic logs in a limestone reservoir;1995年3月14日授予Sinha,B.的题目为“Identification of Stress Induced Anisotropy in Formations”的美国专利申请5,398,215,这些文献所公开的全部内容通过引用在此全文并入]。或者(ii)使用平行频散曲线的特征识别固有各向异性或者裂缝导致的各向异性[1994年的Geophysics中Sinha,B.K.等人的Borehole flexural modes in anisotropic formations;1995年3月14日授予Sinha,B.的题目为“Identificationof Stress Induced Anisotropy in Formations”的美国专利申请5,398,215,这些文献所公开的全部内容通过引用在此全文并入]。然而,当裂缝和应力效应都存在时,或者由于高频率的衰减难于解释频散曲线的分析时[2006年47thAnnu.Logging Symp.,SPWLA中Donald,A.等人的Advancements in acoustic techniques for evaluating natural fractures,此文献所公开的全部内容通过引用在此全文并入],或者当各向异性介质的对称轴线和井眼轴线没有对齐时,所观测的各向异性的解释变得更加具有挑战性。必须提供独立的信息以确认观测和区分不同效应的相对重要性。
当主应力的方向和天然裂缝平面的法向不对齐时,区分不同效应的相对重要性尤其重要。对斯通利波波型反射和衰减的分析允许识别井眼中的张开裂缝,并且估算所述张开裂缝的孔径[1989年9月26日授予Hsu,K.的题目为“Method and apparatus for detecting and evaluating borehole wall”的美国专利4,870,627;1989年的Geophysics中Hornby,B.E等人的Fracture evaluation using reflected Stoneley-wave arrivals;1997年Geophysics中Tezuka,K.等人的Modeling of low-frequency Stoneley-wave propagation in an irregular borehole;1989年5月16日授予Hornby,B.的题目为“BoreholeLogging Method for Fracture Detection and Evaluation”的美国专利4,831,600,这些文献公开的内容通过引用在此全文并入]。此外,井眼图像(电图像和超声图像)的解释可以用于识别张开裂缝或闭合裂缝[2000年的Springer中Luthi,S.M.的Geological well logs:their use in reservoir modeling;1993年9月7日授予Luthi,S.的题目为“Width determination of fracturesintersecting a borehole”的美国专利5,243,521,这些文献公开的内容通过引用在此全文并入]。然后可以计算诸如位置和走向的裂缝特性。然而,不存在实际技术用于正演定量模拟天然裂缝导致的声波各向异性和应力导致的声波各向异性,从而区别它们的相对效应。
发明内容
根据本发明的一个实施例,用于区分天然裂缝导致的声波各向异性和应力导致的声波各向异性的方法包括步骤以下:获得与井眼相关的图像数据和声波数据;使用声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向;使用图像数据识别裂缝类型和走向;计算快横波方位和由于地质断裂引起的快慢横波之间的慢度的差;以及识别具有正被张开的天然裂缝、闭合的天然裂缝和不相等的主应力效应中的至少一个影响时的弹性介质的深度带。
根据本发明的另一方面,用于区分天然裂缝导致的声波各向异性和应力导致的声波各向异性的设备包括:存储器,所述存储器能够操作以存储与井眼相关的图像数据和声波数据;能够操作以使用声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的电路;能够操作以使用数据图像识别裂缝类型和走向的电路;能够操作以计算快横波方位和由地质断裂导致的快慢横波之间的慢度的差的电路;和能够操作以识别具有正被张开的天然裂缝、闭合的天然裂缝和不相等主应力效应中的至少一个所影响时的弹性介质的深度带的电路。
附图说明
图1示出了一种准备用于识别和估算各向异性的声波处理;和
图2示出了具有弹性介质的深度带的识别和评估,所述弹性介质通过天然裂缝、不相等的主应力效应及其任意组合而被改变。
具体实施方式
参照附图1和2,在初始步骤(100)中,利用现有波型(包括但不限于单极波型(mode)、偶极波型和斯通利波波型)中的任一个及其组合获得声波测井图。与各向异性地层相关联的单极P波和S波,单极斯通利波和交叉偶极横波声波数据用于估算一个压缩模量和三个剪切模量。具有垂直对称轴线的正交地层由三个剪切模量c44、c55和c66表征。
在步骤(100)中获取的声波测井图被处理以准备用于识别和估算声波各向异性慢度特性。如步骤(101)中所示的,单极波型纵波慢度被识别。来自单极波型的纵波慢度用作步骤(221)中的输入。如步骤(102)中所示的,由偶极波型数据识别横波各向异性。在垂直井眼中,两个垂直偶极剪切模量(c44和c55)由交叉偶极波形的方位各向异性分析直接估算。快横波方向(或者方位)也被确定。用于确定快横波方向的技术包括但不限于Alford旋转和交叉偶极波形的参数反演。最高剪切模量(c44或c55)用作步骤(221)中的输入。快横波方位和快慢慢度用作步骤(224)中的输入,用于计算预测的声波测量值与观测的声波测量值之间的不吻合值。如步骤(103)中所示,偶极频散(dipole dispersion)曲线可以用来识别各向异性的原因:(i)使用偶极曲线的特征交叉识别应力引起的效应(由于远场不相等的主应力和井眼周围近场压力集中),或者(ii)使用平行频散曲线识别固有各向异性或者由裂缝导致的各向异性。频散曲线用作步骤(300)中的质量控制指示器。如步骤(104)中所示,假设施加校正以消除近井眼变化和测井仪的影响,第三剪切模量c66由垂直井眼中的斯通利波数据估算。当比较c44、c55和c66时,如果来自斯通利波的剪切模量是最高剪切模量,则所述剪切模量用作步骤(221)中的输入。对于具有垂直或者非垂直对称轴线的正交地层中的非垂直井来说,结合步骤(101、102和104)。使用用于沿井眼轴线的Sh波和Sv波速度的综合表达式估算三个弹性模量,而由用于q-P波速的表达式估算弹性模量的第四组合。然后,来自单极的测量的纵波慢度(DTcomeas),和来自偶极的横波慢度(DTsfast meas)或者来自斯通利波的横波慢度(DTsST)中的最低横波慢度用作步骤(221)中的输入。来自斯通利波波型的裂缝分析在步骤(105)中执行。斯通利波波型反应的分析允许识别井眼中的张开裂缝并估算所述张开裂缝的孔径。来自斯通利波波型的裂缝位置的识别用作步骤(300)中的质量控制指示器QC。
现在参照图2,根据在步骤(200)中得到的图像测井图在步骤(201)中执行裂缝分析。井眼图像的解释通常包括层理分析(bedding analysis)、构造分析、不均匀性分析以及裂缝分析。对于电测井仪来说,裂缝的识别基于与主岩相比的电阻率的观测(或者导电裂缝或者具有电阻性的裂缝),而对于超声波测井仪来说,裂缝的识别基于振幅衬度以及传播时间的结合。井眼图像的解释包括区分不同裂缝类型(例如,天然裂缝和由不相等的主应力和井眼周围的应力集中产生的裂缝),并确定裂缝的地质学特性(即,给定深度的位置和走向、倾角和倾斜方位特性)。所述解释还可以包括裂缝孔径的确定[1990年Geophysics中Luthi,S.M.等人的Fracture apertures from electrical borehole scans,该文献所公开的全部内容通过引用在此全文并入],裂缝迹线长度连通性(fracture trace length connectivity)、裂缝孔隙度、储层裂缝间隔和裂缝组的裂缝密度。
地质学分析和地质力学分析允许区分天然裂缝和由不相等的主应力以及井眼周围应力集中导致的裂缝。应力导致的裂缝是在现有平面处的受拉钻进和剪切钻进导致的裂缝、劈开(breakout)和滑移(slip)。受拉钻进导致的裂缝在井眼图像上容易地探测到,因为它们相对于井壁对称,通常接近垂直或者沿垂直于最小主应力的平面产生。剪切钻进产生的裂缝在垂直于最小主应力的方向上被探测到,但是在井壁的相对侧不对称地扩大,且断裂角度(failure angle)小于直角,并与岩性单位的机械强度有关。当一个或者优选地两个无变形弧(即,所述弧的半径接近钻头的半径)被两个受损区带近似180°地间隔开的情况下,产生劈开。在井壁处的切向应力最大的情况下劈开。由于垂直井通常平行于垂直主应力,因此劈开简单地说明最小主应力的方位。倾斜井倾斜于影响最大切应力的三个主应力。在钻井期间或者钻井之后可以重新激活岩石中的滑移面。这种平面通常相对于当前的应力轴线以倾斜角被定位。在图示的技术中,当考虑张开的天然裂缝、闭合的天然裂缝、和应力导致的裂缝时,裂缝的类型(或者裂缝组)的数目可以为三个但不限于三个。步骤(201)中的裂缝分析使得确定了裂缝类型和地质特性(例如,倾斜角度和倾斜方位),所述裂缝类型和地质特性用作用于步骤(221)中的计算的输入。
当存在具有各种走向的若个平坦裂缝并且可以在图像测井图上识别所述裂缝时,张开的裂缝将在弹性介质中产生过柔量(excess ofcompliance),并因此增加声波慢度,即,降低速度[参见示例,1995年的Geophysics中Schoenberg,M.,等人的Seismic anisotropy of fractured rock,该文献所公开的全部通过引用在此全文并入]。然而,愈合的裂缝可以减少介质(例如充满方解石和石英的介质)的整个柔量。假设存在一个或若干个平行裂缝组,长波长的有效介质可以看作为具有不同可能对称性(例如,横向各向同性的、正交的、单斜晶系、甚至三斜晶系)的各向异性弹性岩石。然而,Grechka和Kachanov[2006年的Geophysics中Grechka,V.等人的Seismic characterization of multiple fracture sets:does orthotropy suffice?,该文献所公开的全部内容通过引用在此全文并入]示出不管裂缝组的数量、所述裂缝组裂纹密度和走向,正交对称性通常对干裂缝是足够的。如果裂缝填充有可压缩流体或者如果考虑裂缝表面的粗糙度,则可以得到诸如单斜晶系的和三斜晶系的较低的对称性。基于裂缝类型、结构及其它裂缝特性的解释,可以分析将对弹性介质有影响的所观测的裂缝。由图像测井图和斯通利波方法观测的裂缝张开度典型地大约为10μm-2mm。假设声波频率为1-5kHz并且慢度为100-800μs/ft,典型的声波波长大约为0.25-10ft(0.08-3.05m)。因此,低频声波波长总是大于裂缝开口度。对于裂缝组来说,另一重要参数是相邻裂缝之间的距离,所述相邻裂缝之间的距离由裂缝密度参数(裂缝数量与给定距离的比)来描述。对于非常稀疏的裂缝组来说,典型的裂缝密度值小于0.75(单位是1/米),而对于密集裂缝组来说,所述裂缝密度值可以超过10(1/米)。裂缝间距(裂缝密度的倒数)总是大于裂缝张开度,但是裂缝间距可以小于、等于或者大于声波波长。当裂缝间距比波长小的多时,可以应用有效介质模型,并且该模型的参数仅仅取决于裂缝组的走向和密度、填充裂缝的材料的特性以及主岩石的弹性系数(例如,来自1980年的J.Eng.Mech.Div.Am.Soc.Civ.Eng.中Kachanov,M.的Continuum model of medium with cracks;1980年Journal of the AcousticalSociety of America中Schoenberg,M.的Elastic wave behavior across linear slip interfaces;1980年Mathematical Proceedings of the Cambridge PhilosophicalSociety中Hudson,J.A的Overall properties of a cracked solid,这些文献公开的内容通过引用在此全文并入]的用于非相互作用裂纹的有效介质理论)。当裂缝间距接近或者大于声波波长时,由有效介质理论获得的数据和模拟的定量比较不是严格正确的。然而,当可以忽略第一级散射效应时,仍然可以以定性的方式使用该模型。
Sayers和Kachanov[1995年的J.Geophys.Res.中Sayers.C.M.等人的Microcrack-induced elastic wave anisotropy of brittle rocks和1995年的Geophysics中Schoenberg和Sayers Schoenberg,M.等人的Seismic anisotropy of fractured rock,这些文献所公开的内容通过引用在此全文并入]描述了一种用于包括在通过断裂岩石的地震传播中地质上的实际裂纹和裂缝的效应的简单位移不连续方法。他们考虑任意形状和有限尺寸的纤细裂纹或平坦和平行裂纹组,对于所述纤细裂纹或平坦和平行裂纹组来说,不同裂纹和裂缝是不相互作用的(每一个裂纹被假设为受到不被相邻裂纹干扰的均衡的应力场)。显示裂纹或裂缝相互作用和交叉的非相互作用近似值的数值确认对有效弹性几乎没有影响[2006年的Geophysics中Grechka,V.等人的Effective elasticity of rocks with closely spaced and intersecting cracks,该文献所公开的内容通过引用在此全文并入]。利用此假设,整个弹性柔量sijkl能够通过裂纹或裂缝组的柔量和基底介质(host medium)的柔量的线性相加来分解:
当裂纹或裂缝旋转不变(例如,裂纹或裂缝是平坦的),基底介质的过柔量可以写成:
其中δij是Kronecker符号,αij是二阶张量,而βijkl是四阶张量。对于嵌入在典型的体积V中的单个裂纹来说,Sayers和Kachanov(1995)将αij和βijkl定义为:
对于平坦和平行裂缝组来说,Schoenberg和Sayers(1995)将αij和βijkl定义为:
其中,Z(m) N和Z(m) T是mth裂缝的法向和切向柔量(单位:l/应力),n(m) i是裂缝法向ith分量。可以对特殊情况得到这两个方程之间的关系,例如,在立方体V=L×L×L中具有规则间距s=L/N和相等面积A=L×L的一组N个平行纤细裂纹(裂纹平行于立方体一个面)给出BTN/1=BT/s=ZT和(BN-BT)N/1=(BN-BT)/s=ZN-ZT。
根据由井眼图像解释得出的裂缝特性,如果满足长波长有效介质条件,则该模型可以应用于计算裂缝导致的声波各向异性。主岩石的裂缝柔量sh ijkl和裂缝Z(m) N和Z(m) T中的一个必须被指定(221)。在声波发射器和接收器(由测井仪的几何结构给出)之间探测到的弹性介质限定在其上满足长波长近似值的体积。基底介质柔量sh ijkl涉及无裂缝介质的弹性特性(即,声波慢度)。它们在井中是未知的,因为具有所有混合效应的介质由测量。然而,实际上好的起始点是将基底介质各向异性而考虑。然后,通过由单极测量的纵波慢度(DTco meas)、由偶极测量的最低横波慢度(DTsfast meas)以及密度在离散深度处定义两个各向同性的弹性常数。在正演模拟问题的情况下,法向Z(m) N和切向Z(m) T柔量是未知的,并且对于每个裂缝来说可以不同。目前,在图像测井图上观测的每个裂缝(m)被假定为具有带有柔量[Z(m) N](n)和[Z(m) T](n)的给定类型的裂缝组(n)的一部分。在开始时,每个裂缝根据地质解释属于给定组。在工作流程(220)的后期阶段,每个裂缝可以被重新指定为不同的裂缝组。对于干裂缝或者填充气体的裂缝来说,Z(m) N/Z(m) T~1是良好的临近值,但是通常,对于填充水或填充油的裂缝来说,Z(m) N/Z(m) T<<l(Sayers,2002)。当法向柔量和切向柔量相等时(Z(m) N/Z(m) T=1),βijkl为零,并且弹性张量的最低的可能对称性是斜方晶的。通过对于石英二长石和石灰石的实验室测量,裂缝法向和剪切柔量的估算值在0.01.le-12到0.5.le-12m/pa的范围内变化;[1995年Int.J.Rock Mech.Min.中Pyrak-Nolte,L.J.的The seismic response of fracture and the interrelationships among fracture properties:;2005年Dphil Thesis中LubbeR.的A field and laboratory investigation of the compliance of fractured rocks,这些文献公开的内容通过引用在此全文并入],和通过在主频率为2kHZ下的井间地震实验,所述值从0.25.le-12到3.5.le-12m/pa的范围内变化[2006年Geophys.Prosp.中Lubbe,R.等人的A field investigation of fracture compliance,该文献公开的内容通过引用在此全文并入]。法向和切向柔量Z(m) N和Z(m) T与裂缝的物理特性的关系不在本专利的保护范围内;更多细节参见Pyrak-Nolte(1995)。在本发明的方法中,我们通过不同的裂缝组区分不同的裂缝类型。例如,使用井眼图像解释,三组裂缝可以被定义为张开的天然裂缝、闭合的天然裂缝和应力导致的裂缝。因为裂缝弹性柔量sh ijkl(方程5和6)是各个裂缝贡献的线性叠加的结果,因此总裂缝柔量能够例如以张开的天然裂缝分量sof ijkl、闭合的天然裂缝导致的分量scf ijkl和应力导致的分量sif ijkl的形式被分解成给定数量组的线性叠加:
过张量柔量sof ijkl、scf ijkl、sif ijkl被各自定义如下:
且使用具有相应法向和切向分量的[Z(m) N](of)、[Z(m) T](of)、[Z(m) N](cf)、[Z(m) T](cf)、[Z(m) N](if)和[Z(m) T](if)的方程(5)和方程(6)定义二阶张量αof ij、acf ij和aif ij和四阶张量βof ij、βcf ij和βif ij。在步骤221中,根据定义的组的数量,法向和切向柔量作为[Z(m) N](n)、[Z(m) T](n)被指定给每个组(n)。
在步骤(222)中使用方程(7)来计算由于所有裂缝组的过柔量。
在步骤(223)中执行快慢横波慢度以及方位的预测。一旦已经定义上述参数,则可以在离散深度点处执行计算。在选定的离散深度处,存在于所选择的深度上下方的深度窗(depth window)内的裂缝包括在计算中,并且裂缝深度、倾斜角度和倾斜方位的知识允许为每个裂缝确定n(r) i,并计算所有弹性柔度张量。可以对柔度张量反演以得到刚度张量。要注意的是在此方法中,对弹性介质的任何具体的对称性不做假设。然后使用刚度张量通过Christoffel方程求解用于任意各向异性的运动。假设传播方向沿着井眼轴线,我们对弹性波传播的三种波型的偏振矢量求解一个压缩qP波和两个剪切qS1波和qS2波。在正交于井眼的横波速度的方位变化的分析提供通常在声波各向异性中所观测的两个特性:快横波方位和两个快慢横波慢度。
在步骤(224)中使用步骤(223)的输出,即,预测的快横波方位和快慢横波慢度。
在步骤(224)中,观测值与预测值之间的不吻合值用于在由张开的天然裂缝影响的区带、闭合的天然裂缝影响的区带和应力导致的裂缝影响的区带之间进行区分。在深度z处,观测的快横波方位和慢度差与预测的快横波方位和慢度差之间的预测值误差或不吻合值可以表示为:
然后,下面的不吻合矢量的标准Ln在被定义为:
或者
其中,σi是方程(11)和方程(12)中使用的观测值的标准偏差。
使用方程(11)和方程(12),定义下面的标准:
(方程15)
然后,通过将各个裂缝重新指定在不同类型中,并对每一种裂缝类型重新定义柔量特性,在重复的计算步骤221、222、223、224中执行优化环路(步骤220)。以上方程(分别是独立的方程15、方程16、方程17、方程18,或者用于联合反演(joint inversion)的方程19)可以最小化为反演问题,以对于给定深度间隔为每一个裂缝得到最优类型和柔度参数[Z(m) N](of)、[Z(m) T](of)、[Z(m) N](cf)、[Z(m) T](cf)、[Z(m) N](if)和[Z(m) T](if)。还可以定义最小化方程15、16、17、18和19的反演波型(至少平方或不平方,即,n=2或者n≠2)。当方程15、16、17、18及19中的数量在定义的临界值以下时,停止优化环路。224的输出是已经重新指定给三种裂缝类型(张开的天然裂缝、闭合的天然裂缝和由应力产生的裂缝)中的一个的裂缝的集合,并且对于所述裂缝来说,柔度参数[Z(m)N](of)、[Z(m) T](of)、[Z(m) N](cf)、[Z(m) T](cf)、[Z(m) N](if)和[Z(m) T](if)现在是已知的。
在步骤300中,在深度窗中的离散深度处,对三种类型中的每一种计算裂缝的比例,并显示为颜色表示/阴影的曲线,所述曲线表示被张开的天然裂缝、闭合的天然裂缝和应力产生的裂缝所影响的区带。
可以将使用偶极频散曲线(步骤104)识别的固有各向异性和应力导致的各向异性的区带的结果与使用斯通利波波型(步骤105)已经识别的裂缝的区带的结果进行比较。
虽然已经通过以上示例性实施例说明了本发明,但是本领域的技术人员将理解的是在不背离这里公开的发明构思的情况下可以对图示的实施例进行修改和改变。此外,虽然结合各种示例性结构说明了优选的实施例,但是本领域的技术人员将认识到可以使用各种具体结构实施该系统。因此,除了被所附权利要求的保护范围和精神限制之外,本发明不应该被看作为限制性的。
Claims (14)
1.一种区分天然裂缝导致的声波各向异性和应力导致的声波各向异性的方法,所述方法包括步骤以下:
获得与井眼相关的图像数据和声波数据;
使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向,其中所述使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的步骤确定偶极快横波方向、和快横波慢度与慢横波慢度之间的差;
使用所述图像数据识别裂缝类型和走向;
根据由所述图像数据解释的裂缝特性的知识,通过天然裂缝导致的声波各向异性和应力导致的声波各向异性的正演定量模型计算所述偶极快横波方向和由于地质断裂引起的快横波慢度与慢横波慢度之间的差,其中所述根据由所述图像数据解释的裂缝特性的知识通过天然裂缝导致的声波各向异性和应力导致的声波各向异性的正演定量模型计算所述偶极快横波方向和由于地质断裂引起的快横波慢度与慢横波慢度之间的差的步骤预测偶极快横波方向、和快横波慢度与慢横波慢度之间的差;
通过计算由所述使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的步骤确定的所述偶极快横波方向和快横波慢度与慢横波慢度之间的差、与、由所述根据由所述图像数据解释的裂缝特性的知识通过天然裂缝导致的声波各向异性和应力导致的声波各向异性的正演定量模型计算所述偶极快横波方向和由于地质断裂引起的快横波慢度与慢横波慢度之间的差的步骤预测的所述偶极快横波方向和快横波慢度与慢横波慢度之间的差、之间的不吻合值计算预测值与观测值之间的不吻合值;以及
识别正被张开的天然裂缝、闭合的天然裂缝和应力导致的效应中的至少一个影响的深度带。
2.如权利要求1所述的方法,其中,所述正演定量模型基于力学模型。
3.如权利要求1所述的方法,其中,所述使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的步骤使用单极P波估算纵波慢度。
4.如权利要求3所述的方法,其中,所述使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的步骤使用单极S波估算所述所述快横波慢度和所述慢横波慢度。
5.如权利要求4所述的方法,其中,所述使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的步骤使用交叉偶极横波声波数据估算两个横波慢度以及偶极快横波方向。
6.如权利要求5所述的方法,其中,所述使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的步骤使用单极斯通利波估算所述所述快横波慢度和所述慢横波慢度。
7.如权利要求1所述的方法,还包括以下步骤:
使用所述图像数据以通过地质学分析和地质力学分析识别裂缝特性,所述地质学分析和地质力学分析使得在张开的天然裂缝、闭合的天然裂缝和应力导致的效应之间进行先验区分。
8.一种区分天然裂缝导致的声波各向异性和应力导致的声波各向异性的设备,包括:
存储器,所述存储器能够操作以存储与井眼相关的图像数据和声波数据;
能够操作以使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的电路,其中所述能够操作以使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的电路确定偶极快横波方向、和快横波慢度与慢横波慢度之间的差;
能够操作以使用所述图像数据识别裂缝类型和走向的电路;
能够操作以根据由所述图像数据解释的裂缝特性的知识通过天然裂缝导致的声波各向异性和应力导致的声波各向异性的正演定量模型计算偶极快横波方向和由地质断裂导致的快横波慢度与慢横波慢度之间的差的电路,所述能够操作以根据由所述图像数据解释的裂缝特性的知识通过天然裂缝导致的声波各向异性和应力导致的声波各向异性的正演定量模型计算偶极快横波方向和由地质断裂导致的快横波慢度与慢横波慢度之间的差的电路预测偶极快横波方向、和快横波慢度与慢横波慢度之间的差;
能够操作以通过计算由所述能够操作以使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的电路确定的偶极快横波方向和快横波慢度与慢横波慢度之间的差与由能够操作以根据由所述图像数据解释的裂缝特性的知识通过天然裂缝导致的声波各向异性和应力导致的声波各向异性的正演定量模型计算偶极快横波方向和由地质断裂导致的快横波慢度与慢横波慢度之间的差的电路预测的偶极快横波方向和快横波慢度与慢横波慢度之间的差之间的不吻合值计算预测值与观测值之间的不吻合值的电路;以及
能够操作以识别正被张开的天然裂缝、闭合的天然裂缝和应力导致的效应中的至少一个所影响的深度带的电路。
9.如权利要求8所述的设备,其中,所述正演定量模型基于力学模型。
10.如权利要求8所述的设备,其中,所述能够操作以使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的电路使用单极P波估算纵波慢度。
11.如权利要求10所述的设备,其中,所述能够操作以使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的电路使用单极S波估算所述快横波慢度和所述慢横波慢度。
12.如权利要求11所述的设备,其中,所述能够操作以使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的电路使用交叉偶极横波声波数据估算两个横波慢度以及偶极快横波方向。
13.如权利要求12所述的设备,其中,所述能够操作以使用所述声波数据估算至少一个压缩模量和两个剪切模量以及偶极快横波方向的电路使用单极斯通利波估算所述快横波慢度和所述慢横波慢度。
14.如权利要求8所述的设备,还包括:
能够操作以使用所述图像数据通过地质学分析和地质力学分析识别裂缝特性的电路,所述地质学分析和地质力学分析使得在张开的天然裂缝、闭合的天然裂缝和应力导致的效应之间进行先验区分。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US82529806P | 2006-09-12 | 2006-09-12 | |
US60/825,298 | 2006-09-12 | ||
PCT/US2007/078041 WO2008033770A1 (en) | 2006-09-12 | 2007-09-10 | Discriminating natural fracture- and stress-induced sonic anisotropy using a combination of image and sonic logs |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101553742A CN101553742A (zh) | 2009-10-07 |
CN101553742B true CN101553742B (zh) | 2014-02-12 |
Family
ID=38886775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200780033949.XA Expired - Fee Related CN101553742B (zh) | 2006-09-12 | 2007-09-10 | 使用图像和声波测井图的组合区分天然裂缝导致的声波各向异性和应力导致的声波各向异性 |
Country Status (4)
Country | Link |
---|---|
US (1) | US7457194B2 (zh) |
CN (1) | CN101553742B (zh) |
RU (1) | RU2402045C1 (zh) |
WO (1) | WO2008033770A1 (zh) |
Families Citing this family (51)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7529150B2 (en) * | 2006-02-06 | 2009-05-05 | Precision Energy Services, Ltd. | Borehole apparatus and methods for simultaneous multimode excitation and reception to determine elastic wave velocities, elastic modulii, degree of anisotropy and elastic symmetry configurations |
US8190369B2 (en) * | 2006-09-28 | 2012-05-29 | Baker Hughes Incorporated | System and method for stress field based wellbore steering |
US7777606B2 (en) * | 2007-01-09 | 2010-08-17 | Westerngeco L.L.C. | Fracture cluster mapping |
US7647183B2 (en) * | 2007-08-14 | 2010-01-12 | Schlumberger Technology Corporation | Method for monitoring seismic events |
US7649804B2 (en) * | 2007-10-25 | 2010-01-19 | Schlumberger Technology Corporation | In-situ determination of yield stress state of earth formations |
US9477002B2 (en) * | 2007-12-21 | 2016-10-25 | Schlumberger Technology Corporation | Microhydraulic fracturing with downhole acoustic measurement |
FR2928959B1 (fr) * | 2008-03-21 | 2010-03-12 | Inst Francais Du Petrole | Procede d'estimation de la densite de fractures dans un milieu rocheux |
US8478530B2 (en) * | 2008-07-07 | 2013-07-02 | Baker Hughes Incorporated | Using multicomponent induction data to identify drilling induced fractures while drilling |
US8908473B2 (en) * | 2008-12-23 | 2014-12-09 | Schlumberger Technology Corporation | Method of subsurface imaging using microseismic data |
US8117014B2 (en) * | 2009-03-27 | 2012-02-14 | Schlumberger Technology Corporation | Methods to estimate subsurface deviatoric stress characteristics from borehole sonic log anisotropy directions and image log failure directions |
CN102236101B (zh) * | 2010-05-07 | 2013-03-13 | 中国石油天然气股份有限公司 | 多层裂缝预测方法和装置 |
CN101907725B (zh) * | 2010-06-12 | 2011-12-28 | 中国石油天然气股份有限公司 | 裂缝预测方法和装置 |
US9157318B2 (en) | 2011-01-04 | 2015-10-13 | Schlumberger Technology Corporation | Determining differential stress based on formation curvature and mechanical units using borehole logs |
CN102253415B (zh) * | 2011-04-19 | 2013-03-20 | 中国石油大学(华东) | 基于裂缝等效介质模型的地震响应模式建立方法 |
CN102183788B (zh) * | 2011-05-13 | 2012-11-21 | 上海石油天然气有限公司 | 地层条件下横波测井曲线的合成方法 |
CN102879800B (zh) * | 2011-07-15 | 2015-04-29 | 中国石油天然气集团公司 | 一种横波分裂裂缝检测的方法 |
EP2737343A2 (en) | 2011-07-27 | 2014-06-04 | Schlumberger Technology B.V. | Multi-well anisotropy inversion |
CN102540253B (zh) * | 2011-12-30 | 2014-02-05 | 中国石油天然气股份有限公司 | 一种陡倾角地层及裂缝地震偏移方法和装置 |
US20140078288A1 (en) * | 2012-06-19 | 2014-03-20 | Schlumberger Technology Corporation | Far Field In Situ Maximum Horizontal Stress Direction Estimation Using Multi-Axial Induction And Borehole Image Data |
CN103852798A (zh) * | 2012-11-28 | 2014-06-11 | 中国石油集团长城钻探工程有限公司 | 井孔斯通利波的慢度测量方法 |
US20150309200A1 (en) * | 2012-11-30 | 2015-10-29 | Timur Vyacheslavowich Zharnikov | A method for processing acoustic waveforms |
US20150301213A1 (en) * | 2012-11-30 | 2015-10-22 | Schlumberger Technology Corporation | A method for processing acoustic waveforms |
AU2012396820A1 (en) | 2012-12-11 | 2015-05-14 | Halliburton Energy Services, Inc. | Method and system for direct slowness determination of dispersive waves in a wellbore environment |
US9581716B2 (en) * | 2013-01-21 | 2017-02-28 | Schlumberger Technology Corporation | Methods and apparatus for estimating borehole mud slownesses |
US20140241111A1 (en) * | 2013-02-28 | 2014-08-28 | Weatherford/Lamb, Inc. | Acoustic borehole imaging tool |
US9670770B2 (en) * | 2013-06-24 | 2017-06-06 | Baker Hughes Incorporated | Fracture evaluation through cased boreholes |
US9664816B2 (en) * | 2013-12-06 | 2017-05-30 | Halliburton Energy Services, Inc. | Fracture detection and characterization using resistivity images |
CA2944375C (en) * | 2014-03-31 | 2023-01-24 | John Adam Donald | Subsurface formation modeling with integrated stress profiles |
CN104977606A (zh) * | 2014-04-02 | 2015-10-14 | 中国石油化工股份有限公司 | 一种缝洞型储层地震数值模型建立方法 |
CN103993871B (zh) * | 2014-05-08 | 2017-02-15 | 中国石油化工股份有限公司 | 针对薄互层地层的测井资料标准化处理方法及装置 |
CN105467438B (zh) * | 2014-09-04 | 2017-11-24 | 中国石油化工股份有限公司 | 一种基于三模量的泥页岩地应力三维地震表征方法 |
CN104199097B (zh) * | 2014-09-17 | 2017-12-26 | 西南石油大学 | 用于挤压构造体系潜在断裂定量判识的新方法 |
US10392934B2 (en) | 2014-10-03 | 2019-08-27 | Schlumberger Technology Corporation | Method and apparatus for processing waveforms |
CN104612646B (zh) * | 2014-12-10 | 2018-03-16 | 中国石油大学(华东) | 棋盘格式可视化水力裂缝模拟装置 |
US10724365B2 (en) | 2015-05-19 | 2020-07-28 | Weatherford Technology Holdings, Llc | System and method for stress inversion via image logs and fracturing data |
US10168442B2 (en) * | 2015-05-19 | 2019-01-01 | Schlumberger Technology Corporation | Differential energy analysis for dipole acoustic measurement |
CN105068117B (zh) * | 2015-08-25 | 2016-06-08 | 中国矿业大学(北京) | 用于裂缝性介质的avo反演方法、装置和设备 |
CN106168678B (zh) * | 2016-07-15 | 2018-08-17 | 中国矿业大学(北京) | 一种裂隙中传播横波的分离方法及系统 |
CN109098707B (zh) * | 2017-06-21 | 2022-02-11 | 中国石油化工股份有限公司 | 砂砾岩体油藏的直井缝网压裂的适应性评价方法、砂砾岩体油藏的直井缝网压裂方法 |
CN110295894B (zh) * | 2018-03-22 | 2022-03-01 | 中国石油化工股份有限公司 | 一种建立水平井产能预测模型的方法 |
CN108983312B (zh) * | 2018-05-25 | 2019-09-24 | 中国海洋石油集团有限公司 | 一种裂缝型储层各向异性岩石模量估算方法 |
CN109040047B (zh) * | 2018-07-25 | 2021-03-30 | 中国石油天然气股份有限公司 | 示功图数据压缩、传输方法和装置、计算机设备和存储介质 |
US10977489B2 (en) * | 2018-11-07 | 2021-04-13 | International Business Machines Corporation | Identification of natural fractures in wellbore images using machine learning |
CN111273341B (zh) * | 2018-12-05 | 2022-05-10 | 中国石油天然气股份有限公司 | 根据裂缝空间分布的含裂缝储层岩石物理建模方法 |
CN109902422B (zh) * | 2019-03-08 | 2021-07-27 | 西南石油大学 | 一种井眼轨迹选取方法及装置 |
CN111830563B (zh) * | 2019-04-17 | 2023-06-30 | 中国石油天然气集团有限公司 | 纯横波多层裂缝裂缝方向及快慢波时差确定方法及装置 |
CN110174696B (zh) * | 2019-05-29 | 2021-04-30 | 长安大学 | 一种介质对称轴与观测坐标轴互换的地震波采集方法 |
CN111736218B (zh) * | 2020-05-29 | 2023-10-27 | 中国石油天然气集团有限公司 | 地层各向异性成因定量分析方法、设备及可读存储介质 |
CN112034516B (zh) * | 2020-08-17 | 2021-11-12 | 中国石油大学(华东) | 一种评价井周横波三维各向异性的方法 |
US11525935B1 (en) | 2021-08-31 | 2022-12-13 | Saudi Arabian Oil Company | Determining hydrogen sulfide (H2S) concentration and distribution in carbonate reservoirs using geomechanical properties |
US11921250B2 (en) | 2022-03-09 | 2024-03-05 | Saudi Arabian Oil Company | Geo-mechanical based determination of sweet spot intervals for hydraulic fracturing stimulation |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4831600A (en) * | 1986-12-31 | 1989-05-16 | Schlumberger Technology Corporation | Borehole logging method for fracture detection and evaluation |
US5960371A (en) * | 1997-09-04 | 1999-09-28 | Schlumberger Technology Corporation | Method of determining dips and azimuths of fractures from borehole images |
US6718266B1 (en) * | 2002-10-31 | 2004-04-06 | Schlumberger Technology Corporation | Determination of dipole shear anisotropy of earth formations |
CN1217203C (zh) * | 1999-01-15 | 2005-08-31 | 贝克休斯公司 | 通过野外测量标定确定应力-速度耦合系数的方法 |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4870627A (en) | 1984-12-26 | 1989-09-26 | Schlumberger Technology Corporation | Method and apparatus for detecting and evaluating borehole wall fractures |
US5243521A (en) | 1988-10-03 | 1993-09-07 | Schlumberger Technology Corporation | Width determination of fractures intersecting a borehole |
US5060204A (en) * | 1990-06-27 | 1991-10-22 | Chevron Research And Technology Company | Method of layer stripping to determine fault plane stress build-up |
US5214613A (en) | 1991-03-12 | 1993-05-25 | Schlumberger Technology Corporation | Method and apparatus for determining properties of anisotropicelastic media |
US5398215A (en) | 1993-11-19 | 1995-03-14 | Schlumberger Technology Corporation | Identification of stress induced anisotropy in formations |
US5475650A (en) | 1993-11-19 | 1995-12-12 | Schlumberger Technology Corporation | Measurement of nonlinear properties of formation using sonic borehole tool while changing pressure in borehole |
US5808963A (en) | 1997-01-29 | 1998-09-15 | Schlumberger Technology Corporation | Dipole shear anisotropy logging |
US6944094B1 (en) * | 1997-06-20 | 2005-09-13 | Bp Corporation North America Inc. | High resolution determination of seismic polar anisotropy |
GB2329043B (en) * | 1997-09-05 | 2000-04-26 | Geco As | Method of determining the response caused by model alterations in seismic simulations |
US5999486A (en) * | 1998-07-23 | 1999-12-07 | Colorado School Of Mines | Method for fracture detection using multicomponent seismic data |
US7062072B2 (en) * | 1999-12-22 | 2006-06-13 | Schlumberger Technology Corporation | Methods of producing images of underground formations surrounding a borehole |
US6925031B2 (en) * | 2001-12-13 | 2005-08-02 | Baker Hughes Incorporated | Method of using electrical and acoustic anisotropy measurements for fracture identification |
US6714873B2 (en) | 2001-12-17 | 2004-03-30 | Schlumberger Technology Corporation | System and method for estimating subsurface principal stresses from seismic reflection data |
USH2116H1 (en) * | 2001-12-31 | 2005-04-05 | Southwest Research Institute | Method for detecting, locating, and characterizing single and multiple fluid-filled fractures in fractured formations |
US6714480B2 (en) | 2002-03-06 | 2004-03-30 | Schlumberger Technology Corporation | Determination of anisotropic moduli of earth formations |
US6920082B2 (en) * | 2002-06-27 | 2005-07-19 | Baker Hughes Incorporated | Method and apparatus for determining earth formation shear-wave transverse isotropy from borehole stoneley-wave measurements |
US7035165B2 (en) * | 2003-01-29 | 2006-04-25 | Baker Hughes Incorporated | Imaging near-borehole structure using directional acoustic-wave measurement |
US7516015B2 (en) * | 2005-03-31 | 2009-04-07 | Schlumberger Technology Corporation | System and method for detection of near-wellbore alteration using acoustic data |
US7529150B2 (en) * | 2006-02-06 | 2009-05-05 | Precision Energy Services, Ltd. | Borehole apparatus and methods for simultaneous multimode excitation and reception to determine elastic wave velocities, elastic modulii, degree of anisotropy and elastic symmetry configurations |
US7486589B2 (en) * | 2006-02-09 | 2009-02-03 | Schlumberger Technology Corporation | Methods and apparatus for predicting the hydrocarbon production of a well location |
-
2006
- 2006-10-27 US US11/553,565 patent/US7457194B2/en not_active Expired - Fee Related
-
2007
- 2007-09-10 WO PCT/US2007/078041 patent/WO2008033770A1/en active Application Filing
- 2007-09-10 CN CN200780033949.XA patent/CN101553742B/zh not_active Expired - Fee Related
- 2007-09-10 RU RU2009113531/28A patent/RU2402045C1/ru not_active IP Right Cessation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4831600A (en) * | 1986-12-31 | 1989-05-16 | Schlumberger Technology Corporation | Borehole logging method for fracture detection and evaluation |
US5960371A (en) * | 1997-09-04 | 1999-09-28 | Schlumberger Technology Corporation | Method of determining dips and azimuths of fractures from borehole images |
CN1217203C (zh) * | 1999-01-15 | 2005-08-31 | 贝克休斯公司 | 通过野外测量标定确定应力-速度耦合系数的方法 |
US6718266B1 (en) * | 2002-10-31 | 2004-04-06 | Schlumberger Technology Corporation | Determination of dipole shear anisotropy of earth formations |
Also Published As
Publication number | Publication date |
---|---|
RU2402045C1 (ru) | 2010-10-20 |
CN101553742A (zh) | 2009-10-07 |
US20080062814A1 (en) | 2008-03-13 |
US7457194B2 (en) | 2008-11-25 |
WO2008033770A1 (en) | 2008-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101553742B (zh) | 使用图像和声波测井图的组合区分天然裂缝导致的声波各向异性和应力导致的声波各向异性 | |
CA2323676C (en) | Formation stress identification and estimation using borehole monopole and cross-dipole acoustic measurements | |
Sinha et al. | Dipole dispersion crossover and sonic logs in a limestone reservoir | |
CN109563736B (zh) | 用于表征地层的系统和方法 | |
US6714873B2 (en) | System and method for estimating subsurface principal stresses from seismic reflection data | |
US7872944B2 (en) | Method of monitoring microseismic events | |
US7652950B2 (en) | Radial profiling of formation mobility using horizontal and vertical shear slowness profiles | |
CA2540995C (en) | System and method for detection of near-wellbore alteration using acoustic data | |
US8379483B2 (en) | Method and apparatus for determining radial shear velocity variation from dipole acoustic logging | |
US8553493B2 (en) | Method for permeable zone detection | |
Harrison et al. | Acquisition and analysis of sonic waveforms from a borehole monopole and dipole source for the determination of compressional and shear speeds and their relation to rock mechanical properties and surface seismic data | |
WO2010111077A1 (en) | Methods to estimate subsurface deviatoric stress characteristics from borehole sonic log anisotropy directions and image log failure directions | |
MX2008001681A (es) | Metodo para caracterizador anisotropia de formacion de onda de esfuerzo cortante. | |
Prioul et al. | Forward modeling of fracture-induced sonic anisotropy using a combination of borehole image and sonic logs | |
US8681582B2 (en) | Method for sonic indication of formation porosity and lithology | |
EP0969294A2 (en) | Method of imaging the permeability and fluid content structure within sediment | |
Cheng et al. | Identifying and estimating formation stress from borehole monopole and cross-dipole acoustic measurements | |
US6061300A (en) | Method of imaging the permeability and fluid content structure within sediment | |
Pistre et al. | A new modular sonic tool provides complete acoustic formation characterization | |
Sinha et al. | Hydraulic fracture characterization using cased-hole sonic data | |
Mari et al. | Active and Passive acoustic logging applied to the detection of preferential flow in a sedimentary aquifer | |
Zheng | Dispersion Asymptotic Analysis and its Applications in Acoustic Logging | |
AU2003254730C1 (en) | Estimating formation stress using borehole monopole and cross-dipole acoustic measurements: theory and method | |
Escandón et al. | Detecting and characterizing fractures in sedimentary deposits with stoneley waves | |
Robinson et al. | LWD Sonic Real-Time Leaky Compressional and Quadrupole Measurements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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 |
Granted publication date: 20140212 Termination date: 20140910 |
|
EXPY | Termination of patent right or utility model |