基于水平集方法的复杂地质速度模型表达

张建兴; 杨钦; 孟宪海; 李吉刚

doi:10.13700/j.bh.1001-5965.2015.0106

基于水平集方法的复杂地质速度模型表达

doi: 10.13700/j.bh.1001-5965.2015.0106

北京航空航天大学计算机学院, 北京 100083

基金项目: 国家自然科学基金(61003110);国家科技重大专项(2011ZX05008-004-10)

详细信息

作者简介:
张建兴男,博士研究生。主要研究方向:计算机图形学、三维建模、层析速度建模。Tel.:010-82338474 E-mail:zjx7906@163.com;孟宪海男,博士,讲师。主要研究方向:计算机图形学、计算几何、计算机可视化。Tel.:010-82316283 E-mail:x_h_meng@buaa.edu.cn

通讯作者:
孟宪海,Tel.:010-82316283 E-mail:x_h_meng@buaa.edu.cn

中图分类号: TP391
计量
- 文章访问数: 868
- HTML全文浏览量: 157
- PDF下载量: 576
- 被引次数: 0
出版历程
- 收稿日期: 2015-03-04
- 网络出版日期: 2016-02-20

Representation of complex geological velocity model based on level set method

School of Computer Science and Engineering, Beijing University of Aeronautics and Astronautics, Beijing 100083, China

摘要

摘要: 地质速度模型表达的难点在于对复杂地质现象的精确描述和对正、反演应用的有效支持。提出了一种将速度间断面信息以隐函数的零等值面形式内嵌于速度模型的表达方法,多重拷贝策略使得相邻的速度块体互不干涉,妥善处理了在速度间断面处插值速度信息的矛盾。该表达方法适于描述复杂地质构造下的速度模型,如多值层面和具有非流形性质的断裂层面。在基于射线理论的应用中,利用符号距离场和速度间断面二叉树能迅速获得空间任意点的速度信息,检测射线与速度间断面的相交关系,并使射线在速度间断面处严格执行Snell定律,提高了射线路径的正确性。两组射线追踪的应用实例验证该表达方法的有效性。
- 速度模型 /
- 水平集 /
- 多重拷贝 /
- 速度间断面 /
- 二叉树
Abstract: The difficulties in the expression of geological velocity model are the accurate description of complex geological phenomena and the effective support for the forward and inverse applications. We presented a novel representation method for describing velocity models. Velocity discontinuity was embedded into model in the form of zero-value contour of an implicit function. Multi-copy strategy was to properly handle the contradictory interpolation of velocity value at the boundary of discontinuities. The representation method can be applied to describe very complicated velocity models such as the faulted interface with non-manifold nature and multiple value surfaces. In applications based on the ray theory, signed distance field and velocity discontinuity binary tree enable fast access to velocity information of arbitrary point and easy detection of ray-interface intersection. The model can make the ray propagate positively and honor Snell's law strictly when facing velocity interface and thus improve the accuracy of ray path. Two groups of ray tracing examples verify the effectiveness of the method.
- velocity model /
- level set /
- multi-copy strategy /
- velocity discontinuity /
- binary tree

HTML全文

参考文献(16)

[1]	YILMAZ Ö,DOHERTY S M.Seismic data analysis:Processing,inversion,and interpretation of seismic data[M].Tulsa:Society Exploration Geophysicists,2000:14.
[2]	LEVANDER A R,NOLET G,AMERICAN G U.Seismic earth:Array analysis of broadband seismograms[M].California:American Geophysical Union,2005:49-65.
[3]	孟宪海,杨钦,李吉刚.基于层面结构的三维闭合地质区块构造算法[J].北京航空航天大学学报,2005,31(2):182-186. MENG X H,YANG Q,LI J G.Construction of coherent 3D geological blocks from stratified geological structure[J].Journal of Beijing University of Aeronautics and Astronautics,2005,31(2):182-186(in Chinese).
[4]	ZELT C,SMITH R.Seismic traveltime inversion for 2-D crustal velocity structure[J].Geophysical Journal International,1992,108(1):16-34.
[5]	PEREYRA V.Modeling,ray tracing,and block nonlinear travel-time inversion in 3D[J].Pure and Applied Geophysics,1996,148(3):345-386.
[6]	XU T,ZHANG Z,GAO E,et al.Segmentally iterative ray tracing in complex 2D and 3D heterogeneous block models[J].Bulletin of the Seismological Society of America,2010,100(2):841-850.
[7]	张俊安,刘瑞刚,杨钦,等.复杂地质结构的四维地质层面自动生成算法[J].北京航空航天大学学报,2007,33(9):1094-1098. ZHANG J A,LIU R G,YANG Q,et al.Auto-construction algorithm of four-dimensional(4D) geological stratum on complex geological structure[J].Journal of Beijing University of Aeronautics and Astronautics,2007,33(9):1094-1098(in Chinese).
[8]	OSHER S,SETHIAN J A.Fronts propagating with curvature-dependent speed:Algorithms based on Hamilton-Jacobi formulations[J].Journal of Computational Physics,1988,79(1):12-49.
[9]	OSHER S,FEDKIW R.Level set methods and dynamic implicit surfaces[M].New York:Springer,2003:7.
[10]	ZHAO H K,OSHER S,FEDKIW R.Fast surface reconstruction using the level set method[C]//Proceedings of IEEE Workshop on Variational and Level Set Methods in Computer Vision,2001.Piscataway,NJ:IEEE Press,2001:194-201.
[11]	孟宪海,李吉刚,杨钦,等.复杂限定Delaunay三角化算法[J].中国科学:信息科学,2010,40(3):381-392. MENG X H,LI J G,YANG Q,et al.Algorithm of complex conforming delaunay triangulation[J].Science China:Information Sciences,2010,40(3):381-392(in Chinese).
[12]	ERVENY V.Seismic ray theory[M].Cambridge:Cambridge University Press,2001:vii.
[13]	SAMBRIDGE M,BRAUN J,MCQUEEN H.Geophysical parametrization and interpolation of irregular data using natural neighbours[J].Geophysical Journal International,1995,122(3):837-857.
[14]	VINJE V,IVERSEN E,ÅSTEBØL K,et al.Estimation of multivalued arrivals in 3D models using wavefront construction-Part Ⅰ[J].Geophysical Prospecting,1996,44(5):819-842.
[15]	VINJE V,IVERSEN E,ÅSTEBØL K,et al.Estimation of multivalued arrivals in 3D models using wavefront construction-Part Ⅱ:Tracing and interpolation[J].Geophysical Prospecting,1996,44(5):843-858.
[16]	VINJE V,IVERSEN E,GJOYSTDAL H.Traveltime and amplitude estimation using wavefront construction[J].Geophysics,1993,58(8):1157-1166.