CN112835101A - 裂缝流体因子预测方法及系统 - Google Patents
裂缝流体因子预测方法及系统 Download PDFInfo
- Publication number
- CN112835101A CN112835101A CN202110271299.4A CN202110271299A CN112835101A CN 112835101 A CN112835101 A CN 112835101A CN 202110271299 A CN202110271299 A CN 202110271299A CN 112835101 A CN112835101 A CN 112835101A
- Authority
- CN
- China
- Prior art keywords
- fluid factor
- impedance
- fracture fluid
- azimuth
- zero
- 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
- 239000012530 fluid Substances 0.000 title claims abstract description 133
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000014509 gene expression Effects 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 15
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- IAZDPXIOMUYVGZ-UHFFFAOYSA-N Dimethylsulphoxide Chemical compound CS(C)=O IAZDPXIOMUYVGZ-UHFFFAOYSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 239000011541 reaction mixture Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
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/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)
Abstract
本发明提供一种裂缝流体因子预测方法,属于流体识别技术领域,计算各向异性截距阻抗和梯度阻抗数据体,计算零阶扩展方位傅里叶系数,确定预测流体因子的最佳旋转角,通过线性拟合得到零阶扩展傅里叶系数与裂缝流体因子的关系式,结合零阶扩展傅里叶系数与裂缝流体因子的关系式,得到裂缝流体因子。本发明推导扩展方位弹性阻抗方程及其扩展方位傅里叶系数表达式,获取裂缝流体因子与零阶扩展方位傅里叶系数的相关系数确定最优的旋转角,利用反演得到的各向异性截距阻抗和梯度阻抗及最优的旋转角计算得到零阶扩展方位傅里叶系数,不需要估测裂缝弱度,避免裂缝弱度反演不准确引起的累计误差,即可实现裂缝流体因子的准确估测,流体识别精度高。
Description
技术领域
本发明涉及地震流体识别技术领域,具体涉及一种基于扩展方位傅里叶系数的裂缝流体因子预测方法及系统。
背景技术
流体因子是能够指示储层中填充流体类型的重要参数,通常通过岩石物理分析优选与储层中填充流体类型密切相关的敏感参数作为流体因子。当前常用的流体因子主要为物理意义较为明确的以泊松比及弹性模量参数等为代表的流体因子以及基于各种弹性参数组合运算的流体因子等。对裂缝型储层而言,现有的各向同性介质流体指示因子由于均未考虑到裂缝的影响,因此各向同性流体指示因子并不完全适用于裂缝中充填流体类型的识别。
基于线性滑动理论的流体因子预测,通过详细研究裂缝流体对裂缝参数的影响,提出了能够有效识别裂缝充填流体类型的指示因子KKrN(即法向柔度和切向柔度的比值)。而含气裂缝的流体指示因子数值明显大于含油裂缝的流体指示因子,且含气裂缝流体指示因子的变化特征也比含油裂缝的流体指示因子变化更明显。
传统的求取裂缝流体因子的方法是通过方位叠前地震反演首先得到弹性参数及裂缝弱度参数,然后再间接计算得到裂缝流体因子,这种间接计算方法依赖于方位叠前地震反演的精度。因为裂缝弱度的反演需要裂缝方位的先验信息,特别是在三维地震反演中,整个三维地区的裂缝方位先验信息往往难以获取,这大大降低了裂缝弱度的预测精度,而这种误差累积会进一步影响到裂缝流体因子的准确估测。
发明内容
本发明的目的在于提供一种不需要估测裂缝弱度,可避免裂缝弱度反演不准确引起的累计误差的基于扩展方位傅里叶系数的裂缝流体因子预测方法及系统,以解决上述背景技术中存在的至少一项技术问题。
为了实现上述目的,本发明采取了如下技术方案:
一方面,本发明提供的一种裂缝流体因子预测方法,包括:
计算各向异性截距阻抗和梯度阻抗数据体,计算零阶扩展方位傅里叶系数,确定预测流体因子的最佳χ角,通过线性拟合得到零阶扩展傅里叶系数与裂缝流体因子的关系式,结合零阶扩展傅里叶系数与裂缝流体因子的关系式,得到裂缝流体因子。
优选的,利用方位部分角度叠加地震数据、方位角度地震子波以及测井曲线插值得到各向异性截距阻抗和梯度阻抗低频模型,通过分方位叠前地震贝叶斯反演得到各向异性截距阻抗和梯度阻抗数据体。
优选的,利用各向异性截距阻抗和梯度阻抗反演结果,计算不同方位角下的扩展方位弹性阻抗,对扩展方位弹性阻抗进行傅里叶级数展开,得到零阶扩展方位傅里叶系数。
优选的,计算零阶扩展方位傅里叶系数与裂缝流体因子的相关系数,确定预测流体因子的最佳χ角。
优选的,结合各向异性截距阻抗和梯度阻抗数据体以及最佳χ角,计算得到扩展方位弹性阻抗数据体,对扩展方位弹性阻抗数据体进行傅里叶级数展开,得到零阶扩展方位傅里叶系数数据体,利用零阶扩展方位傅里叶系数与裂缝流体因子的关系式,将零阶扩展方位傅里叶系数数据体转化为裂缝流体因子。
第二方面,本发明提供一种基于如上所述的裂缝流体因子预测方法的裂缝流体因子预测系统,包括:
第一计算模块,用于计算各向异性截距阻抗和梯度阻抗数据体;
第二计算模块,用于计算零阶扩展方位傅里叶系数;
确定模块,用于确定预测裂缝流体因子的最佳χ角;
拟合模块,用于将零阶扩展方位傅里叶系数与裂缝流体因子进行线性拟合,得到零阶扩展傅里叶系数与裂缝流体因子的关系式;
转化模块,用于结合各向异性截距阻抗和梯度阻抗数据体、最佳χ角以及零阶扩展傅里叶系数与裂缝流体因子的关系式,将零阶扩展方位傅里叶系数转化为裂缝流体因子。
第三方面,本发明提供一种裂缝流体因子预测装置,包括如权利要求7所述的裂缝流体因子预测系统。
第四方面,本发明提供一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质包括用于执行如上所述的裂缝流体因子预测方法的指令。
第五方面,本发明提供一种电子设备,包括如上所述的非暂态计算机可读存储介质;以及能够执行所述非暂态计算机可读存储介质的所述指令的一个或多个处理器。
本发明有益效果:推导扩展方位弹性阻抗方程及其扩展方位傅里叶系数表达式,通过获取裂缝流体因子与零阶扩展方位傅里叶系数的相关系数确定最优的χ角,利用反演得到的各向异性截距阻抗和梯度阻抗及最优的χ角计算得到零阶扩展方位傅里叶系数,不需要估测裂缝弱度,可以避免裂缝弱度反演不准确引起的累计误差,即可实现裂缝流体因子的准确估测,该方法计算简单,精度较高。
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例1所述的裂缝流体因子预测系统功能原理框图。
图2为本发明实施例2所述的裂缝流体因子预测方法流程图。
图3为本发明实施例2所述的测井曲线示意图。
图4为本发明实施例2所述的不同χ角下的零阶傅里叶系数示意图。
图5为本发明实施例2所述的不同χ角情况下的零阶扩展方位傅里叶系数与裂缝流体因子的相关系数。
图6(a)为本发明实施例2所述的利用测井曲线合成四个方位的角道集数据示意图。
图6(b)为本发明实施例2所述的利用测井曲线合成四个方位的角道集中添加高斯噪声后的数据示意图。
图7为本发明实施例2所述的无噪声情况下的反演的各向异性截距阻抗和梯度阻抗示意图。
图8为本发明实施例2所述的高斯噪声情况下的反演的各向异性截距阻抗和梯度阻抗示意图。
图9(a)为本发明实施例2所述的无噪声情况下的预测的裂缝流体因子结果示意图。
图9(b)为本发明实施例2所述的高斯噪声情况下的预测的裂缝流体因子结果示意图。
具体实施方式
下面详细叙述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。
还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件和/或它们的组。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
为便于理解本发明,下面结合附图以具体实施例对本发明作进一步解释说明,且具体实施例并不构成对本发明实施例的限定。
本领域技术人员应该理解,附图只是实施例的示意图,附图中的部件并不一定是实施本发明所必须的。
实施例1
如图1所示,本实施例1提出一种裂缝流体因子预测系统,该系统包括:
第一计算模块,用于计算各向异性截距阻抗和梯度阻抗数据体;
第二计算模块,用于计算零阶扩展方位傅里叶系数;
确定模块,用于确定预测流体因子的最佳χ角;
拟合模块,用于将零阶扩展方位傅里叶系数与裂缝流体因子进行线性拟合,得到零阶扩展傅里叶系数与裂缝流体因子的关系式;
转化模块,用于结合各向异性截距阻抗和梯度阻抗数据体、最佳χ角以及零阶扩展傅里叶系数与裂缝流体因子的关系式,将零阶扩展方位傅里叶系数转化为裂缝流体因子。
在本实施例1中,所述第一计算模块被配置为:基于方位部分角度叠加地震数据、方位角度地震子波以及测井曲线插值得到各向异性截距阻抗和梯度阻抗低频模型,通过分方位叠前地震贝叶斯反演得到各向异性截距阻抗和梯度阻抗数据体。
所述第二计算模块被配置为:利用各向异性截距阻抗和梯度阻抗反演结果,计算不同χ角下的扩展方位弹性阻抗,对扩展方位弹性阻抗进行傅里叶级数展开,得到零阶扩展方位傅里叶系数。
所述确定模块被配置为:计算零阶扩展方位傅里叶系数与裂缝流体因子的相关系数,确定预测流体因子的最佳χ角。
所述转化模块被配置为:结合各向异性截距阻抗和梯度阻抗数据体以及最佳χ角,计算得到扩展方位弹性阻抗数据体,对扩展方位弹性阻抗数据体进行傅里叶级数展开,得到零阶扩展方位傅里叶系数数据体,利用零阶扩展方位傅里叶系数与裂缝流体因子的关系式,将零阶扩展方位傅里叶系数数据体转化为裂缝流体因子。
在本实施例1中,第一计算模块还包括有构建单元,用于利用测井曲线插值得到各向异性截距阻抗和梯度阻抗低频模型。
在本实施例1中,利用上述的裂缝流体因子预测系统实现了裂缝流体因子预测方法,该方法包括:
计算各向异性截距阻抗和梯度阻抗数据体,计算零阶扩展方位傅里叶系数,确定预测流体因子的最佳χ角,通过线性拟合得到零阶扩展傅里叶系数与裂缝流体因子的关系式,结合零阶扩展傅里叶系数与裂缝流体因子的关系式,得到裂缝流体因子。
在本实施例1中,计算各向异性截距阻抗和梯度阻抗数据体包括:基于方位部分角度叠加地震数据、方位角度地震子波以及测井曲线插值得到各向异性截距阻抗和梯度阻抗低频模型,通过分方位叠前地震贝叶斯反演得到各向异性截距阻抗和梯度阻抗数据体。
计算零阶扩展方位傅里叶系数包括:利用各向异性截距阻抗和梯度阻抗反演结果,计算不同χ角下的扩展方位弹性阻抗,对扩展方位弹性阻抗进行傅里叶级数展开,得到零阶扩展方位傅里叶系数。
确定预测流体因子的最佳χ角包括:计算零阶扩展方位傅里叶系数与裂缝流体因子的相关系数,确定预测流体因子的最佳χ角。
结合各向异性截距阻抗和梯度阻抗数据体以及最佳χ角,计算得到扩展方位弹性阻抗数据体,对扩展方位弹性阻抗数据体进行傅里叶级数展开,得到零阶扩展方位傅里叶系数数据体,利用零阶扩展方位傅里叶系数与裂缝流体因子的关系式,将零阶扩展方位傅里叶系数数据体转化为裂缝流体因子。
实施例2
如图2所示,本实施例2提出了一种基于扩展方位傅里叶系数的裂缝流体因子预测方法流程,具体包括以下步骤:
(1)利用方位部分角度叠加地震数据、方位角度地震子波以及测井曲线插值得到的各向异性截距阻抗和梯度阻抗低频模型,通过分方位叠前地震贝叶斯反演得到各向异性截距阻抗和梯度阻抗数据体;
(2)利用井旁道的各向异性截距阻抗和梯度阻抗反演结果计算不同χ角情况下的扩展方位弹性阻抗,对其进行傅里叶级数展开得到井旁道的零阶扩展傅里叶系数,并计算其与井中裂缝流体因子的相关系数,确定估测裂缝流体因子的最佳χ角,通过线性拟合得到零阶扩展傅里叶系数与裂缝流体因子的关系式;
(3)结合反演的各向异性截距阻抗和梯度阻抗数据体以及最佳χ角计算得到扩展方位弹性阻抗数据体,对其进行傅里叶级数展开得到零阶扩展方位傅里叶系数数据体,利用零阶扩展方位傅里叶系数与裂缝流体因子的关系式,将零阶扩展方位傅里叶系数数据体转化为裂缝流体因子。
方法原理:
HTI介质中的纵波反射系数方程如下式(1)所示:
其中,aα(θ)=sec2θ,aβ(θ)=-8gsin2θ,aρ(θ)=1-4gsin2θ;
aδT(θ,φ)=gsin2θcos2φ(1-tan2θsin2φ);
其中,α、β、ρ分别表示各向同性背景介质的纵、横波速度及密度,α0、β0、ρ0分别表示各向同性背景介质的纵、横波速度及密度的平均值;δN、δT分别表示法向及切向弱度。ΔδN=δN2-δN1、ΔδT=δT2-δT1分别表示上下地层之间的法向弱度差值和切向弱度差值。g=β2/α2表示横波速度平方与纵波速度平方之比,通常设置为一个常数。θ表示入射角,表示地震测线观测方位与裂缝对称轴方位之间的夹角。表示地震测线观测方位,表示裂缝对称轴方位。
方程(1)可重写为:
RPP(θ,φ)=A+B(φ)sin2θ+C(φ)sin2θtan2θ (4)
当入射角小于30°时,省略掉大角度项(sin2θtan2θ),得到小角度近似式,如式(8)所示:
RPP(θ,φ)=A+B(φ)sin2θ (8)
类比各向同性介质中的扩展弹性阻抗推导过程,将tanχ=sin2θ带入方程(8),左右两边同时乘以cosχ,得到:
RPP(χ,φ)=Rcosχ=Acosχ+B(φ)sinχ (9)
其中,χ表示旋转角,RPP(χ,φ)表示尺度化的方位反射系数。
通过推导方程(9)得到HTI介质中的扩展方位弹性阻抗方程如下:
式中,AIiso=αρ、分别表示各向同性截距阻抗和梯度阻抗。AI、BI(φ)分别表示各向异性截距阻抗和梯度阻抗,AI0=α0ρ0表示纵波阻抗均值。EEI(χ,φ)表示扩展方位弹性阻抗,当χ为0°时,EEI(χ,φ)=AI表示各向异性截距阻抗;当χ为90°时,EEI(χ,φ)=BI(φ)表示各向异性梯度阻抗。
对方程(10)两边同时取对数,得到:
LEEI(χ,φ)=cosχLAI+sinχLBI (13)
对方程(13)进行傅里叶级数展开,得到:
LEEI(χ,φ)=R0(χ)+R2(χ)cos(2φ) (14)
R2(χ)=[g(δT-kδN)]sinχ (16)
式中,Rn(χ),n=(0,2)表示HTI介质中的扩展方位傅里叶系数。二阶扩展方位傅里叶系数表示尺度化的各向异性梯度(Bani=g(δT-kδN))。
方程(14)可表示为另一种形式:
其中,Un(χ)、Vn(χ)分别表示扩展方位傅里叶系数Rn(χ),n=(0,2)的余弦分量和正弦分量,其表达式分别为:
在裂缝密度较小情况下,切向弱度δT不随裂缝填充物变化,而法向弱度δN则受裂缝填充物的影响。为了指示裂缝填充物的变化,基于裂缝流体因子:
当裂缝完全充填流体时,裂缝流体因子为0;当裂缝完全充填气体时,裂缝流体因子接近1;因此裂缝流体因子能够较好区分裂缝中的填充流体。
根据方程(15)可知,零阶扩展方位傅里叶系数R0(χ)与地层系数g、法向弱度和切向弱度均有关。
在本实施例2中,以某工区实际井为例,分析了零阶扩展方位傅里叶系数R0(χ)与裂缝流体因子之间的关系。
如图3展示了测井曲线,包括纵、横波速度,密度,法向弱度及切向弱度。利用方程(15)合成不同旋转角情况下的零阶扩展方位傅里叶系数,并计算其与裂缝流体因子的相关系数。
图4展示了不同旋转角情况下合成的零阶扩展方位傅里叶系数。可以看出零阶扩展方位傅里叶系数随着旋转角χ的变化而变化。图5展示了不同旋转角情况下的零阶扩展方位傅里叶系数与裂缝流体因子的相关系数。可以看出零阶扩展方位傅里叶系数与裂缝流体因子的相关系数(点线)最高达到0.99,因此,零阶扩展方位傅里叶系数可用于估测裂缝流体因子。
在本实施例2中,选取图3所示井曲线验证提出方法的有效性。利用井曲线合成四个方位的角道集数据如图6(a)所示,为了验证方法的抗噪性,向合成的方位角道集中添加信噪比为2的高斯噪声,如图6(b),以模拟观测地震数据。图7和图8分别为无噪声和信噪比为2情况下反演的各向异性截距阻抗和梯度阻抗,其中黑色实线表示真实值,点线表示初始模型,虚线表示反演结果。可以看出无噪声情况下反演的各向异性截距阻抗和梯度阻抗与真实值完全重合;即使在信噪比为2情况下,反演的各向异性截距阻抗和梯度阻抗也与真实值较为一致。图9展示了无噪声和信噪比为2情况下预测的裂缝流体因子。可以看出,无噪声情况下预测的裂缝流体因子与真实值完全重合,信噪比为2情况下预测的裂缝流体因子也与真实值基本一致。
实施例3
本实施例3提供一种非暂态计算机可读存储介质,所述非暂态计算机可读存储介质包括用于执行如上所述的裂缝流体因子预测方法方法的指令,所述的方法包括如下步骤:
计算各向异性截距阻抗和梯度阻抗数据体,计算零阶扩展方位傅里叶系数,确定预测流体因子的最佳方位角,通过线性拟合得到零阶扩展傅里叶系数与裂缝流体因子的关系式,结合零阶扩展傅里叶系数与裂缝流体因子的关系式,得到裂缝流体因子。
实施例4
本实施例4提供一种电子设备,包括如实施例3所述的非暂态计算机可读存储介质;以及能够执行所述非暂态计算机可读存储介质的所述指令的一个或多个处理器。
综上所述,本发明实施例所述的裂缝流体因子预测方法及系统,首先推导扩展方位弹性阻抗方程及其扩展方位傅里叶系数表达式,通过获取裂缝流体因子与零阶扩展方位傅里叶系数的相关系数确定最优的旋转角,利用反演得到的各向异性截距阻抗和梯度阻抗及最优的方位角计算得到零阶扩展方位傅里叶系数,不需要估测裂缝弱度,可以避免裂缝弱度反演不准确引起的累计误差,即可实现裂缝流体因子的准确估测,该方法计算简单,精度较高。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本公开的优选实施例,并不用于限制本公开,对于本领域的技术人员来说,本公开可以有各种更改和变化。凡在本公开的精神和原则之内所作的任何修改、等同替换、改进等,均应包含在本公开的保护范围之内。
上述虽然结合附图对本公开的具体实施方式进行了描述,但并非对本公开保护范围的限制,所属领域技术人员应该明白,在本发明公开的技术方案的基础上,本领域技术人员在不需要付出创造性劳动即可做出的各种修改或变形,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种裂缝流体因子预测方法,其特征在于,包括:
计算各向异性截距阻抗和梯度阻抗数据体,计算零阶扩展方位傅里叶系数,确定预测流体因子的最佳旋转角,通过线性拟合得到零阶扩展傅里叶系数与裂缝流体因子的关系式,结合零阶扩展傅里叶系数与裂缝流体因子的关系式,得到裂缝流体因子。
2.根据权利要求1所述的裂缝流体因子预测方法,其特征在于,基于方位部分角度叠加地震数据、方位角度地震子波以及测井曲线插值得到各向异性截距阻抗和梯度阻抗低频模型,通过分方位叠前地震贝叶斯反演得到各向异性截距阻抗和梯度阻抗数据体。
3.根据权利要求2所述的裂缝流体因子预测方法,其特征在于,利用各向异性截距阻抗和梯度阻抗反演结果,计算不同旋转角下的扩展方位弹性阻抗,对扩展方位弹性阻抗进行傅里叶级数展开,得到零阶扩展方位傅里叶系数。
4.根据权利要求1所述的裂缝流体因子预测方法,其特征在于,计算零阶扩展方位傅里叶系数与裂缝流体因子的相关系数,确定预测流体因子的最佳的旋转角。
5.根据权利要求4所述的裂缝流体因子预测方法,其特征在于,结合各向异性截距阻抗和梯度阻抗数据体以及最佳旋转角,计算得到扩展方位弹性阻抗数据体,对扩展方位弹性阻抗数据体进行傅里叶级数展开,得到零阶扩展方位傅里叶系数数据体,利用零阶扩展方位傅里叶系数与裂缝流体因子的关系式,将零阶扩展方位傅里叶系数数据体转化为裂缝流体因子。
6.根据权利要求2所述的裂缝流体因子预测方法,其特征在于,利用测井曲线插值得到各向异性截距阻抗和梯度阻抗低频模型。
7.一种基于如权利要求1-6任一项所述的裂缝流体因子预测方法的裂缝流体因子预测系统,其特征在于,包括:
第一计算模块,用于计算各向异性截距阻抗和梯度阻抗数据体;
第二计算模块,用于计算零阶扩展方位傅里叶系数;
确定模块,用于确定预测流体因子的最佳旋转角;
拟合模块,用于将零阶扩展方位傅里叶系数与裂缝流体因子进行线性拟合,得到零阶扩展傅里叶系数与裂缝流体因子的关系式;
转化模块,用于结合各向异性截距阻抗和梯度阻抗数据体、最佳旋转角以及零阶扩展傅里叶系数与裂缝流体因子的关系式,将零阶扩展方位傅里叶系数转化为裂缝流体因子。
8.一种裂缝流体因子预测装置,其特征在于,包括如权利要求7所述的裂缝流体因子预测系统。
9.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质包括用于执行如权利要求1-6任一项所述的裂缝流体因子预测方法的指令。
10.一种电子设备,其特征在于,包括如权利要求9所述的非暂态计算机可读存储介质;以及能够执行所述非暂态计算机可读存储介质的所述指令的一个或多个处理器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110271299.4A CN112835101B (zh) | 2021-03-13 | 2021-03-13 | 裂缝流体因子预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110271299.4A CN112835101B (zh) | 2021-03-13 | 2021-03-13 | 裂缝流体因子预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112835101A true CN112835101A (zh) | 2021-05-25 |
CN112835101B CN112835101B (zh) | 2022-04-15 |
Family
ID=75930089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110271299.4A Active CN112835101B (zh) | 2021-03-13 | 2021-03-13 | 裂缝流体因子预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112835101B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113740896A (zh) * | 2021-07-30 | 2021-12-03 | 中国石油大学(华东) | 一种预测正交各向异性储层裂缝的方法及存储介质 |
CN113835119A (zh) * | 2021-09-18 | 2021-12-24 | 中国石油大学(华东) | 基于方位扩展弹性阻抗的页岩储层水平地应力预测方法 |
CN114002738A (zh) * | 2021-11-02 | 2022-02-01 | 中海石油(中国)有限公司 | 高角度断裂约束的叠前裂缝预测方法、装置及存储介质 |
CN114019565A (zh) * | 2021-11-02 | 2022-02-08 | 中海石油(中国)有限公司 | 快速叠前裂缝预测方法、装置、计算机存储介质及设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100042325A1 (en) * | 2008-08-18 | 2010-02-18 | Beasley Craig J | Determining characteristics of a subterranean body using pressure data and seismic data |
CN103675913A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 一种裂缝型储层检测方法 |
CN106290105A (zh) * | 2016-07-20 | 2017-01-04 | 中国石油大学(华东) | 一种碳酸盐岩储层溶蚀孔隙体积含量预测方法 |
CN109425893A (zh) * | 2017-08-23 | 2019-03-05 | 中国石油化工股份有限公司 | 一种碳酸盐缝洞体系发育带预测方法及装置 |
CN111208580A (zh) * | 2020-03-16 | 2020-05-29 | 中国石油大学(北京) | 确定页岩气储层非均质性的方法和装置 |
-
2021
- 2021-03-13 CN CN202110271299.4A patent/CN112835101B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100042325A1 (en) * | 2008-08-18 | 2010-02-18 | Beasley Craig J | Determining characteristics of a subterranean body using pressure data and seismic data |
CN103675913A (zh) * | 2012-09-04 | 2014-03-26 | 中国石油天然气集团公司 | 一种裂缝型储层检测方法 |
CN106290105A (zh) * | 2016-07-20 | 2017-01-04 | 中国石油大学(华东) | 一种碳酸盐岩储层溶蚀孔隙体积含量预测方法 |
CN109425893A (zh) * | 2017-08-23 | 2019-03-05 | 中国石油化工股份有限公司 | 一种碳酸盐缝洞体系发育带预测方法及装置 |
CN111208580A (zh) * | 2020-03-16 | 2020-05-29 | 中国石油大学(北京) | 确定页岩气储层非均质性的方法和装置 |
Non-Patent Citations (1)
Title |
---|
苗青,等: "塔里木轮古西奥陶系碳酸盐岩岩溶孔洞充填特征识别方法", 《中国石油勘探》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113740896A (zh) * | 2021-07-30 | 2021-12-03 | 中国石油大学(华东) | 一种预测正交各向异性储层裂缝的方法及存储介质 |
CN113740896B (zh) * | 2021-07-30 | 2022-06-03 | 中国石油大学(华东) | 一种预测正交各向异性储层裂缝的方法及存储介质 |
CN113835119A (zh) * | 2021-09-18 | 2021-12-24 | 中国石油大学(华东) | 基于方位扩展弹性阻抗的页岩储层水平地应力预测方法 |
CN114002738A (zh) * | 2021-11-02 | 2022-02-01 | 中海石油(中国)有限公司 | 高角度断裂约束的叠前裂缝预测方法、装置及存储介质 |
CN114019565A (zh) * | 2021-11-02 | 2022-02-08 | 中海石油(中国)有限公司 | 快速叠前裂缝预测方法、装置、计算机存储介质及设备 |
CN114002738B (zh) * | 2021-11-02 | 2023-08-15 | 中海石油(中国)有限公司 | 高角度断裂约束的叠前裂缝预测方法、装置及存储介质 |
CN114019565B (zh) * | 2021-11-02 | 2023-10-13 | 中海石油(中国)有限公司 | 快速叠前裂缝预测方法、装置、计算机存储介质及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN112835101B (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112835101B (zh) | 裂缝流体因子预测方法及系统 | |
CN106842313B (zh) | 基于方位叠前地震数据的各向异性参数反演方法 | |
Downton et al. | Azimuthal Fourier coefficients: A simple method to estimate fracture parameters | |
CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
US9244188B2 (en) | System and method for estimating a nuclear magnetic resonance relaxation time cutoff | |
CN104252007B (zh) | 一种相容性岩石物理建模方法 | |
CN111208560B (zh) | 正交介质裂缝型储层水平裂缝及垂直裂缝同步预测方法 | |
CN110133720A (zh) | 一种基于统计岩石物理模型的横波速度预测方法 | |
CN113312583A (zh) | 裂缝性页岩气储层性能参数预测方法及系统 | |
CN112230276B (zh) | 裂缝型致密储层流体识别方法、系统、识别仪、介质及应用 | |
Jones et al. | A High-resolution, Three-dimensional Model of the M2, M4M6, S2, N2, K1and O1Tides in the Eastern Irish Sea | |
CN104101904B (zh) | 一种快速求取地层横波速度的方法 | |
CN109655909B (zh) | 计算页岩全角度各向异性脆性指数的方法及系统 | |
CN103115625A (zh) | 一种浮体横纵荡及升沉运动的测量方法及系统 | |
EP3428655A1 (en) | Meter for ship speed relative to water and method for correcting measurement values of meter for ship speed relative to water | |
CN110320575B (zh) | 基于岩石物理模型的页岩有机质含量确定方法及装置 | |
CN115993649A (zh) | 基于等效方位杨氏模量的裂缝参数预测方法及系统 | |
CN113552624A (zh) | 孔隙度预测方法及装置 | |
CN110954949A (zh) | 一种致密砂岩软孔隙度分布反演方法 | |
CN105988136A (zh) | 一种利用纵横波速度信息进行含气性分析的方法 | |
CN112379437B (zh) | 页岩储层各向异性参数求取方法及装置 | |
CN106249293B (zh) | 一种各向异性流体识别因子反演方法及系统 | |
CN112346130A (zh) | 一种富有机质岩石横波速度预测方法、存储介质及系统 | |
CN116381795B (zh) | 一种含气储层的声波时差校正方法和装置 | |
CN111830562A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |