CN108627871A - 一种tti介质裂缝性质参数的反演方法 - Google Patents
一种tti介质裂缝性质参数的反演方法 Download PDFInfo
- Publication number
- CN108627871A CN108627871A CN201710186472.4A CN201710186472A CN108627871A CN 108627871 A CN108627871 A CN 108627871A CN 201710186472 A CN201710186472 A CN 201710186472A CN 108627871 A CN108627871 A CN 108627871A
- Authority
- CN
- China
- Prior art keywords
- parameter
- character parameter
- cracks character
- cracks
- tti
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- 230000014509 gene expression Effects 0.000 claims abstract description 9
- 239000012530 fluid Substances 0.000 claims description 19
- 238000010606 normalization Methods 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000010008 shearing Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 206010017076 Fracture Diseases 0.000 description 56
- 208000010392 Bone Fractures Diseases 0.000 description 53
- 238000005516 engineering process Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000012267 brine Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000007598 dipping method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
公开了一种TTI介质裂缝性质参数的反演方法。该方法包括以下步骤：建立TTI介质PP波反射系数与纵横波模量、各向异性参数和密度的关系式；基于所述各向异性参数与密度归一化的TTI介质弹性参数矩阵的关系将所述反射系数与所述弹性参数矩阵相关联；将所述反射系数表示为频率以及裂缝性质参数的函数；基于所述反射系数构建反演裂缝性质参数的目标函数；通过对所述目标函数反演得到的裂缝性质参数进行归一化概率密度计算，将归一化概率密度最大值的位置所对应的裂缝性质参数的数值确定为反演的最终结果。该方法增加了预测的裂缝性质参数即裂缝旋转角度，可以更详细的描述裂缝发育情况。
Description
技术领域
本发明涉及油汽地球物理技术领域，更具体地，涉及一种TTI介质裂缝性质参数的反演方法。
背景技术
不同尺度的裂缝具有不同的描述方式：1)对于油藏数模中百米级大尺度裂缝的描述，主要基于构造应力场分析，结合几何型叠后地震属性(如相干、曲率等)进行；2)对于勘探阶段期望的几十米级中小尺度裂缝，主要依靠叠前叠后地震几何属性或物理规律进行。叠后地震数据的相干、曲率、蚂蚁追踪或其衍生技术可以识别出相对大尺度的断裂(大于1/8波长)，但对于更小尺度的裂缝(1/100-1/8波长)预测，则需要通过叠前地震数据和方法来进行；3)至于英尺级的极小尺度断裂，只能通过钻测井资料予以描述。对于中小尺度的裂缝精细预测是目前的研究重点。
目前国内外的裂缝预测方法主要针对具有近水平对称轴的近垂直裂缝(HTI)，主要依靠地震波传播过程中的运动与动力学参数依据HTI反射近似进行裂缝密度和方位的反演(Ruger等，“Using AVO for fracture detection：Analytic basis and practicalsolution”，The Leading Edge，1997年16期)。然而，实际的裂缝发育绝非仅是垂直裂缝，其发育角度是多变的，例如TTI介质(具有倾斜对称轴的横向各向同性介质)，但目前关于TTI介质的弹性波动形式研究十分薄弱，尚无可用于裂缝性质参数反演的方法。此外，实验室观测表明地震各向异性与频率有关，而且当裂缝介质内饱含不同流体时，各向异性强度相对于地震频率的变化梯度(各向异性频散)是有差异的。但目前关于裂缝-孔隙介质参数引起的频散或衰减特征尚不清晰。
综上所述，现有的裂缝预测方法主要还是针对双相HTI介质，Ali提出了基于Chapman理论反演HTI介质裂缝走向、裂缝密度、裂缝开度和裂缝半径并估计渗透率(“Anisotropic permeability in fractured reservoirs from frequency-dependentseismic amplitude versus angle and azimuth data”，Geophysical Prospecting，2013年62期)，在对于双相TTI介质裂缝预测研究方面几乎空白。因此，有必要研究TTI介质裂缝性质参数的反演方法。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
裂缝既可以作为重要的油气储集空间，也可以作为油气运移的渗流通道。国内外无数的勘探与开发经验表明：有利的裂缝发育带往往预示着具有强产能的优质储层集合体，因而裂缝的精细预测与评价对勘探、开发与生产都至关重要，特别是对于我国具有复杂构造与沉积背景的东部砂岩和西部碳酸盐岩储层，以及日益重要的裂缝型非常规储层，例如致密砂岩、非常规碳酸盐岩和泥页岩储层等。AVAZ方法是裂缝储层评价的重要手段，然而它受到静态等效地震岩石物理模型的限制，一般用于反演裂缝密度和裂缝走向。随着裂缝性孔隙介质等效模型的发展，频变AVAZ响应特征被证明可以携带更多的裂缝性质信息。本发明针对现有技术中的以上问题，针对TTI介质提出了基于Chapman理论分析流体类型和裂缝长度对频变AVAZ响应特征的影响并建立反演裂缝密度、裂缝走向、裂缝长度、裂缝充填流体类型和裂缝旋转角度等裂缝性质参数的反演方法。
本发明提出了一种TTI介质裂缝性质参数的反演方法，该方法可以包括以下步骤：
建立TTI介质PP波反射系数与纵横波模量、各向异性参数和密度的关系式；
基于所述各向异性参数与密度归一化的TTI介质弹性参数矩阵的关系将所述反射系数与所述弹性参数矩阵相关联；
将所述反射系数表示为频率以及裂缝性质参数的函数；
基于所述反射系数构建反演裂缝性质参数的目标函数；
通过对所述目标函数反演得到的裂缝性质参数进行归一化概率密度计算，将归一化概率密度最大值的位置所对应的裂缝性质参数的数值确定为反演的最终结果。
优选地，所述裂缝性质参数包括裂缝密度、裂缝走向、裂缝长度、裂缝充填流体类型和裂缝旋转角度。
优选地，所述反射系数与纵横波模量、各向异性参数和密度的关系式为：
其中，表示TTI介质PP波反射系数，是方位角，θ是入射角，α和β是各向同性面上的纵波速度和横波速度，Z是垂向纵波波阻抗，Z＝ρα，G是垂向剪切模量，G＝ρβ2，δ、γ、χ、ε为各向异性参数，ρ为密度，上置符号“—”代表上下两层介质参数的平均值，前置符号“Δ”代表上下两层介质参数的差，为各向同性介质PP波反射系数，其计算公式为：
优选地，所述归一化的TTI介质弹性参数矩阵为：
优选地，所述各性异性参数δ、γ、χ、ε与所述弹性参数矩阵的分量的关系为：
优选地，所述目标函数为：
其中，f为频率，e为裂缝密度，l为裂缝长度，t为裂缝充填流体类型，o为裂缝旋转角度，I为入射角样点数，J为方位角样点数，K为频率样点数，N为裂缝密度样点数，P为裂缝长度样点数，Q为裂缝充填流体类型样点数，Rijk是在入射角θi，方位角和频率fk情况下的反射系数观测数值。
优选地，通过最小二乘法对所述目标函数进行反演。
优选地，裂缝性质参数的归一化概率密度通过以下步骤计算：
确定裂缝性质参数的概率密度函数f(x)：
其中，x是裂缝性质参数的可能值，μ是对所述目标函数反演得到该裂缝性质参数的值，σ是反射系数观测数值与将反演得到的裂缝性质参数值代入公式(1)的差值；
裂缝性质参数的归一化概率密度为函数f(x)与其最大值的比值。
本发明针对TTI介质，通过分析流体类型和TTI反射系数表达式构建了反演的目标函数，对目标函数进行反演并通过计算整个空间步长的概率密度得到裂缝密度、裂缝走向、裂缝长度、裂缝充填流体类型和裂缝旋转角度在各自的空间步长上的归一化概率密度，归一化概率密度最大值对应的步长位置即裂缝密度、裂缝走向、裂缝长度、裂缝充填流体类型和裂缝旋转角度的预测结果。该方法增加了预测的裂缝性质参数即裂缝旋转角度，可以更详细的描述裂缝发育情况。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1为根据本发明的一个实施例的TTI介质裂缝性质参数的反演方法的流程图。
图2a-2e分别为方位角、裂缝密度、裂缝长度、裂缝充填流体类型、裂缝旋转角度的归一化概率密度图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
实施例
图1为根据本发明的一个实施例的TTI介质裂缝性质参数的反演方法的流程图。
该方法可以包括以下步骤：
(1)建立TTI介质PP波反射系数与纵横波模量、各向异性参数和密度的关系式。
基于Ivan Psencik的理论，可以建立TTI介质PP波反射系数与纵横波模量、各向异性参数和密度的关系式如下：
其中，表示TTI介质PP波反射系数，是方位角，θ是入射角，方位角是表征裂缝走向的参数；α和β是各向同性面上的纵波速度和横波速度，Z是垂向纵波波阻抗，Z＝ρα，G是垂向剪切模量，G＝ρβ2，α、β、Z、G属于纵横波模量；δ、γ、χ、ε属于各向异性参数，上置符号“—”代表上下两层介质参数的平均值，前置符号“Δ”代表上下两层参数的差，为各向同性介质PP波反射系数，其计算公式为：
(2)基于所述各向异性参数与密度归一化的TTI介质弹性参数矩阵的关系将所述反射系数与所述弹性参数矩阵相关联。
密度归一化的TTI介质弹性参数矩阵为：
所述各性异性参数δ、γ、χ、ε与所述弹性参数矩阵的分量的关系为：
(3)将所述反射系数表示为频率以及裂缝性质参数的函数。
根据Chapman理论，弹性参数矩阵A是随频率变化并受裂缝密度、裂缝走向、裂缝长度、裂缝充填流体类型和裂缝旋转角度影响。因此，TTI介质PP波反射系数可以认为是随频率、裂缝密度、裂缝长度、裂缝充填流体类型和裂缝旋转角度变化的函数：其中，f是频率，e是裂缝密度，l是裂缝长度，t是流体类型，o是裂缝旋转角度。裂缝充填流体类型包括盐水、油、气体。
(4)基于所述反射系数构建反演裂缝性质参数的目标函数。
所构建的反演裂缝性质参数的目标函数为：
其中，f为频率，e为裂缝密度，l为裂缝长度，t为裂缝充填流体类型，o为裂缝旋转角度，I为入射角样点数，J为方位角样点数，K为频率样点数，N为裂缝密度样点数，P为裂缝长度样点数，Q为裂缝充填流体类型样点数，Rij k是在入射角θi，方位角和频率fk情况下的反射系数观测数值。
(5)通过对所述目标函数进行反演得到裂缝性质参数的归一化概率密度，将归一化概率密度最大值的位置所对应的裂缝性质参数的数值确定为反演的最终结果。
首先，对目标函数(12)进行最小二乘反演，得到一组裂缝密度、裂缝走向、裂缝长度、裂缝充填流体类型和裂缝旋转角度的反演结果作为下一步的初始值。接下来，根据上一步的初始值，利用蒙特卡洛马尔科夫链的方法，在整个空间步长上分别搜索任意一个参数不同空间步长下得到的归一化概率密度，归一化概率密度值为1时所对应的值则是裂缝性质参数的反演结果。
裂缝性质参数的归一化概率密度通过以下步骤计算：
确定裂缝性质参数的概率密度函数f(x)：
其中，x是裂缝性质参数的可能值，μ是对所述目标函数反演得到该裂缝性质参数的值，σ是反射系数观测数值与将反演得到的裂缝性质参数值代入公式(1)的差值；
裂缝性质参数的归一化概率密度为函数f(x)与其最大值的比值。
应用示例
下面对某研究区应用本发明提出的方法对TTI介质裂缝性质参数进行反演。
基于Chapman理论计算裂缝旋转角度为0°的TTI介质(即VTI介质)的等效弹性参数，此时由于考虑了速度频散与衰减的现象，得到的弹性参数是随频率变化的。推导TTI介质的纵波反射系数表达式，将其与TTI介质双相等效弹性参数相结合，得到TTI介质模型随频率变化的反射系数即模型数据。将实际观测到的反射系数在不同频率下的观测值与模型模拟值进行拟合，利用蒙特卡洛的方法进行参数步长空间范围内概率密度计算，最终实现裂缝走向、裂缝密度、裂缝长度、裂缝填充类型和裂缝旋转角度的预测。
对TTI介质进行数值模拟，数值模拟过程中参数值即真实值为：方位角0°，裂缝密度0.08，裂缝长度1m，裂缝类型为盐水饱和，裂缝旋转角度为60°，数值模拟出不同频率状态下纵波反射系数随方位角和入射角的变化。根据图1归纳的反演步骤，可反演得到各个参数的归一化概率密度，归一化概率密度最大值位置为反演的最终结果，如图2所示，方位角反演结果为0°-10°，裂缝密度反演结果为0.08，裂缝长度反演结果为1m-1.2m，裂缝类型反演结果为盐水填充，裂缝旋转角度反演结果为56°。对比真实值与反演结果对比，裂缝密度、裂缝长度和裂缝类型的反演结果与真实值吻合，方位角的反演结果区间包含了真实值，裂缝旋转角度的反演结果与真实值存在6％的误差。与目前HTI介质裂缝性质反演对比，增加了一个描述裂缝性质的参数即裂缝旋转角度，因此，基于本发明的反演结果基本与真实值吻合，验证了该方法的有效性。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。
Claims (8)
1.一种TTI介质裂缝性质参数的反演方法，其特征在于，所述方法包括以下步骤：
建立TTI介质PP波反射系数与纵横波模量、各向异性参数和密度的关系式；
基于所述各向异性参数与密度归一化的TTI介质弹性参数矩阵的关系将所述反射系数与所述弹性参数矩阵相关联；
将所述反射系数表示为频率以及裂缝性质参数的函数；
基于所述反射系数构建反演裂缝性质参数的目标函数；
通过对所述目标函数反演得到的裂缝性质参数进行归一化概率密度计算，将归一化概率密度最大值的位置所对应的裂缝性质参数的数值确定为反演的最终结果。
2.根据权利要求1所述的TTI介质裂缝性质参数的反演方法，其特征在于，所述裂缝性质参数包括裂缝密度、裂缝走向、裂缝长度、裂缝充填流体类型和裂缝旋转角度。
3.根据权利要求2所述的TTI介质裂缝性质参数的反演方法，其特征在于，所述反射系数与纵横波模量、各向异性参数和密度的关系式为：
其中，表示TTI介质PP波反射系数，是方位角，θ是入射角，α和β是各向同性面上的纵波速度和横波速度，Z是垂向纵波波阻抗，Z＝ρα，G是垂向剪切模量，G＝ρβ2，δ、γ、χ、ε为各向异性参数，ρ为密度，上置符号“—”代表上下两层介质参数的平均值，前置符号“Δ”代表上下两层介质参数的差，为各向同性介质PP波反射系数，其计算公式为：
4.根据权利要求3所述的TTI介质裂缝性质参数的反演方法，其特征在于，所述密度归一化的TTI介质弹性参数矩阵为：
5.根据权利要求4所述的TTI介质裂缝性质参数的反演方法，其特征在于，所述各性异性参数δ、γ、χ、ε与所述弹性参数矩阵的分量的关系为：
6.根据权利要求5所述的TTI介质裂缝性质参数的反演方法，其特征在于，所述目标函数为：
其中，f为频率，e为裂缝密度，l为裂缝长度，t为裂缝充填流体类型，o为裂缝旋转角度，I为入射角样点数，J为方位角样点数，K为频率样点数，N为裂缝密度样点数，P为裂缝长度样点数，Q为裂缝充填流体类型样点数，Rijk是在入射角θi，方位角和频率fk情况下的反射系数观测数值。
7.根据权利要求6所述的TTI介质裂缝性质参数的反演方法，其特征在于，通过最小二乘法对所述目标函数进行反演。
8.根据权利要求1所述的TTI介质裂缝性质参数的反演方法，其特征在于，裂缝性质参数的归一化概率密度通过以下步骤计算：
确定裂缝性质参数的概率密度函数f(x)：
其中，x是裂缝性质参数的可能值，μ是对所述目标函数反演得到该裂缝性质参数的值，σ是反射系数观测数值与将反演得到的裂缝性质参数值代入公式(1)的差值；
裂缝性质参数的归一化概率密度为函数f(x)与其最大值的比值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710186472.4A CN108627871B (zh) | 2017-03-24 | 2017-03-24 | 一种tti介质裂缝性质参数的反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710186472.4A CN108627871B (zh) | 2017-03-24 | 2017-03-24 | 一种tti介质裂缝性质参数的反演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108627871A true CN108627871A (zh) | 2018-10-09 |
| CN108627871B CN108627871B (zh) | 2020-02-21 |
Family
ID=63706796
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710186472.4A Active CN108627871B (zh) | 2017-03-24 | 2017-03-24 | 一种tti介质裂缝性质参数的反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108627871B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957526A (zh) * | 2017-05-17 | 2018-12-07 | 中国石油化工股份有限公司 | 获得裂缝融合数据体的方法 |
| CN110161562A (zh) * | 2019-04-29 | 2019-08-23 | 中国石油大学(华东) | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
| CN110988988A (zh) * | 2019-11-25 | 2020-04-10 | 中国矿业大学(北京) | 基于垂直裂缝介质的地震波波场模拟方法及装置 |
| CN113075728A (zh) * | 2021-02-26 | 2021-07-06 | 河海大学 | 一种建立致密砂岩多尺度三维岩石物理图板的方法 |
| CN113568045A (zh) * | 2021-07-26 | 2021-10-29 | 中南大学 | 一种tti介质裂缝倾角的反演方法 |
| CN114063163A (zh) * | 2021-11-11 | 2022-02-18 | 中南大学 | 一种裂缝型储层单斜等效介质地震表征与反演方法及系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104076391A (zh) * | 2014-04-16 | 2014-10-01 | 孙赞东 | 基于tti介质四阶旅行时方程的局部角度域各向异性偏移方法 |
| CN104316965A (zh) * | 2014-10-29 | 2015-01-28 | 中国石油天然气集团公司 | 一种裂缝方位和强度的预测方法及系统 |
| CN104407378A (zh) * | 2014-11-25 | 2015-03-11 | 中国石油天然气股份有限公司 | 一种各向异性参数反演方法及装置 |
| WO2016069851A1 (en) * | 2014-10-30 | 2016-05-06 | Schlumberger Canada Limited | Inversion technique for fracture characterization in highly inclined wells using multiaxial induction measurements |
| US20160266275A1 (en) * | 2015-03-10 | 2016-09-15 | Schlumberger Technology Corporation | Methods for estimating formation parameters |
| CN106324665A (zh) * | 2015-07-01 | 2017-01-11 | 中国石油化工股份有限公司 | 反演裂缝密度的方法及系统 |
-
2017
- 2017-03-24 CN CN201710186472.4A patent/CN108627871B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104076391A (zh) * | 2014-04-16 | 2014-10-01 | 孙赞东 | 基于tti介质四阶旅行时方程的局部角度域各向异性偏移方法 |
| CN104316965A (zh) * | 2014-10-29 | 2015-01-28 | 中国石油天然气集团公司 | 一种裂缝方位和强度的预测方法及系统 |
| WO2016069851A1 (en) * | 2014-10-30 | 2016-05-06 | Schlumberger Canada Limited | Inversion technique for fracture characterization in highly inclined wells using multiaxial induction measurements |
| CN104407378A (zh) * | 2014-11-25 | 2015-03-11 | 中国石油天然气股份有限公司 | 一种各向异性参数反演方法及装置 |
| US20160266275A1 (en) * | 2015-03-10 | 2016-09-15 | Schlumberger Technology Corporation | Methods for estimating formation parameters |
| CN106324665A (zh) * | 2015-07-01 | 2017-01-11 | 中国石油化工股份有限公司 | 反演裂缝密度的方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 刘喜武等: "裂缝性孔隙介质频变AVAZ反演方法研究进展", 《石油物探》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957526A (zh) * | 2017-05-17 | 2018-12-07 | 中国石油化工股份有限公司 | 获得裂缝融合数据体的方法 |
| CN110161562A (zh) * | 2019-04-29 | 2019-08-23 | 中国石油大学(华东) | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
| CN110161562B (zh) * | 2019-04-29 | 2020-09-18 | 中南大学 | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
| CN110988988A (zh) * | 2019-11-25 | 2020-04-10 | 中国矿业大学(北京) | 基于垂直裂缝介质的地震波波场模拟方法及装置 |
| CN113075728A (zh) * | 2021-02-26 | 2021-07-06 | 河海大学 | 一种建立致密砂岩多尺度三维岩石物理图板的方法 |
| CN113075728B (zh) * | 2021-02-26 | 2021-11-09 | 河海大学 | 一种建立致密砂岩多尺度三维岩石物理图板的方法 |
| CN113568045A (zh) * | 2021-07-26 | 2021-10-29 | 中南大学 | 一种tti介质裂缝倾角的反演方法 |
| CN114063163A (zh) * | 2021-11-11 | 2022-02-18 | 中南大学 | 一种裂缝型储层单斜等效介质地震表征与反演方法及系统 |
| CN114063163B (zh) * | 2021-11-11 | 2022-11-11 | 中南大学 | 一种裂缝型储层单斜等效介质地震表征与反演方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108627871B (zh) | 2020-02-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108627871A (zh) | 一种tti介质裂缝性质参数的反演方法 | |
| US20080086287A1 (en) | Integrated Anisotropic Rock Physics Model | |
| RU2661489C1 (ru) | Способ комплексирования исходных данных для уточнения фильтрационного строения неоднородных карбонатных коллекторов | |
| CN103293563B (zh) | 一种确定油气储层岩石裂隙发育度和流体性质的方法 | |
| Xie et al. | Construction of a 3D Geomechanical Model for Development of a shale Gas Reservoir in the Sichuan Basin | |
| Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
| US10073185B2 (en) | Predicting hydraulic fracture propagation | |
| CN106405639B (zh) | 一种叠前地震储层岩性参数的反演方法 | |
| CN110318740A (zh) | 一种随钻声波测井评价地层各向异性的方法 | |
| Xiao et al. | A method to evaluate pore structures of fractured tight sandstone reservoirs using borehole electrical image logging | |
| CN104834003B (zh) | 一种非常规致密储层的相控压缩系数地震预测方法 | |
| CN104007462A (zh) | 基于衰减各向异性的裂缝预测方法 | |
| CN107797141A (zh) | 一种反演裂缝性质的方法 | |
| Varela et al. | Geomechanics: Pressure, stress field, and hydraulic fractures | |
| Gu et al. | Shale fracturing characterization and optimization by using anisotropic acoustic interpretation, 3D fracture modeling, and supervised machine learning | |
| Zeng et al. | Unreliable determination of in situ stress orientation by borehole breakouts in fractured tight reservoirs: A case study of the upper Eocene Hetaoyuan Formation in the Anpeng field, Nanxiang Basin, China | |
| CN110456412B (zh) | 一种基于叠后地震数据识别碳酸盐岩储层流体饱和度的方法 | |
| AU2015275302B2 (en) | Method and apparatus of determining stiffness coefficients of formation | |
| CN108957527A (zh) | 岩层网状裂缝的地震预测方法 | |
| WO2020153972A1 (en) | Evaluating anisotropic effective permeability in rock formations having natural fracture networks | |
| WO2016209822A1 (en) | Predicting hydraulic fracture propagation | |
| Tang et al. | Assessing rock brittleness and fracability from radial variation of elastic wave velocities from borehole acoustic logging | |
| Deng et al. | A new index used to characterize the near-wellbore fracture network in naturally fractured gas reservoirs | |
| CN109339771A (zh) | 一种页岩油气层孔隙压力预测方法及系统 | |
| Li et al. | Heterogeneous reservoir prediction of ultra-deep strike-slip fault-damaged zone constrained with local seismic anomaly data |
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 |