复杂重磁异常多阶段决策的最优化反演

X()(),,,11,2002,37(4):349353,,,,,,ABSTRACTLiuTianyou,ZhuLixin,DengRonglaiandLiQinghao.Optimuminversionofsophisticatedgravitationalandmagneticanomaliesbymulti-stagedecision.OGP,2002,37(4):349353Theunconvergenceoftenoccurredduringtheinversionofsophisticatedgravitationalandmagneticanomaliesproducedbyseveralgeologicbodieswhenconventionaldampedleastsquareinversionorgeneralizedinversematrixinversionhasbeenused.Usingthemethodofmulti2stagedecisionbasedonnon2linearprogrammingtodomulti2scalewaveletdecompositionofobjectivefunctionanddoingsuccessivelyinversionfromlarge2scaletosmall2scale,theredundanceofobjectivefunctionhasbeenovercome.Thepracticeprovedthatthemethodcansolvetheoptimuminversionofsophisticatedanomalieswithmulti2peakgeneratedbymulti2geologicbodies,anditsprocessofsearchingminimumthroughdecomposingtheobjectivefunctionbymulti2scaleisreasonableandconvergence.Themethodhasbeenusedtoprocessandinterpretthedistributionofigneousrockinacertainoilfield.Keywords:gravitationalandmagneticanomalies,optimuminversion,wavelet,decompositionbymulti2scale,dynamicprogramming2070,,()(Johnson,1969;Al2chailabi,1970,1972;Kunaratnan,1972;Byax,1973;,,,,,1974,1980),,,,,20028374XLiuTianyou,DepartmentofGeophysics,ChinaUniversityofGeosciences,WuhanCity,HubeiProvince,430074,China2001711,2002528,2080,Backus2Gilbert,(Pedersen,1977,1979;,1982),30,,,(),,2090,,(,,1997;,1999)(,1998)(,1999),,,,;,,10,,[2],[3]Mallat,,,,,,,=mk=1$Tk-f(x,b1,b2,,bn)2(1):$Tk(k=1,2,,m);f(x,b1,b2,,bn)nb1,b2,,bnfb(0)i,=mk=1$Tk-f(0)(k)+5f(0)(k)5b1D1+++5f(0)(k)5bnDn2(2):bi=b(0)i+Di,b(0)i;Dibi,AD=g(3)Ann=PTPPmn=5f(0)(k)5bign1=PT($T-f(0)),(3),,,,(3)(A+KI)D0=g(4):I;K[1],,,,,,,f(x)=x2(-PxP),,f(x)=x2f(x)=P23-4n=1(-1)n+1cosnxn2=P23-4cosx12-cos2x22+cos3x32-1:n=1,f1(x)x=0;n=2,f2(x)3x=0,x=0532002P4;n=3,f3(x)5x=0,x=P4,x=P6,,,(),,,,,,,,,,,1(n)Bellman(1957,1962),Bellman,,,(1),Jacobi(2)Jacobi,(3)(4),(1)(5),;,1,,;,;L1L2:,2Harr,4,10-4,14,:(Xi,Zi),2bi,2li,isi,Msi,Ai(i=1,2)2;;;;2,,isiMsiis1=60Ms1=1000010-3Aömis2=45Ms2=1000010-3Aömis1=60Ms1=800010-3Aömis2=60Ms2=800010-3Aömis1=5812Ms1=85611610-3Aömis2=4118Ms2=94031810-3Aöm153374:,,,,,,,,200010-3Aöm3,,,,50100m,,,3;;;(1)L2,PTPD=PTY(5):D;Y,(5)PTPPTY,JacobiL1PTD=Y(6)PTMN,M,N(6)JacobiY,[7](2),,,,,,?,?4,,,4(),[1],,1().,1980,5491[2],1.,1994[3]BunksC,SaleckFM,ZaleskiSandChaaventG.Multiscaleseismicwaveforminversion.Geophysics,1995,60(5):145714732532002[4],..,1997,36(2)[5],..,1998,41(5):704711[6]..,1999,38(2):2030[7]..,1999,486517[8],.1,1986,25(4),,,,,,2001,200251521,,3390,,4,,,:,;,(),,,,,,,,;,,,,TACross,,:,,;,,AöS,,,,;,,(),,,;,,6,353374:,1966;1989,,19982001;,1975;1997(),2000;,,,1964;(),1962;1982;;,(),1956;1982,1987,1989,1997;,1941;1964;1987,1991,1993,1964;1986,1989;;,1965;1986;;,;;,,,1945;1968,1982(),;,1956;1980;,1960;1984;,10;;,,1971;(),1965;1991,;CTMPP,;30;,1962;1984;,1964;1989,;20;,1962;1982,1996;,10;,;,1963;1985,2001();;,,1978;1999(),,;10;5;IET,,1962;1982,19891998;,,30;(),20028374