CAJ | 학술논문

针对层析成像模型矩形网格剖分存在的一些问题，提出了复杂结构三角网最小走时射线追踪层析成像方法。以 Delaunay 三角剖分的优化准则，根据模型点、线、面的几何结构关系，进行三角网格剖分；采用三角网波行面扩展法计算声波初至走时，追踪出任意接收点至激发点处的射线路径，射线路径包含了从接收点至源点的坐标、所在三角单元等信息，形成矩阵方程。应用稳定性较好的联合迭代重构技术求解矩阵方程，得到模型的速度分布。数值模拟结果表明，三角网射线层析成像方法分辨率高，成像结果更接近实际结构形态，有效地解决了矩形网剖分对复杂区域网格参数化灵活性差，速度界面描述精度低等问题。
침대층석성상모형구형망격부분존재적일사문제，제출료복잡결구삼각망최소주시사선추종층석성상방법。이 Delaunay 삼각부분적우화준칙，근거모형점、선、면적궤하결구관계，진행삼각망격부분；채용삼각망파행면확전법계산성파초지주시，추종출임의접수점지격발점처적사선로경，사선로경포함료종접수점지원점적좌표、소재삼각단원등신식，형성구진방정。응용은정성교호적연합질대중구기술구해구진방정，득도모형적속도분포。수치모의결과표명，삼각망사선층석성상방법분변솔고，성상결과경접근실제결구형태，유효지해결료구형망부분대복잡구역망격삼수화령활성차，속도계면묘술정도저등문제。
In view of problems existed in tomography on rectangular net subdividing to a model, we proposed a triangular net minimum travel time ray tracing tomography method for complex structure. A triangle grid was subdivided based on the Delaunay triangular subdividing criterion and the relationships among the points, lines, and surfaces of the model. The first travel times was calculated with triangular net plane wave propagation method, followed by tracing the ray paths of any receiving point to the excitation point. The ray paths include the informations of the coordinates from receiving points to source points and all triangle cells, and the matrix equation is formed. The matrix equation is conducted by applying SIRT algorithm and the velocity distribution is obtained. The simulation results show that the use of triangular net ray tomography method achieved higher resolution version, and more precisely reflected the actual morphological characteristics. The problems of inflexibility of subdividing elastic wave tomography for complex structure models, and poor definition of velocity interfaces is effectively resolved.