【技术实现步骤摘要】
一种裂缝介质PP波与分裂PS波AVO联合反演方法
本专利技术涉及地震勘探
，具体涉及一种裂缝介质PP波与分裂PS波AVO联合反演方法。
技术介绍
AVO(amplitudevariationwithoffset)/AVA(amplitudevariationwithangle)反演技术是以弹性波理论为基础，利用叠前地震道集反射振幅随炮检距(或入射角)的变化规律，来估算预测储层的物性参数和岩石特征的重要技术。当地震波入射到反射界面时，Zoeppritz方程描述了纵波入射横波入射与各自激发的反射透射纵波横波之间的关系，为AVO技术提供了理论基础；Aki和Richards[1]简化了Zoeppritz方程，得到了基于纵横波速度和密度的PP波和PS波AVO近似公式；Goodway[2]分析了拉梅常数对碳氢化合物的敏感程度推导了AVO近似公式；Gray[3]基于Aki和Richards的公式上推导了用拉梅常数及密度表达的反射系数近似方程；Rüger[4]推导了用Thomsen各向异性参数表达的AVO反射近似公式，并对其公式的精确性及适用介质进行了对比与分析。裂缝诱导的...
【技术保护点】
1.一种裂缝介质PP波与分裂PS波AVO联合反演方法，其特征在于，包括如下步骤：S1、抽取PP波AVA道集、PS1波AVA道集和PS2波AVA道集，并将PS1波、PS2波的时间压缩到PP波的时间上；S2、初始模型的建立：通过测井数据获得初始的纵波的垂直速度、快横波的垂直速度以及密度模型，并根据岩石物理测井信息或者速度分析，建立初始的各向异性参数模型；S3、反演：S3.1对各向同性弹性参数即纵波的垂直速度α、快横波的垂直速度β和密度ρ进行反演：S3.1.1、根据初始的纵波的垂直速度、快横波的垂直速度以及密度模型的模型数据，用Rüger的反射系数近似公式分别计算PP波、PS1波...
【技术特征摘要】
1.一种裂缝介质PP波与分裂PS波AVO联合反演方法，其特征在于，包括如下步骤：S1、抽取PP波AVA道集、PS1波AVA道集和PS2波AVA道集，并将PS1波、PS2波的时间压缩到PP波的时间上；S2、初始模型的建立：通过测井数据获得初始的纵波的垂直速度、快横波的垂直速度以及密度模型，并根据岩石物理测井信息或者速度分析，建立初始的各向异性参数模型；S3、反演：S3.1对各向同性弹性参数即纵波的垂直速度α、快横波的垂直速度β和密度ρ进行反演：S3.1.1、根据初始的纵波的垂直速度、快横波的垂直速度以及密度模型的模型数据，用Rüger的反射系数近似公式分别计算PP波、PS1波、PS2波的反射系数，此时各向异性参数ε(V)、δ(V)、γ为0，不对各向异性参数进行反演；所述Rüger的反射系数近似公式如下：其中，Rpp、Rps1、Rps2分别表示PP波、PS1波、PS2波的反射系数，α表示纵波的垂直速度，β表示快横波的垂直速度，ε(V)指纵波水平速度和垂直速度的差值与纵波垂直速度的比值，表示纵波的各向异性程度；δ(V)表示纵波和横波各向异性相对大小；γ表示横波速度的各向异性发育程度；Z表示垂直纵波阻抗Z＝ρα，G表示对应的垂直传播的波的剪切模量G＝ρβ2，i为纵波入射角，j为快横波反射角，θ为入射方位角，下标1代表上层介质的物性参数，下标2代表下层介质的参数；S3.1.2、采用理论雷克子波，利用计算得到的PP波、PS1波、PS2波的反射系数分别制作PP波、PS1波、PS2波的合成记录；S3.1.3、联合PP波、PS1波、PS2波计算目标函数如下：Q＝||(H+λI)ΔΜ-JT(Robs-RM)||2；其中，M＝(αβρεδγ)T为参数模型矩阵，根据高斯牛顿法，在真实值R附近进行一阶泰勒展开得到JΔM＝Rk+1-Rk；k代表迭代次数，ΔM＝(ΔαΔβΔρΔεΔδΔγ)T是更新的模型，J为HTI介质的雅可比矩阵；Robs和RM分别为观测反射系数矩阵和模型反射系数矩阵；S3.1.4、当目标函数不满足精度要求，则修改纵波的垂直速度、快横波的垂直速度以及密度模型的模型数据，并重新返回步骤S3.1.1的反射系统计算步骤，再次进行反演；当目标函数满足精度要求，停止迭代；最终得到达到测井分辨率的反演结果纵波的垂直速度α、快横波的垂直速度β和密度ρ；S3.2、以步骤S3.1的反演结果作为约束条件，对各向异性参数进行反演：S3.2.1、利用初始的各向异性参数模型的模型数据，用Rüger的反射系数近似公式分别计算PP波、PS1波、PS2波的反射系数；S3.2.2、采用理论雷克子波，利用计算得到的...
【专利技术属性】
技术研发人员：芦俊，杨帅，王赟，
申请(专利权)人：中国地质大学北京，
类型：发明
国别省市：北京,11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。