Three-dimensional crustal stress numerical simulation method based on deep and shallow dynamic bidirectional coupling

Abstract

The invention relates to the technical field of numerical simulation of earth crust dynamics, and provides a three-dimensional ground stress numerical simulation method based on deep and shallow dynamic bidirectional coupling, which comprises the following steps: constructing a three-dimensional geometric model of a target area and configuring deep rock elastic-plastic rheological parameters and a superficial erosion rate equation; generating a deep unstructured grid and a superficial irregular grid by adopting a differentiated grid strategy, and setting adaptive time step control; a nonlinear solver is configured to achieve staged solution of a thermal-mechanical coupling equation, and cross-scale data transmission is completed through bilinear interpolation; operating a simulation output material attribute field, a pressure field and a strain rate field; and the result is converted into a standardized format, and the stress concentration area is visually marked. Through a deep and shallow dynamic bidirectional coupling mechanism, synchronous simulation of structure driving and surface response is realized, and simulation precision and engineering applicability of earth crust dynamic evolution are remarkably improved.

Images