基于3D Delaunay剖分算法的重力建模与分析.pdfVIP

第55卷第7期 地 球 物 理 学 报 v01．55，N。．7 2012年7月 CHINESEJOURNALOFGEOPHYSICS Jul．，2012 2267，doi：10．6038／j． 李振海，罗志才，钟波．基于3DDelaunay剖分算法的重力建模与分析．地球物理学报，2012，55(7)：2259 issn．0001-5733．2012．07．012． and basedon3D I，iZH。LuoZ analyzing Delaunaytriangulationalgorithm．ChineseJ． C，ZhongB．Gravitymodeling Chinese)，2012，55(7)doi：10．6038／j．issn．000卜5733．2012．07．012． Geophys．(in 基于3D Delaunay剖分算法的重力建模与分析 李振海1，罗志才1’2’3。，钟波1’2 1武汉大学测绘学院，武汉430079 2武汉大学地球空间环境与大地测量教育部重点实验室，武汉430079 3武汉大学测绘遥感信息工程国家重点实验室，武汉430079 摘要地质体的重力建模是正确解释和应用重力资料的关键问题之一．针对非规则形状变密度的三度体，本文 Delaunay剖分算法将三维目标地质体分解为若干变 提出了基于3DDelaunay剖分算法的重力建模方法．采用3D 密度四面体体元，推导了基于四面体体元的重力正演公式，建立了剩余密度值与重力异常值的线性方程组；以变密 度的长方体和倾斜台阶组合体为例，比较分析了常规块体算法和3DDelaunay剖分算法应用于重力正演的有效性， 并采用共轭梯度法加密度约束条件对非规则形状变密度的倾斜台阶组合体进行了密度反演．计算结果验证了本文 方法的正确性和有效性．基于3DDelaunay剖分算法的重力建模可应用于存在褶皱、断层、裂缝等复杂地质体的重 力正反演计算， 关键词3DDelaunay剖分算法，重力建模，三度体，重力异常 中图分类号P312 收稿日期2011-01—24，201205—15收修定稿 doi：10．6038／j．issn．0001—5733．2012．07．012 and basedon3D Gravitymodelinganalyzing Delaunaytriangulationalgorithm Zhi—Cail’2”。ZHONGB01’2 LIZhen—Hail．LU0 1School andGeomatics，Wuhan 430079，China ofGeodesy University，Wuhan 2 Environmentand Education， KeyLaboratoryofGeospace Geodesy，Ministryof Wuhan 430079，China University，Wuhan in andRemote 3State Survey