geolog6.7.1學習材料geolog教程-nmr processing

1、200310200310數(shù)據(jù)處理流常規(guī)處理步相關支持文CMR數(shù)據(jù)處理流常規(guī)處理步相關支持文CMR1.2.2asMRILC*1.2.3Haliburton/Numar CMR200QC孔隙度刻非完全極化校等待時間校2.1.4 信號相位檢3CMR測井資料的處3. 1處理流程說相位反相位角平3.1.3濾數(shù)據(jù)整轉換為孔隙孔隙度刻3.2.1方法3.2.2方法相位角計相位角平相位角翻回波串濾3.6.2Geolog1 數(shù)據(jù)處理流toT2 Geolog1 數(shù)據(jù)處理流toT2 QC QC 1.2.1CMR#CMRQC#Thisfilespecifiesthenameoftheconstantsandlogsre

2、quired#CMRQCand#ThelimitsandstandarddeviationsoftheQCcurvesareusedfor#QCsingoftheCMR#QC Curves -GeometricmeanT2 0-1.2.1CMR#CMRQC#Thisfilespecifiesthenameoftheconstantsandlogsrequired#CMRQCand#ThelimitsandstandarddeviationsoftheQCcurvesareusedfor#QCsingoftheCMR#QC Curves -GeometricmeanT2 0-(CMRP)-0-(

3、CMRP)-0-Polarizationcorrectionapp to porosity(p.u)-1-Polarizationcorrection-CMRsystem-DeltaB0field1-CMRsignal:noise-1-StandarddeviationofCMR-CMRfree-Standarddeviation-Signalphase(checkedinhighporosity-Gammasmoothingofthe0-00-Waitflagplete#Raw Data EchosignalfromRcomponent EchosignalfromX#CMRsedGeolo

4、g_NameCMR porosity Boundfluidvolume CMR free fluid CMR Permeability Geometric mean T2#CMRT2DistributionlitudedistributionStandarddeviationoflitude#CMRinimumboundfluidheequationforpermeability umallowedCMRumvalueforpolarizationusedtotriggertheinsufficientwaittimestackingbytheMAXIS) filteringlevelHydr

5、ogen index of the formation fluid CoefficientAfortheSDRpermeabilityCoefficientBfortheSDRinimumboundfluidheequationforpermeability umallowedCMRumvalueforpolarizationusedtotriggertheinsufficientwaittimestackingbytheMAXIS) filteringlevelHydrogen index of the formation fluid CoefficientAfortheSDRpermeab

6、ilityCoefficientBfortheSDRCoefficientCfortheSDRCoefficientAfortheCoefficientBfortheTimurCoefficientCfortheTimurT1:T2ratiosetforsignalumT1:T2ratiocalculationofCMRinimumT1:T2ratiocalculationofCMRLogarithmicallyevenlydistributedvalues usedtoplottheT2T2LMthresholdtotriggerthebad-umvalueforT2singoftheT2s

7、ingoftheT2inimumvalueforT2T2cutoffT2cutoffT2cutoffT2cutoffT2cutoffT2cutoffT2cutofflowsignalzones T2orboundT2ofthemud-filtrateusedtoformud-filtratesignalatshortT2Timenthe90 xand180ypulseswhichequaltoonehalftheechoNumberofheacquisitioncycle NumberofstacksRealtimehydrogenindexcorrectionforoilPermeabili

8、tyformulaePorositythresholdtotriggerthebad-holeflag Repetition number for each wait-timeUsedtosetthewaittimeofeach#End#MRILQCasMRILC*#Thisfilespecifiesthenameoftheconstantsandlogsrequiredfor # MRIL QC and PRO#ThelimitsandstandarddeviationsoftheQCcurvesareusedfor#QCsingoftheMRIL#QC-erecho22-22erecho-

9、Ringing-Offset-Tool-Minimumhigh-umhighvoltage -Unregulatedhigh-Instrument -erecho22-22erecho-Ringing-Offset-Tool-Minimumhigh-umhighvoltage -Unregulatedhigh-Instrument -Instrument5-Instrument5-Instrument02-Goodnessoffit#Raw Data Qualifiers minimumCPMGchannelleveldataCPMGchannelleveldataCPMGchannellev

10、eldataCPMGchannelleveldataEchocurvesforCPMGdatafromEX,EY#Aliases for Raw Data AliasE1R= AliasE2R= Alias EX = Alias EY =AliasECHO=EC#QualifiersforRawDataQualifiersare heformxyz(e.gE1Rxyz).ZinnotxStandardT2dataorlongerTW(waittime)inDTWmodeshorterTE(echospacing)inDTETheshorterTW(waittime)echodatainDTWT

11、he longer TE (echo spacing data) in DThe 0.6ms TE short echo data in C/TP logging modeDepthbasedcurves(channellevelannotQTimebasedcurves(channellevelFrequency#MRILsedlitudeofEffectiveporosity(MBI+MBVM) Bulk volume irreducibleBulk volume movable VolumeofclayboundwaterApparent#MRILT2DistributionArrays

12、 T2distributionfromECH(longerTW,shorterT2distributionfromECHB(shorterTW,shorter#End#MRILQC#ThisfilelitudeofEffectiveporosity(MBI+MBVM) Bulk volume irreducibleBulk volume movable VolumeofclayboundwaterApparent#MRILT2DistributionArrays T2distributionfromECH(longerTW,shorterT2distributionfromECHB(short

13、erTW,shorter#End#MRILQC#Thisfilespecifiesthenameoftheconstantsandlogsrequiredfor # Haliburton . Numar MRIL QC and PRO#ThelimitsandstandarddeviationsoftheQCcurvesareusedfor#QCsingoftheMRIL#0#Check field AllowablevariationB1magneticfield#Noise Check Standarddeviationofambient00-MeanofambientStandardde

14、viationoferechonoise Meanofambient0808-Standarddeviationofambient0505erecho-Meanofambient#RINGINGQC-#RINGING_QC-#QC Curves -LowvoltagesensorUr15high -Lowvoltagesensor-Lowvoltagesensor-Lowvoltagesensor15LowvoltagesensorUr15-LowvoltagesensorStandarddeviationofambient0505erecho-Meanofambient#RINGINGQC-

15、#RINGING_QC-#QC Curves -LowvoltagesensorUr15high -Lowvoltagesensor-Lowvoltagesensor-Lowvoltagesensor15LowvoltagesensorUr15-Lowvoltagesensorunreg5v -Lowvoltagesensor-Lowvoltagesensor-5 -Lowvoltagesensor+5Minimumhighumhigh-B1-MSIG-0MPHIA- -2-MPHI-#Raw Data ImaginaryXRealechoes,withrunningaverageechoes

16、,withrunningaverageXXer,salinityand Realechoes.AVR afterapplyingphasecorrection,correctedechoes.AVIafterapplyingphasecorrection,correctedfor er,salinityaveragedphaseangle. XEstimatedXdeviationofnoiseestimatedfromXMeannoiseafterrunningaverageandphaseFiterror(qualityofXX#Aliases for Raw Data = REA = I

17、MA =REALC= IMAGC= PHNO= PHER= SNR = CHI =#QualifiersforRawDataQualifiersheformABCDxyzwherexyzthedenotedbyasingleletter(e.gAVRA=AVRqualifierA) # y z not defined for the MRILNoacquisitiondata or short TE data.ABShortTWorlongTECPartialrecoverydataforclay(small#MRILsed0.5ms-CGroup 1.0ms-CGroup 2.0ms-CGr

18、oup 4.0ms-CGroup 4.0ms-BGroup 8.0ms-BGroup 16ms-BGroup 32ms-BGroup 64ms-BGroup Fiterror(qualityofXX#Aliases for Raw Data = REA = IMA =REALC= IMAGC= PHNO= PHER= SNR = CHI =#QualifiersforRawDataQualifiersheformABCDxyzwherexyzthedenotedbyasingleletter(e.gAVRA=AVRqualifierA) # y z not defined for the MR

19、ILNoacquisitiondata or short TE data.ABShortTWorlongTECPartialrecoverydataforclay(small#MRILsed0.5ms-CGroup 1.0ms-CGroup 2.0ms-CGroup 4.0ms-CGroup 4.0ms-BGroup 8.0ms-BGroup 16ms-BGroup 32ms-BGroup 64ms-BGroup 128ms-BGroup 256ms-BGroup 512ms-BGroup 1024ms-BGroup 2048ms-BFiltered2BinFiltered2BinFilter

20、ed2BinFiltered2BinClay-boundWater(PR06)(0.5-ms+1-ms+2-2222222222SumTotalPorosity-PR06(0.5-ms+1-msSumTotalPorosity(P4PR+P2B+.ms+4-msBins)+BGroup(8-ms+.(4-2048msBins) 4.0ms-AGroup 8.0ms-AGroup 16ms-AGroup 32ms-AGroup 64ms-A2Bin2Bin2Bin2Bin2Bin128ms-AGroup 256ms-AGroup 512ms-AGroup 1024ms-AGroup 2048ms

21、-AGroup Filtered2BinFiltered2BinFiltered2BinFiltered2BinFiltered2BinSumTotalPorosity-PR06(0.5-ms+1-msEffectivePorosity(P4PR+P2A+.ms+4-msBins)+A128ms-AGroup 256ms-AGroup 512ms-AGroup 1024ms-AGroup 2048ms-AGroup Filtered2BinFiltered2BinFiltered2BinFiltered2BinFiltered2BinSumTotalPorosity-PR06(0.5-ms+1

22、-msEffectivePorosity(P4PR+P2A+.ms+4-msBins)+AGroup(8-ms+.SplicedAGroup(4-2048-ms) (P4PR+P2A+.+P10A) (P4PR+P2A+.+P10A) Cut-VI-SplicedASpectralBVI-SplicedAumBVIumofCBVITA,SplicedAGroup(P4PR+P2A+.PermeabilityEstimate(CoatesMPHITA)-SplicedAGroup(P4PR+P2A+.PermeabilityEstimateMethod)-SplicedAGroup(P4PR+P

23、2A+.RaycombinedwithMRIL CalibratedFilteredFormationWaterResistivityatMud-FiltrateResistivityatFormationFinalTotalPorosity(dependsonEffectivePorosity(claycorrected) MRIAN MRIL Total PorosityMRIAN MRIL Effective Porosity MRIANMRILBulkVolumeIrreducibleEstimatedWforDualWaterM (TPORxSWT) (CBVWT-CLPOR) Cl

24、ayBoundWater Total Porosity TotalBulkWaterVolumeffective BulkWater VolumeTotalBoundIrreducible(MBVITA+Total Water Saturation EffectiveWaterSaturation(CBVTotalIrreducibleWaterTotalFreeFluidEquivalentConductivityBoundWaterSaturation(afterPermeability Estimate (using MPHI) Permeability to Water (using

25、MPHI) PermeabilitytoHydrocarbons(usingMPHI) Permeability Estimate (using TPOR) Permeability to Water (using TPOR) PermeabilitytoHydrocarbons(usingTPOR)putedDensityPorosity(fromputedNeutronPorosityputedSonicPorosity(fromPHIN - RHOB x-plot porosity Porosity Estimate fas Phase Porosity Estimate for Oil

26、 Phase Porosity Estimate for WatTDA Estimate of Bulk WaterTDAEstimateofBulkOil TDAEstimateofBulkPolarizationandHydrogenIndexPolarizationandHydrogenIndexEffectivePorosityrmeabilityEstimate(Coates DrinsicRelaxationT2Water Limit-Long Te BulkFormationFluidD)(forx-plot) BulkWaterrmeabilityEstimate(Coates

27、 DrinsicRelaxationT2Water Limit-Long Te BulkFormationFluidD)(forx-plot) BulkWaterD/DWInverseWater Saturation (from DIFAN) BulkVolumeWater(fromDIFAN)RelativePermeabilitytoOil(fromDIFAN) Relative Permeability to Water (from#MRILT2DistributionArrays #EndT2distributionfromECHOA(longerTW,shorterTE) T2 di

28、stribution from ECHOB (shorter TW, shorterT2distributionfrompartial2 CMR200 質量2.1CMR200QCCMR_PHI_CONV=(constant_calib*_MC* *(CMR_TEMP+CMR_PHI_CONV=(constant_calib*_MC* *(CMR_TEMP+CMR_PHI_CONV=CMRporosityconstant_claib=constantofcalibrationfor_MC=litudemasterB0_MC=B0masterFREQ_OP=OperatingCMR_TEMP=Do

29、wnholeHI=Hydrogen =Hydrogen_MC TemperatureduringmastercalibrationTestlooplitudeduringmastercalibration(v). litude during master calibration.=InverselitudeduringmastercalibrationSicfieldstrengthduringthemasterTemperatureduringloggingTestlitudeduringSicfieldstrengthduring=Hydrogen_MC Temperatureduring

30、mastercalibrationTestlooplitudeduringmastercalibration(v). litude during master calibration.=InverselitudeduringmastercalibrationSicfieldstrengthduringthemasterTemperatureduringloggingTestlitudeduringSicfieldstrengthduring21 DH BHIMC DH MC =CMRporosity,calibratedandcorrected=litudecomputeddownholeMa

31、ster=litudeof litudeicthetest theloopsignalduringcalibration=loop signalduringduringthemastercalibration. during logging.S=Sic=TemperatureofthewaterbottleduringcalibrationTemperatureduringlogging=Hydrogenindexofformation=NumberofT2heinputT2Input =litudeT1T2_MIN=MinimumT1T2ratio(defaultT1T2_MAXumT1T2

32、ratio(default=Waittimelog(polarizationtime)fortheCMRloggingpassT1T2_MIN=MinimumT1T2ratio(defaultT1T2_MAXumT1T2ratio(default=Waittimelog(polarizationtime)fortheCMRloggingpass=T2valueof binheT2distribution=T2valueof thelastbinheT2distribution=Numberof T2heinputT2CMR 正方法采用改變 T1:T2 比值的方法。非完全極化校正采用如下算法：A

33、(corrected) Ai(measuredi W 1RT2i litudeofthe=binoftheT2distributioncorrectedpletePolarizationlitudeoftheithbinofthe=Waittime =T1:T2T2 relaxation of the ithbin ofthe T2=litudeofthe=binoftheT2distributioncorrectedpletePolarizationlitudeoftheithbinofthe=Waittime =T1:T2T2 relaxation of the ithbin ofthe

34、T2=CMR_PHI_CONV=CMRporosityfactor,calculatedusinghequalitycontrolflowPorosityCalibration T2_CORRECT_MINT1:T2T2_CORRECT_STDT1:T2T2_CORRECT_MAXT1:T2moduleminimumvalueforstandardvalueforumvalueforERRMINUSCMRPCMRP_MIN ERRPLUS = CMRP - CMRP_MAX=umallowableporositydeviation A (measuredA(corrected) i=umall

35、owableporositydeviation A (measuredA(corrected) ii =litudeoftheithbinoftheT2Polarizationlitudeoftheithbinofthe=Waittime(s) T1:T2 ratioT2 relaxation of the ithbin ofthe T2Shalecut-offusedtoitableforcheckingthephase.Wherevshaleisthe checksnthecut-offervalexcluded=Porositycut-offusedtoShalecut-offusedt

36、oitableforcheckingthephase.Wherevshaleisthe checksnthecut-offervalexcluded=Porositycut-offusedtoidentifyforthesignalexcludedfromWhereporosityischecksnthecut-off erval=10o). umallowabledeviationforthesignalphase(defaultsignalphaseexceedDEV_CUTandtheformationiscleanandporouslimitssetbyandPHI_CUT)theSP

37、HASE_THRflagissetto Input =Volumeof=Output =SignalphaseoftheCMR=Signalphasethresholdtthesignalphasehasumallowabledeviationincleanandtan1/M 上述公式中，Mx My tan1/M 上述公式中，Mx My 在多孔介質中，信號的相位應保持相對穩(wěn)定( 100)。在泥質含量較高、孔隙度較低的3 CMR 測井資料的31Xk tan1 k Ri= Phase =ithechooftheechok=numberofechoestobehephasei= Phase =ith

38、echooftheechok=numberofechoestobehephaseangleR&X=InphaseandQuadraturedetectedcomponentofthesignali (Ri cos) (Xi sinnoisei(Ri sin) (Xi cossignali=signaloftheith noisei=noiseoftheRi=inphasecomponentXi=quadraturetheith of theith 用理。 T23.23.2 =Hydrogen=Temperatureduringmaster=Hydrogen=Temperatureduringm

39、astercalibrationlitudeduringmastercalibrationTestlitudeduringmaster=Inverselitudeduringmastercalibration SicfieldstrengthduringthemasterTemperatureduringloggingTestlitudeduringSicfieldstrengthduring2=HydrogenCONSTANT_CALIB=Tooldependentcalibrationconstant(default2=HydrogenCONSTANT_CALIB=Tooldependentcalibrationconstant(default =0.