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

权利要求书

1.一种基于岩石物理的叠前裂缝定量预测方法，其特征在于，所述叠前裂缝定量预测方法包括：

根据待测储层的叠前CMP道集中的方位角及入射角的大小，将所述叠前CMP道集划分为多个方位角道集：方位角道集1，方位角道集2，…，方位角道集n；

对各所述方位角道集进行叠加、偏移处理，计算各方位角道集的纵波反射系数；

根据Ruger方程精确式以及各所述方位角道集及其纵波反射系数，计算待测储层中裂缝的切向各异性系数Δ_T；

根据所述切向各异性系数Δ_T计算待测储层的裂缝密度e；

其中，所述的根据Ruger方程精确式以及各所述方位角道集及其纵波反射系数，计算待测储层中裂缝的切向各异性系数Δ_T，包括：

根据HTI介质中各向异性参数与裂缝岩石物理参数之间的关系，对Ruger方程精确式进行变换，得到纵波反射系数简化方程：

x1+x2sin2i+x3sin2i tan2i+x4[-g(1+sin2φtan2i)+2g]cos2φsin2i+x5g cos2φsin2i[(g sin2φ+g-1)tan2i+2g-1]=r;]]>

将所述方位角道集及其对应的纵波反射系数代入所述纵波反射系数简化方程得到矩阵方程：

Ax＝b；

根据所述矩阵方程建立一用于反演所述x的最优化目标函数：

min E(x)obj=||Ax-b||22;]]>

对所述最优化目标函数进行求解，得到所述切向各异性系数Δ_T；

其中，i和φ分别表示纵波的入射角和方位角；r为与入射角i和方位角φ有关的纵波反射系数；z＝ρα为纵波波阻抗；ρ为介质密度，g/cm³；α为纵波速度，m/s；G＝ρβ²为横波切向模量，β为横波速度；ε、δ和γ为Thomsen的各向异性参数；Δ_T为裂缝岩石物理参数中的切向各异性系数；Δ_N为裂缝岩石物理参数中的法向各异性系数；g＝μ/(λ+2μ)为横纵波速度比的平方，λ、μ为不含裂隙各向同性岩石的拉梅参数；

A=1sin2i1W1[-g(1+R1)+2g]P1P1[(g sin2φ1+g-1)tan2i1+2g-1]g1sin2i2W2[-g(1+R2)+2g]P2P2[(g sin2φ2+g-1)tan2i2+2g-1]g...............1sin2inWn[-g(1+Rn)+2g]PnPn[(g sin2φn+g-1)tan2in+2g-1]g,]]>

R_k、P_k及W_k为中间变量，R_k＝sin²φ_k tan²i_k，P_k＝cos²φ_k sin²i_k，W_k＝sin²i_k tan²i_k，r_k＝r_k(i_k,φ_k)为方位角道集k对应的纵波反射系数，b＝[r₁(i₁,φ₁)；r₂(i₂,φ₂)；...；r_n(i_n,φ_n)]^T，x＝[x₁,x₂,x₃,x₄,x₅]^T，k∈[1，n]，n为方位角道集的数量。