基于Delaunay剖分的二维非规则重力建模及重力计算  

2-D irregular gravity modeling and computation of gravity based on Delaunay triangulation

作　　者：张岭[1] 郝天珧[1] 

机构地区：[1]中国科学院地质与地球物理研究所,北京100029

出　　处：《地球物理学报》2006年第3期877-884,共8页Chinese Journal of Geophysics

基　　金：中国科学院知识创新工程重大项目(KZCX1-SW-18-01);863探索性项目(2004AA615010)联合资助

摘　　要：针对计算截面为任意形状的非均匀密度二度体的重力问题,提出了二维情况下的非规则重力建模方法,即利用二维Delaunay剖分方法,将其截面分割为若干三角形,将该二度体分解成为若干三角棱柱的组合;通过求取变密度的三角棱柱体的重力异常,并将其累加,最后得出了截面为任意形状的水平无限延伸的二度体的重力计算公式.以千米桥油田的一条实测重力剖面为例,用该方法进行了反演计算,并将反演结果与常规方法做了分析比较,进行了相应的分析和讨论.Aiming at calculating gravity anomaly of 2-D complicated and heterogeneous geological body, we propose Irregular Gravity Modeling （IGM） approach based on Delaunay triangulation. As we know, the most fundamental cell is triangle in 2-D space. With Delaunay triangulation, we can divide a two-dimensional body, extending horizontally with arbitrary cross-section, into various triangular prisms with varying density in optimal way. After inversing the density of those nodes, we can figure out the density distribution. Using analytic equations derived, our approach is convenient and efficient in calculating. In comparison with traditional approach, the implementation of Delaunay triangulation helps us to use various triangles to approximate the geometry of geological body, and varying density is more closo to the real density distribution when figuring out 2-D geological body with complicated geometry and heterogeneous density. Here, we implement our method in inversing a gravity profile measured in Qianmiqiao Oil Field, analyze the result, compare our method with the conventional, and make a discussion.

关 键 词：非规则重力建模 DELAUNAY剖分 变密度 二度体 重力异常 

分 类 号：P631[天文地球—地质矿产勘探]