现代试井解释方法.doc

1现代试井解释方法现代试井解释方法现代试井解释时期以 70 年代初雷米发表关于均质油层双对数拟合图版为开始标志其特点为： 建立双对数拟合分析法，可以运用早期试井数据； 给出半对数直线段出现时间，使常规分析更可靠； 采用图版曲线拟合法和数值模拟法，使用计算机，解释模型多； 解释过程是“边解释，边检验边解释，边检验”的过程，保证解释的可靠性 试井解释模型可按照基本模型及边条件划分： 基本模型：1. 均质油藏；2．非均质油藏：多层油藏，渗透率变化；3．双重空隙介质油藏：拟稳态窜流，不稳态窜流4．双孔双渗介质油藏：拟稳态窜流，不稳态窜流 内边界条件：1. 井筒储存； 外边界条件：1. 无限大地层； 2. 表皮系数； 2. 不渗透边界； 3. 裂缝切割井； 3. 恒压边界； 4. 打开不完善； 4 封闭边界； 5．水平井； 由基本模型, 内边界条件和外边界条件，可组合出许多试井解释模型，它们的拟合图版曲线可用计算机 快速计算出来§1 试井使用的无量纲物理量 （1-1） wDrrr2 wtt DrCktt)(ppBqkhpi pDwe eDrrr)(wfi pwDppBqkhp2 wtc DhrCCC其中是单位制换算系数，各单位制的单位及换算系数如下所示：ctp,,单位制基本SI实用SI工程制英矿场制长度 r,h,LMmmft时间 tShhh压力 pPaMPaPsi渗透率 k(=mD)mDmD产量q (油)(地面)(地面)1.84219.033141.210.00360.00026370.8936503121/3.683*10438.06MSCF/d粘度Pa.smPa.s (= cp)cpcp产量qsc (气)1042/cmkg2m23*10msm /3dm /3dm /3dSTB/sm /3dm /3dm /3p21t410484. 3c212121由于无量纲物理量与单位制无关，利用此表可方便地进行单位制换算。

利用上述无量纲表达式，基本微分方程式变成：2（1-2）DDDD D DDrp rprrr     1将边条件及初条件无因次化，与上式一同求解，即得问题的解 （1-1）式给出了问题解的无因次)(DDtp量与有因次量之间的关系 （1-1）式取对数得：, (1-3)BqkhpppDlglglg2lglglgwtt Drcktt上式说明，若将 p-t 关系绘成双对数坐标图，无因次曲线与有因次曲线形状完全相同，解的无因次量坐 标与有因次量坐标之间相差同一常数因此，通过上下左右移动有因次曲线，可以和无因次曲线完全重 合这就是可用双对数坐标曲线拟合的原因§2 流动阶段的识别诊断曲线——双对数图上绘出的曲线特种识别曲线—— 特定坐标下体现其流动特征的曲线 （1） 早期阶段 井筒储存期：压力与时间关系为， ， （2-1）DD DCtptCqBpctp CqBtpctp lglglg这说明在井筒储存期，双对数图上是斜率为一的直线，直角坐标下也是直线双对数图上是斜率为一的直线，直角坐标下也是直线 无限导流垂直裂缝井：压力与时间关系为， （2-2）tLp Ltplglglg21这说明该井早期阶段的双对数图上是斜率为 0.5 的直线，为一直线。

tp_有限导流垂直裂缝井：压力与时间关系为， （2-3）41 Ltp Ltplglg41lg这说明该井诊断曲线为早期阶段的双对数图上是斜率为 0.25 的直线，为一直线4_tp（2） 无限大作用径向流阶段压降方程：（2-4） SrcktkhqBppwttp iwf87. 04lglg15. 12对于实用 SI 制变成：（2-5） 0923. 287. 0lglg12. 22SrcktkhqBppwtiwf3可见—lgt 呈直线是特种识别曲线wfp压恢方程，对于赫诺法有（2-6）   ttt khqBpppp iwslg15. 1对于实用 SI 制变成：（2-7）   ttt khqBpppiwslg12. 2对于 MDH 法有（2-8）  SrcktkhqBppwttp wfws87. 04lglg15. 12对于实用 SI 制变成：（2-9） 0923. 287. 0lglg12. 22SrcktkhqBppwtwfws可见特种识别曲线是—或者—呈直线。

wfptlgwfptttplg（3） 外边界反映阶段 1） 恆压边界：压力最后稳定边平； 2） 不渗透边界或者直线断层附近一口井：半对数图出现两个直线段，后者斜率是前者二倍；3） 封闭油藏：达到拟稳态后，，p—t 图上出现直线consttp（4） 识别油藏类型的重要性 对于不同的油藏类型，同一参数的意义可能不同例如解释出的 kh，对于均质油层是单层 kh 值， S1000 未污染井 5—1000 酸化见效井 0.5—5 压裂见效井 <0.5 解释图版还用红线标出纯井筒储存期结束和半对数直线出现时间Bourdet 给解释图版加入了压力导数曲线，即对于每条—曲线，又增加了相应的)/(DDwDCtpDDCt/—关系曲线，构成复合图版，大大提高了拟合的唯一性)('DD D DD CtpCtDDCt/压力导数曲线的特点及应用在纯井筒储存期有关系式：， DD DCtpD DD D DDpCtpCt'因而，对于不同的值，压力曲线都与压力导数曲线重合，是斜率为 1 的直线S DeC2在半对数直线期有关系式：StStEippDDwDD)4ln(21)41(21 于是有：821)('DD DDDDDD D DD dtdptCtddp CtpCt因而，对于不同的值，压力导数曲线重合于纵坐标为 0.5 的水平线。

在纯井筒储存期和半对S DeC2数直线期之间，压力导数曲线是一组上凸曲线 压降试井的拟合解释过程1）计算每个测点的压差及其数值导数，列出数据表wfippp'p',,ptpt2）拟合过程：取一张和解释图版相同的空白透明双对数坐标纸，绘上实测压力数据和实)_(tp测压力导数数据)_('tpt3）将绘上实测压力数据的透明双对数坐标纸蒙在解释图版纸上，保持坐标线平行，上下左右移动透明双对数坐标纸，使实测压力及导数数据与相应于某个值的压力及导数导数曲线拟S DeC2合最佳，记录此曲线拟合值为MS DeC)(24）在坐标纸上任选一点作为拟合点，读出该点在透明双对数坐标纸上和解释图版纸上的坐标拟合值：，，，Mp)(Mt)(MDp )(M DD Ct)(5）根据已获得的拟合值，由下式计算参数：纵坐标拟合：， （3-11）MMD ptpqBkh )()()(kh hk 横坐标拟合：由 ，Ctkh CtCtDD  得 ， （3-12）M DDMCtCttkhC )()(  2 wtc DhrCCC曲线拟合得： （3-13）DMS D CeCS)(ln212 6）在双对数拟合图上读出半对数直线段开始时间，绘制 的半对数曲线，若半对数直线)_(tp段足够长，则用传统方法解释 k 和 S，结果应和拟合结果接近，否则重新拟合解释，直到二 者结果一致为止。

若半对数直线段太短，可不作半对数分析9§4 压恢试井的拟合解释过程1)利用压降图版进行压恢试井的拟合解释10由叠加原理得（4-1）)()()()(DDDpDDpDDpwsiwfiwfwswstpttptpkhqBppppppp由于总有 ，故有)()(DpDDpDDttptp)()(tptpkhqBpDDp ws因此，关井时间的压力恢复值小于生产同一时间产生的压降值因此，一般说来压力恢复曲线不同t 于压力降落曲线由下式计算压力恢复标准曲线： （4-2）)()()(DDDpDDpDDwsDtpttptpp因为它与有关，因而没有压恢图版，压恢试井的拟合解释利用压降图版进行，原因如下pt若是满足条件，故有： ， pppptttttttmax)()(DpDDpDDttptp)()(tptpkhqBpDDp ws这时关井时间的压力恢复值近似等于生产同一时间产生的压降值，因此整条压力恢复曲线可用t 压降图版进行拟合解释若是条件不满足，这时不能拟合整条压力恢复曲线，但是对于满足的点子，maxttpttp仍有，也就是说，压力恢复曲线早期近似等于早期压降曲线，)()(tptpkhqBpDDp ws可用压降图版拟合，其它部分偏离压降图版曲线不能拟合。

由偏离时间可算出之值，tptt /它应与右边所标之值一致，否则重新拟合 ptt /压力恢复资料的校正处理 利用（4-1）式和大时间的半对数近似式，不难得出：     tttmtptttkhqBtptpttptpkhqBtpppDpDpDpDDDpDDpDDp wslg)()(4ln214ln21)()()()()(于是（4-3）)(lg)(tptttmtppp ws11因此，若关井前的压降过程已达到半对数直线段和可以估计其斜率 m,则可用上式将压恢压差折 算为相应的压降压差，于是可用压降图版拟合解释 压力恢复的导数曲线由（4-2）得)()()()(DD D DDpD D DpD D DD wsDCtpCttpCtpCtp当半对数近似可用时，对求导得：DD Ct于是pDDpDDDDDpDDDDDpD D DDD DDwsDtttCt Ctt CtCttpCtpCtp5 . 05 . 05 . 0)()()('''5 . 0'pDDpDDDwsDtttCtp因此，如果把压降导数图版的纵坐标看作（） ，而不是原来的（） ，把 pDDpDDDwsDtttCtp'DDDCtp'横坐标看作，则压降导数图版就成为压恢导数图版。

由于DD CtppppDDpDDDwsDttttpqBkh tttCtp***.'' 故有因次的压力导数应取 pp ttttp**'2） 压恢试井拟合解释步骤12计算各测压点的压差及导数wfwswsppp pp wsttttp**'利用压降综合图版进行拟合，步骤与压降拟合相似 根据横坐标拟合值，纵坐标拟合值及曲线拟合值，计算参数 kh，S，C，方法与压降试井相同利用计算机，通过计算无因次赫诺曲线，决定原始地层压力 由叠加原理得：，)()(DDDpDDp wsitpttpkhqBpp， DMMp pp khqB)()()(DDDpDD MDMwsitpttppppp计算参数 kh，S，C 之后，即可计算上式右边，正确的原始地层压力应使得左右两边很好拟ip合3） 试井中应用自动拟合法 以上所述试井解释方法，基本属于手工解释范围，它有下列问题： 压恢试井关。