[发明专利]基于岩石物理的叠前裂缝定量预测方法及系统

摘要

本发明提供了一种基于岩石物理的叠前裂缝定量预测方法及系统，所述叠前裂缝定量预测方法包括根据待测储层的叠前CMP道集中的方位角及入射角的大小，将所述叠前CMP道集划分为多个方位角道集方位角道集1，方位角道集2，…，方位角道集n；对各所述方位角道集进行叠加、偏移处理，计算各方位角道集的纵波反射系数；根据Ruger方程精确式以及各所述方位角道集及其纵波反射系数，计算待测储层中裂缝的切向各异性系数ΔT；根据所述切向各异性系数ΔT计算待测储层的裂缝密度e。本发明能够准确地实现裂缝密度的定量估算，进而准确地进行油气藏预测和指导开发井位的部署，达到了降低油气藏预测风险的目的。

主权项

一种基于岩石物理的叠前裂缝定量预测方法，其特征在于，所述叠前裂缝定量预测方法包括：根据待测储层的叠前CMP道集中的方位角及入射角的大小，将所述叠前CMP道集划分为多个方位角道集：方位角道集1，方位角道集2，…，方位角道集n；对各所述方位角道集进行叠加、偏移处理，计算各方位角道集的纵波反射系数；根据Ruger方程精确式以及各所述方位角道集及其纵波反射系数，计算待测储层中裂缝的切向各异性系数ΔT；根据所述切向各异性系数ΔT计算待测储层的裂缝密度e；其中，所述的根据Ruger方程精确式以及各所述方位角道集及其纵波反射系数，计算待测储层中裂缝的切向各异性系数ΔT，包括：根据HTI介质中各向异性参数与裂缝岩石物理参数之间的关系，对Ruger方程精确式进行变换，得到纵波反射系数简化方程：x1+x2sin2i+x3sin2itan2i+x4[-g(1+sin2φtan2i)+2g]cos2φsin2i+x5gcos2φsin2i[(gsin2φ+g-1)tan2i+2g-1]=r;]]>将所述方位角道集及其对应的纵波反射系数代入所述纵波反射系数简化方程得到矩阵方程：Ax＝b；根据所述矩阵方程建立一用于反演所述x的最优化目标函数：minE(x)obj=||Ax-b||22;]]>对所述最优化目标函数进行求解，得到所述切向各异性系数ΔT；其中，i和φ分别表示纵波的入射角和方位角；r为与入射角i和方位角φ有关的纵波反射系数；z＝ρα为纵波波阻抗；ρ为介质密度，g/cm3；α为纵波速度，m/s；G＝ρβ2为横波切向模量，β为横波速度；ε、δ和γ为Thomsen的各向异性参数；ΔT为裂缝岩石物理参数中的切向各异性系数；ΔN为裂缝岩石物理参数中的法向各异性系数；g＝μ/(λ+2μ)为横纵波速度比的平方，λ、μ为不含裂隙各向同性岩石的拉梅参数；A=1sin2i1W1[-g(1+R1)+2g]P1P1[(gsin2φ1+g-1)tan2i1+2g-1]g1sin2i2W2[-g(1+R2)+2g]P2P2[(gsin2φ2+g-1)tan2i2+2g-1]g...............1sin2inWn[-g(1+Rn)+2g]PnPn[(gsin2φn+g-1)tan2in+2g-1]g,]]>Rk、Pk及Wk为中间变量，Rk＝sin2φk tan2ik，Pk＝cos2φk sin2ik，Wk＝sin2ik tan2ik，rk＝rk(ik,φk)为方位角道集k对应的纵波反射系数，b＝[r1(i1,φ1)；r2(i2,φ2)；...；rn(in,φn)]T，x＝[x1,x2,x3,x4,x5]T，k∈[1，n]，n为方位角道集的数量。