气井的现代试井解释方法.

12流体在多孔介质中渗流时，压力变化服从方程：221136prrprpt.假设：流体是弱可压缩的，且其压缩系数为常数。气体：粘度和压缩系数等都是压力的函数；真实气体的偏差系数Z也是压力的函数。31965年，Al-Hussainy和Ramey引进了“真实气体的势函数(Realgaspotential)”，或称为“拟压力(Pseudopressure)”的概念：()PPZdPPPo2气体渗流方程221136rrrt.4§１拟压力的计算可用最简单的“梯形法”计算拟压力：()[()()]()pPZdPPZPZPPPPjnjjjj212220111(P)P，MPa5§２试井解释方法气井的现代试井解释和油井十分相似。一、Gringarten(格林加坦)图版拟合气井无因次压力的定义是：pKhqTTppKhqTpDscfscf002714378489.().()式中：(p)－拟压力差，MPa2/(mPa·s)ppptptpiwfwswf()[()][()]()压降情形恢复情形6q－气井产量，104m3/d；Tf－气层温度，K；Psc－标准状态的压力，Psc=0.101325MPa；Tsc－标准状态的温度，Tsc=293.15K；K－气层渗透率，m2；h－气层厚度，m。tD和CD的定义和油井情形相同。7用格林加坦图版进行气井试井解释的步骤实测压力P拟压力(P)拟压力差(P)()()ptpt关系曲线（压降情形）关系曲线（恢复情形）拟合8KqThppfD78489.()拟合CKhtCtKhtCtDDDD721721.(/).(/)拟合拟合CCChrDtw22SCeCaDSD122ln()拟合Sa－拟表皮系数9拟表皮系数Sa＝真表皮系数S+非达西流造成的无因次附加压降D·q式中：D－惯性-湍流系数，(104m3/d)Saq(104m3/d)S10格林加坦图版是压降图版。压力恢复测试同油井的压力恢复解释当关井前生产时间很长时，压力恢复的双对数曲线才能真正与格林加坦图版中的某一条样板曲线相拟合。在关井前生产时间不够长的情形，将只有前一段实测曲线可以与图版中某一条样板曲线相拟合，而后一段则不能。11例11某砂岩气藏的一口探井压力降落测试。图版拟合：CDe2S=106选拟合点：PD=1.70(P)=104tD/CD=25.4t=1KqThppmfD784894304443757848917117010000241942.()......拟合12CKhtCtmMPaDD7217200024191710021254118423.(/)....././拟合CCChrDtw21842201310014841710153919522......SCeCaDSD121210391953926ln()ln..拟合13二、Bourdet（布德）图版拟合分析KqThptCdpdttfDDD78489.(/)()/'拟合KqThptCdpdtttttfDDDpp78489.(/)()/()/'拟合压降情形恢复情形计算C、CD和Sa的公式与格林加坦图版拟合分析相同。14例１２气井压力恢复测试的导数曲线分析。tCpdpdtttttMPaMPastCthCeDDDppDDDS'.()()/().14910309110422615KqThptCdpdtttttmfDDDpp78489430444375784891711491000021242.(/)()/()/......'拟合CKhtCtmMPaDD721720002121710021309113273.(/)....././拟合16CCChrDtw21327201310014841710152823622......SCeCaDSD121210282364126ln()ln..拟合17三、试井解释步骤１、初拟合划分流动阶段。２、特种识别曲线分析（ptwf)~lg压降情形恢复情形()~lgptttwsp霍纳曲线()~lgptwsMDH曲线18气井的压降方程()().(lg..)pppTqTKhKtCrSwfiscscftwa424209077086862压力恢复的霍纳方程为()().lgpppTqTKhtttwsiscscfp4242MDH方程()().(lg..)pppTqTKhKtCrSwswfscscftwa42420907708686219斜率的绝对值为mpTqTKhscscf4242.可算得：KpTqTmhqTmhscscff4242001467..SpphmKCraiwftw11511090772.()[(()]lg.压降SphpmKCrawswftw11511090772.[(()]()lg.恢复20３、终拟合由中期特种识别曲线直线段的斜率计算压力拟合值：().pmD拟合1151用它对初拟合进行修正，并计算各项参数。21１、调整参数，产生样板曲线，与实测压力曲线进行拟合；２、绘制无因次霍纳曲线，进行解释结果的检验；３、进行压力史拟合，进一步检验解释结果的可靠性。计算机进行解释：22§３拟压力的简化１、(P)简化为P210203040P，MPaPw13.8MPaZ=C0()ppZdpCpp2100223()().(lg..)pppTqTKhKtCrSwfiscscftwa424209077086862ppCpTqTKhKtCrSwfiscscftwa220242420907708686.(lg..)ppCpTqTKhtttwsiscscfp2204242.lgppCpTqTKhKtCrSwswfscscftwa220242420907708686.(lg..)24压力平方法(p)P2压降曲线霍纳曲线ptttwsp2lgMDH曲线ptws2lgptwf2lgptwflgptttwsplgptwslg25KZpqTTmhZqTmhscfscf4242001467..SpphmKCraiwftw1151109077222.()lg.SphpmKCrawswftw1151109077222.()lg.26在图版拟合分析中，无因次压力定义为：pKhqZTTppDscfsc00271432.()KqZThppfD784892.()拟合27压力导数拟合：dpdttdpdttttttpp()()()22布德图版KqZThptCdpdttfDDD784892.(/)()/'拟合KqZThptCdpdtttttfDDDpp784892.(/)()//'拟合28２、(P)简化为P10203040PPw20.7MPaZPC0'()'PpZdpCPp220029ppCpTqTKhKtCrSwfiscscftwa2121090770868602.(lg..)'ppCpTqTKhtttwsiscscfp21210.lg'ppCpTqTKhKtCrSwswfscscftwa2121090770868602.(lg..)'30用压力本身进行试井解释(p)P压降曲线霍纳曲线ptttwsplgMDH曲线ptwslgptwflgptwflgptttwsplgptwslg31KZpqTpTmhZqTpmhiiscfisciifi21217335103..SpphmKCraiwftw11511090772.()lg.SphpmKCrawswftw11511090772.()lg.用半对数曲线分析进行解释：32图版拟合分析pKhqTTpZppKhqTpZpDscfiiiscfiii0054286157.KqZTphppiifiD157()拟合格林加坦图版拟合布德图版拟合KqZTphptCdpdttiifiDDD157'(/)[()/]拟合KqZTphptCdpdtttttiifiDDDpp157'(/)[()/()]()/拟合