A method of reconstructing complex stratigraphic surfaces with multitype fault constraints
Deng Shi-Wu1, Jia Yu1, Yao Xing-Miao2, and Liu Zhi-Ning2
1. School of Nuclear Technology and Automation Engineering, Chengdu University of Technology, Chengdu 610059, China.
2. School of Resources and Environment, University of Electronic and Technology of China, Chengdu 611731, China.
Abstract The construction of complex stratigraphic surfaces is widely employed in many fields, such as petroleum exploration, geological modeling, and geological structure analysis. It also serves as an important foundation for data visualization and visual analysis in these fields. The existing surface construction methods have several deficiencies and face various difficulties, such as the presence of multitype faults and roughness of resulting surfaces. In this paper, a surface modeling method that uses geometric partial differential equations (PDEs) is introduced for the construction of stratigraphic surfaces. It effectively solves the problem of surface roughness caused by the irregularity of stratigraphic data distribution. To cope with the presence of multitype complex faults, a two-way projection algorithm between three-dimensional space and a two-dimensional plane is proposed. Using this algorithm, a unified method based on geometric PDEs is developed for dealing with multitype faults. Moreover, the corresponding geometric PDE is derived, and an algorithm based on an evolutionary solution is developed. The algorithm proposed for constructing spatial surfaces with real data verifies its computational efficiency and its ability to handle irregular data distribution. In particular, it can reconstruct faulty surfaces, especially those with overthrust faults.
This work was financially supported by the National Natural Science foundation of China (No. U1562218).
Cite this article:
. A method of reconstructing complex stratigraphic surfaces with multitype fault constraints[J]. APPLIED GEOPHYSICS, 0, (): 195-204.
[1]
Cai, Q., Yang, Q., and Chen, Q. M., 2004, Conforming delaunay triangulation for geological structure with overlapping domains: Journal of Computer-Aided Design & Computer Graphics, 16(6), 766−771.
[2]
Carmo M. P., 1994, Differential geometry of surfaces: Differential Forms and Applications, 77−98.
[3]
Du, H., and Qin, H., 2007, Free-form geometric modeling by integrating parametric and implicit PDEs: IEEE Trans: The Visual Computer, 13(3), 549−561.
[4]
Aurenhammer, F, 1991, Voronoi diagrams—a survey of a fundamental geometric data structure: ACM Computing Surveys (CSUR), 23(3), 345-405.
[5]
Frank, T., Tertois, A. L., and Mallet, J. L., 2007, 3D-reconstruction of complex geological interfaces from irregularly distributed and noisy point data: Computers & Geosciences, 33(7), 932−943.
[6]
Farshbaf, P. S., Khatib, M. M., and Nazari, H., 2016, Solid meshing of 3D geological model in finite element analysis: a case study of east azerbaijan, NW Iran: Modeling Earth Systems and Environment, 2(1), 1−7.
[7]
Bloor, M., and Wilson, M., 1989, Generating blend surfaces using partial differential equations: Computer Aided Design, 21(3), 165−171.
[8]
Frey, P. J., Borouchaki, H., and George, P. L., 1998, 3D Delaunay mesh generation coupled with an advancing-front approach: Computer Methods in Applied Mechanics and Engineering, 157(1-2), 115−131.
[9]
Cai, Q., Yang, Q., and Chen, Q. M., 2004, Conforming delaunay triangulation for geological structure with overlapping domains: Journal of Computer-Aided Design & Computer Graphics, 16(6), 766−771.
[10]
Carmo M. P., 1994, Differential geometry of surfaces: Differential Forms and Applications, 77−98.
[11]
Hillier, M. J, Schetselaar, E. M., de Kemp, E. A., et al. 2014, Three-dimensional modelling of geological surfaces using generalized interpolation with radial basis functions: Mathematical Geosciences, 46(8), 931−953.
[12]
Jia, Y., Deng, S. W., and Yao, X., 2015, Kriging interpolation algorithm based on constraint particle swarm optimization: Journal of Chengdu University of Technology, 42(1), 104−109.
[13]
Kuwert, E., and Schätzle, R., 2001, The Willmore flow with small initial energy: J. Diff. Geom., 57(3), 409−441.
[14]
Du, H., and Qin, H., 2007, Free-form geometric modeling by integrating parametric and implicit PDEs: IEEE Trans: The Visual Computer, 13(3), 549−561.
[15]
Frank, T., Tertois, A. L., and Mallet, J. L., 2007, 3D-reconstruction of complex geological interfaces from irregularly distributed and noisy point data: Computers & Geosciences, 33(7), 932−943.
[16]
Farshbaf, P. S., Khatib, M. M., and Nazari, H., 2016, Solid meshing of 3D geological model in finite element analysis: a case study of east azerbaijan, NW Iran: Modeling Earth Systems and Environment, 2(1), 1−7.
[17]
Frey, P. J., Borouchaki, H., and George, P. L., 1998, 3D Delaunay mesh generation coupled with an advancing-front approach: Computer Methods in Applied Mechanics and Engineering, 157(1-2), 115−131.
[18]
Li, M. C., and Miu, Z., 2011, 3D interpolation-approximation fitting construction method for complex geological surfaces: Engineering Science, 13(12), 103−107.
[19]
Hillier, M. J, Schetselaar, E. M., de Kemp, E. A., et al. 2014, Three-dimensional modelling of geological surfaces using generalized interpolation with radial basis functions: Mathematical Geosciences, 46(8), 931−953.
[20]
Jia, Y., Deng, S. W., and Yao, X., 2015, Kriging interpolation algorithm based on constraint particle swarm optimization: Journal of Chengdu University of Technology, 42(1), 104−109.
[21]
Kuwert, E., and Schätzle, R., 2001, The Willmore flow with small initial energy: J. Diff. Geom., 57(3), 409−441.
[22]
Liu, S., Hu, X., Xi, Y., 2015, 2D inverse modeling for potential fields on rugged observation surface using constrained delaunay triangulation: Computers & Geosciences, 76, 18−30.
[23]
Liu, D., Xu, G., and Zhang, Q., 2008, A discrete scheme of Laplace-Beltrami operator and its convergence over quadrilateral meshes: Computers and Mathematics with Applications, 55(6), 1081−1093
[24]
Li, M. C., and Miu, Z., 2011, 3D interpolation-approximation fitting construction method for complex geological surfaces: Engineering Science, 13(12), 103−107.
[25]
Liu, S., Hu, X., Xi, Y., 2015, 2D inverse modeling for potential fields on rugged observation surface using constrained delaunay triangulation: Computers & Geosciences, 76, 18−30.
[26]
Mallet, J. L., 1989, Discrete smooth interpolation: ACM Transactions on Graphics, 8(2), 121−144.
[27]
Marghany, M., 2012, Fuzzy B-spline algorithm for 3-D lineament reconstruction: Int. J. Phys. Sci., 7(15), 2294−301.
[28]
Liu, D., Xu, G., and Zhang, Q., 2008, A discrete scheme of Laplace-Beltrami operator and its convergence over quadrilateral meshes: Computers and Mathematics with Applications, 55(6), 1081−1093
[29]
Meyer, M., and Desbrun, M., Schroder, P., 2002, Discrete differential-geometry operators for triangulated 2-manifolds: Proceedings of Visual Mathematics’02. Berlin, Germany.
[30]
Mallet, J. L., 1989, Discrete smooth interpolation: ACM Transactions on Graphics, 8(2), 121−144.
[31]
Marghany, M., 2012, Fuzzy B-spline algorithm for 3-D lineament reconstruction: Int. J. Phys. Sci., 7(15), 2294−301.
[32]
Meyer, M., and Desbrun, M., Schroder, P., 2002, Discrete differential-geometry operators for triangulated 2-manifolds: Proceedings of Visual Mathematics’02. Berlin, Germany.
[33]
Ma, L., and Zhu, X., 1998, Application of PDE method to free form surface design: Chinese Journal of Computers 21(3), 357−362.
[34]
Ma, L., and Zhu, X., 1998, Application of PDE method to free form surface design: Chinese Journal of Computers 21(3), 357−362.
[35]
Pellerin, J., Lévy, B., and Caumon, G., 2014, Automatic surface remeshing of 3D structural models at specified resolution: A method based on Voronoi diagrams: Computers & Geosciences, 62, 103−116.
[36]
Pellerin, J., Lévy, B., and Caumon, G., 2014, Automatic surface remeshing of 3D structural models at specified resolution: A method based on Voronoi diagrams: Computers & Geosciences, 62, 103−116.
[37]
Plummer, C. C., McGeary, D., and Carlson D. H., 1991, Physical geology: Wm. C. Brown.
[38]
Plummer, C. C., McGeary, D., and Carlson D. H., 1991, Physical geology: Wm. C. Brown.
[39]
Schneider, R., and Kobbelt, L., 2001, Geometric fairing of irregular meshes for free-form surface design: Computer Aided Geometric Design, 18(4), 359−379.
[40]
Schneider, R., and Kobbelt, L., 2001, Geometric fairing of irregular meshes for free-form surface design: Computer Aided Geometric Design, 18(4), 359−379.
[41]
White, B., 2002, Evolution of curves and surfaces by mean curvature. In: Proceedings of the International Congress of Mathematicians, vol. I, Beijing. 525−538.
Yu, Z. W., 1987, A New Method for interpolating geological surface: Journal of China University of Mining & Technology, 16(4), 69−76.
[44]
Yan, H. W., and Wu, J. P., 2012, Research on genetic algorithm kriging optimized by weight least square: Computer Technology and Development, 22(3), 92−95.
[45]
Yao, X., Wang, Q., and Liu, Z., 2015, Fast 3D kriging interpolation using delaunay tetrahedron with CUDA-enabled GPU: 85th Annual International Meeting, SEG, Expanded Abstracts, 1739−1743.
[46]
Yao, X., Luo, C., Hu, G., et al., 2014, A method of geological surface reconstruction with reverse fault, 84th Annual International Meeting, SEG, Expanded Abstracts, 1559−1564.
[47]
White, B., 2002, Evolution of curves and surfaces by mean curvature. In: Proceedings of the International Congress of Mathematicians, vol. I, Beijing. 525−538.
[48]
Zhong, D. H., Li, M. C., and Song, L. G., 2006, Enhanced NURBS modeling and visualization for large 3D geoengineering applications: an example from the Jinping first-level hydropower engineering project, China: Computers & Geosciences, 32(9), 1270−1282.
Zeng, J., 2005, Some Questions of partial differential equation method in surface modeling: Master Thesis, Dalian University of Technology.
[51]
Yu, Z. W., 1987, A New Method for interpolating geological surface: Journal of China University of Mining & Technology, 16(4), 69−76.
[52]
Yan, H. W., and Wu, J. P., 2012, Research on genetic algorithm kriging optimized by weight least square: Computer Technology and Development, 22(3), 92−95.
[53]
Yao, X., Wang, Q., and Liu, Z., 2015, Fast 3D kriging interpolation using delaunay tetrahedron with CUDA-enabled GPU: 85th Annual International Meeting, SEG, Expanded Abstracts, 1739−1743.
[54]
Yao, X., Luo, C., Hu, G., et al., 2014, A method of geological surface reconstruction with reverse fault, 84th Annual International Meeting, SEG, Expanded Abstracts, 1559−1564.
[55]
Zhong, D. H., Li, M. C., and Song, L. G., 2006, Enhanced NURBS modeling and visualization for large 3D geoengineering applications: an example from the Jinping first-level hydropower engineering project, China: Computers & Geosciences, 32(9), 1270−1282.
[56]
Zeng, J., 2005, Some Questions of partial differential equation method in surface modeling: Master Thesis, Dalian University of Technology.
[57]
Zheng, C., 2005, Some research on the generating and manipulating of PDE transition surface: Master Thesis, Dalian University of Technology.
[58]
Zhao, H. X., and Xu, G., 2006, Triangular surface mesh fairing via Gaussian curvature flow: Journal of Computational and Applied Mathematics, 195(1), 300−311.
[59]
Zinck, G, Donias, M, Daniel, J., et al., 2013, Fast seismic horizon reconstruction based on local dip transformation: Journal of Applied Geophysics, 96, 11−18.
[60]
Zheng, C., 2005, Some research on the generating and manipulating of PDE transition surface: Master Thesis, Dalian University of Technology.
[61]
Zinck, G., Donias, M., and Guillon, S., 2014, Discontinuous seismic horizon reconstruction based on local dip transformation: 2014 IEEE International Conference on Image Processing (ICIP), 5856−5860.
[62]
Zhao, H. X., and Xu, G., 2006, Triangular surface mesh fairing via Gaussian curvature flow: Journal of Computational and Applied Mathematics, 195(1), 300−311.
[63]
Zinck, G, Donias, M, Daniel, J., et al., 2013, Fast seismic horizon reconstruction based on local dip transformation: Journal of Applied Geophysics, 96, 11−18.
[64]
Zinck, G., Donias, M., and Guillon, S., 2014, Discontinuous seismic horizon reconstruction based on local dip transformation: 2014 IEEE International Conference on Image Processing (ICIP), 5856−5860.
0, Vol. 
