基于弛豫-扩散的二维核磁共振流体识别方法-胡法龙

552201210PETROLEUMEXPLORATIONANDDEVELOPMENTVol.39No.5文章编号：1000-0747(2012)05-0552-07基于弛豫-扩散的二维核磁共振流体识别方法胡法龙1,2，周灿灿1,2，李潮流1,2，徐红军1,2，周凤鸣3，司兆伟3（1.中国石油勘探开发研究院；2.中国石油天然气股份有限公司测井重点实验室；3.中国石油冀东油田公司）基金项目：中国石油勘探开发研究院院级项目“二维核磁共振测井流体识别新技术研究”（115201C90001）摘要：基于MRIL-Prime核磁共振测井仪器现有采集模式，将不同采集模式测井信息进行组合后获得二维核磁共振信号，利用多回波串联合反演技术获得孔隙流体弛豫-扩散的二维核磁共振信息分布，用以识别复杂储集层流体性质。相对一维核磁共振测井，该流体性质识别方法增加了扩散域流体信息，可以在二维空间内将油、气、水信号分离，提高核磁共振测井流体性质识别能力。利用MRIL-Prime仪器对南堡凹陷A井油层和B井水层进行多回波间隔的二维核磁共振测井试验，解释结果与试油结果相吻合，说明二维核磁共振测井在轻质油识别和大孔隙储集层流体识别方面相对一维核磁共振测井技术有明显优势。图6参20关键词：二维核磁共振测井；采集模式；多回波串联合反演；扩散系数中图分类号：P631.811文献标识码：AFluididentificationmethodbasedon2Ddiffusion-relaxationnuclearmagneticresonance(NMR)HuFalong1,2,ZhouCancan1,2,LiChaoliu1,2,XuHongjun1,2,ZhouFengming3,SiZhaowei3(1.PetroChinaResearchInstituteofPetroleumExplorationandDevelopment,Beijing100083,China;2.PetroChinaResearchLaboratoryofWellLogging,Beijing100083,China;3.PetroChinaJidongOilfieldCompany,Tangshan063004,China)Abstract:BasedoncurrentacquisitionmodesofMRIL-PrimeNMRloggingtool,2DNMRsignalscouldbeobtainedbythecombinationofloggingdatafromdifferentmodes,thenthefluidpropertiesincomplicatedreservoirscouldbedistinguishedby2Ddiffusion-relaxationNMRloggingdatadistributionofporefluids,generatedbymulti-echotrainjointinversion.Incomparisonwith1DNMRlogging,thismethodcouldincreasefluidinformationindiffusionregime,separateoil,gasandwatersignalsin2DspaceandenhancetheidentificationcapacityoffluidpropertiesfromNMRlogging.The2DNMRlogginginthemulti-echowaveintervalwasappliedintheoilpaysinWellAandthewaterlayersinWellBintheNanpuSagbyMRIL-Primetool,andtheinterpretationmatchesthewelltestingresult.Itindicatesthat2DNMRlogginghasadvantagesontheidentificationoflightoil,andfluidsinmacroporereservoirsthan1DNMRlogging.Keywords:2DNMRlogging;acquisitionmode;jointmulti-echotraininversion;diffusioncoefficient0[1]T2DSMSSM(TDA)(DIFAN)(EDM)[2-5]SunHürlimann[6-7][8-9]MRIL-Prime201210-553MRIL-Prime-T2DDT2[10]1T222B2S2D1111TTTT1292D1()1012DGTETG1/T2D02int2B2S111TTT22922int2D2int1111()+=+1012DGTETTTT3“”13T21a31b-T2D1c-T2D[11-12]MRIL-PrimeTWTET2图1孔隙介质中不同流体组分的二维核磁共振信息分布图2图2（T2，D）二维核磁共振测井不同回波间隔采集模式示意图554Vol.39No.5CPMG1,2int,2D,/(1/1/)11(1e)ejijlpmTWTtTTijlijlgf1,2,,in424TW2CPMGT12(,,,)nggggT2Djlf2int,2D,2(1/1/)1111emin()min2ijlpmntTTjlijlifgf2112pmjljlf52int,2D,2int,2D,(1/1/)(1/1/)111ee0ijlijlpmntTTtTTjliijljlifgf(1,2,,;in1,2,,;jm1,2,,)lp6jlfT2D2int,2D,(1/1/)1max0,eijlntTTjliifc(1,2,,;jm1,2,,)lp7672int,2D,2int,2D,(1/1/)(1/1/)11eeijlijlpmtTTtTTiiijlccg(1,2,,)in88icT1()2cMIccg99T2DjlfnT11122(,,)nngggg[13-15]3[16-20]MRIL-PrimeAMRIL-PrimeMRIL-Prime0.91.23.66.0ms12.998sT231.2ms6.0msMRIL-PrimeDTE4120.9ms3.6msMRIL-PrimeD9TWE32955～2970mT2TET2T24mPasT2444-42.5×105cm2/s62%99t201210-555图3MRIL-Prime仪器在A井采集的不同回波间隔核磁共振T2分布图4A井不同深度点的（T2，D）二维核磁共振信息分布图MRIL-Prime556Vol.39No.5A5BT2D9TWE3510564～869图5MRIL-Prime仪器在B井采集的不同回波间隔核磁共振T2分布图6B井（T2，D）二维核磁共振信息分布106T2D201210-5573164.0～3175.6m63170mT274m34MRIL-PrimeT2DT2DT2T2D2.5×105cm2/s符号注释：T2msT2BmsT2SmsT2DmsDcm2/sHz/GGG/cmTEmsT2intms%TWmsgii%nmT2pDT1,jjmsfjljT2lDT2D%tiimsT2int,jjT2msT2D,llmsεii%gααiiciicMISPmVGRAPIdhcmRLLDΩmRLLSΩmRMSFLΩmCNL%ρg/cm3tus/mEDSPECms参考文献：[1]CoatesGR,XiaoLZ,PrammerMG.NMRloggingprinciplesandapplications[M].Texas:GulfPublishingCompany,1999.[2]AkkurtR,PrammerMG,MooreMA.SelectionofoptimalacquisitionparametersforMRILlogs[C]//SPWLA37thAnnualLoggingSymposium.NewOrleans:SPWLA,1996.[3]AkkurtR.Determinationofresidualoilsaturationusingenhanceddiffusion[R].SPE49014,1998.[4],,.[J].,2001,25(5):365-368.WangZhongdong,WangHao,XiangTiande.IntegrationofNMRDTWlogsandDTElogstoimprovetheoilbedinterpretationaccuracy[J].WellLoggingTechology,2001,25(5):365-368.[5],,,.[J].,2010,37(6):703-708.LiuZhonghua,LiXia,ZhaoWenzhi,etal.Enhanceddiffusiontheoryofnuclearmagneticresonance(NMR)anditsapplicationtofluididentificationofcomplexreservoirs[J].PetroleumExplorationandDevelopment,2010,37(6):703-708.[6]SunB,DunnKJ.CoreanalysiswithtwodimensionalNMR[C]//InternationalSymposiumoftheSocietyofCoreAnalysts.Monterey:SCA,2002.[7]HürlimannMD,VenkataramananL.Quantitativemeasurementoftwo-dimensionaldistributionfunctionsofdiffusionandrelaxationingrosslyinhomogeneousfields[J].JournalofMagneticResonance,2002,157:31-42.[8],,,.[J].,2005,29(5):430-434.XieRanhong,XiaoLizhi,DengKejun,etal.Two-dimensionalNMRlogging[J].WellLoggingTechology,2005,29(5):430-434.[9],.(T2,D)[J].,2009,52(9):2410-2418.XieRanhong,XiaoLizhi.The(T2,D)NMRloggingmethodforfluidscharacterization[J].ChineseJournalofGeophysics,2009,52(9):2410-2418.[10].[J].,2007,31(5):401-407.XiaoLizhi.SomeimportantissuesforNMRloggingapplicationsinChina[J].WellLoggingTechology,2007,31(5):401-407.[11],.[J].,2006,26(5):557-560.HuFalong,XiaoLizhi.AmethodfordeterminingNMRlogsignalintensity[J].NuclearElectronics&DetectionTechnology,2006,26(5):557-560.[12],.[J].,2008,23(1):173-177.HuFalong,XiaoLizhi.NumericalsimulationofstaticmagneticfielddistributionofNMRloggingtoo