CN104407378B - 一种各向异性参数反演方法及装置 - Google Patents
一种各向异性参数反演方法及装置 Download PDFInfo
- Publication number
- CN104407378B CN104407378B CN201410686071.1A CN201410686071A CN104407378B CN 104407378 B CN104407378 B CN 104407378B CN 201410686071 A CN201410686071 A CN 201410686071A CN 104407378 B CN104407378 B CN 104407378B
- Authority
- CN
- China
- Prior art keywords
- delta
- overbar
- anisotropic parameters
- wave
- anisotropic
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 208000010392 Bone Fractures Diseases 0.000 claims abstract description 75
- 206010017076 Fracture Diseases 0.000 claims abstract description 75
- 238000013508 migration Methods 0.000 claims abstract description 9
- 230000005012 migration Effects 0.000 claims abstract description 9
- 238000010008 shearing Methods 0.000 claims description 24
- 238000010276 construction Methods 0.000 claims description 20
- 238000001615 p wave Methods 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000013016 damping Methods 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 238000010304 firing Methods 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 229920000535 Tan II Polymers 0.000 claims 4
- 238000009825 accumulation Methods 0.000 claims 2
- 239000007789 gas Substances 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 101100234408 Danio rerio kif7 gene Proteins 0.000 description 1
- 101100221620 Drosophila melanogaster cos gene Proteins 0.000 description 1
- 101100398237 Xenopus tropicalis kif11 gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 101150118300 cos gene Proteins 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请公开了一种各向异性参数反演方法及装置,该反演方法包括以下步骤:将获取的地震纵波资料数据进行各向异性偏移处理,生成分方位的角道集;利用所生成的角道集,考虑大于30度的入射角对反射振幅的影响,构建裂缝介质的各向异性参数的反演模型;通过对所构建的反演模型进行计算来获取地下储层的各向异性参数;利用所获取的各向异性参数来分析目标储层的裂缝发育特征。本申请公开的反演装置相应的包括角道集生成单元、反演模型构建单元、各向异性参数获取单元以及裂缝发育特征分析单元。通过本申请中的各向异性参数反演方法及装置能提高各向异性参数的反演精度,实现准确预测地下介质裂缝发育特征的目的。
Description
技术领域
本申请涉及石油勘探领域,尤其涉及一种各向异性参数反演方法及装置。
背景技术
油气在裂缝型圈闭中聚集而成的油气藏称为裂缝型油气藏,裂缝型油气藏在世界石油和天然气的产量、储量中占有十分重要的地位,裂缝型油气藏的产量占全世界石油天然气总产量的一半以上。在国内的油气勘探中,裂缝型油气藏是一个重要的勘探领域,在勘探裂缝型油气藏时,最重要的是分析和认识裂缝发育特征。储层内发育裂缝时对外表现为各向异性特征,所以用各向异性参数来分析裂缝发育特征。
随着采集处理技术的发展和各向异性理论的完善,裂缝预测技术呈现多样化,由叠后走到叠前,预测精度逐步提高。对于垂直排列的平行裂缝,可以将其等效为HTI介质(Horizontal Transverse Isotropy,即具有水平对称轴的横向各向同性介质)模型来研究。若入射角为θ,测线方位角为φ,Ruger认为HTI介质分界面处的纵波反射系数是各向异性参数、入射角和方位角的函数:
上式中:φs为裂缝平面的对称轴方向,即裂缝面法向,(φ-φs)为测线与裂缝面法向的夹角;α,β,Z,G分别是纵波速度、横波速度、纵波阻抗和剪切模量;εV,δV,γ为各向异性参数。分别表示上下两层介质间纵波速度的平均值、横波速度的平均值、纵波阻抗的平均值和剪切模量的平均值,Δα,ΔZ,ΔG,ΔεV,ΔδV,Δγ分别表示上下两层介质间纵波速度的差值、纵波阻抗的差值、剪切模量的差值以及三个各向异性参数的差值。
在石油工业领域,目前常用下面公式(2)进行AVAZ反演;用各向异性AVO梯度Bani(也被称作拟各向异性参数)描述裂缝发育强度;用参数φs表示裂缝面法向。
R(θ,φ)≈R0+[Biso+Bani cos2(φ-φs)]sin2θ (2)
其中:
将公式(2)表示裂缝面法向的φs从余弦项中展开,其转化为:
R(θ,φ)=a1+(sin2θ)a2+(cos2φsin2θ)a3+(sin2φsin2θ)a4 (3)
其中:a1=R0,
依据公式(3)构建反演方程组,解方程组得XI=(a1,a2,a3,a4)T,再由公式(4)计算各向异性梯度Bani和裂缝面法向φs等参数:
R0=a1 R0=a1
或
开平方运算过程中正负号的选取,使得由公式(4)计算的各向异性梯度Bani存在二义性,进而也影响裂缝面法向φs的取值。由此计算的各向异性组合参数存在较大误差,不能准确地刻画裂缝发育特征。
发明内容
本申请通过提供一种各向异性参数反演方法及装置来克服常规各向异性参数反演的不足,以准确预测地下介质裂缝发育特征。
为达到上述技术目的,本申请通过以下技术方案实现:
本申请提供了一种各向异性参数反演方法,包括如下步骤:
将获取的地震纵波资料数据进行各向异性偏移处理,生成分方位的角道集,所述角道集包含有地震纵波的入射角、方位角和反射振幅的相关信息;
利用所述角道集,利用含有入射角三角函数四次项的Ruger公式,考虑大于30度的入射角对反射振幅的影响,构建裂缝介质的各向异性参数的反演模型;
对所构建的反演模型进行计算,同时获取地下储层的各向异性参数εV,δV,γ,其中εV表示裂缝介质的纵波各向异性程度,δV表示纵波在横向和垂向之间的各向异性变化的快慢程度,γ表示横波各向异性程度;
利用所获取的各向异性参数εV,δV,γ来分析目标储层的裂缝发育特征。
优选的,所述构建裂缝介质的各向异性参数的反演模型包括:
获取裂缝介质的纵横波速度比以及裂缝平面的对称轴方向φs,并对纵波阻抗和剪切模量进行近似分解;
利用所获取的纵横波速度比裂缝平面的对称轴方向φs以及分解后的纵波阻抗和剪切模量来简化地震纵波的反射振幅R与入射角θ、方位角φ以及各向异性参数之间的关系式,简化后的关系式如下:
其中A,B,C,D,E,F是入射角θ和方位角φ的组合:
φ'=φ-φs;
利用简化后的地震纵波的反射振幅与入射角、方位角以及各向异性参数之间的关系式,构建各向异性参数的反演方程组。
优选的,所述对纵波阻抗和剪切模量进行近似分解包括:
利用下式对纵波阻抗的相对比值和剪切模量的相对比值进行近似分解:
其中,分别表示上下两层介质间纵波速度的平均值、横波速度的平均值以及密度的平均值,Δα,Δβ,Δρ分别表示上下两层介质间纵波速度的差值、横波速度的差值以及密度的差值。
优选的,所述各向异性参数的反演方程组包括下述超定方程组:
Mn×6X6×1=Rn×1
上式中n为覆盖次数,M,X和R分别为:
优选的,所述对所构建的反演模型进行计算来获取地下储层的各向异性参数包括:
设定阻尼因子,然后利用奇异值分解法求解所述超定方程组,求解的结果为:
X=M+R
其中M+为矩阵M的广义逆;
通过求解得到的X来获取地下储层的各向异性参数。
本申请还提供了一种各向异性参数的反演装置,包括:
角道集生成单元,所述角道集生成单元用于将获取的地震纵波资料数据进行各向异性偏移处理,生成分方位的角道集,所述角道集包含有地震纵波的入射角、方位角和反射振幅的相关信息;
反演模型构建单元,所述反演模型构建单元用于利用所述角道集,考虑大于30度的入射角对反射振幅的影响,构建裂缝介质的各向异性参数的反演模型;
各向异性参数获取单元,所述各向异性参数获取单元用于对所构建的反演模型进行计算,同时获取地下储层的各向异性参数εV,δV,γ,其中εV表示裂缝介质的纵波各向异性程度,δV表示纵波在横向和垂向之间的各向异性变化的快慢程度,γ表示横波各向异性程度;
裂缝发育特征分析单元,裂缝发育特征分析单元用于利用获取的各向异性参数εV,δV,γ来分析目标储层的裂缝发育特征。
优选的,所述反演模型构建单元包括:
预处理单元,所述预处理单元用于获取裂缝介质的纵横波速度比以及裂缝平面的对称轴方向φs,并对纵波阻抗和剪切模量进行近似分解;
关系式简化单元,所述关系式简化单元用于利用所获取的来纵横波速度比裂缝平面的对称轴方向φs以及分解后的纵波阻抗和剪切模量来简化地震纵波的反射振幅与入射角、方位角以及各向异性参数之间的关系式,简化后的关系式如下:
其中A,B,C,D,E,F是入射角θ和方位角φ的组合:
φ'=φ-φs;
反演方程组构建单元,所述反演方程组构建单元用于利用简化后的地震纵波的反射振幅与入射角、方位角以及各向异性参数之间的关系式,构建各向异性参数的反演方程组。
优选的,所述预处理单元包括近似分解单元,所述近似分解单元用于利用下式对纵波阻抗的相对比值和剪切模量的相对比值进行近似分解:
其中,分别表示上下两层介质间纵波速度的平均值、横波速度的平均值以及密度的平均值,Δα,Δβ,Δρ分别表示上下两层介质间纵波速度的差值、横波速度的差值以及密度的差值。
优选的,所述反演方程组构建单元包括超定方程组构建单元,所述超定方程组构建单元用于构建下述超定方程组:
Mn×6X6×1=Rn×1
上式中n为覆盖次数,M,X和R分别为:
优选的,各向异性参数获取单元包括求解单元,所述求解单元用于通过设定阻尼因子并利用奇异值分解法来求解所述超定方程组,求解的结果为:
X=M+R
其中M+为矩阵M的广义逆;通过求解得到的X来获取地下储层的各向异性参数。
与现有技术相比,本申请具有下有益效果:
通过考虑公式(1)中有关入射角θ的四次项sin2θtan2θ对反射振幅的影响以及利用高精度的反演来减小各向异性参数的相对误差,从而实现准确预测地下介质裂缝发育特征的目的。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请中各向异性参数反演方法的流程示意图。
图2为地震波射线在裂缝介质中的传播示意图。
图3为本申请实施例一中各向异性参数反演方法的流程示意图。
图4为反射振幅的等值线图。
图5为sin2θ和sin2θtan2θ与入射角θ的关系曲线。
图6为不同入射角范围的各向异性参数反演结果的相对误差柱状图。
图7为不同方位角范围的各向异性参数反演结果相对误差柱状图。
图8为不同信噪比的各向异性参数反演结果相对误差柱状图。
图9为各向异性参数与裂缝密度的关系曲线。
图10为本申请实施例二中各向异性参数反演装置的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
本申请提供一种各向异性参数反演方法,如图1所示,该反演方法包括如下步骤:
步骤S1,将获取的地震纵波资料数据进行各向异性偏移处理,生成分方位的角道集数据,所述角道集数据包含有地震纵波的入射角、方位角和反射振幅的相关信息;
步骤S2,利用所述角道集数据,考虑大于30度的入射角对反射振幅的影响,构建裂缝介质的各向异性参数的反演模型;
步骤S3,对所构建的反演模型进行计算来获取地下储层的各向异性参数;
步骤S4,利用获取的各向异性参数来分析目标储层的裂缝发育特征。
与现有的各向异性参数反演方法相比,通过本申请中的技术方案能提高各向异性参数的反演精度,以达到准确分析裂缝发育特征的目的。
下面结合附图,对本发明的具体实施方式作进一步的详细说明。
实施例一
图2为地震波射线在裂缝介质中的传播示意图,图中入射角θ为地震波射线与反射面法向的夹角,方位角φ为地震波的激发点和接收点的连线(即测线)与x轴的夹角。储层内发育裂缝时表现为各向异性特征,其反射系数(即反射振幅)不仅随入射角变化,也随方位角变化。本申请利用反射振幅R与入射角θ、方位角φ以及三个各向异性参数之间的关系,通过高精度的AVAZ(Amplitude variation with Angle and Azimuth振幅随入射角、方位角变化)反演来获取地下储层的各向异性参数,以准确分析目标储层的裂缝发育特征。
如图3所示,本申请实施例所提供的各向异性参数反演方法包括以下步骤:
步骤S01,将获取的地震纵波资料数据进行各向异性偏移处理,生成分方位的角道集。
对所获取的整个研究区域的三维地震纵波资料数据进行各向异性偏移处理,生成分方位的角道集,例如角道集1,角道集2……角道集n,这些角道集中含有入射角、方位角以及反射振幅的相关信息,其中入射角θ、方位角φ以及反射振幅R三者之间具有一定的关系,如图4所示,该图为反射振幅的等值线图,横轴表示方位角的大小,纵轴表示入射角的大小,图中曲线上的数值表示反射振幅的大小,当地震波沿着不同的射线传播时,振幅不仅随着入射角变化,还随着方位角而变化。例如,当方位角φ为0°时,随着入射角θ的增大,反射振幅R增大;当方位角φ为90°时,随着入射角θ的增大,反射振幅R减小;对于同一个入射角θ,反射振幅R关于90°的方位角对称。
步骤S102,利用反射振幅R与入射角θ、方位角φ以及三个各向异性参数之间的关系,考虑大偏移距对反射振幅R的影响,构建裂缝介质的各向异性参数的反演模型。
各向异性裂缝介质模型中含有一系列垂直排列且相互平行的裂缝。对于垂直排列的平行裂缝,可以将其等效为HTI介质(Horizontal Transverse Isotropy,即具有水平对称轴的横向各向同性介质)模型来研究。Ruger认为HTI介质分界面处的纵波反射系数(即反射振幅)是各向异性参数、入射角和方位角的函数:
上式中:R为反射振幅,θ为入射角、φ为方位角,φs为裂缝平面的对称轴方向,即裂缝面法向,(φ-φs)为测线与裂缝面法向的夹角;α,β,Z,G分别是纵波速度、横波速度、纵波阻抗和剪切模量;εV,δV,γ为各向异性参数。分别表示上下两层介质间纵波速度的平均值、横波速度的平均值、纵波阻抗的平均值和剪切模量的平均值,Δα,ΔZ,ΔG,ΔεV,ΔδV,Δγ分别表示上下两层介质间纵波速度的差值、纵波阻抗的差值、剪切模量的差值以及三个各向异性参数的差值。
Ruger反射系数公式(1)包含入射角三角函数的二次项sin2θ和四次项sin2θtan2θ,两者取值大小随入射角θ变化,如图5所示。从图中可以看出,当θ≤30°时,sin2θtan2θ≤1/12<<1,此时可忽略该四次项对反射振幅R的影响;而当30°≤θ≤45°时,sin2θtan2θ随着入射角θ的增大而迅速增大,此时不能忽略该四次项对反射振幅R的影响;当θ>45°时,sin2θtan2θ>sin2θ>1/2,该四次项对反射振幅的影响大过二次项,更不能忽略。
所以为了提高反演精度,考虑大偏移距(对应于大于30度的入射角)对反射振幅R的影响,本申请采用保留Ruger公式中四次项sin2θtan2θ的高精度的AVAZ算法进行各向异性参数的反演,构建裂缝介质的各向异性参数的反演模型,具体步骤如下:
(1)利用地震波叠后波阻抗或者叠前AVO(Amplitude Versus Offset振幅随偏移距的变化)反演,计算裂缝介质的纵横波速度比
其中利用叠后波阻抗反演主要是利用道积分方法来进行反演,得到纵横波阻抗,再计算裂缝介质的纵横波速度比叠前AVO反演主要是利用Shuey近似方程进行反演,该方程的具体形式如下:
其中:
通过对式(5)进行反演可以得到参数A,B,C的值,再通过对下述方程组进行逆转换并求解便可以得到纵横波速度α,β,进而便可以得到裂缝介质的纵横波速度比
其中:
上述求解中一般是先给f设定一初值,然后通过迭代求解出纵横波速度α,β。上述是通过AVO三参数反演来计算出纵横波速度比但不限于这一种方式,也可以通过裂缝介质的泊松比或拉梅常数等来计算出裂缝介质的纵横波速度比
(2)利用不同方位的地震属性,通过AVD(Amplitude Versus Direction振幅随方向变化)方法反演裂缝平面的对称轴方向φs。
在各向异性介质中,入射角固定时,振幅随方位角呈余弦变化,动校正速度、走时和AVO梯度等属性随测线方位角的变化也可近似为余弦函数
F(φ)=CAVD+GAVD cos2(φ-φs) (7)
上式中:F为地震属性,CAVD为截距,GAVD为梯度,φs为裂缝法向的方位角,即裂缝平面的对称轴方向,φ为测线方位角,(φ-φs)为测线与裂缝面法向的夹角。
利用式(7)便可以反演计算得到裂缝平面的对称轴方向φs。
(3)将纵波阻抗和剪切模量进行近似分解,主要是将纵波阻抗的相对比值和剪切模量的相对比值近用纵波速度、横波速度以及密度这三个物理量来近似表达,即:
(8)
(4)利用式(8)以及式(6)和式(7)中计算得到的纵横波速度比裂缝平面的对称轴方向φs对式(1)进行化简,将反射振幅R化简为参数向量的函数,化简后的Ruger公式如下:
其中A,B,C,D,E,F是入射角θ和方位角φ的组合:
(10)
(5)利用式(9)构建各向异性参数的反演方程组,在本实施例中所构建的反演方程组为超定方程组,如下式所示:
Mn×6X6×1=Rn×1 (11)
上式中n为覆盖次数,M,X和R分别为:
步骤S103,利用步骤S102中所构建的反演模型,来获取地下储层的各向异性参数。
在本实施例中,该步骤主要是指利用超定方程组(11)来求解地下储层的各向异性参数。具体求解过程如下:
当M是列满秩矩阵时,超定方程组MX=R存在唯一解X=M+R,其中M+为矩阵M的广义逆。此时,可通过奇异值分解算法计算广义逆M+,计算过程中设置阻尼因子μ等于矩阵M奇异值总和的10%,该值只是为经验值,阻尼因子μ并不限于为此值,其取值随奇异值总和的改变而改变。
若M的奇异值分解为:
其中U为n×n阶矩阵;VT为矩阵V的共轭转置,为6×6阶矩阵;r为奇异值的个数,Λr为对角线矩阵,对角线上的元素即为M的奇异值,其从大到小排列为奇异值σ1,σ2,...,σp,σp+1,...σr,那么广义逆M+为:
计算中当出现奇异值σp+1<μ时,可将式中σp+1以后的倒数项置为零,这样得到的广义逆M+如式(14)所示:
最终可求得超定方程组(11)的解为:
由于M+MX=M+R,因此称M+M为分辨率矩阵,该分辨率矩阵越接近单位矩阵,所求解就越准确。
从式(16)中便可以得到与地下储层的各向异性参数εV,δV,γ相关的三个参数ΔεV,ΔδV,Δγ以及其他三个参数的取值,反演的求解结果亦可以用于各向同性梯度Biso和各向异性梯度Bani的计算。
上述反演求解的结果为各向异性参数的差值,本申请中利用以下转换关系式,对各向异性参数的差值这三个反演参数进行累计求和便可以得到地下储层的各向异性参数εV,δV,γ。这三个各向异性参数都具有明确的物理意义,均可描述裂缝发育特征,εV表示裂缝介质的纵波各向异性程度,δV表示纵波在横向和垂向之间的各向异性变化的快慢程度,γ表示横波各向异性程度。
其中j=1,2,...i;i=1,2,...imax
以下以不同入射角范围、方位角范围和信噪比的反演效果为例来说明本申请所提供的反演方法的有效性。
图6为不同入射角范围的各向异性参数反演结果的相对误差柱状图。其中,最小入射角固定为5°,最大入射角分别为30°,35°,40°和45°。横轴为三个各向异性参数的相关参数(ΔεV,ΔδV,Δγ),纵轴为相对误差,从图中可以看出,对于这四个入射角范围,三个各向异性参数的相关参数的相对误差很小,接近为0。
图7为不同方位角范围的各向异性参数反演结果相对误差柱状图。图中所示的是用方位角范围为30°的数据来模拟窄方位角道集,角道集的方位角范围分别为[0°,30°],[30°,60°],[60°,90°]。图中横轴为三个反演参数(ΔεV,ΔδV,Δγ),纵轴为相对误差。从图中可以看出,每个窄方位所得的参数Δγ的相对误差均小于5%,比ΔεV和ΔδV稳定。窄方位[30°,60°]所得到的三个反演参数的相对误差均小于9%,比另外两个窄方位稳定。
考虑到信噪比是描述角道集数据品质的重要参数,现考虑信噪比为60db、40db以及没有噪音的情形,各向异性参数的反演结果相对误差如图8所示。该图的横轴为三个反演参数(ΔεV,ΔδV,Δγ),纵轴为相对误差,参数ΔεV和Δγ的相对误差均小于1%,参数ΔδV的相对误差小于6%。
综上所述,通过本申请所提供的各向异性参数反演方法能提高各向异性参数的反演精度,从而能准确的预测地下介质中裂缝的发育特征。
步骤S104,利用步骤S103中所获得的各向异性参数,来分析研究区域中目标储层的裂缝发育特征。
图9为各向异性参数与裂缝密度的关系曲线,从图中可以看出,各向异性参数与裂缝密度呈线性关系,但该曲线的斜率对于不同流体和背景介质会有所差异。针对所研究区域,如果所获得的各向异性参数越大,该地区的裂缝就发育得越好。
实施例二
本申请还提供了一种各向异性参数的反演装置,如图10所示,该装置包括角道集生成单元901、反演模型构建单元902、各向异性参数获取单元903以及裂缝发育特征分析单元904。其中角道集生成单元901用于将获取的地震纵波资料数据进行各向异性偏移处理,生成分方位的角道集,该角道集包含有地震纵波的入射角、方位角和反射振幅的相关信息。
反演模型构建单元902用于利用所生成的角道集,并考虑大于30度的入射角对反射振幅的影响,构建裂缝介质的各向异性参数的反演模型,反演模型构建单元902包括三个子单元,分别为预处理单元9021、关系式简化单元9022以及反演方程组构建单元9023。预处理单元9021用于获取裂缝介质的纵横波速度比以及裂缝平面的对称轴方向φs,并对纵波阻抗和剪切模量进行近似分解,其包含有近似分解单元;关系式简化单元9022用于利用从预处理单元9021中所获取的来纵横波速度比裂缝平面的对称轴方向φs以及分解后的纵波阻抗和剪切模量来简化地震纵波的反射振幅与入射角、方位角以及各向异性参数之间的关系式,简化后的关系式为式(9);反演方程组构建单元9023用于利用式(9),构建各向异性参数的反演方程组,在本实施例中为超定方程组(11)。
各向异性参数获取单元903用于利用反演模型构建单元902中所构建的反演模型,通过设定阻尼因子并对该反演模型进行计算来获取地下储层的各向异性参数;
裂缝发育特征分析单元904用于利用所获取的各向异性参数来分析目标储层的裂缝发育特征。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (10)
1.一种各向异性参数反演方法,其特征在于,包括如下步骤:
将获取的地震纵波资料数据进行各向异性偏移处理,生成分方位的角道集,所述角道集包含有地震纵波的入射角、方位角和反射振幅的相关信息;
利用所述角道集,考虑大于30度的入射角对反射振幅的影响,构建裂缝介质的各向异性参数的反演模型;
对所构建的反演模型进行反演计算以获取地下储层中上下两层介质间的各向异性参数εV,δV,γ的差值;
对所获取的所述各向异性参数εV,δV,γ的差值分别进行累积求和计算,得到所述各向异性参数εV,δV,γ,其中εV表示裂缝介质的纵波各向异性程度,δV表示纵波在横向和垂向之间的各向异性变化的快慢程度,γ表示横波各向异性程度;
利用所得到的各向异性参数εV,δV,γ来分析目标储层的裂缝发育特征。
2.根据权利要求1所述的各向异性参数反演方法,其特征在于,所述构建裂缝介质的各向异性参数的反演模型包括:
获取裂缝介质的纵横波速度比以及裂缝平面的对称轴方向φs,并对纵波阻抗和剪切模量进行近似分解;
利用所获取的纵横波速度比裂缝平面的对称轴方向φs以及分解后的纵波阻抗和剪切模量来简化地震纵波的反射振幅R与入射角θ、方位角φ以及各向异性参数之间的关系式,简化后的关系式如下:
其中A,B,C,D,E,F是入射角θ和方位角φ的组合:
φ'=φ-φs;
利用简化后的地震纵波的反射振幅与入射角、方位角以及各向异性参数之间的关系式,构建各向异性参数的反演方程组。
3.根据权利要求2所述的各向异性参数反演方法,其特征在于,所述对纵波阻抗和剪切模量进行近似分解包括:
利用下式对纵波阻抗的相对比值和剪切模量的相对比值进行近似分解:
其中,分别表示上下两层介质间纵波速度的平均值、横波速度的平均值以及密度的平均值,Δα,Δβ,Δρ分别表示上下两层介质间纵波速度的差值、横波速度的差值以及密度的差值。
4.根据权利要求2或3所述的各向异性参数反演方法,其特征在于,所述各向异性参数的反演方程组包括下述超定方程组:
Mn×6X6×1=Rn×1
上式中n为覆盖次数,M,X和R分别为:
5.根据权利要求4所述的各向异性参数反演方法,其特征在于,所述对所构建的反演模型进行反演计算以获取地下储层中上下两层介质间的各向异性参数εV,δV,γ的差值包括:
设定阻尼因子,然后利用奇异值分解法求解所述超定方程组,求解的结果为:
X=M+R
其中M+为矩阵M的广义逆;
从求解得到的X中获取所述各向异性参数εV,δV,γ的差值。
6.一种各向异性参数的反演装置,其特征在于,包括:
角道集生成单元,所述角道集生成单元用于将获取的地震纵波资料数据进行各向异性偏移处理,生成分方位的角道集,所述角道集包含有地震纵波的入射角、方位角和反射振幅的相关信息;
反演模型构建单元,所述反演模型构建单元用于利用所述角道集,考虑大于30度的入射角对反射振幅的影响,构建裂缝介质的各向异性参数的反演模型;
各向异性参数获取单元,所述各向异性参数获取单元用于对所构建的反演模型进行反演计算以获取地下储层中上下两层介质间的各向异性参数εV,δV,γ的差值,对所获取的所述各向异性参数εV,δV,γ的差值分别进行累积求和计算,得到所述各向异性参数εV,δV,γ,其中εV表示裂缝介质的纵波各向异性程度,δV表示纵波在横向和垂向之间的各向异性变化的快慢程度,γ表示横波各向异性程度;
裂缝发育特征分析单元,裂缝发育特征分析单元用于利用所得到的各向异性参数εV,δV,γ来分析目标储层的裂缝发育特征。
7.根据权利要求6所述的各向异性参数的反演装置,其特征在于,所述反演模型构建单元包括:
预处理单元,所述预处理单元用于获取裂缝介质的纵横波速度比以及裂缝平面的对称轴方向φs,并对纵波阻抗和剪切模量进行近似分解;
关系式简化单元,所述关系式简化单元用于利用所获取的来纵横波速度比裂缝平面的对称轴方向φs以及分解后的纵波阻抗和剪切模量来简化地震纵波的反射振幅与入射角、方位角以及各向异性参数之间的关系式,简化后的关系式如下:
其中A,B,C,D,E,F是入射角θ和方位角φ的组合:
φ'=φ-φs;
反演方程组构建单元,所述反演方程组构建单元用于利用简化后的地震纵波的反射振幅与入射角、方位角以及各向异性参数之间的关系式,构建各向异性参数的反演方程组。
8.根据权利要求7所述的各向异性参数的反演装置,其特征在于,所述预处理单元包括近似分解单元,所述近似分解单元用于利用下式对纵波阻抗的相对比值和剪切模量的相对比值进行近似分解:
其中,分别表示上下两层介质间纵波速度的平均值、横波速度的平均值以及密度的平均值,Δα,Δβ,Δρ分别表示上下两层介质间纵波速度的差值、横波速度的差值以及密度的差值。
9.根据权利要求7或8所述的各向异性参数的反演装置,其特征在于,所述反演方程组构建单元包括超定方程组构建单元,所述超定方程组构建单元用于构建下述超定方程组:
Mn×6X6×1=Rn×1
上式中n为覆盖次数,M,X和R分别为:
10.根据权利要求9所述的各向异性参数的反演装置,其特征在于,各向异性参数获取单元包括求解单元,所述求解单元用于通过设定阻尼因子并利用奇异值分解法来求解所述超定方程组,求解的结果为:
X=M+R
其中M+为矩阵M的广义逆;从求解得到的X中获取所述各向异性参数εV,δV,γ的差值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410686071.1A CN104407378B (zh) | 2014-11-25 | 2014-11-25 | 一种各向异性参数反演方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410686071.1A CN104407378B (zh) | 2014-11-25 | 2014-11-25 | 一种各向异性参数反演方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104407378A CN104407378A (zh) | 2015-03-11 |
CN104407378B true CN104407378B (zh) | 2017-05-10 |
Family
ID=52645022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410686071.1A Active CN104407378B (zh) | 2014-11-25 | 2014-11-25 | 一种各向异性参数反演方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104407378B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104820765A (zh) * | 2015-05-26 | 2015-08-05 | 中国科学院声学研究所东海研究站 | 一种海底声学特性的单参数模型的建立方法 |
GB2555365B (en) * | 2015-07-28 | 2021-07-28 | Geoquest Systems Bv | Seismic constrained discrete fracture network |
CN108260359A (zh) * | 2015-07-30 | 2018-07-06 | 中国科学院地质与地球物理研究所 | 一种地震信号处理方法、装置和系统 |
CN107436449A (zh) * | 2016-05-26 | 2017-12-05 | 中国石油化工股份有限公司 | 一种提取水平层状介质各向异性参数的方法 |
CN106528955A (zh) * | 2016-10-18 | 2017-03-22 | 西安交通大学 | 基于阻尼奇异值分解的运行工况传递路径分析方法 |
CN108627871B (zh) * | 2017-03-24 | 2020-02-21 | 中国石油化工股份有限公司 | 一种tti介质裂缝性质参数的反演方法 |
CN108732623B (zh) * | 2017-04-17 | 2020-05-19 | 中国石油化工股份有限公司 | Tti介质层析反演成像方法及计算机可读存储介质 |
CN109143351B (zh) * | 2017-06-28 | 2020-07-24 | 中国石油化工股份有限公司 | 叠前各向异性特征参数反演方法及计算机可读存储介质 |
CN109655889B (zh) * | 2017-10-11 | 2021-04-02 | 中国石油化工股份有限公司 | 一种各向异性参数联合反演方法及系统 |
CN109683197B (zh) * | 2018-12-10 | 2020-05-01 | 中国石油天然气集团有限公司 | 水平裂缝预测方法、装置和设备 |
CN111399050B (zh) * | 2019-01-02 | 2021-04-30 | 中国石油天然气股份有限公司 | 高角度裂缝预测方法及装置 |
CN112444872B (zh) * | 2019-08-30 | 2023-07-18 | 中国石油化工股份有限公司 | 基于正交介质的叠前裂缝预测方法和存储介质 |
CN113126153A (zh) * | 2019-12-30 | 2021-07-16 | 中国石油天然气集团有限公司 | 基于数据组合的叠前深度偏移方法及装置 |
CN113075727B (zh) * | 2021-02-20 | 2023-09-26 | 中国石油天然气股份有限公司 | 基于振幅方位导数的储层裂缝分布特征确定方法及系统 |
CN113120213B (zh) * | 2021-03-31 | 2023-11-28 | 中国飞机强度研究所 | 一种可变形乘波体耐高温柔性蒙皮及其设计方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5583825A (en) * | 1994-09-02 | 1996-12-10 | Exxon Production Research Company | Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data |
US6928367B2 (en) * | 2002-09-27 | 2005-08-09 | Veritas Dgc Inc. | Reservoir fracture characterization |
CN102033242B (zh) * | 2010-10-22 | 2012-09-26 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
-
2014
- 2014-11-25 CN CN201410686071.1A patent/CN104407378B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104407378A (zh) | 2015-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104407378B (zh) | 一种各向异性参数反演方法及装置 | |
Zhubayev et al. | Ultrasonic velocity and attenuation anisotropy of shales, Whitby, United Kingdom | |
Shapiro | Fluid-induced seismicity | |
Guo et al. | Effective elastic properties of rocks with transversely isotropic background permeated by aligned penny‐shaped cracks | |
Hatchell et al. | Rocks under strain: Strain-induced time-lapse time shifts are observed for depleting reservoirs | |
Dutta et al. | Attenuation and dispersion of compressional waves in fluid-filled porous rocks with partial gas saturation (White model)—Part I: Biot theory | |
CN102692645B (zh) | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 | |
Barbosa et al. | Estimation of fracture compliance from attenuation and velocity analysis of full‐waveform sonic log data | |
Picotti et al. | Anisotropy and crystalline fabric of Whillans Ice Stream (West Antarctica) inferred from multicomponent seismic data | |
Huang et al. | 3D-printed rock models: Elastic properties and the effects of penny-shaped inclusions with fluid substitution | |
CN102788994A (zh) | 一种储层裂缝的确定方法 | |
CN103487831B (zh) | Avo地震正演计算方法 | |
Far et al. | Seismic characterization of reservoirs with asymmetric fractures | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
Kang et al. | Sequential approach to joint flow‐seismic inversion for improved characterization of fractured media | |
Guo et al. | Effects of coupling between wave‐induced fluid flow and elastic scattering on P‐wave dispersion and attenuation in rocks with aligned fractures | |
Gao et al. | Comparison of full and quasi‐static seismoelectric analytically based modeling | |
CN106501872B (zh) | 一种裂缝储层地应力特征的计算方法及装置 | |
Ball et al. | A joint M onte C arlo analysis of seafloor compliance, R ayleigh wave dispersion and receiver functions at ocean bottom seismic stations offshore N ew Z ealand | |
Korneev et al. | Seismic velocity changes caused by an overburden stress | |
CN105510965B (zh) | 流体识别方法和装置 | |
CN107728205A (zh) | 一种地层压力预测方法 | |
Robinson et al. | A probabilistic approach for estimating the separation between a pair of earthquakes directly from their coda waves | |
CN104422955B (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
Barbosa et al. | Estimates of individual fracture compliances along boreholes from full‐waveform sonic log data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |