[发明专利]一种裂缝介质PP波与分裂PS波AVO联合反演方法有效

申请号：	201810894632.5	申请日：	2018-08-08
公开（公告）号：	CN108828660B	公开（公告）日：	2019-10-22
发明（设计）人：	芦俊;杨帅;王赟	申请（专利权）人：	中国地质大学（北京）
主分类号：	G01V1/28	分类号：	G01V1/28
代理公司：	北京国昊天诚知识产权代理有限公司 11315	代理人：	许志勇
地址：	100083***	国省代码：	北京;11
权利要求书：	查看更多	说明书：	查看更多
摘要：
搜索关键词：	裂缝反演反射系数近似各向同性介质各向异性参数牛顿迭代算法初始模型弹性参数时间压缩有效信息分裂联合高斯流体岩性地层匹配
钻瓜网技术展会专利词库专利权人专利榜在售专利公布日期热门专利

【权利要求书】：

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的反射系数近似公式如下：

其中，R_pp、R_ps1、R_ps2分别表示PP波、PS1波、PS2波的反射系数，α表示纵波的垂直速度，β表示快横波的垂直速度，ε^(V)指纵波水平速度和垂直速度的差值与纵波垂直速度的比值，表示纵波的各向异性程度；δ^(V)表示纵波和横波各向异性相对大小；γ表示横波速度的各向异性发育程度；Z表示垂直纵波阻抗Z＝ρα，G表示对应的垂直传播的波的剪切模量G＝ρβ²，i为纵波入射角，j为快横波反射角，θ为入射方位角，下标1代表上层介质的物性参数，下标2代表下层介质的参数；

S3.1.2、采用理论雷克子波，利用计算得到的PP波、PS1波、PS2波的反射系数分别制作PP波、PS1波、PS2波的合成记录；

S3.1.3、联合PP波、PS1波、PS2波计算目标函数如下：

Q＝||(H+λI)ΔΜ-J^T(R_obs-R_M)||²；

其中，M＝(α β ρ ε δ γ)^T为参数模型矩阵，根据高斯牛顿法，在真实值R附近进行一阶泰勒展开得到JΔM＝R_k+1-R_k；k代表迭代次数，ΔM＝(Δα Δβ Δρ Δε Δδ Δγ)^T是更新的模型，J为HTI介质的雅可比矩阵；R_obs和R_M分别为观测反射系数矩阵和模型反射系数矩阵；λ为一个标量，I为单位矩阵；

S3.1.4、当目标函数不满足精度要求，则修改纵波的垂直速度、快横波的垂直速度以及密度模型的模型数据，并重新返回步骤S3.1.1的反射系统计算步骤，再次进行反演；当目标函数满足精度要求，停止迭代；最终得到达到测井分辨率的反演结果纵波的垂直速度α、快横波的垂直速度β和密度ρ；

S3.2、以步骤S3.1的反演结果作为约束条件，对各向异性参数进行反演：

S3.2.1、利用初始的各向异性参数模型的模型数据，用Rüger的反射系数近似公式分别计算PP波、PS1波、PS2波的反射系数，所述Rüger的反射系数近似公式中：

S3.2.2、采用理论雷克子波，利用计算得到的PP波、PS1波、PS2波的反射系数分别制作PP波、PS1波、PS2波的合成记录；

S3.2.3、按照步骤S3.1.3的方法联合PP波、PS1波、PS2波计算目标函数；

S3.2.4、当目标函数不满足精度要求，则修改纵波的垂直速度、快横波的垂直速度以及密度模型的模型数据，并重新返回步骤S3.2.1的反射系统计算步骤，再次进行反演；当目标函数满足精度要求，停止迭代；最终得到达到测井分辨率的反演结果各向异性参数ε^(V)、δ^(V)、γ；

其中，

所述目标函数中，雅可比矩阵J在AVO里的扩展表达式为：

反射系数关于弹性参数和各向异性参数的一阶偏导，用中心差分方法表示：

运用PP波、PS1波、PS2波AVA联合反演，则雅可比矩阵J为：

且H＝J^TJ，则有：

其中，e_PP、e_PS1和e_PS2分别是PP波，PS1波、PS2波AVO的残留误差，θ₁…θ_i代表了纵波的每一个入射角，R₁…R_i代表了每个角度所对应的PP波、PS1波、PS2波的反射系数。