CN101957454A - 基于入射角的avo近似公式及属性提取方法 - Google Patents
基于入射角的avo近似公式及属性提取方法 Download PDFInfo
- Publication number
- CN101957454A CN101957454A CN 201010285247 CN201010285247A CN101957454A CN 101957454 A CN101957454 A CN 101957454A CN 201010285247 CN201010285247 CN 201010285247 CN 201010285247 A CN201010285247 A CN 201010285247A CN 101957454 A CN101957454 A CN 101957454A
- Authority
- CN
- China
- Prior art keywords
- delta
- theta
- rho
- sin
- formula
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
所属技术领域
本发明涉及地震资料处理领域,是一种利用地震数据进行储层识别的方法。
背景技术
AVO技术是根据振幅随炮检距的变化规律所反映出的地下岩性及其孔隙流体的性质来直接预测油气和估计地壳岩性参数的一项技术。其理论基础是描述平面波在水平分界面上反射和折射的Zoeppritz方程。完整的Zoeppritz方程形式复杂,物理意义不明确,无法直接应用于地震数据的处理和参数反演。尽管人们早就认识到了这一点,但在早期的研究中,人们仅从数值计算的角度对Zoeppritz描述的反射系数进行了各种讨论和近似。Muskat系统地计算了弹性介质中平面波的反射和透射系数;Koefoed利用17组不同弹性参数的模型,详细地研究了泊松比对两个各向同性介质之间的反射界面所产生的产生的反射系数的影响,并第一次给出了将泊松比与反射系数直接联系起来的Zoeppritz近似方程。Bortfeld第一个给出了物理意义比较明确的反射系数近似计算公式,使得基于Zoeppritz方程的反射和透射问题真正应用到了勘探地球物理领域。在以上工作的基础上,不少人做了大量的工作,分别从不同的方面对Zoeppritz方程进行了详细的讨论和简化,提出了不同假设条件下的纵波反射振幅的近似表达式。但是所有的近似公式中角度都是应用的上层入射角和下层透射角的平均,给进一步的参数估算带来一定的误差。
发明内容
本发明的目的在于提供基于入射角的AVO近似公式及属性提取方法。
本发明所采用的技术方案是:基于入射角的AVO近似公式及属性提取方法,1)基于入射角的AVO近似公式的提出,2)基于入射角的AVO近似公式的理论推导,3)应用常规方法提取P、G,Rs的属性;4)将属性归一化;5)由井数据或者研究区域的经验公式确定b值,然后应用公式(14)和公式(15)对常规方法提取的属性进行修正;6)将修正结果与井数据进行对比分析,根据分析结果对修正结果进行校正,得到最终的属性结果。
具体的求取方法为:
1)由Aki和Shuey近似可得:
其中,
Vp=(Vp1+Vp2)/2,Vs=(Vs1+Vs2)/2,ρ=(ρ1+ρ2)/2
ΔVp=Vp2-Vp1,ΔVs=Vs2-Vs1,Δρ=ρ2-ρ1,θ=(θ1+θ2)/2
Vp1和Vp2分别是反射界面上下上、下地层的纵波速度,Vs1和Vs2分别是上、下地层的横波速度,θ1和θ2分别是地层的入射角和透射角
2)根据Snell定律有:
对公式(1)左边项进行泰勒展开,及
将公式(3)代入公式(2)可得:
又因为:
所以可得:
当采用一阶近似时,上式可以写成:
令
则在此基础上有
而当θ1和θ2相差不大时,
sinθ1 sinθ2=sin2θ
3)这样可以得到基于入射角的AVO近似公式:
4)Gardner等人于1974提出的密度与纵波速度之间的经验公式的形式为[6]:
其中,a和b为系数(一般取a=0.31,b=0.25)
由上式可得
lnρ=lna+blnVp
则
Δlnρ=bΔlnVp(9)
又因为
代入上式得
所以
即:
将上式代入公式(8)得:
即:
其中,
5)一般由Shuey近似计算截距和梯度属性P和G(进一步可计算横波阻抗反射系数Rs等属性)的公式为:
R(θ)=P+G0sin2θ(12)
其中,G0是根据平均角计算出的梯度
6)在实际应用公式(11)计算属性时,由于求取角度θ2不方便,常常用入射角θ1代替θ进行估算,即应用下式计算P和Rs:
R(θ)=P+G0sin2θ1(13)
比较式(11)和式(13)可得:
所以,计算横波阻抗反射系数Rs的公式为:
这里的P和Rs0由公式(13)计算得到。
本发明的优点是:计算方便、精确,与实际的参数误差小。
按照上面的方法原理,首先应用模型进行了分析,下面给出具体分析结果。
a、模型简介
本技术应用了三个模型,包括Ostrander所构造的含气砂岩模型[7]、Goodway根据实测资料给出的含气砂岩模型和含气砂岩的四种典型AVO模型[8],以及Marmousi II中的气层模型[9]。其参数分别见表1-3。
表1Ostrander(1984)的三层含油气砂岩与页岩模型
表2Goodway(1997)三层含油气砂岩与页岩模型
表3Marmousi II含气砂岩与页岩模型
b、修改前后P-P波反射系数曲线对比分析
图1-3中,实线为精确的zoeppritz方程计算结果,点线是将精确的shuey近似表示成随入射角变换的函数得到的结果,“+”线是应用修改之后的公式计算的结果。
从图1-3中,可以看出当上、下储层差值比较大时,如果直接应用入射角代替平均角存在比较大的误差。这样就为下一步的角道集提取,以及属性的计算带来很大的误差。而修改之后的公式在角度比较小时((0°<θ≤30°)),结果是非常准确的,与精确的zoeppritz方程计算结果之间差别非常的小。修改前后,误差定量分析如图所示。
在图4中,“1”代表Ostrander模型;“2”代表Goodway模型;“3”Marmousi II中的气层模型。黑色为修正之前的相对误差绝对值,灰色为修正之后的相对误差绝对值。
从图4中可以看出,经过修正之后误差大大降低。特别是气层引起上、下储层物性差别比较大的时候。如图中的“1”和“3”所示。
c、Marmousi II模型分析
从MarmousiII模型中抽取一道进行方法验证,所取的道集中包含了气、油层以及水层。
图6中黑线是用阻抗差计算的结果,星号是利用常规的方法计算Rs属性,菱形是用修正公式计算的结果。从图中可以看到在气层出现的地方,常规方法提取的Rs属性与用阻抗差公式计算得到结果是完全相反的。而应用修正公式修正后(如图中,菱形所示),结果能够比较好的与精确值吻合,特别是在气层出现的位置。图7是应用各种方法计算Rs然后与40Hz雷克子波褶积得到的合成地震记录。图7(a)是由模型的横波速度和密度直接计算Rs,然后再与40Hz雷克子波褶积得到的合成地震记录。图7(b)是应用常规方法计算的Rs与40Hz雷克子波褶积得到的合成地震记录。图7(c)是应用修正公式计算的Rs与40Hz雷克子波褶积得到的合成地震记录。从图7的对比中可以看出,常规方法计算的Rs在气层位置会产生很大的误差,出现了数值反转(如图7(b)所示),经过修正公式修正之后气层处的误差消失,并且数值与精确计算结果一致,如图7(c)所示。
附图说明:
图1为Ostrander模型。
图2为Goodway模型。
图3为Marmousi II模型中的气层。
图4为修改前后误差对比分析。
图5为从mar2模型中抽取的一道数据。
图6为Marmousi II模型的属性对比。
图7Marmousi II模型横波合成地震记录的比较。其中(a)为垂直入射计算的Rs与子波褶积得到的合成地震记录。(b)为常规方法计算的Rs与子波褶积得到的合成地震记录。(c)为修正之后的Rs与子波褶积得到的合成地震记录。
图8为实施例的叠前角度道集。其中(a)为某地区8°角度叠加道集。(b)为某地区16°角度叠加道集。(c)为某地区24°角度叠加道集。
图9为实施例基于修正公式计算的Rs属性。
具体实施方式
基于入射角的AVO近似公式及属性提取方法,实现过程包括:1)基于入射角的AVO近似公式的提出,2)基于入射角的AVO近似公式的理论推导,3)应用常规方法提取P、G,Rs的属性;4)将属性归一化;5)由井数据或者研究区域的经验公式确定b值,然后应用公式(14)和公式(15)对常规方法提取的属性进行修正;6)将修正结果与井数据进行对比分析,根据分析结果对修正结果进行校正,得到最终的属性结果。
具体的求取方法为:
1)由Aki和Shuey近似可得:
其中,
Vp=(Vp1+Vp2)/2,Vs=(Vs1+Vs2)/2,ρ=(ρ1+ρ2)/2
ΔVp=Vp2-Vp1,ΔVs=Vs2-Vs1,Δρ=ρ2-ρ1,θ=(θ1+θ2)/2
Vp1和Vp2分别是反射界面上下上、下地层的纵波速度,Vs1和Vs2分别是上、下地层的横波速度,θ1和θ2分别是地层的入射角和透射角
2)根据Snell定律有:
对公式(1)左边项进行泰勒展开,及
将公式(3)代入公式(2)可得:
又因为:
所以可得:
(6)
当采用一阶近似时,上式可以写成:
令
则在此基础上有
而当θ1和θ2相差不大时,
sinθ1sinθ2=sin2θ
3)这样可以得到基于入射角的AVO近似公式:
4)Gardner等人于1974提出的密度与纵波速度之间的经验公式的形式为[6]:
其中,a和b为系数(一般取a=0.31,b=0.25)
由上式可得
lnρ=lna+blnVp
则
Δlnρ=bΔlnVp(9)
又因为
代入上式得
所以
即:
将上式代入公式(8)得:
即:
其中,
5)一般由Shuey近似计算截距和梯度属性P和G(进一步可计算横波阻抗反射系数Rs等属性)的公式为:
R(θ)=P+G0sin2θ(12)
其中,G0是根据平均角计算出的梯度
6)在实际应用公式(11)计算属性时,由于求取角度θ2不方便,常常用入射角θ1代替θ进行估算,即应用下式计算P和Rs:
R(θ)=P+G0sin2θ1(13)
比较式(11)和式(13)可得:
所以,计算横波阻抗反射系数Rs的公式为:
这里的P和Rs0由公式(13)计算得到。
根据上述理论,对实际地震数据进行了叠前地震属性提取及分析。
Claims (2)
2.根据权利要求1所述的基于入射角的AVO近似公式及属性提取方法,其特征在于具体的求取方法为:
1)由Aki和Shuey近似可得:
其中,
Vp=(Vp1+Vp2)/2,Vs=(Vs1+Vs2)/2,ρ=(ρ1+ρ2)/2
ΔVp=Vp2-Vp1,ΔVs=Vs2-Vs1,Δρ=ρ2-ρ1,θ=(θ1+θ2)/2
Vp1和Vp2分别是反射界面上下上、下地层的纵波速度,Vs1和Vs2分别是上、下地层的横波速度,θ1和θ2分别是地层的入射角和透射角
2)根据Snell定律有:
对公式(1)左边项进行泰勒展开,及
将公式(3)代入公式(2)可得:
又因为:
(5)
所以可得:
当采用一阶近似时,上式可以写成:
令
则在此基础上有
而当θ1和θ2相差不大时,
sinθ1 sinθ2=sin2θ
3)这样可以得到基于入射角的AVO近似公式:
4)Gardner等人于1974提出的密度与纵波速度之间的经验公式的形式为[6]:
其中,a和b为系数(一般取a=0.31,b=0.25)
由上式可得
lnρ=lna+blnVp
则
Δlnρ=bΔlnVp(9)
又因为
代入上式得
所以
即:
将上式代入公式(8)得:
即:
其中,
5)一般由Shuey近似计算截距和梯度属性P和G(进一步可计算横波阻抗反射系数Rs等属性)的公式为:
R(θ)=P+G0sin2θ(12)
其中,G0是根据平均角计算出的梯度
6)在实际应用公式(11)计算属性时,由于求取角度θ2不方便,常常用入射角θ1代替θ进行估算,即应用下式计算P和Rs:
R(θ)=P+G0sin2θ1(13)
比较式(11)和式(13)可得:
所以,计算横波阻抗反射系数Rs的公式为:
这里的P和Rs0由公式(13)计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010285247 CN101957454A (zh) | 2010-09-17 | 2010-09-17 | 基于入射角的avo近似公式及属性提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010285247 CN101957454A (zh) | 2010-09-17 | 2010-09-17 | 基于入射角的avo近似公式及属性提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101957454A true CN101957454A (zh) | 2011-01-26 |
Family
ID=43484893
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010285247 Pending CN101957454A (zh) | 2010-09-17 | 2010-09-17 | 基于入射角的avo近似公式及属性提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101957454A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013123624A1 (zh) * | 2012-02-24 | 2013-08-29 | Chen Xingping | 一个探测煤层气富集高渗区的碳氢检测因子与使用方法 |
CN103675911A (zh) * | 2014-01-06 | 2014-03-26 | 中国石油大学(华东) | 一种基于纵波和转换横波联合反演截距和梯度的方法 |
CN104820239A (zh) * | 2015-05-13 | 2015-08-05 | 中国石油大学(华东) | 一种方位叠前地震属性解耦提取方法 |
CN104965226A (zh) * | 2015-06-17 | 2015-10-07 | 成都高普石油工程技术有限公司 | 一种岩体中油气信息的分析方法 |
CN106556862A (zh) * | 2015-09-29 | 2017-04-05 | 中国石油化工股份有限公司 | 用于页岩气勘探avo技术的pp波反射系数计算方法 |
CN106772615A (zh) * | 2016-11-09 | 2017-05-31 | 中国石油大学(华东) | 一种宽频地震多域联合avo反演方法 |
CN106908839A (zh) * | 2017-05-05 | 2017-06-30 | 中国石油大学(华东) | 一种地震波精确反射系数解析方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030046006A1 (en) * | 2001-08-31 | 2003-03-06 | Exxonmobil Upstream Research Company | Method for classifying AVO data using an interpreter-trained neural network |
US20040199330A1 (en) * | 2003-04-01 | 2004-10-07 | Conocophillips Company | Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data |
US20060004521A1 (en) * | 2004-07-02 | 2006-01-05 | Veritas Dgc Inc. | Seismic data calibration using crossplotted AVO attributes |
CN101281253A (zh) * | 2007-04-06 | 2008-10-08 | 中国石油集团东方地球物理勘探有限责任公司 | 一种利用振幅随偏移距变化特征提高油气检测精度的方法 |
CN101551466A (zh) * | 2008-04-03 | 2009-10-07 | 中国石油天然气集团公司 | 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法 |
CN101598805A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种多分量地震资料层位对比和标定方法 |
-
2010
- 2010-09-17 CN CN 201010285247 patent/CN101957454A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030046006A1 (en) * | 2001-08-31 | 2003-03-06 | Exxonmobil Upstream Research Company | Method for classifying AVO data using an interpreter-trained neural network |
US20040199330A1 (en) * | 2003-04-01 | 2004-10-07 | Conocophillips Company | Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data |
US20060004521A1 (en) * | 2004-07-02 | 2006-01-05 | Veritas Dgc Inc. | Seismic data calibration using crossplotted AVO attributes |
CN101281253A (zh) * | 2007-04-06 | 2008-10-08 | 中国石油集团东方地球物理勘探有限责任公司 | 一种利用振幅随偏移距变化特征提高油气检测精度的方法 |
CN101551466A (zh) * | 2008-04-03 | 2009-10-07 | 中国石油天然气集团公司 | 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法 |
CN101598805A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种多分量地震资料层位对比和标定方法 |
Non-Patent Citations (3)
Title |
---|
《地球学报》 20060228 孙鹏远等 P-SV波反射系数近似及其AVO属性特征 85-89页 1-2 , 第01期 2 * |
《石油地球物理勘探》 20070415 王振国等 AVO方法检测油气应用实例分析 194-197页 1-2 , 第02期 2 * |
《西北地震学报》 20080930 刘前坤等 基于等效Thomsen参数的P-SV波AVO属性研究 214-221页 1-2 , 第03期 2 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013123624A1 (zh) * | 2012-02-24 | 2013-08-29 | Chen Xingping | 一个探测煤层气富集高渗区的碳氢检测因子与使用方法 |
CN103675911A (zh) * | 2014-01-06 | 2014-03-26 | 中国石油大学(华东) | 一种基于纵波和转换横波联合反演截距和梯度的方法 |
CN104820239A (zh) * | 2015-05-13 | 2015-08-05 | 中国石油大学(华东) | 一种方位叠前地震属性解耦提取方法 |
CN104820239B (zh) * | 2015-05-13 | 2018-07-06 | 中国石油大学(华东) | 一种方位叠前地震属性解耦提取方法 |
CN104965226A (zh) * | 2015-06-17 | 2015-10-07 | 成都高普石油工程技术有限公司 | 一种岩体中油气信息的分析方法 |
CN106556862A (zh) * | 2015-09-29 | 2017-04-05 | 中国石油化工股份有限公司 | 用于页岩气勘探avo技术的pp波反射系数计算方法 |
CN106772615A (zh) * | 2016-11-09 | 2017-05-31 | 中国石油大学(华东) | 一种宽频地震多域联合avo反演方法 |
CN106772615B (zh) * | 2016-11-09 | 2018-11-13 | 中国石油大学(华东) | 一种宽频地震多域联合avo反演方法 |
CN106908839A (zh) * | 2017-05-05 | 2017-06-30 | 中国石油大学(华东) | 一种地震波精确反射系数解析方法及装置 |
CN106908839B (zh) * | 2017-05-05 | 2018-09-28 | 中国石油大学(华东) | 一种地震波精确反射系数解析方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11016211B2 (en) | 4D time shift and amplitude joint inversion for obtaining quantitative saturation and pressure separation | |
Maraschini et al. | A new misfit function for multimodal inversion of surface waves | |
CN101957454A (zh) | 基于入射角的avo近似公式及属性提取方法 | |
CN105759310B (zh) | 复杂非均质储层介质中地震波衰减及速度频散预测方法 | |
CN103733089B (zh) | 用于包括不确定性估计的地下表征的系统和方法 | |
CN103293552A (zh) | 一种叠前地震资料的反演方法及系统 | |
CN101446645A (zh) | 一种利用地震流体阻抗进行流体确定的方法 | |
Hutton | Geophysical modeling and structural interpretation of a 3D reflection seismic survey in Farnsworth Unit, TX | |
CN102928870A (zh) | 基于正则化的非线性地震叠前弹性参数反演方法 | |
CN103675911A (zh) | 一种基于纵波和转换横波联合反演截距和梯度的方法 | |
CN107390270A (zh) | 一种基于弹性波逆时偏移ADCIGs的AVA分析方法 | |
Leaney | Microseismic source inversion in anisotropic media | |
Mallick et al. | Azimuthal anisotropy analysis of P-wave seismic data and estimation of the orientation of the in situ stress fields: An example from the Rock-Springs uplift, Wyoming, USA | |
CN104199088B (zh) | 一种提取入射角道集的方法及系统 | |
CN107728205A (zh) | 一种地层压力预测方法 | |
CN104977610A (zh) | 基于入射角的avo近似公式进行属性提取的方法 | |
CN104422955B (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
CN103984018A (zh) | 多波联合振幅随入射角变化反演方法 | |
EP3301485B1 (en) | Crosswell microseismic system | |
US20120099396A1 (en) | System and method for characterization with non-unique solutions of anisotropic velocities | |
CN104280767A (zh) | 一种基于柯西分布的稀疏脉冲反演方法 | |
Babasafari et al. | Practical workflows for monitoring saturation and pressure changes from 4D seismic data: A case study of Malay Basin | |
CN105527648A (zh) | 用于各向异性参数反演的敏感度矩阵的计算方法及系统 | |
CN104965226A (zh) | 一种岩体中油气信息的分析方法 | |
CN101290356A (zh) | 一种炸药震源勘探的辅助方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20110126 |