CN109655908A - 基于HTI介质的方位AVAz的正演模拟方法及系统 - Google Patents
基于HTI介质的方位AVAz的正演模拟方法及系统 Download PDFInfo
- Publication number
- CN109655908A CN109655908A CN201710942969.4A CN201710942969A CN109655908A CN 109655908 A CN109655908 A CN 109655908A CN 201710942969 A CN201710942969 A CN 201710942969A CN 109655908 A CN109655908 A CN 109655908A
- Authority
- CN
- China
- Prior art keywords
- orientation
- avaz
- crack
- data
- seismic
- 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
- 230000008030 elimination Effects 0.000 claims abstract description 30
- 238000003379 elimination reaction Methods 0.000 claims abstract description 30
- 238000004458 analytical method Methods 0.000 claims abstract description 25
- 238000012482 interaction analysis Methods 0.000 claims abstract description 17
- 238000004088 simulation Methods 0.000 claims description 14
- 230000003993 interaction Effects 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 11
- 238000000034 method Methods 0.000 abstract description 11
- 230000008901 benefit Effects 0.000 abstract description 5
- 230000004044 response Effects 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 3
- 238000012360 testing method Methods 0.000 abstract description 3
- 230000005540 biological transmission Effects 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000000284 extract Substances 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm 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
- 230000037361 pathway Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000003079 shale oil Substances 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
- G01V2210/632—Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
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
本发明公开了一种基于HTI介质的方位AVAz的正演模拟方法及系统,方法包括:建立裂缝型储层的岩石物理模型;基于方位角道集的AVAz特征分析计算;建立裂缝椭圆交互分析视图;基于裂缝椭圆交互分析视图,分析裂缝引起的地震振幅属性的方位AVAz特征。其优点在于:利用地震属性数据实时分析地下裂缝型储层的地震响应规律;解决地震振幅数据方位各向异性特征的不确定性导致的各向异性椭圆解释成果错误问题,确定椭圆的正确方位指示;为叠前裂缝检测提供正演分析工具,完善并扩充检测流程。在实际生产应用的效果表明,正演的解释结论能够为叠前裂缝检测提供有效的理论依据,增强了反演裂缝方位数据的准确性与可信度。
Description
技术领域
本发明涉及油气地球物理勘探领域,更具体地,涉及一种基于HTI介质的方位AVAz的正演模拟方法及系统。
背景技术
世界油气勘探开发实践表明,几乎所有的油气藏都和裂缝有关,超过60%的油气藏可以归结为裂缝型油气藏,裂缝为储层提供了基本的存储空间和运移通道,在大多数碳酸盐岩和非常规致密油气、页岩油气储层中具有至关重要的作用,所以裂缝储层一直是非常规油气藏的研究热点。通过长期的研究发现:地震波在HTI介质传播时,其反射振幅、衰减、频率和速度与传播的方位相关(AVAz),即与裂缝的走向和密度有关。根据这个原理,利用叠前地震资料提取方位地震属性如振幅、衰减、频率、方位速度等可以检测HTI或近似HTI型的裂缝。一些研究结论表明:1、纵波垂直于裂缝带传播时会产生比较明显的旅行时延迟和衰减,并有反射强度降低和频率变低等现象;2、地震频率的衰减与裂缝密度场的空间变化有关,在沿裂缝走向方向随偏移距(入射角)的增加衰减慢,而垂直裂缝走向方向随偏移距(入射角)的增加衰减快,裂缝密度越大衰减幅度越快;3、振幅随偏移距的变化(Amplitude variation with offset,AVO)随方位角的变化关系(即AVO梯度)反映了岩石硬度的变化。这些研究表明,利用叠前地震资料提取方位地震属性如振幅、速度、主频、衰减等检测裂缝型储层是完全可行的,可以实现对裂缝方位、裂缝密度等参数信息的描述。通过对与方位角有关的地震属性研究结论发现,部分地震属性数据在方位各向异性上都有比较固定的特征,但地震振幅属性却存在不确定性(表1和表2)。
表1地震属性的各向异性特征
表2方位椭圆指示关系
因此,有必要开发一种基于HTI介质的方位AVAz的正演模拟方法及系统,能够为叠前裂缝检测提供有效的理论依据,增强了反演裂缝方位数据的准确性与可信度。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种基于HTI介质的方位AVAz的正演模拟方法及系统,其能够通过利用地质、测井资料和岩石物理参数,根据岩石物理理论建立裂缝储层的理论地质模型,模拟地震波在已知裂缝模型中传播时各个方位的地震响应,并分析任一区段由裂缝引起的地震振幅属性的方位AVO(AVAZ)特征,增强了反演裂缝方位数据的准确性与可信度。
根据本发明的一方面,提出了一种基于HTI介质的方位AVAz的正演模拟方法。所述正演模拟方法包括:
建立裂缝型储层的岩石物理模型;
基于方位角道集的AVAz特征分析计算;
建立裂缝椭圆交互分析视图;
基于所述裂缝椭圆交互分析视图,分析裂缝引起的地震振幅属性的方位AVAz特征。
优选地,所述建立裂缝型储层的岩石物理模型包括:通过拾取实际井的横波、纵波、密度数据和子波数据,结合各向异性参数和分析层段参数,计算任一层段反射系数与子波的褶积,建立所述裂缝型储层的岩石物理模型。
优选地,所述基于方位角道集的AVAz特征分析计算包括:
通过对正演模拟所述方位角道集裂缝型储层的交互拾取,获取各个方位的地震属性数据;
基于所述地震属性数据,获取一套分方位、分角度的地震属性数据交汇图;
基于所述地震属性数据交汇图,实时显示方位角道集的振幅随偏移距的变化特征曲线;
通过在所述振幅随偏移距的变化特征曲线的交互拾取,提取各方位任一入射角的所述地震属性数据,拟合在椭圆方程上。
优选地,对任一平面上的一组样本点,以所述椭圆方程为模型进行拟合,并计算出所述椭圆方程的椭圆长轴长度、短轴长度、长轴与正北方向顺时针的夹角的参数值。
优选地,在HTI垂直裂缝型介质中,所述正演模拟时与正北方向的夹角为0°为裂缝走向,与正北方向的夹角为90°时代表裂缝法向。
优选地,所述建立裂缝椭圆交互分析视图包括:
获取所述各个方位的任一时间点振幅数据,投影至所述地震属性数据交汇图;
基于所述地震属性数据交汇图,获取任一入射角的所述各个方位的振幅数据;
将所述各个方位的振幅数据投影至极坐标;
基于所述椭圆方程的拟合,获取裂缝的长短轴指示方向。
优选地,基于所述HTI垂直裂缝型介质的方位定义,分析各向异性椭圆方程视图,得到裂缝方向指示为长轴,即该层位所述各向异性椭圆方程的长轴平行于所述裂缝方向。
优选地,所述地震振幅属性的方位AVAz特征包括实时获取所述各向异性椭圆方程长轴的方向与裂缝方向的关系和裂缝模型的密度。
根据本发明的另一方面,提出了一种基于HTI介质的方位AVAz的正演模拟系统,所述正演模拟系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
建立裂缝型储层的岩石物理模型;
基于方位角道集的AVAz特征分析计算;
建立裂缝椭圆交互分析视图;
基于所述裂缝椭圆交互分析视图,分析裂缝引起的地震振幅属性的方位AVAz特征。
优选地,所述建立裂缝型储层的岩石物理模型包括:通过拾取实际井的横波、纵波、密度数据和子波数据,结合各向异性参数和分析层段参数,计算任一层段反射系数与子波的褶积,建立所述裂缝型储层的岩石物理模型。
根据本发明的一种基于HTI介质的方位AVAz的正演模拟方法及系统,其优点在于:利用地震属性数据实时分析地下裂缝型储层的地震响应规律;解决地震振幅数据方位各向异性特征的不确定性导致的各向异性椭圆解释成果错误问题,确定椭圆的正确方位指示;为叠前裂缝检测提供正演分析工具,完善并扩充检测流程。在实际生产应用的效果表明,正演的解释结论能够为叠前裂缝检测提供有效的理论依据,增强了反演裂缝方位数据的准确性与可信度。
本发明的方法及系统具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的,或者将在并入本文中的附图和随后的具体实施例中进行详细陈述,这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的附图标记通常代表相同部件。
图1示出了根据本发明的一个示例性实施例的一种基于HTI介质的方位AVAz的正演模拟方法的步骤的流程图。
图2示出了根据本发明的一个示例性实施例的裂缝型储层的地质模型。
图3示出了本发明图2的裂缝型储层的地质模型定义的示意图。
图4a、图4b和图4c分别示出了根据本发明的一个示例性实施例的裂缝椭圆交互分析视图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
本发明提出了一种基于HTI介质的方位AVAz的正演模拟方法,该正演模拟方法包括:
建立裂缝型储层的岩石物理模型;
基于方位角道集的AVAz特征分析计算;
建立裂缝椭圆交互分析视图;
基于裂缝椭圆交互分析视图,分析裂缝引起的地震振幅属性的方位AVAz特征。
其中,HTI介质反射、透射系数方程完整地表述了该类型介质中平面波反射、透射系数与入射角的关系,但其方程组解析解的表达式十分复杂,很难直接用来分析介质参数对振幅变化的影响。为了明确地表达反射系数、透射系数与弹性常数的关系,利用近似解的方式简化HTI介质P波入射时的振幅随偏移距的变化(Amplitude variation with offset,AVO)方程:
式中,为垂直方向纵波的阻抗;
为横波的模量;
为界面上下两侧垂向纵波速度的平均值;
为横波速度的平均值;
Δε为上下两层各向异性参数差值;
Δγ为上下两层各向异性参数差值;
Δδ为上下两层各向异性参数差值;
ρ为密度。
作为优选方案,建立裂缝型储层的岩石物理模型包括:通过拾取实际井的横波(S波)、纵波(P波)、密度数据和子波数据,结合各向异性参数和分析层段参数,计算任一层段反射系数与子波的褶积,建立裂缝型储层的岩石物理模型,即裂缝型储层的地质模型,如图2所示。
利用地震属性数据实时分析地下裂缝型储层的地震响应规律。
作为优选方案,基于方位角道集的AVAz特征分析计算包括:
通过对正演模拟方位角道集裂缝型储层的交互拾取,获取各个方位的地震属性数据;
基于地震属性数据,获取一套分方位、分角度的地震属性数据交汇图;
基于地震属性数据交汇图,实时显示方位角道集的AVO特征曲线;
通过在AVO特征曲线的交互拾取,提取各方位任一入射角的地震属性数据,拟合在椭圆方程上。
其中,AVAz为振幅随方位角、入射角的变化(Amplitude variation with Angleand Azimuth)。
解决地震振幅数据方位各向异性特征的不确定性导致的各向异性椭圆解释成果错误问题,确定椭圆的正确方位指示。
根据叠前各向异性裂缝检测的原理,在HTI介质中,纵波地震资料的AVO梯度和动校正(Normal Moveout,NMO)速度等多个参数都随着方位角呈椭圆形变化,这一特点是利用纵波地震资料检测裂缝的技术基础。
为叠前裂缝检测提供正演分析工具,完善并扩充检测流程。
作为优选方案,对任一平面上的一组样本点,以所述椭圆方程为模型进行拟合,并计算出所述椭圆方程的椭圆长轴长度、短轴长度、长轴与正北方向顺时针的夹角的参数值。
方位椭圆拟合中,椭圆方程为:
Ax2+Bxy+Cy2+Dx+Ey+F=0
可以定义一点(x,y)到椭圆的几何距离为:
F(a,x)=a·x=Ax2+Bxy+Cy2+Dx+Ey+F
其中,a=[A B C D E F]T,x=[x2 xy y2 x y 1]T;
当参与拟合的各点与椭圆的几何距离的平方和最小时,即可得到最小二乘意义下的椭圆拟合参数最优解,得到:
写成矩阵形式如下:
D·a=0
展开即为:
该方程可写成广义特征值问题的形式:
DTDa=λCa
其中,λ为特征值,C为约束矩阵,由于椭圆方程存在以下关系:
则C有如下形式:
求解以上方程,可得六个特征值λ,其中仅有唯一的λ∈R+,及其对应的特征值向量a为方程的解,此时a中的元素即为最优椭圆参数值,通过解析几何中椭圆一般方程与椭圆参数的关系,可按以下公式计算椭圆几何参数:
上述各式中,x,y为椭圆中心相对坐标;
a,b分别为椭圆的长轴和短轴;
θ为椭圆长轴与正北方向的夹角。
作为优选方案,在HTI垂直裂缝型介质中,正演模拟时与正北方向的夹角为0°为裂缝走向,与正北方向的夹角为90°时代表裂缝法向。
其中,建立裂缝椭圆交互分析视图包括:
获取各个方位的任一时间点振幅数据,投影至地震属性数据交汇图;
基于地震属性数据交汇图,获取任一入射角的各个方位的振幅数据;
将各个方位的振幅数据投影至极坐标;
基于椭圆方程的拟合,获取裂缝的长短轴指示方向。
其中,获取如图2所示的各个方位的某时间点振幅数据,投影到AVAz交汇视图,再通过AVAz交汇视图获取某入射角(图中范围是0-30°)各个方位的振幅数据,再将各个方位的振幅数据投影到极坐标,利用椭圆拟合方法,求取裂缝的长短轴指示方向,沿裂缝方向,振幅随入射角增大而增大;垂直于裂缝方向,振幅随入射角增大而减小,根据HTI介质裂缝模型的方位定义,与正北方向的夹角为0度且代表裂缝走向,再分析各向异性椭圆视图,裂缝方向指示为长轴,代表该层位各向异性椭圆的长轴平行于裂缝方向,这个结论将应用于对实际数据的成果解释。
作为优选方案,地震振幅属性的方位AVAz特征包括实时获取各向异性椭圆方程长轴的方向与裂缝方向的关系和裂缝模型的密度。
在实际生产应用的效果表明,正演的解释结论能够为叠前裂缝检测提供有效的理论依据,增强了反演裂缝方位数据的准确性与可信度。
实施例
图1示出了根据本发明的一个示例性实施例的一种基于HTI介质的方位AVAz的正演模拟方法的步骤的流程图。
如图1所示,本实施例的一种基于HTI介质的方位AVAz的正演模拟方法,包括:
建立裂缝型储层的岩石物理模型;
基于方位角道集的AVAz特征分析计算;
建立裂缝椭圆交互分析视图;
基于裂缝椭圆交互分析视图,分析裂缝引起的地震振幅属性的方位AVAz特征。
其中,HTI介质反射、透射系数方程完整地表述了该类型介质中平面波反射、透射系数与入射角的关系,但其方程组解析解的表达式十分复杂,很难直接用来分析介质参数对振幅变化的影响。为了明确地表达反射系数、透射系数与弹性常数的关系,利用近似解的方式简化HTI介质P波入射时的AVO方程:
式中,为垂直方向纵波的阻抗;
为横波的模量;
为界面上下两侧垂向纵波速度的平均值;
为横波速度的平均值;
Δε、Δγ、Δδ均为上下两层各向异性参数差值
ρ为密度。
图2示出了根据本发明的一个示例性实施例的裂缝型储层的地质模型。
如图2所示,建立裂缝型储层的岩石物理模型包括:通过拾取实际井的横波(S波)、纵波(P波)、密度数据和子波数据,结合各向异性参数和分析层段参数,计算任一层段反射系数与子波的褶积,建立裂缝型储层的岩石物理模型,即裂缝型储层的地质模型。
其中,通过对正演模拟方位角道集裂缝型储层的交互拾取,获取各个方位的地震属性数据,开发一套分方位、分角度的地震属性数据交汇图,实时显示方位角道集的AVO特征曲线,通过用户在AVO特征曲线的交互拾取,提取各方位某入射角的属性数据,在拟合椭圆上展布。根据叠前各向异性裂缝检测的原理,在HTI介质中,纵波地震资料的AVO梯度和NMO速度等多个参数都随着方位角呈椭圆形变化,这一特点是利用纵波地震资料检测裂缝的技术基础。对给定平面上的一组样本点,以椭圆方程为模型进行拟合,使某一椭圆方程尽量满足这些数据点,并求出该椭圆方程的各个参数。
方位椭圆拟合中,椭圆方程为:
Ax2+Bxy+Cy2+Dx+Ey+F=0
可以定义一点(x,y)到椭圆的几何距离为:
F(a,x)=a·x=Ax2+Bxy+Cy2+Dx+Ey+F
其中,a=[A B C D E F]T,x=[x2 xy y2 x y 1]T;
当参与拟合的各点与椭圆的几何距离的平方和最小时,即可得到最小二乘意义下的椭圆拟合参数最优解,得到:
写成矩阵形式如下:
D·a=0
展开即为:
该方程可写成广义特征值问题的形式:
DTDa=λCa
其中,λ为特征值,C为约束矩阵,由于椭圆方程存在以下关系:
则C有如下形式:
求解以上方程,可得六个特征值λ,其中仅有唯一的λ∈R+(此处引用了一篇文献的结论”Andrew Fitzgibbon,Pilu,Robert B.Fisher.Direct Least Square Fitting ofEllipses[J].IEEE Transaction on Pattern Analysis and Machine Intelligence,1999,21(5):476-480”中已经证明该广义特性系统只有1个正的特征值,并且这个特征值对应的特征向量为所求取的椭圆的参数向量),及其对应的特征值向量a为方程的解,此时a中的元素即为最优椭圆参数值,通过解析几何中椭圆一般方程与椭圆参数的关系,可按以下公式计算椭圆几何参数:
上述各式中,x,y为椭圆中心相对坐标;
a,b分别为椭圆的长轴和短轴;
θ为椭圆长轴与正北方向的夹角。
图3示出了本发明图2的裂缝型储层的地质模型定义的示意图。
如图3所示,在HTI垂直裂缝型介质中,正演模拟时与正北方向的夹角为0°为裂缝走向,与正北方向的夹角为90°时代表裂缝法向。
获取上述图2中各个方位的某时间点振幅数据,投影到AVAZ交汇视图,再通过AVAZ交汇视图获取某入射角(图中范围是0-30°)各个方位的振幅数据,再将各个方位的振幅数据投影到极坐标,利用椭圆拟合方法,求取裂缝的长短轴指示方向。
图4a、图4b和图4c分别示出了根据本发明的一个示例性实施例的裂缝椭圆交互分析视图。
如图4a-图4c所示,沿裂缝方向,振幅随入射角增大而增大;垂直于裂缝方向,振幅随入射角增大而减小,根据HTI介质裂缝模型的方位定义,与正北方向的夹角为0度且代表裂缝走向,再分析各向异性椭圆视图,裂缝方向指示为长轴,代表该层位各向异性椭圆的长轴平行于裂缝方向,这个结论将应用于对实际数据的成果解释。另外,10°入射角示例中,椭圆比率1.092代表了模型定义的裂缝密度,随着入射角的增大,椭圆比率变大,裂缝密度也逐渐增大。
其中,包括井旁地震道合成记录视图、裂缝型储层的岩石物理模型合成记录视图、裂缝模型分方位地震响应模拟视图、以及各向异性椭圆分析视图等,视图之间建立联动机制,通过鼠标在层位上的拾取操作,实时获取各向异性椭圆长轴的方向与裂缝方向的关系以及裂缝模型的密度。
以上已经描述了本发明的实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的实施例。在不偏离所说明的实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释实施例的原理、实际应用或对市场中的技术的改进,或者使本技术领域的其它普通技术人员能理解本文披露的实施例。
Claims (10)
1.一种基于HTI介质的方位AVAz的正演模拟方法,其特征在于,所述正演模拟方法包括:
建立裂缝型储层的岩石物理模型;
基于方位角道集的AVAz特征分析计算;
建立裂缝椭圆交互分析视图;
基于所述裂缝椭圆交互分析视图,分析裂缝引起的地震振幅属性的方位AVAz特征。
2.根据权利要求1所述的基于HTI介质的方位AVAz的正演模拟方法,其中,所述建立裂缝型储层的岩石物理模型包括:通过拾取实际井的横波、纵波、密度数据和子波数据,结合各向异性参数和分析层段参数,计算任一层段反射系数与子波的褶积,建立所述裂缝型储层的岩石物理模型。
3.根据权利要求1所述的基于HTI介质的方位AVAz的正演模拟方法,其中,所述基于方位角道集的AVAz特征分析计算包括:
通过对正演模拟所述方位角道集裂缝型储层的交互拾取,获取各个方位的地震属性数据;
基于所述地震属性数据,获取一套分方位、分角度的地震属性数据交汇图;
基于所述地震属性数据交汇图,实时显示方位角道集的振幅随偏移距的变化特征曲线;
通过在所述振幅随偏移距的变化特征曲线的交互拾取,提取各方位任一入射角的所述地震属性数据,拟合在椭圆方程上。
4.根据权利要求3所述的基于HTI介质的方位AVAz的正演模拟方法,其中,对任一平面上的一组样本点,以所述椭圆方程为模型进行拟合,并计算出所述椭圆方程的椭圆长轴长度、短轴长度、长轴与正北方向顺时针的夹角的参数值。
5.根据权利要求3所述的基于HTI介质的方位AVAz的正演模拟方法,其中,在HTI垂直裂缝型介质中,所述正演模拟时与正北方向的夹角为0°为裂缝走向,与正北方向的夹角为90°时代表裂缝法向。
6.根据权利要求5所述的基于HTI介质的方位AVAz的正演模拟方法,其中,所述建立裂缝椭圆交互分析视图包括:
获取所述各个方位的任一时间点振幅数据,投影至所述地震属性数据交汇图;
基于所述地震属性数据交汇图,获取任一入射角的所述各个方位的振幅数据;
将所述各个方位的振幅数据投影至极坐标;
基于所述椭圆方程的拟合,获取裂缝的长短轴指示方向。
7.根据权利要求6所述的基于HTI介质的方位AVAz的正演模拟方法,其中,基于所述HTI垂直裂缝型介质的方位定义,分析各向异性椭圆方程视图,得到裂缝方向指示为长轴,即该层位所述各向异性椭圆方程的长轴平行于所述裂缝方向。
8.根据权利要求7所述的基于HTI介质的方位AVAz的正演模拟方法,其中,所述地震振幅属性的方位AVAz特征包括实时获取所述各向异性椭圆方程长轴的方向与裂缝方向的关系和裂缝模型的密度。
9.一种基于HTI介质的方位AVAz的正演模拟系统,其特征在于,所述正演模拟系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
建立裂缝型储层的岩石物理模型;
基于方位角道集的AVAz特征分析计算;
建立裂缝椭圆交互分析视图;
基于所述裂缝椭圆交互分析视图,分析裂缝引起的地震振幅属性的方位AVAz特征。
10.根据权利要求9所述的基于HTI介质的方位AVAz的正演模拟系统,其中,所述建立裂缝型储层的岩石物理模型包括:通过拾取实际井的横波、纵波、密度数据和子波数据,结合各向异性参数和分析层段参数,计算任一层段反射系数与子波的褶积,建立所述裂缝型储层的岩石物理模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710942969.4A CN109655908A (zh) | 2017-10-11 | 2017-10-11 | 基于HTI介质的方位AVAz的正演模拟方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710942969.4A CN109655908A (zh) | 2017-10-11 | 2017-10-11 | 基于HTI介质的方位AVAz的正演模拟方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109655908A true CN109655908A (zh) | 2019-04-19 |
Family
ID=66109642
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710942969.4A Pending CN109655908A (zh) | 2017-10-11 | 2017-10-11 | 基于HTI介质的方位AVAz的正演模拟方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109655908A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112198549A (zh) * | 2019-07-08 | 2021-01-08 | 中国石油天然气集团有限公司 | 一种基于地震正演模板的叠前裂缝确定方法及系统 |
CN112305594A (zh) * | 2019-07-23 | 2021-02-02 | 中国石油天然气股份有限公司 | 非均质储层的油气分布确定方法及系统 |
CN112462421A (zh) * | 2020-10-30 | 2021-03-09 | 中国石油天然气集团有限公司 | 储层信息预测方法、装置、电子设备及存储介质 |
CN113075727A (zh) * | 2021-02-20 | 2021-07-06 | 中国石油天然气股份有限公司 | 基于振幅方位导数的储层裂缝分布特征确定方法及系统 |
CN113419280A (zh) * | 2021-06-21 | 2021-09-21 | 大庆油田有限责任公司 | 基于改进椭圆拟合的叠前裂缝密度估算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0756182A1 (fr) * | 1995-07-26 | 1997-01-29 | Institut Français du Pétrole | Méthode et dispositif de filtrage d'ondes elliptiques se propageant dans un milieu |
CN103869359A (zh) * | 2014-02-25 | 2014-06-18 | 中国石油天然气股份有限公司 | 地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 |
CN104155701A (zh) * | 2014-07-03 | 2014-11-19 | 中国石油天然气集团公司 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
CN107045143A (zh) * | 2017-02-22 | 2017-08-15 | 中国石油天然气股份有限公司 | 一种预测裂缝发育的方法及装置 |
-
2017
- 2017-10-11 CN CN201710942969.4A patent/CN109655908A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0756182A1 (fr) * | 1995-07-26 | 1997-01-29 | Institut Français du Pétrole | Méthode et dispositif de filtrage d'ondes elliptiques se propageant dans un milieu |
CN103869359A (zh) * | 2014-02-25 | 2014-06-18 | 中国石油天然气股份有限公司 | 地震纵波多方位属性椭圆拟合预测裂缝的方法及装置 |
CN104155701A (zh) * | 2014-07-03 | 2014-11-19 | 中国石油天然气集团公司 | 一种利用叠前地震资料及井信息的多尺度裂缝预测方法 |
CN107045143A (zh) * | 2017-02-22 | 2017-08-15 | 中国石油天然气股份有限公司 | 一种预测裂缝发育的方法及装置 |
Non-Patent Citations (3)
Title |
---|
孙炜,等: "基于HTI介质各向异性正演的裂缝预测属性优选", 《石油物探》 * |
陈楠,等: "HTI介质方位AVO(AVAZ)正演模拟软件研制与应用", 《地球物理学进展》 * |
陈楠,等: "基于NEWS平台的叠前叠后联合解释系统设计与开发", 《华北地震科学》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112198549A (zh) * | 2019-07-08 | 2021-01-08 | 中国石油天然气集团有限公司 | 一种基于地震正演模板的叠前裂缝确定方法及系统 |
CN112198549B (zh) * | 2019-07-08 | 2024-05-28 | 中国石油天然气集团有限公司 | 一种基于地震正演模板的叠前裂缝确定方法及系统 |
CN112305594A (zh) * | 2019-07-23 | 2021-02-02 | 中国石油天然气股份有限公司 | 非均质储层的油气分布确定方法及系统 |
CN112305594B (zh) * | 2019-07-23 | 2023-10-27 | 中国石油天然气股份有限公司 | 非均质储层的油气分布确定方法及系统 |
CN112462421A (zh) * | 2020-10-30 | 2021-03-09 | 中国石油天然气集团有限公司 | 储层信息预测方法、装置、电子设备及存储介质 |
CN113075727A (zh) * | 2021-02-20 | 2021-07-06 | 中国石油天然气股份有限公司 | 基于振幅方位导数的储层裂缝分布特征确定方法及系统 |
CN113075727B (zh) * | 2021-02-20 | 2023-09-26 | 中国石油天然气股份有限公司 | 基于振幅方位导数的储层裂缝分布特征确定方法及系统 |
CN113419280A (zh) * | 2021-06-21 | 2021-09-21 | 大庆油田有限责任公司 | 基于改进椭圆拟合的叠前裂缝密度估算方法 |
CN113419280B (zh) * | 2021-06-21 | 2022-03-22 | 大庆油田有限责任公司 | 基于改进椭圆拟合的叠前裂缝密度估算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109655908A (zh) | 基于HTI介质的方位AVAz的正演模拟方法及系统 | |
US10816686B2 (en) | Seismic constrained discrete fracture network | |
EP3488267B1 (en) | Seismic spectral balancing | |
CN107728204B (zh) | 基于叠前纵波各向异性的裂缝预测方法及系统 | |
CN105005080B (zh) | 一种利用振幅比率属性识别地层圈闭尖灭线的方法 | |
CN103376469B (zh) | 一种基于超声成像测井的裂缝定量评价方法 | |
CN106443774A (zh) | 一种提高起伏地表叠前深度偏移成像精度的方法 | |
CN104656136A (zh) | 基于实际模型地震模拟指导下的油气储层低频阴影识别技术 | |
CN107656308B (zh) | 一种基于时间深度扫描的共散射点叠前时间偏移成像方法 | |
CN109164485A (zh) | 一种影响低序级断层识别精度的定量分析方法 | |
Muerdter et al. | Understanding subsalt illumination through ray-trace modeling, Part 2: Dipping salt bodies, salt peaks, and nonreciprocity of subsalt amplitude response | |
CN105938203A (zh) | 一种储层特性的检测方法及装置 | |
CN109212592A (zh) | 一种叠前方位p波各向异性裂缝检测方法及装置 | |
CN107390270A (zh) | 一种基于弹性波逆时偏移ADCIGs的AVA分析方法 | |
CN103399345B (zh) | 一种潜山裂缝分布的勘测方法与装置 | |
CN106597547A (zh) | 薄储层地震准确描述方法 | |
CN109507729A (zh) | 依赖频率的avo最优频散流体因子的反演方法 | |
CN111308558B (zh) | 页岩气水平井纵波时差校正方法 | |
CN106501851B (zh) | 一种地震属性优化方法及装置 | |
CN109100803A (zh) | 微断裂的确定方法和装置 | |
CN107942373B (zh) | 基于裂缝性油气储层断裂系统检测的相干算法 | |
Ningkai et al. | Stepped and detailed seismic prediction of shallow-thin reservoirs in Chunfeng oilfield of Junggar Basin, NW China | |
Zhang et al. | Identifying minor faults on top of coalfield Ordovician limestone stratum using seismic attributes derived from azimuthally stacked data | |
CN112394398A (zh) | 一种基于岩相与地震相的地层尖灭线识别方法 | |
Yao et al. | Using seismic methods to detect connectivity of fracture networks controlled by strike-slip faults in ultra-deep carbonate reservoirs: A case study in northern tarim basin, China |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190419 |