inspiredmkxk.com

mkxk.com  时间:2021-04-10  阅读:()
ComputGeosci(2009)13:235–244DOI10.
1007/s10596-008-9087-9ORIGINALPAPERAniterativeensembleKalmanlterforreservoirengineeringapplicationsM.
V.
Krymskaya·R.
G.
Hanea·M.
VerlaanReceived:29November2007/Accepted:14April2008/Publishedonline:14June2008TheAuthor(s)2008AbstractThestudyhasbeenfocusedonexaminingtheusageandtheapplicabilityofensembleKalmanlteringtechniquestothehistorymatchingprocedures.
TheensembleKalmanlter(EnKF)isoftenappliednowadaystosolvingsuchaproblem.
Meanwhile,tradi-tionalEnKFrequiresassumptionofthedistribution'snormality.
Besides,itisbasedonthelinearupdateoftheanalysisequations.
Thesefactsmaycauseprob-lemswhenlterisusedinreservoirapplicationsandresultinsamplingerror.
Thesituationbecomesmoreproblematiciftheaprioriinformationonthereservoirstructureispoorandinitialguessaboutthe,e.
g.
,per-meabilityeldisfarfromtheactualone.
TheabovecircumstanceexplainsareasontoperformsomefurtherresearchconcernedwithanalyzingspecicmodicationoftheEnKF-basedapproach,namely,theiterativeEnKF(IEnKF)scheme,whichallowsrestartingtheM.
V.
Krymskaya(B)·M.
VerlaanFacultyofElectricalEngineering,MathematicsandComputerScience,DelftUniversityofTechnology,Mekelweg4,2628CDDelft,TheNetherlandse-mail:m.
v.
krymskaya@tudelft.
nlM.
Verlaane-mail:verlaan@dutita2.
twi.
tudelft.
nlR.
G.
HaneaTNOBuiltEnvironmentandGeosciences,BusinessUnitGeoEnergyandGeoInformation,TNO,Princetonlaan6,3584CBUtrecht,TheNetherlandse-mail:remus.
hanea@tno.
nlR.
G.
HaneaFacultyofCivilEngineeringandGeosciences,DelftUniversityofTechnology,Stevinweg1,2628CNDelft,TheNetherlandsprocedurewithanewinitialguessthatisclosertotheactualsolutionand,hence,requireslessimprovementbythealgorithmwhileprovidingbetterestimationoftheparameters.
ThepaperpresentssomeexamplesforwhichtheIEnKFalgorithmworksbetterthantra-ditionalEnKF.
Thealgorithmsarecomparedwhileestimatingthepermeabilityeldinrelationtothetwo-phase,two-dimensionaluidowmodel.
KeywordsReservoirengineering·Historymatching·Permeability·EnsembleKalmantler·IterativeensembleKalmanlter1IntroductionThemathematicalmodelingapproachtotheanalysisofreservoirperformancehasgainedpopularitythrough-outtheyears.
However,themodelcanbeusedtoforecastreservoirbehavioronlyifithasbeencali-bratedbeforehand.
Thecalibrationstage,called"his-torymatching"inthereservoirengineeringcontext,aimsatadjustingtheparametersofthereservoirsimu-lationmodelinsuchawaythatthecomputedvaluesofobservablevariablesatindividualwellsareconsistentwithavailablemeasurementsofthosequantities.
Asthemodelsbecomemorecomplicatedandlargerscaled,thereincreasesaneedofautomatichistorymatchingtechniques.
Themainproblemthathastobesolvedviaauto-matichistorymatchingissearchingforthecombina-tionofreservoirparametersforwhichanerrorfunction(objectivefunction)attainsitsminimum.
Thisfunc-tionrepresentsasumofsquareddifferencesbetween236ComputGeosci(2009)13:235–244theobservedreservoirperformanceandtheresultsofsimulationduringthehistoricalperiod[2]:fE=npari=1wi(xioxis)2,(1)wherefEdenoteserrorfunction,nparisthenumberofreservoirparameters,wiistheweightingcoefcient,andxioandxiscorrespondinglyrefertoobservedandsimulateddatathathavetobematched.
Historymatch-ingisusuallyanill-posedproblemsincetherearemoreunknownsthanconstraintstoresolvealltheunde-terminedquantities.
Then,theGauss–NewtonmethodappliedtominimizecostfunctionfailsbecausetheHessianisill-conditioned.
Suchaproblemcanbeover-comebyapplyingsomeregularizationstrategy(e.
g.
,byuseofthepriorgeostatisticalmodel[8]).
Thereareseveralapproachestoautomatichistorymatching,whichdifferinthewaytheyobtainthesetofparametervaluesminimizinganobjectivefunction.
Thechoiceofminimizationtechniqueismainlybasedonwhethertheerrorfunctionhaslinearornonlinearform.
However,mostofthesetraditionalhistorymatchingapproachesareeitherlimitedtothesmall-scaledandsimplereservoirmodelsorinefcientintermsofcom-putationalcosts.
Ingeneral,thesemethodsalsoper-formthetreatmentofuncertaintyviarepeatedhistorymatchingprocessesforvariousinitialmodels,whichre-sultsinevengreatercomputationalefforts.
Moreover,traditionalhistorymatchingdoesnotallowcontinu-ousmodelupdating.
Namely,asthenewdatabecomeavailableforbeingincludedintothematch,thewholehistorymatchingprocesshastoberepeatedusingallobserveddata.
Atthesametime,theamountofde-ployedsensorsforpermanentmonitoringofpressure,temperatureorowratesincreases.
Thisfactyieldstheincreaseofdataoutputfrequencyandlightensupaproblemofincorporatingobtaineddatainthemodelassoonasitbecomeavailablesothatthemodelisalwaysup-to-date.
TheKalmanlteringtechniquesareknownasthemostpopularmethodologyforassimilatingthenewmeasurementstocontinuouslyupdatethestateofthesystem.
Originally,theKalmanlterwasdevelopedforoperatingonthelinearmodels,whilenonlinearityrequiresusingsomefurthermodications,e.
g.
,theextendedKalmanlter.
However,whenthemodelishighlynonlinearorthescaleofthespacevectoristoolarge,applicationoftheextendedKalmanlteralsomeetsdifculties.
ThesedifcultiesareovercomebyapplyingtheensembleKalmanltering(EnKF)algo-rithmbasedontheMonte-Carloapproach.
Thegreatmajorityoftheproblemsinreservoiren-gineeringarehighlynonlinearandcharacterizedbyalargenumberofvariables;thus,theideatouseEnKFinreservoirsimulationseemstobenatural.
Inparticular,itispresentedinthepublications[5,10].
Otherpapers[4,7]reporttheresultsofusingtheEnKFapproachinhistorymatchingprocesses.
Theyconsidertheapplica-tionofEnKFtoaPUNQ-S3model.
AlthoughthesestudiesclearlyshowthatEnKFissuccessfulinassim-ilatingproductiondatatoupdatetheinitialreservoirmodelanditsapplicationallowsreducingcomputa-tionalcostsforhistorymatching,thereisstillenoughspaceforfurtherinvestigationandimprovement.
Specically,theresearchdescribedin[11]hasshownthat,forsomenonlinearmodels,theEnKFdoesnotprovidecompletelyacceptablecharacterizationsoftheuncertainties.
Thesituationbecomesmoreproblematiciftheaprioriinformationonthereservoirstructureispoorandtheinitialguessaboutthesystemstateisfarfromtheactualone.
ThisleadstotheideaofusingimprovedEnKFmodications,namely,iterativeEnKF(IEnKF)schemes,whichallowimprovingtheinitialensembleusedforsimulationand,hence,theresultingestimatedstatevector.
Inthispaper,weapplyaniterativealgorithmbasedonasequentialdataassimilationschemetoestimatingthepermeabilityeldinreservoirmodels.
Thisalgo-rithmisinspiredbyinvestigationsofA.
H.
JazwinskionaniterativeextendedKalmanlteringapproach[6].
TheimplementationandtheimprovementofthismethodagainsttheclassicalEnKFarepresentedwitha2Dreservoirsimulationorientedtowardsestimatingthepermeabilityeld.
2TheensembleKalmanlterKalmanlteringisapowerfultechniquedesignedforsolvingdataassimilationproblems.
ThissectionpresentsthegeneralideaofKalmanlteringinaman-nersimilarto[9]andofensembleKalmanlteringasgivenin[3].
Letusrestrictourselvestothecaseofthefollowinglinearsystem:xk+1=Fkxk+Bkuk+Gkwk,(2)zk=Mkxk+vk,(3)whereFk,Bk,Gk,Mkarematrices,kisthetimeindex,xkdenotesthestateofthesystem,ukisthesysteminput,zkisthevectorofmeasurements,wkisGaussianwhitesystemnoiseprocesswithzeromeanandcovari-ancematrixQk,andvkisGaussianwhitemeasurementnoiseprocesswithzeromeanandcovariancematrixComputGeosci(2009)13:235–244237Rk.
Theinitialstatex0isassumedtobeGaussianwithmeanx0andcovariancematrixP0.
Moreover,processesx0,wk,andvkareassumedtobeindependentfromeachother.
Vectorxk,whichcontainsinformationonthecurrentsystemstate,cannotbedirectlyobserved.
However,itispossibletomeasurezk,whichissomefunctionofxkaffectedbynoiseprocessvk.
Theideaistousetheavailablemeasurementszkforestimatingthestateofthesystemxk.
Tosolvethelteringproblem,Eqs.
2and3,wehavetodeterminetheprobabilitydensityofthestatexkconditionedonthehistoryofavailablemeasurementsz1,zl.
ItturnsoutthatthisconditionaldensityfunctionisGaussian;hence,itcanbecharacterizedbymeansofacovariancematrix.
However,fornonlinearmodeloperatorFk(whichispreciselythecaseofreser-voirengineeringapplications),suchconditionaldensityfunctioncanberepresentedbythersttwomomentsonlyapproximately.
TheEnKFhasbeenexaminedandappliedinanumberofstudiessinceitwasrstintroducedbyGeirEvensenin1994.
Thislteringapproachisrelativelyeasytoimplementandhasaffordablecomputationalcosts.
TheEnKFisbasedonarepresentationoftheprobabilitydensityofthestateestimateattimekbyanitenumberNofrandomlygeneratedsystemstatesxk,i,i=1,N.
EquationstoobtainthemeanxkandcovariancematrixPkofprobabilitydensityofstatexkattimekconditionedonthehistoryofthemeasure-mentsz1,zlviaEnKFalgorithmcanbeformulatedasfollows[3](wheresuperscriptspandustandforpredictedandupdatedsystemstates):Initializationstep:xu0,iN(x0,P0),i=1,N.
(4)Forwardstep:xpk,i=Fkxuk1,i+Bkuk+Gkwk,i,i=1,N,(5)xpk=1NNi=1xpk,i,(6)Lpk=xpk,1xpkxpk,NxpkT,(7)whereLpkdenesanapproximationofthecovari-ancematrixPpkwithrankN:Ppk=1N1LpkLpTk.
(8)Assimilationstep:Kk=LpkLpTkMTkxMkLpkLpTkMTk+(N1)Rk1,(9)whereKkistheso-called"Kalmangain"matrixdeterminingtheweightswithwhichthemeasure-mentshavetobeincorporatedintothemodelupdateoutcome,xuk,i=xpk,i+KkzkMkxpk,i+vk,i,i=1,N,(10)xk=1NNi=1xuk,i.
(11)NotethatEq.
10involvesgeneratingadditionalnoisevk,iwhileconstructingthemeasurementsetcorre-spondingtotheensemble.
Thisnoisevk,ihasthesamestatisticsasassumedfortheobservationerrors.
Theperturbedmeasurementsarenecessaryduetothefactthattheabsenceofperturbationleadstotheupdatedensemble,whichhastoolowvarianceandcausesthedivergenceofthealgorithm[1].
Actually,theforwardmodelintegrationstepwithinreservoirengineeringframeworkcanbeperformedbymakingaforwardrunofthereservoirsimulatorxpk,i=Fxuk1,i,(12)whichcanbedevelopedseparatelyandusedasablackboxinEnKFanalysis.
ItturnsoutthatparameterestimationviaEnKFisalsopossible.
Thiscanbedonebyconstructingthefollowingstatevector:x=my,whereyconsistsofdynamicvariableschangingwithtimeandmisavectorofstaticmodelparameters,whichareconstantintimeandhavetobeestimated.
Now,theKalmanlteranalysisisperformedontheaugmentedstatevector.
Theforwardstepofthealgo-rithmresultsinupdatingonlythedynamicvariableswithtimeandconservingthevaluesofstaticparame-ters.
However,attheassimilationstep,thevariablesofbothtypesaresimultaneouslyupdated,providingcorrectedestimationsofthestatevectorand,hence,modelparameters.
Themodeldescribingmultiphaseuidowinreser-voirishighlynonlinearandthenumberofvariablesincludedintostatespacevectorisverylarge,normallyatleasttwopergridblock.
AlthoughEnKFperformsfairlygoodforthiskindofproblem,itsometimesfails238ComputGeosci(2009)13:235–244toprovideappropriatecharacterizationofuncertainty.
Anexampleisgivenin[11]inrelationtothecasewhentheconditionalpdfforthereservoirmodelismulti-modal.
SuchaphenomenonresultsfromthefactthatmodelnonlinearitydestroysthenormalityofapriorandaposteriordistributionswithinKalmanlteringanalysis.
Moreover,thetaskofobtainingaccurateestimationofthestatevectorbecomesharderiftheaprioriin-formationaboutthereservoirstructureispoorandtheinitialguessabout,e.
g.
,permeabilityeldisfarfromactualone.
TheKalmanltertechniquesaredesignedinsuchawaythatinitialconditionstendtobeforgottenbythelteralgorithmasmoredataareassimilated[9].
However,itisimportantforreservoirengineeringcasesthatreasonableestimationsofthemodelparametersareobtainedbasedonthedatacollectedbeforethewa-terbreakthroughevents.
Inturn,thesedatamaynotbeenoughtoimproveapoorinitialguessabouttheeldstructure.
Insuchasituation,forsomeapplications,theideaofiteratingltergloballycanbeofhelpsinceitallowsrestartingtheprocedurewithanewinitialguessthatisclosertotheactualsolutionand,hence,requireslessimprovementbythealgorithm.
WearegoingtoconsiderthehistorymatchingviaEnKFalgorithmasthestartingpointforfurtherin-vestigations.
WecontinuewithsomeintroductionintoalternativeEnKFtechniques.
3IterativeKalmanlteringThecurrentsectionpresentstheideasofKalmanl-teringalgorithmsthat,inouropinion,canbealterna-tivelyappliedtosolvingthehistorymatchingproblem.
IterativeformsoftheKalmanlterarenotcompletelynewwithinthescopeofreservoirengineeringapplica-tions.
Thesemethodsaimatobtaininganyensemblethatprovidesimprovingtherepresentationofthestatedistribution.
Thereexistseveralapproachesinpetro-leumengineeringliterature,e.
g.
,thead-hocconrmingEnKFmethodproposedby[10].
WewouldliketoexploittheideaofJazwinskitoiteratethelterglobally[6].
AlthoughitwasoriginallysuggestedtoiteratetheextendedKalmanlter,wemodifytheapproachforthecaseofparameterestima-tionviaEnKFtechnique.
Thealgorithmlooksasfollows.
IncorporatingallavailabledataviaEnKFstartingwithx0andP0,weobtaintheestimatedvaluesofxtendandPtend,wheretenddenotestheendtimepointofdataassimilationperiod.
Ifthenumberofavailablemeasurementsissuf-cientlylarge,wecanexpectthattheestimatedmodelparametervaluemtendisclosertothe"true"onethaninitialm0.
Theestimatedmodelparametermtendnowreplacesm0andbecomesanewinitialguessforthenextglobaliteration,whichisdonebyrerunningtheEnKFbasedonthesamebunchofobservations.
Af-terwards,thisprocedurecanberepeateduntilnosuf-cientchangeinestimatedmodelparameterisobtained.
Notethat,whenrerunningthelter,wechangeonlythemeanestimatorofinitialguessaboutthemodelparameterandnotthestatisticsy0andP0.
Thisyieldsthat,inthecaseofGaussianinitialensemble,thenewinitialdistributioncanbegeneratedbyresamplingwithupdatedmeanmtend,y0Tandthesameinitialcovari-anceP0.
TheowchartofsuchanIEnKFispresentedinFig.
1,wherethedashedblockscorrespondtothestepsthatactuallyarethepartsoftheEnKFalgorithm.
Althoughthemoreeducatedchoiceofinitialguessnaturallyshouldresultinbetterestimation,thereisnoguaranteethatiterationwillconverge.
Thus,theplanistoinvestigatethefeaturesoftheaboveIEnKFtech-niqueandtocheckwhetheritindeedallowsimprovingthestatevectorestimations.
ThenextsectionoutlinesthesettingsoftheexperimentusedtotestEnKFandIEnKFperformances.
Fig.
1IEnKFalgorithmowchart00000ComputGeosci(2009)13:235–2442394ExperimentalenvironmentThestudyisaccomplishedonthebasisofatwo-dimensional,two-phaseuidowmodel.
Themodelimplementationisprovidedbytheforwardreservoirsimulatorusedasablackboxtoperformthetimeupdateinlteringalgorithm.
Themodelisappliedtoatwo-dimensionalsquaredpetroleumreservoirwithasizeof700*700mequippedwithuniformcartesiangridconsistingof21gridcellsineachdirection.
Thereservoiristakentobe2minheight;however,weassumethatallquantitiesareverticallyhomogeneous,whichallowsconsideringuidowprocessesonlyintwodimensions.
Weconsiderthewateroodingstageoftherecoveryprocess,whichisperformedthroughtheexploitationoftheinjectionwelllocatedatthecenterofthereservoirandfourproductionwellsestablishedatthecornersoftheeld.
Theinjectionwellisconstrainedbyaprescribedinjectionrateof0.
002m3/sandproductionwells—bybottomholepressureof2.
5*107Pa.
4.
1StatespacerepresentationInthecaseofreservoirengineeringapplications,thestatevectornormallyconsistsofvectorsofpressures(p)andsaturations(s)correspondingtoeachgridblock.
Toperformparameterestimation,wehavetoincludetheparameterofinterestintothespacevector.
Thestudyisfocusedonestimatingthepermeabilityeld.
Itturnsoutthatthenaturallogarithmofpermeabilityisnormallydistributed,hence,wewouldliketoaug-mentthestatevectorbythevectoroflog-permeability(logk)correspondingtoeachgridcell.
Moreover,whileoperatingonaeld,onecanmeasurethefollowingparametersatthewells:bottomholepressures(pwell),oil(qwell,o),andwater(qwell,w)owrates.
Wecanalsoincludeintoavailablemeasurementsthepressureandsaturationquantitiesatthewells.
Finally,themegastatespacevectortakesthefollowingform:x=logkpspwellqwell,oqwell,w.
Themodelparameterisconsideredasstatic,i.
e.
,time-invariant.
Meanwhile,thevalueofthestaticparameteriscorrectedwithinthedataassimilationstep.
Notethatthestatevectorconsistsof441permeabil-ityvalues,441pressurevalues,441watersaturationvalues,1observedbottomholepressureattheinjectionwell,4observedoilowrates,and4observedwaterowratesattheproductionwells,orsimply,x∈R1332.
TherelationEq.
3betweenthemodelvariablesandthemeasurementsforourexampleisassumedtobegiventhroughtrivialmeasurementmatrixM∈R19*1332withonlyzerosandonesasitselements,arrangedintheblockform:M=0M10000M10000M2withblocksM1∈R5*441andM2∈R9*9ofthefollow-ingform:ElementsofmatrixM1correspondingtotheob-servationsatthewellgridblocks[i.
e.
,elementsindexedas(1,1),(2,21),(3,221),(4,421),and(5,421)aresettoone,therestofthematrixiflledinwithzeros].
M2is,infact,anidentitymatrix.
Letusnotethat,althoughthemeasurementoperatorhasalinearform,theactualrelationbetweenthemodelstatevectorandobservablevariablesisnonlinear.
Thederivednotationonlyshiftsthesourceofnonlinearityanddoesnotvanishitseffects.
WeconsiderthemodelEq.
12asbeingperfect,whichmightseemtobenotveryrealistic.
However,suchanassumptionspeciesbetterenvironmentforinvestigatingaparticularIEnKFmethod.
Weexpectthat,inthecasewhentheensemblespreadisnotin-uencedbymodelnoise,theiterativetechniqueshavetodemonstratetheirspecicfeatures.
Onthecontrary,thevaluesofobservablevariablesareassumedtobeimprecise.
Toinitializethelter,oneneedsgeneratinginitialensemblesofonlystaticanddynamicvariablesbecausethereisnoproductiondataavailableatthestartingtime.
Sincethereservoiristypicallyinastateofequi-libriumatthetimewhenproductionstarts,theinitialdynamicvariables(i.
e.
,initialpressuresandwatersatu-rationscorrespondingtoeachgridblock)areassumedtobeperfectlyknown(withoutuncertainty).
Therefore,theyarethesameforeachensemblememberandequaltotheinitialconditionofthe"true"model(i.
e.
,p=3107PaandSw=0.
2).
Thus,attheinitialmoment,theonlypermutationscontainedintheinitialensemblearecausedbyinitialensembleofpermeabilitymodels.
Withinthestudy,wearegoingtousetheinitialpermeabilityensembleconsistingof999members.
Theensemblemeanand240ComputGeosci(2009)13:235–244Meanofpermeabilityfieldsensemble(log(m2))5101520246810121416182029.
12928.
928.
828.
728.
628.
5(a)Varianceofpermeabilityfieldsensemble((log(m2))2)5101520246810121416182000.
20.
40.
60.
81(b)Fig.
2Mean(a)andvariance(b)ofpermeabilityeldsensemblevariancearevisualizedinFig.
2,wherethetoppicturecorrespondstotheensemblemeanandthebottomimagetothevariance,respectively.
4.
2SyntheticmeasurementsgenerationTotesttheperformanceofEnKFalgorithms,wearegoingtodoaso-calledtwinexperiment.
Itrequiresthatthe"true"valuesofobservablevariablesaregen-eratedsyntheticallybyapreliminaryrunofthemodelitselfandthenoisyobservationsarethencreatedbypermutatingthetruevalueswiththemeasurementer-rornoise.
Thisprocedureensuresthatthemodelisindeedabletomatchthedata.
Thereafter,thesyntheticdataareusedintheassimilationexperiments.
Theimplementationofourin-housesimulatorprovidesthe"true"permeabilityeld(seeFig.
3),whichoriginatesfromthetrainingimageofmeanderingchannels.
Nowitispossibletogeneratesyntheticdatainitial-izingthesimulatorwithtruepermeabilityeld,gridblockpressuresp=3107Pa,andwatersaturationsSw=0.
2.
Theerrorineachobservablevariableistakentobe5%ofitsactualscale.
Thesamecovariancematrixisthenusedtorepresentthemeasurementsnoisewithindataassimilationanalysis.
4.
3MeasuresoflterperformanceWequantifythequalityofestimatingatruepermeabil-ityeldintermsofthefollowingspaceaveragedrootmeansquare(RMS)errorsattimek:RMS(logk)k=(logk)k(logk)true22dim(logk),(13)orRMS(logk)k,i=(logk)k,i(logk)true22dim(logk),(14)wheredim(logk)statesforthesizeofvectorofesti-matedparameters(i.
e.
,dim(logk)=441inourstudy),(logk)kistheestimatedvectoroflog-permeabilityaf-terthekthassimilationstep,vector(logk)k,idenotestheithupdatedensemblememberafterthekthassim-ilationstep,andvector(logk)truerepresentsthetruepermeabilityeld.
Truepermeabilityfield(log(m2))510152024681012141618203231.
53130.
53029.
52928.
528Fig.
3"True"permeabilityeldComputGeosci(2009)13:235–2442415ResultsanddiscussionAsequenceofsimulationshasbeenaccomplishedtotesttheperformanceofEnKFandIEnKFassimilationalgorithmsintheframeworkofestimatingthemodelparametersforatwo-phasetwo-dimensionaluidowmodel.
Weproceedbydescribingparticularinstancesanddiscussingtheobtainedresults.
Thestudyofeachparticularalgorithmincludessolv-ingthehistorymatchingproblemandobtainingtheestimateofmodelstaticparameter(i.
e.
,permeability).
Thelteranalysisisdonefromtimet0=0daystilltend=510days,whichensuresthatthewaterbreak-througheventoccursinnoneoftheproductionwells.
TheuseofclassicalEnKFinreservoirengineeringframeworkmeetsanimportantobstacleconcernedwithobtainingphysicallyunreasonablevaluesofthestatevariables.
Itoriginatesfromperformingdataassimila-tiononthestatevectorwithoutanyconstraintcomingfromthephysicalnatureoftheparameters.
Hence,theupdateddynamicvariablesmaybecomeunfeasibleandinconsistentwithestimatedstaticvariables.
Theauthorsof[10]proposedtoincludeoneadditionalso-called"conrming"stepintotheEnKFalgorithminordertoensurethattheupdatedstateisphysicallyplau-sible.
Theideaoftheconrmationstepisthefollowing:Startingattimemomentk1,we,atrst,performaforwardsimulatorrunuptotimekandthenadataassimilationstep.
Then,wetakeonlyrecentlyupdatedstaticmodelparametersandrunagaintheowsimu-latorfromcurrenttimektothenexttimemomentk.
ThedynamicvariablesobtainedreplacethosegotafterthemeasurementupdatestageofEnKFandbecomeaninitialguessforthenexttimeupdatestep.
Thisprocedureavoidsnonphysicalresultsofthemodeling.
Theinclusionoftheconrmationstepintotheal-gorithmresultsinalmostdoublingthecomputationaltimeduetoadditionalforwardmodelrunperensemblememberateachtimestep.
Infact,weusethecon-rmingEnKFtechniqueinsteadofclassicalensembleKalmanlteringforourinvestigations.
So,fromnowon,wemeanconrmingEnKFtechniqueundertheabbreviationofEnKF.
5.
1HistorymatchingviaEnKFItturnsoutthatEnKFfacesanimportantpracticalproblem,namely,standarddeviationoftheerrorsinthestateestimateconvergesveryslowlywiththenumberofensembles.
ThismakestheEnKFquitesensitivetothenumberofensemblemembersusedforsimulation.
Thus,itseemsreasonabletostartbyinvestigatinghowlterperformancedependsonanumberofensemble01002003004005006000.
40.
60.
811.
21.
41.
61.
82Time(days)RMSerror(log(m2))EnsemblemeanEnsemblemembersFig.
4EnKF:RMSerrorinestimatedpermeabilityvstimemembersandndsomeoptimumensemblesize.
Pre-liminaryanalysisshowsthat,inourcase,perform-ingfurtherEnKFrunsonN=60ensemblemembersseemstobeoptimal.
SinceiterativemodicationofEnKFactuallyhasthesameorigin,wenditappro-priatetousethesameensemblesizeasbeingoptimalalsoforIEnKFruns.
Letusnowpresenttheoutcomeofthedataassimi-lationprocedureaccomplishedviatheEnKFalgorithmwithrespecttotheoptimalnumberofensemblemem-bers.
Weconsiderthequalityofestimatingthemodelparameter.
Forthatpurpose,space-averagedRMSer-rorsEqs.
13and14areplottedintime(seeFig.
4).
Thesequantitiesarerelatedtothepartofensemblemeanandensemblemembers'valuescorrespondingtoevaluatedpermeability.
Thegraphdemonstratesim-provementoftheparameterestimationintherstfewdataassimilationstepsfollowedbystabilizationoftheerror,andreductionoftheuncertaintyforestimatedvalue(sincetheensemblespreadbecomesnarrower).
Thismeansthat,atthelatertimes,assimilateddatacarrylessusefulinformationonreservoirstructurethanattheearlytimes.
Indeed,weobtainapermeabilityeldresemblingthetrueone,althoughsomeunderes-timationofthevaluesintheupperrightandoveresti-mationofthevaluesinthebottomleftcornerareaoftheeldoccur(compareFig.
3andleftbottomchartonFig.
5).
Thevarianceeldisactuallyobtainedasthediagonaltermsofcovariancematrixcomputedfromthestatisticalpropertiesoftheensemble.
Thediffer-encebetweenthetoprightandthebottomrightsub-plotsinFig.
5indicatesreductionofthevarianceand,therefore,uncertaintyintheestimation.
242ComputGeosci(2009)13:235–244Meanofinitialensemble(log(m2))510152051015203231302928Varianceofinitialensemble((log(m2))2)5101520510152000.
20.
40.
60.
81Varianceofestimatedensemble((log(m2))2)5101520510152000.
20.
40.
60.
81Estimatedpermeabilityfield(log(m2))510152051015203231302928Fig.
5EnKF:Initialandestimatedpermeabilityeldsandcorre-spondingvariancesAlthoughhistorymatchingonthebasisofEnKFtechniquehasdemonstrateditsefciencyforproperestimatingmodelparameters,thereisstillspaceforimprovement.
Wemayaimatobtainingbetterrepre-sentationofreservoirheterogeneousstructures,whichinturnwillresultinincreasingthequalityofforecasts.
5.
2HistorymatchingviaIEnKFWeproceedbyrunningtheIEnKFalgorithmforthetrialexample.
Infact,weaccomplishthesecondglobaliterationoftheEnKFmethod.
SpaceaveragedRMSerrorsEqs.
13and14areplottedintime(seeFig.
6)toevaluatethequalityofestimatingthemodelparame-ter.
Thegraphdemonstratesimprovementforneitherparameterestimationnoruncertaintycharacterization,01002003004005006000.
40.
60.
811.
21.
41.
61.
82Time(days)RMSerror(log(m2))EnsemblemeanEnsemblemembersFig.
6IEnKF:RMSerrorinestimatedpermeabilityvstimewhichcanbeexpectedsincetherstEnKFiterationdoesnotprovidereducingtheparameterestimationerrorinlatertimesandactuallygivesusarelativelyaccurateestimate.
Indeed,thereisalmostnovisualdifferencebetweenpermeabilityeldsobtainedwithEnKFandIEnKFalgorithms(compareFigs.
5and7).
Considernowasituationwhenaprioriinforma-tiononthevaluesofmodelparametersisfarfromreal.
Forthatpurpose,wetaketheinitialensem-bleoflog-permeabilityeldsandshifteachmemberofitbyaddingthevector0.
5Ishift,whereshiftingvectorIshiftconsistsofonesandIshift∈R1*441.
Notethatsuchashiftdoesnotaffectthevariancestatis-tics,hence,thestructureoftheinitialensembleiskept.
Thedataassimilationisperformedfromtimet0=0daysuptotimemomenttend=510days,andthecovarianceofmeasurementerrorisscaledbythefactorof102topreventlterdivergence.
Suchpara-metersfordataassimilationallowsomereducingtheerrorinestimationofpermeabilityvaluesperformedviaIEnKF(seeRMSerrorsEqs.
13and14plottedintimeonFig.
8).
Indeed,weobtainapermeabilityeldwithastructureresemblingthetrueone,al-thoughsomeoverestimationofthevaluescorrespond-ingtolow-permeabilityareasoftheeldoccurs(seeFig.
9).
Theparametervaluescorrespondingtotheseareasare,inparticular,improvedafterglobaliteration(compareFig.
9aandb).
ThedifferencebetweenthebottomrightchartsinFig.
9aandbindicatesreduc-tionofthevarianceand,therefore,uncertaintyintheestimation.
Notethat,althoughregularlyprovidingoveresti-matedvalues,theltertendstocapturethestructureMeanofinitialensemble(log(m2))510152051015203231302928Varianceofinitialensemble((log(m2))2)5101520510152000.
20.
40.
60.
81Varianceofestimatedensemble((log(m2))2)5101520510152000.
20.
40.
60.
81Estimatedpermeabilityfield(log(m2))510152051015203231302928Fig.
7IEnKF:Initialandestimatedpermeabilityeldsandcor-respondingvariancesComputGeosci(2009)13:235–24424301002003004005000.
60.
811.
21.
41.
61.
822.
2Time(days)RMSerror(log(m2))EnsemblemeanEnsemblemembers(a)01002003004005000.
60.
811.
21.
41.
61.
822.
2Time(days)RMSerror(log(m2))EnsemblemeanEnsemblemembers(b)Fig.
8IEnKF:RMSerrorforestimatedpermeabilityvstime(shiftedinitialensemblewith0.
5Ishiftandmeasurementerrorcovariancematrix102Rareusedinexperiment).
(a)Firstiteration.
(b)Seconditerationoftruepermeabilityeld.
Thishappensbecausetheensembleofpermeabilityelds,usedforthecurrenttest,isonlytheshiftedversionofanalogousensemblepreviouslyusedforinvestigations.
Suchanensem-blecontainssomeinformationontheeldstructure,whichisnotinuencedbysimpleshiftingsinceashiftchangestheensemblemeanandnotthecovariance.
Thegiveninitialstatisticscannotbechangedbecauseitcomesfromthestatisticsofensemblepopulation.
Thus,thepossibilitiesofimprovingparameterestimationbyvaryingstatisticsofinitiallyguessedvaluesofmodelparameterare,inacertainsense,restricted.
AlthoughdemonstratingausageofIEnKFapproachforestimatingpermeabilityvalues,theperformedtestsraiseadditionalproblemstobesolved.
Thelistofsuchproblemsincludesndingcriteriatoevaluatewhethergloballteriterationisneededintherealcasewhennotruepermeabilityvaluesareavailable.
WesupposethatonemayconsidertheRMSdifferencesbetweentheparametervaluesobtainedattwosequentialdataassimilationsteps.
Anotherissueisconcernedwithde-terminingappropriateerrorstatisticsthatcanhaveagreatimpactonimprovementoftheestimationsandthenumberofglobalEnKFiterationsrequiredforthatpurpose.
Summarizing,wemayassertthathistoryMeanofinitialensemble(log(m2))510152051015203231302928Varianceofinitialensemble((log(m2))2)5101520510152000.
20.
40.
60.
81Varianceofestimatedensemble((log(m2))2)5101520510152000.
20.
40.
60.
81Estimatedpermeabilityfield(log(m2))510152051015203231302928(a)Meanofinitialensemble(log(m2))510152051015203231302928Varianceofinitialensemble((log(m2))2)5101520510152000.
20.
40.
60.
81Varianceofestimatedensemble((log(m2))2)5101520510152000.
20.
40.
60.
81Estimatedpermeabilityfield(log(m2))510152051015203231302928(b)Fig.
9IEnKF:Initialandestimatedpermeabilityeldsandcorrespondingvariances(shiftedinitialensemblewith5Ishiftandmeasurementerrorcovariancematrix102Rareusedinexperiment).
(a)Firstiteration.
(b)Seconditeration244ComputGeosci(2009)13:235–244matchingonthebasisoftheIEnKFtechniquehasdemonstrateditsefciencyforimprovingmodelpara-meterestimation.
6ConclusionThestudyhasbeenfocusedontheanalysisoftheusageandtheapplicabilityofensembleKalmanlter-ingtechniqueswithrespecttohistorymatchingprob-lems.
Followingtheideapresentedin[6],aniterativemodicationofEnKFisproposed.
TheaccomplishedcasestudyhasconrmedtheusefulnessoftheEnKFtechniqueforsolvingthehistorymatchingproblemandestimatingreservoirmodelparameter.
TheremightoccurproblemsatwhichtheEnKFalgorithmdoesnotprovideresultsofsufcientaccuracy.
AnappropriateuseoftheIEnKFmethodinsuchacasecanimprovetheestimations.
AcknowledgementsTheauthorswouldliketothankProf.
J.
D.
Jansenforprovidingthereservoirsimulatorandsyntheticreservoirstructure.
WearegratefultoProf.
A.
W.
Heeminkforusefulremarks.
Thecommentsofthetwoanonymousreviewershelpedtoclarifytheobjectiveofthepaper.
OpenAccessThisarticleisdistributedunderthetermsoftheCreativeCommonsAttributionNoncommercialLicensewhichpermitsanynoncommercialuse,distribution,andreproductioninanymedium,providedtheoriginalauthor(s)andsourcearecredited.
References1.
Burgers,G.
,Leeuwen,P.
,Evensen,G.
:AnalysisschemeintheensembleKalmanlter.
Mon.
WeatherRev.
126,1719–1724(1998)2.
Ertekin,T.
,Abou-Kassen,J.
H.
,King,G.
R.
:BasicAp-pliedReservoirSimulation.
SocietyofPetroleumEngineers,Richardson(2001)3.
Evensen,G.
:TheensembleKalmanlter:theoreticalformu-lationandpracticalimplementation.
OceanDyn.
53,343–367(2003)4.
Gu,Y.
,Oliver,D.
S.
:HistorymatchingofthePUNQ-S3reser-voirmodelusingtheensembleKalmanlter.
In:SPE89942,SPEAnnualTechnicalConferenceandExhibition,Houston,26–29September20045.
Gu,Y.
,Oliver,D.
S.
:TheensembleKalmanlterforcon-tinuousupdatingofreservoirsimulationmodels.
J.
EnergyResour.
Technol.
128,79–87(2006)6.
Jazwinski,A.
H.
:StochasticProcessesandFilteringTheory.
Academic,NewYork(1970)7.
Lorentzen,R.
J.
,Nvdal,G.
,Vallès,B.
,Berg,A.
M.
:AnalysisoftheensembleKalmanlterforestimationofpermeabilityandporosityinreservoirmodels.
In:SPE96375,SPEAnnualTechnicalConferenceandExhibition,Dallas,9–12October20058.
Ruijian,L.
,Reynolds,A.
C.
,Oliver,D.
S.
:Historymatchingofthree-phaseowproductiondata.
SPEJ.
8(4),328–340(2003)9.
Simon,D.
:Kalmanltering.
EmbeddedSyst.
Program.
14,72–79(2001)10.
Wen,X.
-H.
,Chen,W.
C.
:Real-timereservoirmodelupdat-ingusingensembleKalmanlter.
SPE92991,SPEReservoirSimulationSymposium,Houston,31January–2February200511.
Zafari,M.
,Reynolds,A.
C.
:AssessingtheuncertaintyinreservoirdescriptionandperformancepredicitionswiththeensembleKalmanlter.
SPE95750,SPEAnnualTechnicalConferenceandExhibition,Dallas,9–12October2005

易探云2核2G5M仅330元/年起,国内挂机宝云服务器,独立ip

易探云怎么样?易探云是国内一家云计算服务商家,致力香港服务器、国内外服务器租用及托管等互联网业务,目前主要地区为运作香港BGP、香港CN2、广东、北京、深圳等地区。目前,易探云推出深圳或北京地区的适合挂机和建站的云服务器,国内挂机宝云服务器(可选深圳或北京地区),独立ip;2核2G5M挂机云服务器仅330元/年起!点击进入:易探云官方网站地址易探云国内挂机宝云服务器推荐:1、国内入门型挂机云服务器...

IonSwitch:$1.75/月KVM-1GB/10G SSD/1TB/爱达荷州

IonSwitch是一家2016年成立的国外VPS主机商,部落上一次分享的信息还停留在2019年,主机商提供基于KVM架构的VPS产品,数据中心之前在美国西雅图,目前是美国爱达荷州科德阿伦(美国西北部,西接华盛顿州和俄勒冈州),为新建的自营数据中心。商家针对新数据中心运行及4号独立日提供了一个5折优惠码,优惠后最低1GB内存套餐每月仅1.75美元起。下面列出部分套餐配置信息。CPU:1core内存...

Virmach($7.2/年)特价机器发放

在八月份的时候有分享到 Virmach 暑期的促销活动有低至年付12美元的便宜VPS主机,这不开学季商家又发布五款年付VPS主机方案,而且是有可以选择七个数据中心。如果我们有需要低价年付便宜VPS主机的可以选择,且最低年付7.2美元(这款目前已经缺货)。这里需要注意的,这次发布的几款便宜年付方案,会在2021年9月30日或者2022年4月39日,分两个时间段会将INTEL CPU迁移至AMD CP...

mkxk.com为你推荐
小度商城小度智能屏Air哪里可以买?大家都怎么入手的?12306崩溃12306网站显示异常,什么原因啊梦之队官网梦之队是哪个国家的?access数据库ACCESS数据库和SQL有什么区别?杰景新特杰德特这个英雄怎么样www.119mm.comwww.kb119.com 这个网站你们能打开不?www.55125.cn如何登录www.jbjy.cn5xoy.comhttp www.05eee.comwww.bbb551.com100bbb网站怎样上不去了partnersonlinecashfiesta 该怎么使用啊~~
厦门虚拟主机 域名备案网站 site5 softlayer 国外服务器网站 isatap 腾讯云数据库 payoneer directadmin java虚拟主机 服务器硬件防火墙 申请免费空间和域名 电信宽带测速软件 石家庄服务器 美国代理服务器 服务器机柜 server2008 电信测速器在线测网速 天鹰抗ddos防火墙 主机游戏 更多