砂岩储层非均匀损害深度的计算方法

摘要

本发明涉及一种砂岩储层非均匀损害深度的计算方法，属油气勘探开发储层保护技术领域，其特征在于具体步骤如下：（1）收集井身参数、边界参数、地层参数、流体参数、射孔参数、生产参数；（2）将总表皮系数S进行分解，求得地层损害表皮系数S_d；（3）设定给定地层渗透率k与井壁渗透率k_w的比值的初值；（4）用牛顿迭代法求解下列方程$S_d=\frac{h}{h_p}\{-\ln \frac{r_s}{r_w}+{\left(\frac{k}{k_w}\right)}^{\frac{r_s}{r_s-r_w}}[-E_i(-\frac{r_w}{r_s-r_w}\ln \frac{k}{k_w})+E_i(-\frac{r_s}{r_s-r_w}\ln \frac{k}{k_w})\left]\right\};$求得地层损害深度为l_s＝r_s‑r_w。本发明可计算砂岩储层非均匀损害的损害深度，使其更加符合实际，对储层保护效果评价具有重要意义。

主权项

砂岩储层非均匀损害深度的计算方法，其特征在于，具体步骤如下：(1)资料数据的收集收集井身参数、边界参数、地层参数、流体参数、射孔参数和生产参数；(2)地层损害表皮系数的求得将总表皮系数S进行分解，求得地层损害表皮系数S_d，即S_d＝S‑S_pt‑S_pf‑S_nD‑S_CA‑S_θ‑S_CP‑S_an‑S_b‑S_vq‑S_p式中：S‑总表皮系数；S_pt‑部分打开表皮系数；S_pf‑射孔表皮系数；S_nD‑非达西表皮系数；S_CA‑油藏形状表皮系数；S_θ‑井斜表皮系数；S_cp‑相变表皮系数；S_an‑各向异性表皮系数；S_b‑边界表皮系数；S_vq‑变产量表皮系数；S_p‑压敏表皮系数；(3)地层渗透率k与井壁渗透率k_w的比值初值的给定根据实验认识与地层损害表皮系数S_d之间的大小关系，给定地层渗透率k与井壁渗透率k_w的比值的初值；(4)用牛顿迭代法求解下列方程(1)$S_d=\frac{h}{h_p}\{-ln\frac{r_s}{r_w}+{\left(\frac{k}{k_w}\right)}^{\frac{r_s}{r_s-r_w}}[-E_i(-\frac{r_w}{r_s-r_w}ln\frac{k}{k_w})+E_i(-\frac{r_s}{r_s-r_w}ln\frac{k}{k_w})\left]\right\}---\left(1\right)$式中：r_s—损害半径，m；r_w—井半径，m；h‑地层厚度，m；h_p‑打开的地层厚度，m；令$f\left(r_s\right)=-E_i\left(-\frac{r_w}{r_s-r_w}\ln \frac{k}{k_w}\right)+E_i\left(-\frac{r_s}{r_s-r_w}\ln \frac{k}{k_w}\right)-{\left(\frac{k}{k_w}\right)}^{-\frac{r_s}{r_s-r_w}}\left(\ln \frac{r_s}{r_w}+\frac{h_p}{h}{S}_d\right)$$f^{\prime }\left(r_s\right)=\{{\left(\frac{k}{k_w}\right)}^{-\frac{r_a}{r_s-r_w}}-{\left(\frac{k}{k_w}\right)}^{-\frac{r_s}{r_s-r_w}}[1+\frac{r_w}{r_s-r_w}(ln\frac{r_s}{r_w}+\frac{h_p}{h}{S}_d)ln\frac{k}{k_w}\left]\right\}/(r_s-r_w)$r_sn+1＝r_sn‑f(r_sn)/f′(r_sn)(n＝0，1，2，……)  (2)损害半径初值设定为12r_w，如果|r_sn+1‑r_sn|≤0.001，停止迭代，则r_sn+1为方程(2)的解，即：地层损害深度为l_s＝r_s‑r_w；(5)根据砂岩储层地层损害深度不超过1.5米的约束条件来计算非均匀损害深度当迭代求出的地层损害深度大于1.5米时，修正值，将值增加一个步长重新进行迭代求解，直至得到符合要求的地层损害深度。