用常规的阻尼最小二乘法和广义逆矩阵法反演由多个地质体产生的复杂重磁异常, 往往不收敛。采用非线性规划的多阶段决策方法对目标函数作小波多尺度分解, 由大尺度到小尺度逐次反演, 消除了目标函数的冗余度, 并克服了目标函数的复杂性。

2002年8月　　　　　　　　　　　　　　石油地球物理勘探　　　　　　　　　　　　　第37卷　第4期

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

刘天佑

Ξ

　朱立新

(中国地质大学地球物理系 武汉

　　)

邓荣来　李庆浩

(河南石油勘探局)

摘　　　要

刘天佑,朱立新,邓荣来,李庆浩1复杂重磁异常多阶段决策的最优化反演1石油地球物理勘探,2002,37(4):

349～353

用常规的阻尼最小二乘法和广义逆矩阵法反演由多个地质体产生的复杂重磁异常,往往不收敛。采用非线性规划的多阶段决策方法对目标函数作小波多尺度分解,由大尺度到小尺度逐次反演,消除了目标函数的冗余度,并克服了目标函数的复杂性。实践表明,多阶段决策方法能解决由多个地质体产生的复杂的多峰值异常的最优化反演问题,并且通过多尺度分解目标函数来搜索极小的过程是合理和收敛的。运用这一方法成功地处理和解释了某油田火成岩体分布问题。

关键词　重磁异常　最优化反演　小波　多尺度分解　动态规划

LiuTianyou,ZhuLixin,i.muminversionofsophisticated

～353gravitalstagedecision.OGP,2002,37(4):349

Theoftenoccurredduringtheinversionofsophisticatedgravitationaland

ingthemethodofmulti2stagedecisionbasedonnon2linearprogrammingtodomulti2scalewaveletdecompositionofobjectivefunctionanddoingsuccessivelyinversionfromlarge2scaletosmall2scale,theredundanceofobjectivefunctionhasbeenovercome.Thepracticeprovedthatthemethodcansolvetheoptimuminversionofsophisticatedanomalieswithmulti2peakgeneratedbymulti2geologicbodies,anditsprocessofsearchingminimumthroughdecomposingtheobjectivefunctionbymulti2scaleisreasonableandconvergence.Themethodhasbeenusedtoprocessandinterpretthedistributionofigneousrockinacertainoilfield.

Keywords:gravitationalandmagneticanomalies,optimuminversion,wavelet,decomposition

bymulti2scale,dynamicprogramming

Byлax,1973;顾学新,陈亚浙,周熙襄,王家林,沈薇英等,1974,1980),地质模型常常采用二度板状体、

引　　　言

20世纪70年代初,是电子计算机在地学领域应用的早期阶段,人们应用非线性反演方法(如高斯法、最速下降法、阻尼最小二乘法)来反演和解释一

些简单孤立的地质体引起的重磁异常(Johnson,1969;Al2chailabi,1970,1972;Kunaratnan,1972;

三度直立棱柱体及其组合模型。由于模型体几何参数和物性参数与重磁异常之间是非线性函数关系,因此必须采用线性化方法,再通过多次迭代求解反问题。这类方法尽管有完善的最优化方法作基础,但实际应用中的地质效果不理想。如初值选择不好,迭代过程可能停顿在目标函数的次极值上而不收敛。

ΞLiuTianyou,DepartmentofGeophysics,ChinaUniversityofGeosciences,WuhanCity,HubeiProvince,430074,China

本文于2001年7月11日收到,修改稿于2002年5月28日收到。