油田儲層參數(shù)反演研究現(xiàn)狀文獻綜述20世紀60年代，美國物理學家Backus和應用數(shù)學家Gilbert提出的Backus-Gilbert理論奠定了反演理論的基礎(chǔ)。1977年，Matsu′uraADDINCNKISM.Ref.{B014ED8AB173447e808CED8193CD7F22}[42]首次提出“大地測量反演”的概念。1997年MossopADDINCNKISM.Ref.{CD6C3F37BA284785868A129F8AADD576}[43]等通過GPS測量獲取了間歇泉地熱田在1977年至1996年地面沉降，首次利用Mogi源模型反演地熱儲層源位置，并探討了增加點源數(shù)量對改善反演精度的效果。1998年，Vasco等人ADDINCNKISM.Ref.{1E5DDD9787AB4cc3B3ACB2981B06F920}[44]利用格林函數(shù)給出的點源響應函數(shù)，解釋了采油引起的地表和地下變形之間的相互關(guān)系，地面對多孔彈性儲層體積變化的響應得到了很好的證實。1999年，ClaudieCarnecADDINCNKISM.Ref.{2B1B0F4E82BA44a5BC568F3D9A67117B}[45]等利用InSAR技術(shù)處理了ERS-1/2在1993年至1997年期間獲得的圖像，在CerroPrieto地熱田發(fā)現(xiàn)了緩慢的局部沉降，并假設(shè)地面沉降是有單一點源在半空間內(nèi)產(chǎn)生彈性變形的情況下產(chǎn)生的，利用Mogi模型反演到油田儲層的深度和位置。2001年，XuADDINCNKISM.Ref.{338F0CE16FD540b58F88FF15315696E7}[46]等使用InSAR技術(shù)獲取勒Belridge和LostHills油田的地表形變數(shù)據(jù)，然后使用彈性變形理論解來模擬這種力學效應，結(jié)合研究區(qū)的實測數(shù)據(jù)，發(fā)現(xiàn)模擬的表面變形與InSAR的實測值吻合，證實了可以使用衛(wèi)星測量地球變形來監(jiān)測油氣生產(chǎn)。2005年，Vasco等ADDINCNKISM.Ref.{ED34145E1E89455bB8C2BD608992F44F}[47]提出了一種從準靜態(tài)變形觀測中推斷儲層體積變化和滲透率等流動性質(zhì)的方法，接著將此方法應用到加州威爾明頓油田，證明了利用PSInSAR技術(shù)獲得的地表形變對于理解地下流體流動的適用性。2006年，F(xiàn)okker等ADDINCNKISM.Ref.{D28E74A62B24467fB63520F13F81A379}[48]結(jié)合了分析和數(shù)值方法的元素，提出了一種新的油氣開采沉陷預測半解析模型。它適用于均勻的層餅狀地層學，每層粘彈性參數(shù)變化，以及任意衰竭模式；利用沉降資料增加對于油藏水平壓實的約束，大大減少了計算量。2007年，Muntendam-Bos等ADDINCNKISM.Ref.{B813FF2F579A485a8812FC28D98D6977}[49]引入了一種新反演方案，可以通過地面沉降觀測在單一過程中求解油藏壓力降。該理論既能適應在一個或多個時期沒有對地面沉降進行估計的情況，也能適應在任意時期引入新地點的情況。并以自由水位下滲透率變化較大的氣藏為例，對該方法進行了驗證。2008年，Vasco等ADDINCNKISM.Ref.{9EBFB439B6764c38A9CA47FAF9A09E93}[4]通過對阿爾及利亞Krechba油田注入二氧化碳過程中采集的InSAR數(shù)據(jù)的應用，揭示了沿兩個西北走向走廊的壓力，并探究出了儲層一定深度的壓力變化與地表形變之間的關(guān)系。2012年，Khakim等ADDINCNKISM.Ref.{87716334C6CB4355881935A9552A0FFA}[50]測試了一種基于張性矩形位錯模型的兩步反演方法,可以從InSAR反演的地表變形中準確地估計儲層變形的分布以及儲層的長寬和走向等信息。2013年，Rocca等ADDINCNKISM.Ref.{8CB484146F5146a69A320411EA2461FF}[51]利用通過InSAR技術(shù)獲取的地表變形信息成功反演了油氣藏(油氣生產(chǎn)以及二氧化碳封存)深度的地球物理參數(shù)。2015年，YangADDINCNKISM.Ref.{9495EA324EBF4f39A6A3C3BEC70712C1}[52]等采用SBAS技術(shù)發(fā)現(xiàn)美國西德克薩斯州某油田在2007年1月至2011年3月期間有一個明顯的地表抬升區(qū)；然后利用真實生產(chǎn)數(shù)據(jù)，利用假定的孔隙度和滲透率平均值，對儲層壓力和地表位移進行了分析模擬，發(fā)現(xiàn)地面隆起的主要原因是CO2注入。2016年，Ji等ADDINCNKISM.Ref.{BE38C31A53404ef48F424FBE46797A2A}[53]使用SBAS方法研究了新疆克拉瑪依油田黑103井區(qū)的時空演化規(guī)律，并基于Okada模型反演了其儲層幾何參數(shù)。同年，劉一霖ADDINCNKISM.Ref.{E6A85363473C4ab3A7FFE066404C5ADD}[6]利用時序InSAR技術(shù)獲取了黃河三角洲東營油田的沉降，對油田地區(qū)沉降采用Mogi模型和橢球模型進行了反演分析。2017年，孫赫等ADDINCNKISM.Ref.{38F2036B0B8644399B8841691AD52432}[54]利用永久散射體干涉(PSI)技術(shù)監(jiān)測了遼河下游平原2007年至2011年的沉降速率圖。通過考慮重疊區(qū)域上的偏移量，定量分析了油田、海岸帶和鹽區(qū)三種地面沉降類型；并以油區(qū)為研究對象，基于Okada模型成功反演了油田儲層的幾何參數(shù)。2020年，Staniewicz等ADDINCNKISM.Ref.{49DD2D16DEA0477e9D850B88FE2DF5CA}[55]利用InSAR技術(shù)繪制了某二疊紀盆地8萬平方公里的產(chǎn)油區(qū)盆地表面沉降速率圖，在活躍的生產(chǎn)井和處理井附近觀察到大量的沉降和隆起特征；通過Okada模型表明，流體提取和沿正斷層的傾滑是觀測到的地震活動和地表形變的潛在原因。同年，Métois等ADDINCNKISM.Ref.{1C92DA8EADA241d49D1D91E0B57024F9}[56]通過覆蓋阿爾巴尼亞Patos-Marinza油田若干景Sentinel-1數(shù)據(jù)的時間序列分析中監(jiān)測到該地存在一個沉降漏斗，然后利用平面負拉伸位錯模擬盆地和儲層壓實作用的彈性模型可以很好地模擬了這種變形?？傮w上，國內(nèi)外對油田儲層參數(shù)反演的研究較少，并缺乏對不同反演模型的比較分析；現(xiàn)有反演模型主要為Mogi模型、橢球模型和Okada模型等，這些模型目前存在兩個問題：(1)旋轉(zhuǎn)自由度不高，不能充分反應儲層真實參數(shù)；(2)數(shù)值計算不穩(wěn)定。因此，本文以史南油田作為研究對象，首先利用時序InSAR技術(shù)分析該地區(qū)地表形變的時空分布與演化特征，并針對油田儲層參數(shù)反演研究中存在的問題，進一步深入研究。參考文獻ADDINMendeleyBibliographyCSL_BIBLIOGRAPHYADDINCNKISM.Bib[1]胡明城.現(xiàn)代大地測量學的理論及其應用[M].北京：測繪出版社，2003年10月[2]獨知行.基于力學模式的大地測量反演理論及應用[J].測繪學報,2002,31(1):94-94.[3]張朝玉.大地測量反演若干理論問題研究[D].武漢:武漢大學,2006.[4]VascoDW,FerrettiA,NovaliF.Reservoirmonitoringandcharacterizationusingsatellitegeodeticdata:InterferometricsyntheticapertureradarobservationsfromtheKrechbafield,Algeria[J].Geophysics,2008,73(6):WA113-WA122.[5]DusseaultMB.GeomechanicalChallengesinPetroleumReservoirExploitation[J].KsceJournalofCivilEngineering,2011,15(4):669-678.[6]劉一霖.黃河三角洲地面沉降時序InSAR技術(shù)監(jiān)測與地下流體開采相關(guān)性分析[D].中國科學院研究生院(海洋研究所),2016.[7]楊崇.遼河油田地表沉降InSAR監(jiān)測及儲層參數(shù)反演[D].西南交通大學,2019.[8]宋波.(2013).黃河三角洲(東營)油氣資源開發(fā)區(qū)地面沉降監(jiān)測成果分析.城市地質(zhì),8(4):40-43.[9]劉國祥.InSAR基本原理[J].四川測繪,2004(04):187-190.[10]劉國祥,丁曉利,陳永奇,李志林,李志偉.使用衛(wèi)星雷達差分干涉技術(shù)測量香港赤臘角機場沉降場[J].科學通報,2001(14):1224-1228.[11]SherwinCW,RuinaJ,RawcliffeR.SomeEarlyDevelopmentsinSyntheticApertureRadarSystems[J].IreTransactionsonMilitaryElectronics,1962(2):111-115.[12]GrahamLC.SyntheticInterferometerRadarforTopographicMapping[J].ProceedingsoftheIeee,1974,62(6):763-768.[13]GrahamL.TopographicMappingFromInterferometricSarObservations[J].Proc.Ieee,1974,62(3):763-768.[14]GoldsteinRM,ZebkerHA,WernerCL.SatelliteRadarInterferometry:Two-dimensionalPhaseUnwrapping[J].RadioScience,1988,23(4):713-720.[15]GabrielAK,GoldsteinRM,ZebkerHA.MappingSmallElevationChangesOverLargeAreas:DifferentialRadarInterferometry[J].JournalofGeophysicalResearch:SolidEarth,1989,94(B7):9183-9191.[16]ZebkerHA,RosenPA,GoldsteinRM,etal.OntheDerivationofCoseismicDisplacementFieldsUsingDifferentialRadarInterferometry:theLandersEarthquake[J].JournalofGeophysicalResearch:SolidEarth,1994,99(B10):19617-19634.[17]MassonnetD,BrioleP,ArnaudA.DeflationofMountEtnaMonitoredBySpaceborneRadarInterferometry[J].Nature,1995,375(6532):567-570.[18]GeD,WangY,GuoX,etal.LandSubsidenceInvestigationAlongRailwayUsingPermanentScatterersSarInterferometry[C]//IGARSS2008-2008IEEEInternationalGeoscienceandRemoteSensingSymposium.IEEE,2008,2:II-1235-II-1238.[19]唐嘉,劉國祥,宋云帆,等.PALSAR和ASARPSI顯著地表沉降探測與分析[J].遙感學報,2015,19(6):1019-1029.[20]WangX,LiuL,ZhaoL,etal.MappingandInventoryingActiveRockGlaciersintheNorthernTienShanofChinaUsingSatelliteSarInterferometry[J].TheCryosphere,2017,11(2):997-1014.[21]WangX,LiuG,YuB,etal.AnIntegratedMethodBasedonDinsar,MaiandDisplacementGradientTensorforMappingthe3dCoseismicDeformationFieldRelatedtothe2011TarlayEarthquake(myanmar)[J].RemoteSensingofEnvironment,2015,170(2):388-404.[22]DaiK,LiZ,TomásR,etal.MonitoringActivityattheDaguangbaoMega-landslide(china)UsingSentinel-1TopsTimeSeriesInterferometry[J].RemoteSensingofEnvironment,2016,186(5):501-513.[23]SandwellDT,PriceEJ.PhaseGradientApproachtoStackingInterferograms[J].JournalofGeophysicalResearch:SolidEarth,1998,103(B12):30183-30204.[24]FerrettiA,PratiC,RoccaF.PermanentscatterersinSARinterferometry[C]//IEEE1999InternationalGeoscienceandRemoteSensingSymposium.IGARSS'99(Cat.No.99CH36293).IEEE,1999,3:1528-1530.[25]FerrettiA,PratiC,RoccaF.PermanentScatterersinSarInterferometry[J].IeeeTransactionsonGeoscienceandRemoteSensing,2001,39(1):8-20.[26]DehlsJF,BasilicoM,ColesantiC.GrounddeformationmonitoringintheRanafjordareaofNorwaybymeansofthepermanentscattererstechnique[C]//IEEEInternationalGeoscienceandRemoteSensingSymposium.IEEE,2002,1:203-207.[27]BerardinoP,FornaroG,LanariR,etal.ANewAlgorithmforSurfaceDeformationMonitoringBasedonSmallBaselineDifferentialSarInterferograms[J].IeeeTransactionsonGeoscienceandRemoteSensing,2002,40(11):2375-2383.[28]HooperA,ZebkerH,SegallP,etal.AnewmethodformeasuringdeformationonvolcanoesandothernaturalterrainsusingInSARpersistentscatterers[J].Geophysicalresearchletters,2004,31(23).[29]李德仁,廖明生,王艷.永久散射體雷達干涉測量技術(shù)[J].武漢大學學報(信息科學版),2004(08):664-668.[30]LauknesTR,DehlsJ,LarsenY,etal.AcomparisonofSBASandPSERSInSARforsubsidencemonitoringinOslo,Norway[C]//Fringe2005Workshop.2006,610.[31]劉國祥,陳強,丁曉利.基于雷達干涉永久散射體網(wǎng)絡(luò)探測地表形變的算法與實驗結(jié)果[J].測繪學報,2007(1):13-18.[32]CasuF.TheSmallBaselineSubsetTechnique:PerformanceAssessmentandNewDevelopmentsfoSurfaceDeformationAnalysisofVeryExtendAreas[J].XXICycle,2009.[33]FerrettiA,FumagalliA,NovaliF,etal.ANewAlgorithmforProcessingInterferometricData-stacks:Squeesar[J].IeeeTransactionsonGeoscience&RemoteSensing,2011,49(9):3460-3470.[34]LeiZ,DingX,ZhongL.GroundSettlementMonitoringBasedonTemporarilyCoherentPointsBetweenTwoSarAcquisitions[J].IsprsJournalofPhotogrammetry&RemoteSensing,2011,66(1):146-152.[35]劉國祥,張瑞,李陶,等.基于多衛(wèi)星平臺永久散射體雷達干涉提取三維地表形變速度場[J].地球物理學報,2012,55(8):2598-2610.[36]ZhouW,ChenG,LiS,etal.InSARApplicationinDetectionofOilfieldSubsidenceonAlaskaNorthSlope[C]//GoldenRocks2006,The41stUSSymposiumonRockMechanics(USRMS).AmericanRockMechanicsAssociation,2006.[37]FilatovA,YevtyushkinA.DetectionofgroundsurfacedisplacementsinareaofintensiveoilandgasproductionbyInSARdata[C]//ProceedingsoftheESALivingPlanetSymposium.Bergen.Norway.ESASP-686.2010.[38]張金芝.基于InSAR時序分析技術(shù)的現(xiàn)代黃河三角洲地面沉降監(jiān)測及典型影響因子分析[D].中國科學院研究生院(海洋研究所),2015.[39]ShiJ,YangH,PengJ,etal.InsarMonitoringandAnalysisofGroundDeformationDuetoFluidOrGasInjectioninFengchengOilField,Xinjiang,China[J].JournaloftheIndianSocietyofRemoteSensing,2019,47(3):455-466.[40]李乃一,伍吉倉.基于時序InSAR技術(shù)監(jiān)測勝利油田地表沉降[J].工程勘察,2018,46(5):50-54.[41]MirzaiiZ,HasanlouM,Samieie-EsfahanyS,etal.LandSubsidenceMonitoringinAzarOilFieldBasedonTimeSeries